Cloning and expression analysis of the mouse stroma marker Snep encoding a novel nidogen domain protein

Computational and experimental characterization of the novel ECM glycoprotein SNED1 and prediction of its interactome

The extracellular matrix (ECM) is a complex meshwork of proteins and an essential component of multicellular life. We have recently reported the characterization of a novel ECM protein, SNED1, and showed that it promotes breast cancer metastasis and regulates craniofacial development. However, the mechanisms by which it does so remain unknown. ECM proteins exert their functions by binding to cell surface receptors and interacting with other ECM proteins, actions that we can predict using knowledge of protein's sequence, structure, and post-translational modifications. Here, we combined in-silico and in-vitro approaches to characterize the physico-chemical properties of SNED1 and infer its putative functions. To do so, we established a mammalian cell system to produce and purify SNED1 and its N-terminal fragment, which contains a NIDO domain, and demonstrated experimentally SNED1's potential to be glycosylated, phosphorylated, and incorporated into an insoluble ECM. We also determined the secondary and tertiary structures of SNED1 and its N-terminal fragment and obtained a model for its NIDO domain. Using computational predictions, we identified 114 proteins as putative SNED1 interactors, including the ECM protein fibronectin. Pathway analysis of the predicted SNED1 interactome further revealed that it may contribute to signaling through cell surface receptors, such as integrins, and participate in the regulation of ECM organization and developmental processes. Last, using fluorescence microscopy, we showed that SNED1 forms microfibrils within the ECM and partially colocalizes with fibronectin. Altogether, we provide a wealth of information on an understudied yet important ECM protein with the potential to decipher its pathophysiological functions.

PCB126 Exposure Revealed Alterations in m6A RNA Modifications in Transcripts Associated With AHR Activation

Chemical modifications of proteins, DNA, and RNA moieties play critical roles in regulating gene expression. Emerging evidence suggests the RNA modifications (epitranscriptomics) have substantive roles in basic biological processes. One of the most common modifications in mRNA and noncoding RNAs is N6-methyladenosine (m6A). In a subset of mRNAs, m6A sites are preferentially enriched near stop codons, in 3' UTRs, and within exons, suggesting an important role in the regulation of mRNA processing and function including alternative splicing and gene expression. Very little is known about the effect of environmental chemical exposure on m6A modifications. As many of the commonly occurring environmental contaminants alter gene expression profiles and have detrimental effects on physiological processes, it is important to understand the effects of exposure on this important layer of gene regulation. Hence, the objective of this study was to characterize the acute effects of developmental exposure to PCB126, an environmentally relevant dioxin-like PCB, on m6A methylation patterns. We exposed zebrafish embryos to PCB126 for 6 h starting from 72 h post fertilization and profiled m6A RNA using methylated RNA immunoprecipitation followed by sequencing (MeRIP-seq). Our analysis revealed 117 and 217 m6A peaks in the DMSO and PCB126 samples (false discovery rate 5%), respectively. The majority of the peaks were preferentially located around the 3' UTR and stop codons. Statistical analysis revealed 15 m6A marked transcripts to be differentially methylated by PCB126 exposure. These include transcripts that are known to be activated by AHR agonists (eg, ahrra, tiparp, nfe2l2b) as well as others that are important for normal development (vgf, cebpd, sned1). These results suggest that environmental chemicals such as dioxin-like PCBs could affect developmental gene expression patterns by altering m6A levels. Further studies are necessary to understand the functional consequences of exposure-associated alterations in m6A levels.

Knockout of the gene encoding the extracellular matrix protein SNED1 results in early neonatal lethality and craniofacial malformations

BACKGROUND

The extracellular matrix (ECM) is a fundamental component of multicellular organisms that orchestrates developmental processes and controls cell and tissue organization. We previously identified the novel ECM protein SNED1 as a promoter of breast cancer metastasis and showed that its level of expression negatively correlated with breast cancer patient survival. Here, we sought to identify the roles of SNED1 during murine development.

RESULTS

We generated two novel Sned1 knockout mouse strains and showed that Sned1 is essential since homozygous ablation of the gene led to early neonatal lethality. Phenotypic analysis of the surviving knockout mice revealed a role for SNED1 in the development of craniofacial and skeletal structures since Sned1 knockout resulted in growth defects, nasal cavity occlusion, and craniofacial malformations. Sned1 is widely expressed in embryos, notably by cell populations undergoing epithelial-to-mesenchymal transition, such as the neural crest cells. We further show that mice with a neural-crest-cell-specific deletion of Sned1 survive, but display facial anomalies partly phenocopying the global knockout mice.

