Gene expression programs of human smooth muscle cells: tissue-specific differentiation and prognostic significance in breast cancers

Depot-specific mRNA expression programs in human adipocytes suggest physiological specialization via distinct developmental programs

Adipose tissue is distributed in diverse locations throughout the human body. Not much is known about the extent to which anatomically distinct adipose depots are functionally distinct, specialized organs, nor whether depot-specific characteristics result from intrinsic developmental programs, as opposed to reversible physiological responses to differences in tissue microenvironment. We used DNA microarrays to compare mRNA expression patterns of isolated human adipocytes and cultured adipose stem cells, before and after ex vivo adipocyte differentiation, from seven anatomically diverse adipose tissue depots. Adipocytes from different depots display distinct gene expression programs, which are most closely shared with anatomically related depots. mRNAs whose expression differs between anatomically diverse groups of depots (e.g., subcutaneous vs. internal) suggest important functional specializations. These depot-specific differences in gene expression were recapitulated when adipocyte progenitor cells from each site were differentiated ex vivo, suggesting that progenitor cells from specific anatomic sites are deterministically programmed to differentiate into depot-specific adipocytes. Many developmental transcription factors show striking depot-specific patterns of expression, suggesting that adipocytes in each anatomic depot are programmed during early development in concert with anatomically related tissues and organs. Our results support the hypothesis that adipocytes from different depots are functionally distinct and that their depot-specific specialization reflects distinct developmental programs.

Single Nucleus Sequencing of Human Colon Myenteric Plexus-Associated Visceral Smooth Muscle Cells, Platelet Derived Growth Factor Receptor Alpha Cells, and Interstitial Cells of Cajal

BACKGROUND AND AIMS

Smooth muscle cells (SMCs), interstitial cells of Cajal (ICCs), and platelet-derived growth factor receptor alpha (PDGFRα+) cells (PαCs) form a functional syncytium in the bowel known as the "SIP syncytium." The SIP syncytium works in concert with the enteric nervous system (ENS) to coordinate bowel motility. However, our understanding of individual cell types that form this syncytium and how they interact with each other remains limited, with no prior single-cell RNAseq analyses focused on human SIP syncytium cells.

METHODS

We analyzed single-nucleus RNA sequencing data from 10,749 human colon SIP syncytium cells (5572 SMC, 372 ICC, and 4805 PαC nuclei) derived from 15 individuals.

RESULTS

Consistent with critical contractile and pacemaker functions and with known enteric nervous system interactions, SIP syncytium cell types express many ion channels, including mechanosensitive channels in ICCs and PαCs. PαCs also prominently express extracellular matrix-associated genes and the inhibitory neurotransmitter receptor for vasoactive intestinal peptide (VIPR2), a novel finding. We identified 2 PαC clusters that differ in the expression of many ion channels and transcriptional regulators. Interestingly, SIP syncytium cells co-express 6 transcription factors (FOS, MEIS1, MEIS2, PBX1, SCMH1, and ZBTB16) that may be part of a combinatorial signature that specifies these cells. Bowel region-specific differences in SIP syncytium gene expression may correlate with regional differences in function, with right (ascending) colon SMCs and PαCs expressing more transcriptional regulators and ion channels than SMCs and PαCs in left (sigmoid) colon.

CONCLUSION

These studies provide new insights into SIP syncytium biology that may be valuable for understanding bowel motility disorders and lead to future investigation of highlighted genes and pathways.

The Molecular Signature of Human Testicular Peritubular Cells Revealed by Single-Cell Analysis

Peritubular cells of the human testis form a small compartment surrounding the seminiferous tubules. They are crucial for sperm transport, and they emerge as contributors to the spermatogonial stem cell niche. They are among the least known cell types of the human body. We employed single-cell RNA sequencing of cultured human testicular peritubular cells (HTPCs), which had been isolated from testicular samples of donors with normal spermatogenesis. The significant overlap between our results and recently published ex vivo data indicates that HTPCs are a highly adequate cellular model to define and study these cells. Thus, based on the expression of several markers, HTPCs can be classified as testicular smooth muscle cells. Small differences between the in vivo/in vitro expressed genes may be due to cellular plasticity. Plasticity was also shown upon addition of FCS to the culture medium. Based on transcriptome similarities, four cellular states were identified. Further analyses confirmed the presence of known stem cell niche-relevant factors (e.g., GDNF) and identified unknown functions, e.g., the ability to produce retinoic acid. Therefore, HTPCs allow us to define the signature(s) and delineate the functions of human testicular peritubular cells. The data may also serve as a resource for future studies to better understand male (in)fertility.

Identification of Hub Genes of Keloid Fibroblasts by Coexpression Network Analysis and Degree Algorithm

Background

Keloid is a benign dermal tumor characterized by abnormal proliferation and invasion of fibroblasts. The establishment of biomarkers is essential for the diagnosis and treatment of keloids.

Methods

We systematically identified coexpression modules using the weighted gene coexpression network analysis method (WGCNA). Differential expressed genes (DEGs) in GSE145725 and genes in significant modules were integrated to identify overlapping key genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were then performed, as well as protein-protein interaction (PPI) network construction for hub gene screening.

Results

Using the R package of WGCNA, 22 coexpression modules consisting of different genes were identified from the top 5,000 genes with maximum mean absolute deviation in 19 human fibroblast samples. Blue-green and yellow modules were identified as the most important modules, where genes overlapping with DEGs were identified as key genes. We identified the most critical functions and pathways as extracellular structure organization, vascular smooth muscle contraction, and the cGMP-PKG signaling pathway. Hub genes from key genes as BMP4, MSX1, HAND2, TBX2, SIX1, IRX1, EDN1, DLX5, MEF2C, and DLX2 were identified.

Conclusion

The blue-green and yellow modules may play an important role in the pathogenesis of keloid. 10 hub genes were identified as potential biomarkers and therapeutic targets for keloid.

Reviving Cav1.2 as an attractive drug target to treat bladder dysfunction

Inhibition of bladder contraction with antimuscarinics is a common approach to treat bladder hyperactivity, and the L-type voltage-gated calcium channel α_1C (Cav1.2) is crucial for bladder contractility. Therefore, strategies aimed at inhibiting Cav1.2 appear warranted. However, multiple clinical trials that attempted to treat bladder overactivity with calcium channel blockers (CCBs) have been unsuccessful, creating an unsolved mystery. In contrast, cardiologists and epidemiologists have reported strong associations between CCB use and bladder hyperactivity, opposing expectations of urologists. Recent findings from our lab offer a potential explanation. We have demonstrated that ketamine which can cause cystitis, functions, like nifedipine, as a Cav1.2 antagonist. We also show that a Cav1.2 agonist which potentiates muscle contraction, rather than antagonizing it, can increase the volume of voids and reduce voiding frequency. This perspective will discuss in detail the unsuccessful urological trials of CCBs and the promise of Cav1.2 agonists as potential novel therapies for bladder dysfunctions.

Mechanical programming of arterial smooth muscle cells in health and ageing

Arterial smooth muscle cells (ASMCs), the predominant cell type within the arterial wall, detect and respond to external mechanical forces. These forces can be derived from blood flow (i.e. pressure and stretch) or from the supporting extracellular matrix (i.e. stiffness and topography). The healthy arterial wall is elastic, allowing the artery to change shape in response to changes in blood pressure, a property known as arterial compliance. As we age, the mechanical forces applied to ASMCs change; blood pressure and arterial wall rigidity increase and result in a reduction in arterial compliance. These changes in mechanical environment enhance ASMC contractility and promote disease-associated changes in ASMC phenotype. For mechanical stimuli to programme ASMCs, forces must influence the cell's load-bearing apparatus, the cytoskeleton. Comprised of an interconnected network of actin filaments, microtubules and intermediate filaments, each cytoskeletal component has distinct mechanical properties that enable ASMCs to respond to changes within the mechanical environment whilst maintaining cell integrity. In this review, we discuss how mechanically driven cytoskeletal reorganisation programmes ASMC function and phenotypic switching.