CONCLUSIONS

Our results demonstrate requisite roles for SNED1 during development and neonatal survival. Importantly, the deletion of 2q37.3 in humans, a region that includes the SNED1 locus, has been associated with facial dysmorphism and short stature.

Epigenome-wide association study reveals CpG sites related to COG of neuroblastoma

Background. Neuroblastoma (NB) is the most common extracranial solid tumor in infants and children. Its variable location and complex pathogenesis make NB hard for early diagnosis and risk classification. Methodology. We analyzed the methylation data of 236 samples from patients with NB in Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. Kaplan–Meier survival analysis was used for comparing overall survival of NB patients in different groups. Epigenome-wide association study (EWAS) was conducted to screen CpGs significantly associated with NB patients’ Children’s Oncology Group (COG). Logistic regression method was used for constructing a model to predict NB patients’ COG. Results. NB patients in low COG showed significantly superior prognosis than those in high COG. A total of seven CpG sites were found closely related to COG. Logistic regression model based on those CpGs showed superior performance in separating NB patients in different COGs. Conclusions. The present study highlights the important role of DNA methylation in NB development, which might provide evidence for treatment decisions for children NB.

Overactivity or blockade of transforming growth factor-β each generate a specific ureter malformation

Transforming growth factor-β (TGFβ) has been reported to be dysregulated in malformed ureters. There exists, however, little information on whether altered TGFβ levels actually perturb ureter development. We therefore hypothesised that TGFβ has functional effects on ureter morphogenesis. Tgfb1, Tgfb2 and Tgfb3 transcripts coding for TGFβ ligands, as well as Tgfbr1 and Tgfbr2 coding for TGFβ receptors, were detected by quantitative polymerase chain reaction in embryonic mouse ureters collected over a wide range of stages. As assessed by in situ hybridisation and immunohistochemistry, the two receptors were detected in embryonic urothelia. Next, TGFβ1 was added to serum-free cultures of embryonic day 15 mouse ureters. These organs contain immature smooth muscle and urothelial layers and their in vivo potential to grow and acquire peristaltic function can be replicated in serum-free organ culture. Such organs therefore constitute a suitable developmental stage with which to define roles of factors that affect ureter growth and functional differentiation. Exogenous TGFβ1 inhibited growth of the ureter tube and generated cocoon-like dysmorphogenesis. RNA sequencing suggested that altered levels of transcripts encoding certain fibroblast growth factors (FGFs) followed exposure to TGFβ. In serum-free organ culture exogenous FGF10 but not FGF18 abrogated certain dysmorphic effects mediated by exogenous TGFβ1. To assess whether an endogenous TGFβ axis functions in developing ureters, embryonic day 15 explants were exposed to TGFβ receptor chemical blockade; growth of the ureter was enhanced, and aberrant bud-like structures arose from the urothelial tube. The muscle layer was attenuated around these buds, and peristalsis was compromised. To determine whether TGFβ effects were limited to one stage, explants of mouse embryonic day 13 ureters, more primitive organs, were exposed to exogenous TGFβ1, again generating cocoon-like structures, and to TGFβ receptor blockade, again generating ectopic buds. As for the mouse studies, immunostaining of normal embryonic human ureters detected TGFβRI and TGFβRII in urothelia. Collectively, these observations reveal unsuspected regulatory roles for endogenous TGFβ in embryonic ureters, fine-tuning morphogenesis and functional differentiation. Our results also support the hypothesis that the TGFβ up-regulation reported in ureter malformations impacts on pathobiology. Further experiments are needed to unravel the intracellular signalling mechanisms involved in these dysmorphic responses. © 2019 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Molecular patterning of the embryonic cranial mesenchyme revealed by genome-wide transcriptional profiling

In the head of an embryo, a layer of mesenchyme surrounds the brain underneath the surface ectoderm. This cranial mesenchyme gives rise to the meninges, the calvaria (top part of the skull), and the dermis of the scalp. Abnormal development of these structures, especially the meninges and the calvaria, is linked to significant congenital defects in humans. It has been known that different areas of the cranial mesenchyme have different fates. For example, the calvarial bone develops from the cranial mesenchyme on the baso-lateral side of the head just above the eye (supraorbital mesenchyme, SOM), but not from the mesenchyme apical to SOM (early migrating mesenchyme, EMM). However, the molecular basis of this difference is not fully understood. To answer this question, we compared the transcriptomes of EMM and SOM using high-throughput sequencing (RNA-seq). This experiment identified a large number of genes that were differentially expressed in EMM and SOM, and gene ontology analyses found very different terms enriched in each region. We verified the expression of about 40 genes in the head by RNA in situ hybridization, and the expression patterns were annotated to make a map of molecular markers for 6 subdivisions of the cranial mesenchyme. Our data also provided insights into potential novel regulators of cranial mesenchyme development, including several axon guidance pathways, lectin complement pathway, cyclic-adenosine monophosphate (cAMP) signaling pathway, and ZIC family transcription factors. Together, information in this paper will serve as a unique resource to guide future research on cranial mesenchyme development.

Construction of Novel DNA Methylation-Based Prognostic Model to Predict Survival in Glioblastoma

Glioblastoma (GBM) is a most aggressive primary cancer in brain with poor prognosis. This study aimed to identify novel tumor biomarkers with independent prognostic values in GBMs. The DNA methylation profiles were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus database. Differential methylated genes (DMGs) were screened from recurrent GBM samples using limma package in R software. Functional enrichment analysis was performed to identify major biological processes and signaling pathways. Furthermore, critical DMGs associated with the prognosis of GBM were screened according to univariate and multivariate cox regression analysis. A risk score-based prognostic model was constructed for these DMGs and prediction ability of this model was validated in training and validation data set. In total, 495 DMGs were identified between recurrent samples and disease-free samples, including 356 significantly hypermethylated and 139 hypomethylated genes. Functional and pathway items for these DMGs were mainly related to sensory organ development, neuroactive ligand-receptor interaction, pathways in cancer, etc. Five genes with abnormal methylation level were significantly correlated with prognosis according to survival analysis, such as ALX1, KANK1, NUDT12, SNED1, and SVOP. Finally, the risk model provided an effective ability for prognosis prediction both in training and validation data set. We constructed a novel prognostic model for survival prediction of GBMs. In addition, we identified five DMGs as critical prognostic biomarkers in GBM progression.

Structure and Function of Basement Membranes

Basement membranes (BMs) are present in every tissue of the human body. All epithelium and endothelium is in direct association with BMs. BMs are a composite of several large glycoproteins and form an organized scaffold to provide structural support to the tissue and also offer functional input to modulate cellular function. While collagen I is the most abundant protein in the human body, type IV collagen is the most abundant protein in BMs. Matrigel is commonly used as surrogate for BMs in many experiments, but this is a tumor-derived BM–like material and does not contain all of the components that natural BMs possess. The structure of BMs and their functional role in tissues are unique and unlike any other class of proteins in the human body. Increasing evidence suggests that BMs are unique signal input devices that likely fine tune cellular function. Additionally, the resulting endothelial and epithelial heterogeneity in human body is a direct contribution of cell-matrix interaction facilitated by the diverse compositions of BMs.

Genome-wide association study identifies common variants associated with pharmacokinetics of psychotropic drugs

Individual variation in pharmacokinetics of psychotropic drugs, particularly metabolism, is an important factor to consider in pharmacological treatment in psychiatry. A large proportion of this variance is still not accounted for, but evidence so far suggests the involvement of genetic factors. We performed a genome-wide association study (GWAS) with concentration dose ratio (CDR) as sub-phenotype to assess metabolism rate of psychotropic drugs in a homogenous Norwegian sample of 1334 individuals diagnosed with a severe mental disorder. The GWAS revealed one genome-wide significant marker (rs16935279, p-value=3.95×10(-10), pperm=7.5×10(-4)) located in an intronic region of the lncRNA LOC100505718. Carriers of the minor allele have a lower metabolism rate of antiepileptic drugs compared to major allele carriers. In addition, several nominally significant associations between single nucleotide polymorphisms (SNPs) and CDR for antipsychotic, antidepressant and antiepileptic drugs were disclosed. We consider standardised CDR to be a useful measure of the metabolism rate of a drug. The present findings indicate that common gene variants could affect the metabolism of psychotropic drugs. This warrants further investigations into the functional mechanisms involved as it may lead to identification of predictive markers as well as novel drug targets.

Associations of chemo- and radio-resistant phenotypes with the gap junction, adhesion and extracellular matrix in a three-dimensional culture model of soft sarcoma

Background

Three-dimensional (3D) culture models are considered to recapitulate the cell microenvironment in solid tumors, including the extracellular matrix (ECM), cell-cell interactions, and signal transduction. These functions are highly correlated with cellular behaviors and contribute to resistances against chemo- and radio-therapies. However, the biochemical effects and mechanisms remain unknown in soft sarcoma. Therefore, we developed an in vitro 3D model of sarcoma to analyze the reasons of the chemo- and radio-resistance in therapies.