Tissue-specific smooth muscle cell subtypes identified by transcriptional profiling

Smooth muscle cells (SMCs) are specialized cells present in many organs where they serve diverse tissue-specific functions. Using the Tabula Muris compendium of single-cell RNA sequencing data, we extracted individual SMC transcriptomes from eight mouse tissues to investigate the transcriptomic landscape of tissue-specific SMCs. We identified marker genes, signaling pathways, and biological processes enriched in tissue-specific SMCs, and inferred potential ligand-receptor interaction between SMC and other cell types. Our analysis also identified sex differences in SMC gene expression in different tissues. Lastly, we used unsupervised clustering to identify novel SMC subtypes based on their downstream targets of transcription factors. Our results highlight the variable SMC phenotypes and underscore this cell's remarkable adaptability to contribute to diverse tissue function.

Combination of PD98059 and TGF-β1 Efficiently Differentiates Human Urine-Derived Stem Cells into Smooth Muscle Cells

Pluripotent adult stem cells have potential applications in cell therapy and tissue engineering. Urine-derived stem cells (UDSCs) differentiate into various cell types. Here, we attempted to differentiate human UDSCs (hUDSCs) into smooth muscle cells (SMCs) using transforming growth factor-beta 1 (TGF-β1) and/or PD98059, an extracellular signal-regulated kinase (ERK) inhibitor. Both quantitative polymerase chain reaction (qPCR) and Western blot analysis showed that the expression of messenger ribonucleic acid (mRNA) and proteins for alpha-smooth muscle actin (α-SMA), calponin (CNN1), and smooth muscle myosin heavy chain (SM-MHC), which are specific markers for SMCs, increased on day 9 after differentiation and again on day 14. The differentiated cells from human UDSCs (hUDSCs) with a combination of TGF-β1 and PD98059 showed the highest expression of SMC marker proteins. Immunocytochemical staining performed to assess the molecular expression revealed CNN and α-SMA colocalizing in the cytoplasm. The cells that differentiated from hUDSCs with a combination of TGF-β1 and PD98059 showed the strongest expression for CNN1, α-SMA, and SM-MHC. Functional testing of the differentiated cells revealed a stronger contractile capacity for the cells differentiated with a combination of PD98059 and TGF-β1 than those differentiated with a single factor. These results suggest the combination of PD98059 and TGF-β1 to be a more effective differentiation method and that differentiated SMCs could be used for restoring the functions of the sphincter muscle or bladder.

Hox-Positive Adult Mesenchymal Stromal Cells: Beyond Positional Identity

Homeotic genes (Hox) are universal regulators of the body patterning process in embryogenesis of metazoans. The Hox gene expression pattern (Hox code) retains in adult tissues and serves as a cellular positional identity marker. Despite previously existing notions that the Hox code is inherent in all stroma mesenchymal cells as a whole, recent studies have shown that the Hox code may be an attribute of a distinct subpopulation of adult resident mesenchymal stromal cells (MSC). Recent evidence allows suggesting a "non-canonical" role for Hox gene expression which is associated with renewal and regeneration in postnatal organs after damage. In tissues with high regenerative capacity, it has been shown that a special cell population is critical for these processes, a distinctive feature of which is the persistent expression of tissue-specific Hox genes. We believe that in the postnatal period Hox-positive subpopulation of resident MSC may serve as a unique regenerative reserve. These cells coordinate creation and maintenance of the correct structure of the stroma through a tissue-specific combination of mechanisms. In this article, we summarize data on the role of resident MSC with a tissue-specific pattern of Hox gene expression as regulators of correct tissue reconstruction after injury.

Non-invasive evaluation of the equine gastrointestinal mucosal transcriptome

Evaluating the health and function of the gastrointestinal tract can be challenging in all species, but is especially difficult in horses due to their size and length of the gastrointestinal (GI) tract. Isolation of mRNA of cells exfoliated from the GI mucosa into feces (i.e., the exfoliome) offers a novel means of non-invasively examining the gene expression profile of the GI mucosa. This approach has been utilized in people with colorectal cancer. Moreover, we have utilized this approach in a murine model of GI inflammation and demonstrated that the exfoliome reflects the tissue transcriptome. The ability of the equine exfoliome to provide non-invasive information regarding the health and function of the GI tract is not known. The objective of this study was to characterize the gene expression profile found in exfoliated intestinal epithelial cells from normal horses and compare the exfoliome data with the tissue mucosal transcriptome. Mucosal samples were collected from standardized locations along the GI tract (i.e. ileum, cecum, right dorsal colon, and rectum) from four healthy horses immediately following euthanasia. Voided feces were also collected. RNA isolation, library preparation, and RNA sequencing was performed on fecal and intestinal mucosal samples. Comparison of gene expression profiles from the tissue and exfoliome revealed correlation of gene expression. Moreover, the exfoliome contained reads representing the diverse array of cell types found in the GI mucosa suggesting the equine exfoliome serves as a non-invasive means of examining the global gene expression pattern of the equine GI tract.

Single Cell Transcriptomes Derived from Human Cervical and Uterine Tissue during Pregnancy

This work presents the workflow for generating single cell transcriptomes derived from primary human uterine and cervical tissue obtained during planned cesarean hysterectomies. In total, a catalogue of 310 single cell transcriptomes are obtained, cell types present in these biopsies are inferred, and specific genes defining each of the cellular types present in the tissue are identified. Further validation of the inferred cell identity is also demonstrated via meta-analysis of independent repositories in literature generated by bulk sequenced data of fluorescence-activated cell sorting sorted cells.

Comparative analysis of myometrial and vascular smooth muscle cells to determine optimal cells for use in drug discovery