Methods

Four soft sarcoma cell lines, HT1080, RD, SW872, and human osteosarcoma cell line 1 (HOSS1), a cell line established from a patient-derived xenograft, were applied to 3D culture and treated with growth factors in methylcellulose-containing medium. Spheroids were examined morphologically and by western blotting, RT-qPCR, and immunofluorescence staining to analyze cell adhesion, gap junctions, ECM genes, and related factors. Proliferation and colony formation assays were performed to assess chemo- and radio-resistances between 3D and two-dimensional (2D) cell cultures. Annexin V and Propidium Iodide staining was used to detect early apoptotic sarcoma cells treated with Doxorubicin, Gemcitabine, and Docetaxel in the 3D model.

Results

The four soft sarcoma cell lines formed spheres in vitro by culture in modified condition medium. Compared with 2D cell culture, expression of ECM genes and proteins, including COL1A1, LOX, SED1, FN1, and LAMA4, was significantly increased in 3D culture. Analysis of cadherin and gap junction molecules showed significant changes in the gene and protein expression profiles under 3D conditions. These changes affected cell–cell communication and were mainly associated with biological processes such as cell proliferation and apoptosis related to chemo- and radio-resistances.

Conclusions

Our findings revealed significant differences between 3D and 2D cell culture systems, and indicated that cellular responsiveness to external stress such as radiation and chemotherapeutics is influenced by differential expression of genes and proteins involved in regulation of the ECM, cell adhesion, and gap junction signaling.

Electronic supplementary material

The online version of this article (doi:10.1186/s13046-015-0175-0) contains supplementary material, which is available to authorized users.

Serum protein S100A9, SOD3, and MMP9 as new diagnostic biomarkers for pulmonary tuberculosis by iTRAQ-coupled two-dimensional LC-MS/MS

This study aimed to discover the novel noninvasive biomarkers for the diagnosis of pulmonary tuberculosis (TB). We applied iTRAQ 2D LC-MS/MS technique to investigate protein profiles in patients with pulmonary TB and other lung diseases. A total of 34 differentially expressed proteins (24 upregulated proteins and ten downregulated proteins) were identified in the serum of pulmonary TB patients. Significant differences in protein S100-A9 (S100A9), extracellular superoxide dismutase [Cu-Zn] (SOD3), and matrix metalloproteinase 9 (MMP9) were found between pulmonary TB and other lung diseases by ELISA. Correlations analysis revealed that the serum concentration of MMP9 in the pulmonary TB was in moderate correlation with SOD3 (r = 0.581) and S100A9 (r = 0.471), while SOD3 was in weak correlation with S100A9 (r = 0.287). The combination of serum S100A9, SOD3, and MMP9 levels could achieve 92.5% sensitivity and 95% specificity to discriminate between pulmonary TB and healthy controls, 90% sensitivity and 87.5% specificity to discriminate between pulmonary TB and pneumonia, and 85% sensitivity and 92.5% specificity to discriminate between pulmonary TB and lung cancer, respectively. The results showed that S100A9, SOD3, and MMP9 may be potential diagnostic biomarkers for pulmonary TB, and provided experimental basis for the diagnosis of pulmonary TB.

Extracellular matrix signatures of human mammary carcinoma identify novel metastasis promoters

The extracellular matrix (ECM) is a major component of tumors and a significant contributor to cancer progression. In this study, we use proteomics to investigate the ECM of human mammary carcinoma xenografts and show that primary tumors of differing metastatic potential differ in ECM composition. Both tumor cells and stromal cells contribute to the tumor matrix and tumors of differing metastatic ability differ in both tumor- and stroma-derived ECM components. We define ECM signatures of poorly and highly metastatic mammary carcinomas and these signatures reveal up-regulation of signaling pathways including TGFβ and VEGF. We further demonstrate that several proteins characteristic of highly metastatic tumors (LTBP3, SNED1, EGLN1, and S100A2) play causal roles in metastasis, albeit at different steps. Finally we show that high expression of LTBP3 and SNED1 correlates with poor outcome for ER⁻/PR⁻breast cancer patients. This study thus identifies novel biomarkers that may serve as prognostic and diagnostic tools.

DOI:http://dx.doi.org/10.7554/eLife.01308.001

Metastasis is the process whereby tumor cells spread within the body and is the cause of most deaths from cancer. This complex process involves several steps: first the cancer cells invade the tissues that surround the tumor; second, the cancer cells enter the blood stream and travel throughout the body; and third, the cancer cells seed the growth of new tumors in distant organs.

Within tissues, the extracellular matrix forms a complex scaffold of proteins that surrounds cells, to support and organize them: it also provides signals that control how much cells can multiply, how likely cells are to stick together or migrate, and even a cell’s chances of survival. Pathologists have used an accumulation of extracellular matrix proteins in tumors as a sign that the outcome of the disease will likely be unfavorable for a patient, and that treatment will be challenging. However, we still do not have a clear picture of the composition of the tumor extracellular matrix and we do not know all the details of how it affects tumor growth and metastasis.

Now, Naba et al. have explored these questions by injecting different types of human breast tumor cells into mice. Some of the cells were capable of spreading throughout the body and were said to have a high ‘metastatic potential’; others were less capable of spreading and were said to have a low metastatic potential. Naba et al. then analyzed the proteins that made up the extracellular matrix of the tumors that grew in the mice. Some proteins were found in both types of tumor; whereas some proteins were only found in the tumors with low metastatic potential and some were only found in the highly metastatic tumors. Naba et al. also demonstrated that both cancer cells and non-cancer cells—which are also found within the tumors—contributed to the production of the extracellular matrix in the tumor. Moreover, and somewhat surprisingly, the contributions from the non-cancer cells in the two types of tumors were also different.

Computational analysis predicted that the production of several extracellular matrix proteins in the highly metastatic tumors was under the control of signaling pathways that are involved in cancer progression. Furthermore, Naba et al. also demonstrated that several of the extracellular matrix proteins specific to highly metastatic tumors were required for the cancer to spread. These proteins are involved in different stages of the metastatic process, and some of them are commonly over-produced in tumors from patients with some of the worst chances of recovery.

If similar results are consistently observed in clinical samples from humans, the work of Naba et al. could help doctors to discriminate between tumors that will spread and those that will not, which should lead to improved patient care. The proteins and pathways associated with the highly metastatic tumors could be also investigated as potential drug targets.

DOI:http://dx.doi.org/10.7554/eLife.01308.002

3D pancreatic carcinoma spheroids induce a matrix-rich, chemoresistant phenotype offering a better model for drug testing

Background

Pancreatic ductal adenocarcinoma (PDAC) is the fourth most common cause of cancer related death. It is lethal in nearly all patients, due to an almost complete chemoresistance. Most if not all drugs that pass preclinical tests successfully, fail miserably in the patient. This raises the question whether traditional 2D cell culture is the correct tool for drug screening. The objective of this study is to develop a simple, high-throughput 3D model of human PDAC cell lines, and to explore mechanisms underlying the transition from 2D to 3D that might be responsible for chemoresistance.

Methods

Several established human PDAC and a KPC mouse cell lines were tested, whereby Panc-1 was studied in more detail. 3D spheroid formation was facilitated with methylcellulose. Spheroids were studied morphologically, electron microscopically and by qRT-PCR for selected matrix genes, related factors and miRNA. Metabolic studies were performed, and a panel of novel drugs was tested against gemcitabine.

Results

Comparing 3D to 2D cell culture, matrix proteins were significantly increased as were lumican, SNED1, DARP32, and miR-146a. Cell metabolism in 3D was shifted towards glycolysis. All drugs tested were less effective in 3D, except for allicin, MT100 and AX, which demonstrated effect.

Conclusions

We developed a high-throughput 3D cell culture drug screening system for pancreatic cancer, which displays a strongly increased chemoresistance. Features associated to the 3D cell model are increased expression of matrix proteins and miRNA as well as stromal markers such as PPP1R1B and SNED1. This is supporting the concept of cell adhesion mediated drug resistance.

Leveraging existing biological knowledge in the identification of candidate genes for facial dysmorphology

Background

In response to the frequently overwhelming output of high-throughput microarray experiments, we propose a methodology to facilitate interpretation of biological data in the context of existing knowledge. Through the probabilistic integration of explicit and implicit data sources a functional interaction network can be constructed. Each edge connecting two proteins is weighted by a confidence value capturing the strength and reliability of support for that interaction given the combined data sources. The resulting network is examined in conjunction with expression data to identify groups of genes with significant temporal or tissue specific patterns. In contrast to unstructured gene lists, these networks often represent coherent functional groupings.

Results

By linking from shared functional categorizations to primary biological resources we apply this method to craniofacial microarray data, generating biologically testable hypotheses and identifying candidate genes for craniofacial development.

Conclusion

The novel methodology presented here illustrates how the effective integration of pre-existing biological knowledge and high-throughput experimental data drives biological discovery and hypothesis generation.