Novel therapeutic regulators of uterine contractility are needed to manage preterm labor, induce labor and control postpartum hemorrhage. Therefore, we previously developed a high-throughput assay for large-scale screening of small molecular compounds to regulate calcium-mobilization in primary mouse uterine myometrial cells. The goal of this study was to select the optimal myometrial cells for our high-throughput drug discovery assay, as well as determine the similarity or differences of myometrial cells to vascular smooth muscle cells (VSMCs)-the most common off-target of current myometrial therapeutics. Molecular and pharmacological assays were used to compare myometrial cells from four sources: primary cells isolated from term pregnant human and murine myometrium, immortalized pregnant human myometrial (PHM-1) cells and immortalized non-pregnant human myometrial (hTERT-HM) cells. In addition, myometrial cells were compared to vascular SMCs. We found that the transcriptome profiles of hTERT-HM and PHM1 cells were most similar (r = 0.93 and 0.90, respectively) to human primary myometrial cells. Comparative transcriptome profiling of primary human myometrial transcriptome and VSMCs revealed 498 upregulated (p ≤ 0.01, log2FC≥1) genes, of which 142 can serve as uterine-selective druggable targets. In the high-throughput Ca²⁺-assay, PHM1 cells had the most similar response to primary human myometrial cells in OT-induced Ca²⁺-release (E_max = 195% and 143%, EC₅₀ = 30 nM and 120 nM, respectively), while all sources of myometrial cells showed excellent and similar robustness and reproducibility (Z' = 0.52 to 0.77). After testing a panel of 61 compounds, we found that the stimulatory and inhibitory responses of hTERT-HM cells were highly-correlated (r = 0.94 and 0.95, respectively) to human primary cells. Moreover, ten compounds were identified that displayed uterine-selectivity (≥5-fold E_max or EC₅₀ compared to VSMCs). Collectively, this study found that hTERT-HM cells exhibited the most similarity to primary human myometrial cells and, therefore, is an optimal substitute for large-scale screening to identify novel therapeutic regulators of myometrial contractility. Moreover, VSMCs can serve as an important counter-screening tool to assess uterine-selectivity of targets and drugs given the similarity observed in the transcriptome and response to compounds.

Platelet-derived growth factor receptor α gene is regulated by multiple first exons

Transcription of the platelet-derived growth factor receptor α (PDGFRA/Pdgfra) gene is considered to be precisely regulated. We have previously reported that the PDGFRA/Pdgfra gene is regulated by a dual promoter system in human and mouse, in which a novel PDGFRA/Pdgfra transcript has a first exon (exon 1β) different from that of the canonical PDGFRA/Pdgfra transcript (exon 1α). To elucidate the function of each transcript, we first investigated the contribution of different PDGFRA transcripts to final protein levels. Notably, knockdown experiments suggested the existence of other PDGFRA transcripts, and we identified five additional first exons (exons 1γ, 1δ, 1ε, 1ζ, and 1η) in intron 1 in both the human and mouse genes. The first exons of the mouse Pdgfra gene showed unique expression patterns: exon 1α was broadly expressed; exon 1β was highly expressed in embryos; exon 1γ was observed at relatively high levels in the adult central nervous system (CNS); and exon 1δ was expressed at relatively high levels in the developing CNS. Furthermore, in silico analysis of common putative transcription factor binding sites in the upstream regions of the first exons of both human and mouse PDGFRA/Pdgfra genes predicted common (such as Sry, Mzf1, and Cdx) and unique (such as Sox5, Lmo2, and GATA) transcription factors. Our findings show the diversity of the transcriptional regulation of the PDGFRA/Pdgfra gene.

A network diffusion approach to inferring sample-specific function reveals functional changes associated with breast cancer

Guilt-by-association codifies the empirical observation that a gene's function is informed by its neighborhood in a biological network. This would imply that when a gene's network context is altered, for instance in disease condition, so could be the gene's function. Although context-specific changes in biological networks have been explored, the potential changes they may induce on the functional roles of genes are yet to be characterized. Here we analyze, for the first time, the network-induced potential functional changes in breast cancer. Using transcriptomic samples for 1047 breast tumors and 110 healthy breast tissues from TCGA, we derive sample-specific protein interaction networks and assign sample-specific functions to genes via a diffusion strategy. Testing for significant changes in the inferred functions between normal and cancer samples, we find several functions to have significantly gained or lost genes in cancer, not due to differential expression of genes known to perform the function, but rather due to changes in the network topology. Our predicted functional changes are supported by mutational and copy number profiles in breast cancers. Our diffusion-based functional assignment provides a novel characterization of a tumor that is complementary to the standard approach based on functional annotation alone. Importantly, this characterization is effective in predicting patient survival, as well as in predicting several known histopathological subtypes of breast cancer.

The non-invasive exfoliated transcriptome (exfoliome) reflects the tissue-level transcriptome in a mouse model of NSAID enteropathy

Non-steroidal anti-inflammatory drugs (NSAIDs) are among the most frequently used classes of medications in the world, yet they induce an enteropathy that is associated with high morbidity and mortality. A major limitation to better understanding the pathophysiology and diagnosis of this enteropathy is the difficulty of obtaining information about the primary site of injury, namely the distal small intestine. We investigated the utility of using mRNA from exfoliated cells in stool as a means to surveil the distal small intestine in a murine model of NSAID enteropathy. Specifically, we performed RNA-Seq on exfoliated cells found in feces and compared these data to RNA-Seq from both the small intestinal mucosa and colonic mucosa of healthy control mice or those exhibiting NSAID-induced enteropathy. Global gene expression analysis, data intersection, pathway analysis, and computational approaches including linear discriminant analysis (LDA) and sparse canonical correlation analysis (CCA) were used to assess the inter-relatedness of tissue (invasive) and stool (noninvasive) datasets. These analyses revealed that the exfoliated cell transcriptome closely mirrored the transcriptome of the small intestinal mucosa. Thus, the exfoliome may serve as a non-invasive means of detecting and monitoring NSAID enteropathy (and possibly other gastrointestinal mucosal inflammatory diseases).

Vascular aging: Molecular mechanisms and potential treatments for vascular rejuvenation

Aging is the main risk factor contributing to vascular dysfunction and the progression of vascular diseases. In this review, we discuss the causes and mechanisms of vascular aging at the tissue and cellular level. We focus on Endothelial Cell (EC) and Smooth Muscle Cell (SMC) aging due to their critical role in mediating the defective vascular phenotype. We elaborate on two categories that contribute to cellular dysfunction: cell extrinsic and intrinsic factors. Extrinsic factors reflect systemic or environmental changes which alter EC and SMC homeostasis compromising vascular function. Intrinsic factors induce EC and SMC transformation resulting in cellular senescence. Replenishing or rejuvenating the aged/dysfunctional vascular cells is critical to the effective repair of the vasculature. As such, this review also elaborates on recent findings which indicate that stem cell and gene therapies may restore the impaired vascular cell function, reverse vascular aging, and prolong lifespan.

Meta-analysis of 375,000 individuals identifies 38 susceptibility loci for migraine

Migraine is a debilitating neurological disorder affecting around one in seven people worldwide, but its molecular mechanisms remain poorly understood. There is some debate about whether migraine is a disease of vascular dysfunction or a result of neuronal dysfunction with secondary vascular changes. Genome-wide association (GWA) studies have thus far identified 13 independent loci associated with migraine. To identify new susceptibility loci, we carried out a genetic study of migraine on 59,674 affected subjects and 316,078 controls from 22 GWA studies. We identified 44 independent single-nucleotide polymorphisms (SNPs) significantly associated with migraine risk (P < 5 × 10(-8)) that mapped to 38 distinct genomic loci, including 28 loci not previously reported and a locus that to our knowledge is the first to be identified on chromosome X. In subsequent computational analyses, the identified loci showed enrichment for genes expressed in vascular and smooth muscle tissues, consistent with a predominant theory of migraine that highlights vascular etiologies.

Transcriptome of the inner circular smooth muscle of the developing mouse intestine: Evidence for regulation of visceral smooth muscle genes by the hedgehog target gene, cJun

BACKGROUND

Digestion is facilitated by coordinated contractions of the intestinal muscularis externa, a bilayered smooth muscle structure that is composed of inner circular muscles (ICM) and outer longitudinal muscles (OLM). We performed transcriptome analysis of intestinal mesenchyme tissue at E14.5, when the ICM, but not the OLM, is present, to investigate the transcriptional program of the ICM.

RESULTS

We identified 3967 genes enriched in E14.5 intestinal mesenchyme. The gene expression profiles were clustered and annotated to known muscle genes, identifying a muscle-enriched subcluster. Using publically available in situ data, 127 genes were verified as expressed in ICM. Examination of the promoter and regulatory regions for these co-expressed genes revealed enrichment for cJUN transcription factor binding sites, and cJUN protein was enriched in ICM. cJUN ChIP-seq, performed at E14.5, revealed that cJUN regulatory regions contain characteristics of muscle enhancers. Finally, we show that cJun is a target of Hedgehog (Hh), a signaling pathway known to be important in smooth muscle development, and identify a cJun genomic enhancer that is responsive to Hh.

CONCLUSIONS

This work provides the first transcriptional catalog for the developing ICM and suggests that cJun regulates gene expression in the ICM downstream of Hh signaling. Developmental Dynamics 245:614-626, 2016. © 2016 Wiley Periodicals, Inc.

Drug repositioning for non-small cell lung cancer by using machine learning algorithms and topological graph theory

Background

Non-small cell lung cancer (NSCLC) is one of the leading causes of death globally, and research into NSCLC has been accumulating steadily over several years. Drug repositioning is the current trend in the pharmaceutical industry for identifying potential new uses for existing drugs and accelerating the development process of drugs, as well as reducing side effects.

Results

This work integrates two approaches - machine learning algorithms and topological parameter-based classification - to develop a novel pipeline of drug repositioning to analyze four lung cancer microarray datasets, enriched biological processes, potential therapeutic drugs and targeted genes for NSCLC treatments. A total of 7 (8) and 11 (12) promising drugs (targeted genes) were discovered for treating early- and late-stage NSCLC, respectively. The effectiveness of these drugs is supported by the literature, experimentally determined in-vitro IC₅₀ and clinical trials. This work provides better drug prediction accuracy than competitive research according to IC₅₀ measurements.

Conclusions

With the novel pipeline of drug repositioning, the discovery of enriched pathways and potential drugs related to NSCLC can provide insight into the key regulators of tumorigenesis and the treatment of NSCLC. Based on the verified effectiveness of the targeted drugs predicted by this pipeline, we suggest that our drug-finding pipeline is effective for repositioning drugs.

Role of Hox genes in stem cell differentiation

Hox genes are an evolutionary highly conserved gene family. They determine the anterior-posterior body axis in bilateral organisms and influence the developmental fate of cells. Embryonic stem cells are usually devoid of any Hox gene expression, but these transcription factors are activated in varying spatial and temporal patterns defining the development of various body regions. In the adult body, Hox genes are among others responsible for driving the differentiation of tissue stem cells towards their respective lineages in order to repair and maintain the correct function of tissues and organs. Due to their involvement in the embryonic and adult body, they have been suggested to be useable for improving stem cell differentiations in vitro and in vivo. In many studies Hox genes have been found as driving factors in stem cell differentiation towards adipogenesis, in lineages involved in bone and joint formation, mainly chondrogenesis and osteogenesis, in cardiovascular lineages including endothelial and smooth muscle cell differentiations, and in neurogenesis. As life expectancy is rising, the demand for tissue reconstruction continues to increase. Stem cells have become an increasingly popular choice for creating therapies in regenerative medicine due to their self-renewal and differentiation potential. Especially mesenchymal stem cells are used more and more frequently due to their easy handling and accessibility, combined with a low tumorgenicity and little ethical concerns. This review therefore intends to summarize to date known correlations between natural Hox gene expression patterns in body tissues and during the differentiation of various stem cells towards their respective lineages with a major focus on mesenchymal stem cell differentiations. This overview shall help to understand the complex interactions of Hox genes and differentiation processes all over the body as well as in vitro for further improvement of stem cell treatments in future regenerative medicine approaches.

Notch3 is necessary for blood vessel integrity in the central nervous system

OBJECTIVE

Vascular smooth muscle cells (VSMC) are important for contraction, blood flow distribution, and regulation of blood vessel diameter, but to what extent they contribute to the integrity of blood vessels and blood-brain barrier function is less well understood. In this report, we explored the impact of the loss of VSMC in the Notch3(-/-) mouse on blood vessel integrity in the central nervous system.

APPROACH AND RESULTS

Notch3(-/-) mice showed focal disruptions of the blood-brain barrier demonstrated by extravasation of tracers accompanied by fibrin deposition in the retinal vasculature. This blood-brain barrier leakage was accompanied by a regionalized and patchy loss of VSMC, with VSMC gaps predominantly in arterial resistance vessels of larger caliber. The loss of VSMC appeared to be caused by progressive degeneration of VSMC resulting in a gradual loss of VSMC marker expression and a progressive acquisition of an aberrant VSMC phenotype closer to the gaps, followed by enhanced apoptosis and cellular disintegration in the gaps. Arterial VSMC were the only mural cell type that was morphologically affected, despite Notch3 also being expressed in pericytes. Transcriptome analysis of isolated brain microvessels revealed gene expression changes in Notch3(-/-) mice consistent with loss of arterial VSMC and presumably secondary transcriptional changes were observed in endothelial genes, which may explain the compromised vascular integrity.

CONCLUSIONS

We demonstrate that Notch3 is important for survival of VSMC, and reveal a critical role for Notch3 and VSMC in blood vessel integrity and blood-brain barrier function in the mammalian vasculature.

Modeling cerebrovascular pathophysiology in amyloid-β metabolism using neural-crest-derived smooth muscle cells

There is growing recognition of cerebrovascular contributions to neurodegenerative diseases. In the walls of cerebral arteries, amyloid-beta (Aβ) accumulation is evident in a majority of aged people and patients with cerebral amyloid angiopathy. Here, we leverage human pluripotent stem cells to generate vascular smooth muscle cells (SMCs) from neural crest progenitors, recapitulating brain-vasculature-specific attributes of Aβ metabolism. We confirm that the lipoprotein receptor, LRP1, functions in our neural-crest-derived SMCs to mediate Aβ uptake and intracellular lysosomal degradation. Hypoxia significantly compromises the contribution of SMCs to Aβ clearance by suppressing LRP1 expression. This enabled us to develop an assay of Aβ uptake by using the neural crest-derived SMCs with hypoxia as a stress paradigm. We then tested several vascular protective compounds in a high-throughput format, demonstrating the value of stem-cell-based phenotypic screening for novel therapeutics and drug repurposing, aimed at alleviating amyloid burden.

Sulfur dioxide inhibits vascular smooth muscle cell proliferation via suppressing the Erk/MAP kinase pathway mediated by cAMP/PKA signaling

The present study was designed to investigate the role of endogenous sulfur dioxide (SO2) in vascular smooth muscle cell (VSMC) proliferation, and explore the possible role of cross-talk between cyclic adenosine monophosphate (cAMP)/protein kinase A (PKA) and extracellular signal-regulated kinase (Erk)/mitogen-activated protein kinase (MAPK) pathways in this action. By cell counting, growth curve depict, flow cytometry and bromodeoxyuridine (BrdU) labeling assays, we found that SO2 inhibited VSMC proliferation by preventing cell cycle progression from G1 to S phase and by reducing DNA synthesis. SO2 synthase aspartate aminotransferase (AAT1 and AAT2) overexpression significantly inhibited serum-induced proliferating cell nuclear antigen (PCNA) protein expression in VSMCs, demonstrated by western blot analysis. Moreover, overexpression of AAT1 or AAT2 markedly reduced incorporation of BrdU in serum-treated VSMCs. By contrast, either AAT1 or AAT2 knockdown significantly exacerbated serum-stimulated VSMC proliferation. Thus, both exogenous- and endogenous-derived SO2 suppressed serum-induced VSMC proliferation. However, annexin V-propidium iodide (PI) staining and cell cycle analysis demonstrated that SO2 did not influence VSMC apoptosis in the serum-induced proliferation model. In a platelet-derived growth factor (PDGF)-BB-stimulated VSMC proliferation model, SO2 dephosphorylated the active sites of Erk1/2, MAPK kinase 1/2 and RAF proto-oncogene serine/threonine-protein kinase (c-Raf) induced by PDGF-BB. However, the inactivation of the three kinases of the Erk/MAPK pathway was not due to the separate interferences on them by SO2 simultaneously, but a consequence of the influence on the upstream activity of the c-Raf molecule. Hence, we examined the cAMP/PKA pathway, which could inhibit Erk/MAPK transduction in VSMCs. The results showed that SO2 could stimulate the cAMP/PKA pathway to block c-Raf activation, whereas the Ser259 site on c-Raf had an important role in SO2-induced suppression of Erk/MAPK pathway. The present study firstly demonstrated that SO2 exerted a negative regulation of VSMC proliferation via suppressing the Erk/MAPK pathway mediated by cAMP/PKA signaling.

Glimpse into Hox and tale regulation of cell differentiation and reprogramming

During embryonic development, cells become gradually restricted in their developmental potential and start elaborating lineage-specific transcriptional networks to ultimately acquire a unique differentiated state. Hox genes play a central role in specifying regional identities, thereby providing the cell with critical information on positional value along its differentiation path. The exquisite DNA-binding specificity of the Hox proteins is frequently dependent upon their interaction with members of the TALE family of homeodomain proteins. In addition to their function as Hox-cofactors, TALE homeoproteins control multiple crucial developmental processes through Hox-independent mechanisms. Here, we will review recent findings on the function of both Hox and TALE proteins in cell differentiation, referring mostly to vertebrate species. In addition, we will discuss the direct implications of this knowledge on cell plasticity and cell reprogramming.

β1-integrin is essential for vasoregulation and smooth muscle survival in vivo

OBJECTIVE

Integrins contribute to vascular morphogenesis through regulation of adhesion and assembly of the extracellular matrix. However, the role of β1-integrin in the mature vascular wall is less clear.

APPROACH AND RESULTS

We sought to determine the function of β1-integrin in mature smooth muscle cells in vivo using a loss of function approach by crossing a tamoxifen-inducible sm22αCre line to a floxed β1-integrin transgenic line. Adult mice lacking smooth muscle β1-integrin survived only 10 weeks post induction. The deletion of β1-integrin resulted in profound loss of vasomotor control. Histological analysis revealed progressive fibrosis in arteries with associated apoptosis of smooth muscle cells, which was not rescued by adventitial stem cells. Smooth muscle cell apoptosis was detected in arteries with dead cells replaced primarily by collagen. Despite the catastrophic effects on vascular smooth muscle, the deleted visceral smooth muscle remained viable with the exception of a short portion of the colon, indicating that vascular but not visceral smooth muscle is particularly sensitive to changes in β1-integrin.

CONCLUSIONS

This study reveals an essential function of β1-integrin in the maintenance of vasomotor control and highlights a critical role for β1-integrin in vascular, but not visceral, smooth muscle survival.

Embryological-origin-dependent differences in homeobox expression in adult aorta: role in regional phenotypic variability and regulation of NF-κB activity

OBJECTIVE

Different vascular beds show differing susceptibility to the development of atherosclerosis, but the molecular mechanisms underlying these differences are incompletely understood. This study aims to identify factors that contribute to the phenotypic heterogeneity of distinct regions of the adult vasculature.

APPROACH AND RESULTS

High-throughput mRNA profiling in adult mice reveals higher expression of the homeobox paralogous genes 6 to 10 (Hox6-10) in the athero-resistant thoracic aorta (TA) than in the athero-susceptible aortic arch (AA). Higher homeobox gene expression also occurs in rat and porcine TA, and is maintained in primary smooth muscle cells isolated from TA (TA-SMCs) compared with cells from AA (AA-SMCs). This region-specific homeobox gene expression pattern is also observed in human embryonic stem cells differentiated into neuroectoderm-SMCs and paraxial mesoderm-SMCs, which give rise to AA-SMCs and TA-SMCs, respectively. We also find that, compared with AA and AA-SMCs, TA and TA-SMCs have lower activity of the proinflammatory and proatherogenic nuclear factor-κB (NF-κB) and lower expression of NF-κB target genes, at least in part attributable to HOXA9-dependent inhibition. Conversely, NF-κB inhibits HOXA9 promoter activity and mRNA expression in SMCs.

CONCLUSION

Our findings support a model of Hox6-10-specified positional identity in the adult vasculature that is established by embryonic cues independently of environmental factors and is conserved in different mammalian species. Differential homeobox gene expression contributes to maintaining phenotypic differences between SMCs from athero-resistant and athero-susceptible regions, at least in part through feedback regulatory mechanisms involving inflammatory mediators, for example, reciprocal inhibition between HOXA9 and NF-κB.

Glutamine synthetase is a genetic determinant of cell type-specific glutamine independence in breast epithelia

Although significant variations in the metabolic profiles exist among different cells, little is understood in terms of genetic regulations of such cell type-specific metabolic phenotypes and nutrient requirements. While many cancer cells depend on exogenous glutamine for survival to justify the therapeutic targeting of glutamine metabolism, the mechanisms of glutamine dependence and likely response and resistance of such glutamine-targeting strategies among cancers are largely unknown. In this study, we have found a systematic variation in the glutamine dependence among breast tumor subtypes associated with mammary differentiation: basal- but not luminal-type breast cells are more glutamine-dependent and may be susceptible to glutamine-targeting therapeutics. Glutamine independence of luminal-type cells is associated mechanistically with lineage-specific expression of glutamine synthetase (GS). Luminal cells can also rescue basal cells in co-culture without glutamine, indicating a potential for glutamine symbiosis within breast ducts. The luminal-specific expression of GS is directly induced by GATA3 and represses glutaminase expression. Such distinct glutamine dependency and metabolic symbiosis is coupled with the acquisition of the GS and glutamine independence during the mammary differentiation program. Understanding the genetic circuitry governing distinct metabolic patterns is relevant to many symbiotic relationships among different cells and organisms. In addition, the ability of GS to predict patterns of glutamine metabolism and dependency among tumors is also crucial in the rational design and application of glutamine and other metabolic pathway targeted therapies.

Analysis of tumor environmental response and oncogenic pathway activation identifies distinct basal and luminal features in HER2-related breast tumor subtypes

Introduction

Breast cancer heterogeneity occurs as a consequence of the dysregulation of numerous oncogenic pathways as well as many non-genetic factors, including tumor microenvironmental stresses such as hypoxia, lactic acidosis, and glucose deprivation. Although the importance of these non-genetic factors is well recognized, it is not clear how to integrate these factors within the genetic framework of cancer as the next logical step in understanding tumor heterogeneity.

Methods

We report here the development of a series of gene expression signatures to measure the influences of microenvironmental stresses. The pathway activities of hypoxia, lactic acidosis, acidosis and glucose deprivation were investigated in a collection of 1,143 breast tumors, which have been separated into 17 breast tumor subgroups defined by their distinct patterns of oncogenic pathways. A validation dataset comprised of 547 breast tumors was also used to confirm the major findings, and representative breast cancer cell lines were utilized to validate in silico results and mechanistic studies.

Results

Through the integrative pathway analysis of microenvironmental stresses and oncogenic events in breast tumors, we identified many known and novel correlations between these two sources of tumor heterogeneity. Focusing on differences between two human epidermal growth factor receptor 2 (HER2)-related subgroups, previously identified based on patterns of oncogenic pathway activity, we determined that these subgroups differ with regards to tumor microenvironmental signatures, including hypoxia. We further demonstrate that each of these subgroups have features consistent with basal and luminal breast tumors including patterns of oncogenic signaling pathways, expression of subtype specific genes, and cellular mechanisms that regulate the hypoxia response. Importantly, we also demonstrate that the correlated pattern of hypoxia-related gene expression and basal-associated gene expression are consistent across HER2-related tumors whether we analyze the tumors as a function of our pathway-based classification scheme, using the intrinsic gene list (ERBB2+), or based on HER2 IHC status. Our results demonstrate a cell lineage-specific phenomenon in which basal-like tumors, HER2-related tumors with high hypoxia, as well as normal basal epithelial cells express increased mRNA levels of HIF-1α compared to luminal types and silencing of HIF-1α results in decreased expression of hypoxia-induced genes.

Conclusions

This study demonstrates differences in microenvironmental conditions in HER2-related subgroups defined by distinct oncogenic pathway activities, and provides a mechanistic explanation for differences in the observed hypoxia response between these subgroups. Collectively, these data demonstrate the potential of a pathway-based classification strategy as a framework to integrate genetic and non-genetic factors to investigate the basis of tumor heterogeneity.

Regional expression of HOXA4 along the aorta and its potential role in human abdominal aortic aneurysms

BACKGROUND

The infrarenal abdominal aorta exhibits increased disease susceptibility relative to other aortic regions. Allograft studies exchanging thoracic and abdominal segments showed that regional susceptibility is maintained regardless of location, suggesting substantial roles for embryological origin, tissue composition and site-specific gene expression.

RESULTS

We analyzed gene expression with microarrays in baboon aortas, and found that members of the HOX gene family exhibited spatial expression differences. HOXA4 was chosen for further study, since it had decreased expression in the abdominal compared to the thoracic aorta. Western blot analysis from 24 human aortas demonstrated significantly higher HOXA4 protein levels in thoracic compared to abdominal tissues (P < 0.001). Immunohistochemical staining for HOXA4 showed nuclear and perinuclear staining in endothelial and smooth muscle cells in aorta. The HOXA4 transcript levels were significantly decreased in human abdominal aortic aneurysms (AAAs) compared to age-matched non-aneurysmal controls (P < 0.00004). Cultured human aortic endothelial and smooth muscle cells stimulated with INF-γ (an important inflammatory cytokine in AAA pathogenesis) showed decreased levels of HOXA4 protein (P < 0.0007).

CONCLUSIONS

Our results demonstrated spatial variation in expression of HOXA4 in human aortas that persisted into adulthood and that downregulation of HOXA4 expression was associated with AAAs, an important aortic disease of the ageing population.

Comparison of genomics and functional imaging from canine sarcomas treated with thermoradiotherapy predicts therapeutic response and identifies combination therapeutics

PURPOSE

While hyperthermia is an effective adjuvant treatment to radiotherapy, we do not completely understand the nature of the response heterogeneity.

EXPERIMENTAL DESIGN

We performed gene expression analysis of 22 spontaneous canine sarcomas before and after the first hyperthermia treatment administered as an adjuvant to radiotherapy. In parallel, diffusion-weighted MRI (DWI) was done prior to the treatment course and at the end of therapy.

RESULTS

From the integrative analysis of gene expression and DWI, we identified significant correlation between tumor responses with genes involved in VEGF signaling, telomerase, DNA repair, and inflammation. The treatment-induced changes in gene expression identified 2 distinct tumor subtypes with significant differences in their gene expression and treatment response, as defined by changes in DWI. The 2 tumor subtypes could also be readily identified by pretreatment gene expression. The tumor subtypes, with stronger expression response and DWI increase, had higher levels of HSP70, POT1, and centrosomal proteins, and lower levels of CD31, vWF, and transferrin. Such differential gene expression between the 2 subtypes was used to interrogate connectivity map and identify linkages to an HSP90 inhibitor, geldanamycin. We further validated the ability of geldanamycin to enhance cell killing of human tumor cells with hyperthermia and radiotherapy in clonogenic assays.

CONCLUSIONS

To our knowledge, this is one of the first successful attempts to link changes in gene expression and functional imaging to understand the response heterogeneity and identify compounds enhancing thermoradiotherapy. This study also demonstrates the value of canine tumors to provide information generalizable to human tumors.

Identification and optimization of classifier genes from multi-class earthworm microarray dataset

Monitoring, assessment and prediction of environmental risks that chemicals pose demand rapid and accurate diagnostic assays. A variety of toxicological effects have been associated with explosive compounds TNT and RDX. One important goal of microarray experiments is to discover novel biomarkers for toxicity evaluation. We have developed an earthworm microarray containing 15,208 unique oligo probes and have used it to profile gene expression in 248 earthworms exposed to TNT, RDX or neither. We assembled a new machine learning pipeline consisting of several well-established feature filtering/selection and classification techniques to analyze the 248-array dataset in order to construct classifier models that can separate earthworm samples into three groups: control, TNT-treated, and RDX-treated. First, a total of 869 genes differentially expressed in response to TNT or RDX exposure were identified using a univariate statistical algorithm of class comparison. Then, decision tree-based algorithms were applied to select a subset of 354 classifier genes, which were ranked by their overall weight of significance. A multiclass support vector machine (MC-SVM) method and an unsupervised K-mean clustering method were applied to independently refine the classifier, producing a smaller subset of 39 and 30 classifier genes, separately, with 11 common genes being potential biomarkers. The combined 58 genes were considered the refined subset and used to build MC-SVM and clustering models with classification accuracy of 83.5% and 56.9%, respectively. This study demonstrates that the machine learning approach can be used to identify and optimize a small subset of classifier/biomarker genes from high dimensional datasets and generate classification models of acceptable precision for multiple classes.

Vascular bed-specific regulation of the von Willebrand factor promoter in the heart and skeletal muscle

A region of the human von Willebrand factor (VWF) gene between -2812 and the end of the first intron (termed vWF2) was previously shown to direct expression in the endothelium of capillaries and a subset of larger blood vessels in the heart and skeletal muscle. Here, our goal was to delineate the DNA sequences responsible for this effect. A series of constructs containing deletions or mutations of vWF2 coupled to LacZ were targeted to the Hprt locus of mice, and the resulting animals were analyzed for reporter gene expression. The findings demonstrate that DNA sequences between -843 and -620 are necessary for expression in capillary but not large vessel endothelium in heart and skeletal muscle. Further, expression of VWF in capillaries and larger vessels of both tissues required the presence of a native or heterologous intron. In vitro assays implicated a role for ERG-binding ETS motif at -56 in mediating basal expression of VWF. In Hprt-targeted mice, mutation of the ETS consensus motif resulted in loss of LacZ expression in the endothelium of the heart and skeletal muscle. Together, these data indicate that distinct DNA modules regulate vascular bed-specific expression of VWF.

Lactic acidosis triggers starvation response with paradoxical induction of TXNIP through MondoA

Although lactic acidosis is a prominent feature of solid tumors, we still have limited understanding of the mechanisms by which lactic acidosis influences metabolic phenotypes of cancer cells. We compared global transcriptional responses of breast cancer cells in response to three distinct tumor microenvironmental stresses: lactic acidosis, glucose deprivation, and hypoxia. We found that lactic acidosis and glucose deprivation trigger highly similar transcriptional responses, each inducing features of starvation response. In contrast to their comparable effects on gene expression, lactic acidosis and glucose deprivation have opposing effects on glucose uptake. This divergence of metabolic responses in the context of highly similar transcriptional responses allows the identification of a small subset of genes that are regulated in opposite directions by these two conditions. Among these selected genes, TXNIP and its paralogue ARRDC4 are both induced under lactic acidosis and repressed with glucose deprivation. This induction of TXNIP under lactic acidosis is caused by the activation of the glucose-sensing helix-loop-helix transcriptional complex MondoA:Mlx, which is usually triggered upon glucose exposure. Therefore, the upregulation of TXNIP significantly contributes to inhibition of tumor glycolytic phenotypes under lactic acidosis. Expression levels of TXNIP and ARRDC4 in human cancers are also highly correlated with predicted lactic acidosis pathway activities and associated with favorable clinical outcomes. Lactic acidosis triggers features of starvation response while activating the glucose-sensing MondoA-TXNIP pathways and contributing to the “anti-Warburg” metabolic effects and anti-tumor properties of cancer cells. These results stem from integrative analysis of transcriptome and metabolic response data under various tumor microenvironmental stresses and open new paths to explore how these stresses influence phenotypic and metabolic adaptations in human cancers.

Solid tumors usually have many differences in their chemical environments, such as low oxygen, depletion of glucose, high acidity (low pH), and accumulation of lactate, from normal tissues. These changes are usually called tumor microenvironmental stresses. In this study, we have used microarrays to compare the transcriptional response and metabolic adaptation in response to these different stresses seen in the tumor microenvironments. Through these comparisons, we have found that lactic acidosis triggers a starvation response, highly similar to glucose deprivation, even in the presence of abundant nutrients and oxygen. Even the cells seem to be starved; cells under lactic acidosis have decreased glucose uptake. We found this unexpected biological behavior was due to the paradoxical induction of a glucose-sensing Mondo-TXNIP pathway. The activation of this novel anti-tumor pathway under lactic acidosis contributes to the anti-Warburg effect and the restriction of cell growth in tumorigenesis by limiting nutrient availability and its inactivation may be required for tumor progression under these microenvironmental stresses.

Latent factor analysis to discover pathway-associated putative segmental aneuploidies in human cancers

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of “trans”-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1α protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.

Gene signatures are a powerful tool to investigate biological processes in human cancer. However, it is clear that these gene signatures do not fully reflect the complexity of human cancer. Here we demonstrate how a latent factor model can improve the in vivo relevance of these pathway-associated gene signatures by dissecting them into co-regulated transcriptional components which better represent the structure in human cancer. We use this approach to analyze hypoxia and lactic acidosis gene signatures to identify latent factors that represent distinct, interacting components of the various biological processes which are in the initial gene signatures but poorly dissected. Some factors are clustered in small chromosomal regions and their expression values are highly correlated with their DNA copy number in both cancer cell lines and human tumors. Therefore, the gene dosage at the DNA levels may explain the differences in gene expression. Several factors contain genes which are known to directly modulate the hypoxia response and allow us to generate testable hypotheses regarding particular copy number changes and hypoxia signatures. Therefore, the use of latent factor analysis is a powerful means to identify pathway-associated changes in the DNA copy number and gene dosage.

Gene-expression profiles, tumor microenvironment, and cancer stem cells in breast cancer: latest advances towards an integrated approach

During the past decade, several high throughput analytical methods for gene-expression profiling have been developed. DNA microarrays and multiplex RT-PCR have been applied in the field of breast-cancer research to establish new molecular taxonomic classifications, or a selected group of genes able to predict the prognosis of the patients and/or their response to chemotherapy. This technology provides an opportunity to refine the anti-neoplastic treatment and avoid the currently observed under- and over-treatment of breast-cancer patients. In parallel, high throughput technologies for gene-expression analysis have been applied to research on cancer stem cells (CSCs) and the tumor microenvironment, offering a wider vision of the molecular processes that influence carcinogenic events, disease development, and the response to the treatment of breast-cancer patients. In this report, we briefly revisit the most relevant genomic studies on breast-cancer prognosis and prediction to introduce the latest advances in tumor dormancy, its implications in the clinical outcome of disease-free patients and its connection with CSCs biology and microenvironment influence in the metastatic process. Finally, we have discussed the contribution of the results of these studies to the design of new experimental strategies oriented towards personalized medicine.

Anatomic Demarcation of Cells: Genes to Patterns

An organizing principle of the diverse cell types in multicellular organisms is their anatomic location. In turn, anatomic location is patterned by the positional identities of cells along developmental axes. Recent progress in functional genomics and chromatin biology illustrates how cells use specific gene expression programs to encode location. Dynamic chromatin states of key genes, notably the Hox loci, serve as the internal representation in cells of their positional identity within the animal.

Integrative genomics identifies DSCR1 (RCAN1) as a novel NFAT-dependent mediator of phenotypic modulation in vascular smooth muscle cells

Vascular smooth muscle cells (SMCs) display remarkable phenotypic plasticity in response to environmental cues. The nuclear factor of activated T-cells (NFAT) family of transcription factors plays a critical role in vascular pathology. However, known functional NFAT gene targets in vascular SMCs are currently limited. Publicly available whole-genome expression array data sets were analyzed to identify differentially expressed genes in human, mouse and rat SMCs. Comparison between vehicle and phenotypic modulatory stimuli identified 63 species-conserved, upregulated genes. Integration of the 63 upregulated genes with an in silico NFAT-ome (a species-conserved list of gene promoters containing at least one NFAT binding site) identified 18 putative NFAT-dependent genes. Further intersection of these 18 potential NFAT target genes with a mouse in vivo vascular injury microarray identified four putative NFAT-dependent, injury-responsive genes. In vitro validations substantiated the NFAT-dependent role of Cyclooxygenase 2 (COX2/PTGS2) in SMC phenotypic modulation and uncovered Down Syndrome Candidate Region 1 (DSCR1/RCAN1) as a novel NFAT target gene in SMCs. We show that induction of DSCR1 inhibits calcineurin/NFAT signaling through a negative feedback mechanism; DSCR1 overexpression attenuates NFAT transcriptional activity and COX2 protein expression, whereas knockdown of endogenous DSCR1 enhances NFAT transcriptional activity. Our integrative genomics approach illustrates how the combination of publicly available gene expression arrays, computational databases and empirical research methods can answer specific questions in any cell type for a transcriptional network of interest. Herein, we report DSCR1 as a novel NFAT-dependent, injury-inducible, early gene that may serve to negatively regulate SMC phenotypic switching.

Forkhead box M1 transcriptional factor is required for smooth muscle cells during embryonic development of blood vessels and esophagus

The forkhead box m1 (Foxm1 or Foxm1b) transcription factor (previously called HFH-11B, Trident, Win, or MPP2) is expressed in a variety of tissues during embryogenesis, including vascular, airway, and intestinal smooth muscle cells (SMCs). Although global deletion of Foxm1 in Foxm1(-/-) mice is lethal in the embryonic period due to multiple abnormalities in the liver, heart, and lung, the specific role of Foxm1 in SMC remains unknown. In the present study, Foxm1 was deleted conditionally in the developing SMC (smFoxm1(-/-) mice). The majority of smFoxm1(-/-) mice died immediately after birth due to severe pulmonary hemorrhage and structural defects in arterial wall and esophagus. Although Foxm1 deletion did not influence SMC differentiation, decreased proliferation of SMC was found in smFoxm1(-/-) blood vessels and esophagus. Depletion of Foxm1 in cultured SMC caused G(2) arrest and decreased numbers of cells undergoing mitosis. Foxm1-deficiency in vitro and in vivo was associated with reduced expression of cell cycle regulatory genes, including cyclin B1, Cdk1-activator Cdc25b phosphatase, Polo-like 1 and JNK1 kinases, and cMyc transcription factor. Foxm1 is critical for proliferation of smooth muscle cells and is required for proper embryonic development of blood vessels and esophagus.

Myocardin-dependent activation of the CArG box-rich smooth muscle gamma-actin gene: preferential utilization of a single CArG element through functional association with the NKX3.1 homeodomain protein

Serum response factor (SRF) is a ubiquitously expressed transcription factor that binds a 10-bp element known as the CArG box, located in the proximal regulatory region of hundreds of target genes. SRF activates target genes in a cell- and context-dependent manner by assembling unique combinations of cofactors over CArG elements. One particularly strong SRF cofactor, myocardin (MYOCD), acts as a component of a molecular switch for smooth muscle cell (SMC) differentiation by activating cytoskeletal and contractile genes harboring SRF-binding CArG elements. Here we report that the human ACTG2 promoter, containing four conserved CArG elements, displays SMC-specific basal activity and is highly induced in the presence of MYOCD. Stable transfection of a non-SMC cell type with Myocd elicits elevations in endogenous Actg2 mRNA. Gel shift and luciferase assays reveal a strong bias for MYOCD-dependent transactivation through CArG2 of the human ACTG2 promoter. Substitution of CArG2 with other CArGs, including a consensus CArG element, fails to reconstitute full MYOCD-dependent ACTG2 promoter stimulation. Mutation of an adjacent binding site for NKX3.1 reduces MYOCD-dependent transactivation of the ACTG2 promoter. Co-immunoprecipitation, glutathione S-transferase pulldown, and luciferase assays show a physical and functional association between MYOCD and NKX3.1; no such functional relationship is evident with the related NKX2.5 transcription factor despite its interaction with MYOCD. These results demonstrate the ability of MYOCD to discriminate among several juxtaposed CArG elements, presumably through its novel partnership with NKX3.1, to optimally transactivate the human ACTG2 promoter.

Transcriptional profiling of human cavernosal endothelial cells reveals distinctive cell adhesion phenotype and role for claudin 11 in vascular barrier function

To determine specific molecular features of endothelial cells (ECs) relevant to the physiological process of penile erection we compared gene expression of human EC derived from corpus cavernosum of men with and without erectile dysfunction (HCCEC) to coronary artery (HCAEC) and umbilical vein (HUVEC) using Affymetrix GeneChip microarrays and GeneSifter software. Genes differentially expressed across samples were partitioned around medoids to identify genes with highest expression in HCCEC. A total of 190 genes/transcripts were highly expressed only in HCCEC. Gene Ontology classification indicated cavernosal enrichment in genes related to cell adhesion, extracellular matrix, pattern specification and organogenesis. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed high expression of genes relating to ECM-receptor interaction, focal adhesions, and cytokine-cytokine receptor interaction. Real-time PCR confirmed expression differences in cadherins 2 and 11, claudin 11 (CLDN11), desmoplakin, and versican. CLDN11, a component of tight junctions not previously described in ECs, was highly expressed only in HCCEC and its knockdown by siRNA significantly reduced transendothelial electrical resistance in HCCEC. Overall, cavernosal ECs exhibited a transcriptional profile encoding matrix and adhesion proteins that regulate structural and functional characteristics of blood vessels. Contribution of the tight junction protein CLDN11 to barrier function in endothelial cells is novel and may reflect hemodynamic requirements of the corpus cavernosum.

Regeneration, repair and remembering identity: the three Rs of Hox gene expression

Hox genes encode transcription factors that specify embryonic positional identity in cells and guide tissue differentiation. Recent advances have greatly increased our understanding of the epigenetic mechanisms that ensure the faithful expression of Hox genes in adult cells and which involve the interplay of histone methylation, demethylation and intergenic transcription of long non-coding RNAs. The transcriptional memory of Hox genes poses both an opportunity and a challenge for regenerative medicine. Matching the positional identity of transplanted stem cells with that of the host environment, as reflected by their respective Hox profiles, is likely to be required to achieve regenerative healing. Strategies to manipulate the plasticity of Hox gene expression will probably become a major focus in regenerative medicine.

The genomic analysis of lactic acidosis and acidosis response in human cancers

The tumor microenvironment has a significant impact on tumor development. Two important determinants in this environment are hypoxia and lactic acidosis. Although lactic acidosis has long been recognized as an important factor in cancer, relatively little is known about how cells respond to lactic acidosis and how that response relates to cancer phenotypes. We develop genome-scale gene expression studies to dissect transcriptional responses of primary human mammary epithelial cells to lactic acidosis and hypoxia in vitro and to explore how they are linked to clinical tumor phenotypes in vivo. The resulting experimental signatures of responses to lactic acidosis and hypoxia are evaluated in a heterogeneous set of breast cancer datasets. A strong lactic acidosis response signature identifies a subgroup of low-risk breast cancer patients having distinct metabolic profiles suggestive of a preference for aerobic respiration. The association of lactic acidosis response with good survival outcomes may relate to the role of lactic acidosis in directing energy generation toward aerobic respiration and utilization of other energy sources via inhibition of glycolysis. This “inhibition of glycolysis” phenotype in tumors is likely caused by the repression of glycolysis gene expression and Akt inhibition. Our study presents a genomic evaluation of the prognostic information of a lactic acidosis response independent of the hypoxic response. Our results identify causal roles of lactic acidosis in metabolic reprogramming, and the direct functional consequence of lactic acidosis pathway activity on cellular responses and tumor development. The study also demonstrates the utility of genomic analysis that maps expression-based findings from in vitro experiments to human samples to assess links to in vivo clinical phenotypes.

It is well recognized that tumor microenvironments play an important role in modulating tumor progression in human cancers. Although previous studies have highlighted the importance of hypoxia, there is limited knowledge on the effects of other components in tumor microenvironments. Therefore, we use gene expression to compare and analyze how cells respond to lactate, acidity, and hypoxia, as well as how these responses can be utilized to predict the clinical outcomes of patients with breast cancers. We uncover an unexpected association with better clinical outcome of the strong lactic acidosis and acidosis response in breast cancers as a result of their abilities to inhibit glycolysis and favor oxidative phosphorylation for energy generation. This effect is caused by not only the repression of the gene expression of glycolysis genes but also the inhibition of Akt activation of cells exposed to lactic acidosis and acidosis. In conclusion, we propose that lactic acidosis and acidosis to be considered as independent prognostic factors for human cancers.

Evidence for Hox-specified positional identities in adult vasculature

BACKGROUND

The concept of specifying positional information in the adult cardiovascular system is largely unexplored. While the Hox transcriptional regulators have to be viewed as excellent candidates for assuming such a role, little is known about their presumptive cardiovascular control functions and in vivo expression patterns.

RESULTS

We demonstrate that conventional reporter gene analysis in transgenic mice is a useful approach for defining highly complex Hox expression patterns in the adult vascular network as exemplified by our lacZ reporter gene models for Hoxa3 and Hoxc11. These mice revealed expression in subsets of vascular smooth muscle cells (VSMCs) and endothelial cells (ECs) located in distinct regions of the vasculature that roughly correspond to the embryonic expression domains of the two genes. These reporter gene patterns were validated as authentic indicators of endogenous gene expression by immunolabeling and PCR analysis. Furthermore, we show that persistent reporter gene expression in cultured cells derived from vessel explants facilitates in vitro characterization of phenotypic properties as exemplified by the differential response of Hoxc11-lacZ-positive versus-negative cells in migration assays and to serum.

CONCLUSION

The data support a conceptual model of Hox-specified positional identities in adult blood vessels, which is of likely relevance for understanding the mechanisms underlying regional physiological diversities in the cardiovascular system. The data also demonstrate that conventional Hox reporter gene mice are useful tools for visualizing complex Hox expression patterns in the vascular network that might be unattainable otherwise. Finally, these mice are a resource for the isolation and phenotypic characterization of specific subpopulations of vascular cells marked by distinct Hox expression profiles.