Research | Open | Published:
Unexpected evolutionarily conserved rapid effects of viral infection on oxytocin receptor and TGF-β/pSmad3
Skeletal Musclevolume 7, Article number: 7 (2017)
shRNA lentiviral vectors are extensively used for gene knockdowns in mammalian cells, and non-target shRNAs typically are considered the proper experimental control for general changes caused by RNAi. However, the effects of non-target lentivirus controls on the modulation of cell signaling pathways remain largely unknown. In this study, we evaluated the effect of control lentiviral transduction on oxytocin receptor (OXTR) expression through the ERK/MAPK pathway in mouse and human skeletal muscle cells, on myogenic activity, and in vivo on mouse muscle regeneration. Furthermore, we mined published data for the influence of viral infections on OXTR levels in human populations and found that unrelated viral pathologies have a common consequence: diminished levels of OXTR.
We examined the change in OXTR mRNA expression upon transduction with control and Smad3-targeting viral vectors through real time RT-PCR and Western blotting, and confirmed with immunofluorescence. Changes in Smad3 and OXTR expression were examined both in vitro with mouse and human myoblasts and in vivo in mouse satellite cells. The general effects of viral infections on OXTR downregulation in humans were also examined by analyzing published Gene Expression Omnibus (GEO) datasets. The change in myoblast myogenic activity caused by the viral transduction (the percent of Pax7 + Ki67+ cells) was examined by immunofluorescence.
Results shown in this work establish that lentiviral control vectors significantly downregulate OXTR expression at mRNA and protein levels and diminish key downstream effectors of OXTR, ERK signaling, reducing the myogenic proliferation of infected cells. This effect is evolutionarily conserved between mouse and human myogenic cells, and it manifests in satellite cells after control lentiviral transduction of mice in vivo. Furthermore, an examination of published datasets uncovered similar OXTR downregulation in humans that are afflicted with different viral infections. Additionally, cells transduced with Smad3-targeting shRNA downregulate OXTR even more than cells transduced with control viruses.
Our work suggests that experimental cohorts transduced with control viruses may not behave the same as un-transduced cells and animals, specifically that control viral vectors significantly change the intensity of key cell-signaling pathways, such as OXTR/ERK. Our results further demonstrate that lentiviral transduction significantly decreases myogenic proliferation and suggest that viral infections in general may play a role in decreasing muscle health and regeneration, a decline in metabolic health, and a lower sense of well-being, as these rely on effective OXTR signaling. Additionally, our data suggest pathway crosstalk between TGF-β/pSmad3 and OXTR, implying that sustained attenuation of the TGF-β/pSmad3 pathway will reduce pro-regenerative OXTR/pERK signaling.
RNA interference (RNAi) is a form of post-transcriptional gene silencing that is mediated by short double-strand RNAs known as small interfering RNA (siRNA). Experimental down-modulation of the canonical signal transduction by RNAi is a common practice used to access changes in cell behavior and hence, significance of a particular pathway. Control transduction viral particles typically are either empty vector backbones that do not activate the RISC pathway, or encode non-target, scrambled shRNA or shRNA to non-mammalian genes (GFP, for instance). The use of such negative control viruses presumes that there will be no cellular or organismal changes in key conserved cell-signaling pathways as compared to the viruses encoding shRNA that target specific pathway genes (e.g., Smad3, Notch, beta-catenin, etc.) . Even when pathways are not themselves targeted, changes in their signaling intensity can impact the interpretation of the results. While lentiviral-mediated RNAi has proven to be efficient in gene silencing and knockdown, and control viruses serve well in this regard , little has been studied on how non-target lentiviral controls might affect the intensity of cell regulatory networks.
Previous studies have shown that canonical evolutionarily conserved signaling pathways, including Notch, TGF-β, MAPK, and Wnt, are involved in regulating embryonic organogenesis and adult tissue maintenance and repair [2,3,4,5,6,7,8,9,10,11,12]; skewing any of these pathways leads to a decline in tissue regenerative capacity [13,14,15,16,17,18]. For example, transforming growth factor beta (TGF-β) and the downstream pSmad3, have been shown to become elevated with age and inflammation, thereby suppressing muscle regeneration and adult neurogenesis [16, 19]. Oxytocin (OXT), on the other hand, a hormone that acts through the OXTR/pERK pathway, has been established as an indispensable determinant of muscle maintenance and regeneration that declines with age . OXT signaling has been also implicated in bone maintenance, reduction of obesity [21,22,23,24,25] in addition to its better-known role in formation of trust and bonding, and enhancement of OXT signaling is used in clinical trials for combating autism and age-related dementia [26,27,28,29,30].
As part of a study whether there is a correlation between the TGF-β pathway signaling and oxytocin receptor (OXTR) expression, we used shRNA vectors targeting Smad3 as a means to downregulate the TGF-β pathway and examined the corresponding change in OXTR expression. We were interested to find that reducing Smad3 appears to negatively affect OXTR expression, indicating TGF-β signaling opposes OXT. Perhaps not surprisingly, control viral vectors did not reduce TGF-β receptors or Smad3. But, surprisingly, the controls for these studies: lentiviral vectors bearing a non-target, non-mammalian shRNA, or even just the viral backbone, significantly downregulated OXTR expression levels in primary myogenic precursors, which in itself would be predicted to reduce muscle maintenance and regeneration , and in vivo in muscle satellite cells.
This suggests that not only experimental cohorts but also the controls transduced with control viruses (in vitro and in vivo) are likely to differ from the un-transduced groups in the relative intensity of key regulatory signaling pathways. And it further implicates that viral infections may play a role in decreasing OXTR and thus broadly interfering with regeneration and maintenance of multiple tissues.
To follow this lead, a thorough, unbiased analysis of published GEO datasets revealed that humans afflicted with a number of different viral infections manifest a specific downregulation of OXTR expression, but not their cell surface receptors in general. In this regard, our data suggest that viral infections will decrease muscle and bone maintenance and repair, as well as the sense of well-being and psychological health, all of which rely on OXTR signaling.
Additionally, our study establishes that experimental attenuation of the TGFβ/pSmad3 pathway with Smad3-targeted shRNA, further downregulates OXTR, and correspondingly diminishes pERK, as compared to the non-target shRNA control vector. These data suggest a TGF-β/pSmad3 and OXTR pathway-pathway cross-talk and implies that a sustained inhibition of the TGF-β/pSmad3 pathway (that becomes elevated with age) will attenuate oxytocin function that is already diminished with age.
22–24-month-old C57BL/6 mice were obtained from NIH (Bethesda, MD, USA).
Mouse primary myoblasts were cultured on Matrigel-coated dishes at 37 °C, 5% CO2 in growth medium (Ham’s F-10, Mediatech), penicillin/streptomycin antibiotics (500 IU/ml, 0.1 mg/ml; MP Biomedicals), and 20% Bovine Growth Serum (Life Technologies/Hyclone), supplemented with FGF-2 (6 ng/ml).
Human primary myoblasts were cultured on dilute Matrigel coated dishes (applied at 5 ug/cm^2 in phosphate-buffered saline (PBS)), at 37 °C, 5% CO2 in growth medium (Ham’s F-10; Mediatech), penicillin/streptomycin antibiotics (500 IU/ml, 0.1 mg/ml; MP Biomedicals), and 5% Bovine Growth Serum (Life Technologies/Hyclone), supplemented with FGF-2 (12 ng/ml).
shRNA delivery by lentiviral transduction
In vivo transduction:
All procedures were performed in accordance with the administrative panel of the Office of Laboratory Animal Care, UC Berkeley. Old tibias anterior and gastrocnemius muscles were infected, in vivo, with non-target shRNA (Sigma SHC002V) control transduction lentiviral particles. Mouse Smad3 shRNA-producing lentiviral particles (from Sigma) were used for in vivo transduction experiments (target-set generated from accession number NM_016769.2):
Smad3-targeted shRNA1: CCGGCCCATGTTTCTGCATGGATTTCTCGAGAAATCCATGCAGAAACATGGGTTTTTG
Smad3-targeted shRNA2: CCGGCCTTACCACTATCAGAGAGTACTCGAGTACTCTCTGATAGTGGTAAGGTTTTTG
Smad3-targeted shRNA3: CCGGCTGTCCAATGTCAACCGGAATCTCGAGATTCCGGTTGACATTGGACAGTTTTTG.
Lot: 11191506MN; 7.6 × 106 −1.4 × 107 TU/ml stocks were used a MOI 0.3 for each administration of each shRNA-encoding virus. Lentiviral particles were delivered into skeletal muscle by intramuscular injections for two consecutive days, and muscle was injured with cardiotoxin on the second day. 3 days after the injury, the in vivo infected myofiber-associated satellite cells were isolated, cultured overnight, and studied by the RT-PCR.
In vitro mouse and human primary myoblast transduction
Primary mouse and human myoblasts were transduced in vitro for qRT-PCR, WB, and IF analysis. Primary myoblast cultures were transduced with different MOIs (see below) with Opti-MEM (Gibco by Life Technologies). Medium for control groups without any transduction was changed into Opti-MEM. 24 h later, medium containing lentiviral particles were removed from wells as well as the non-transducer control group, and fresh Opti-MEM were added to each group. Medium were changed every 24 h until sample collection.
Primary myoblasts were transduced with Smad3-targeted shRNA at MOI = 0.5, empty vector (SHC001V) at MOI = 0.005, 0.05, and 0.5, and with non-target shRNA (SHC002V) at MOI = 0.002, 0.02, 0.2, 0.5, 1, and 2. Viral titer for lentiviral particles are SHC001V (empty vector control): 7.7E7 TU/ml; SHC002V (non-target shRNA control): 2.6E7 TU/ml
RNA extraction, RT-PCR, and real-time PCR
Total RNA was extracted from mouse and human primary myoblast cell culture using RNeasy Mini Kit (QIAGEN) according to manufacturer’s instructions. Total RNA was extracted from mouse gastrocnemius and tibialis anterior muscle using QIAshredder (QIAGEN) and RNeasy Mini Kit (QIAGEN) according to manufacturer’s instructions. Reverse transcription was performed with Invitrogen Superscript III First-Strand Synthesis System for RT-PCR according to manufacturer’s instructions. For real-time PCR amplification and quantification of gene of interest, 1 ug of total cDNA was used for the initial amplification using specific primers to each gene of interest; amplification was performed with a denaturation step at 95 °C for 10 min, followed by 40 cycles of denaturation at 95 °C for 15 s, and primer extension at 60 °C for 30 s. Real-time PCR was performed using PowerSYBR Green PCR Mastermix (from Applied Biosystems) under QuantStudio3 (Applied Biosystems) and CFX Connect Real-Time System (Bio-Rad). Reactions were run in triplicates. Housekeeping gene GAPDH was used as an internal control to normalize the variability in expression levels, and results were analyzed using the 2-ΔΔCT method described .
Primers used in real-time PCR:
Smad3 5′-F: CTGGCTACCTGAGTGAAGATGGAGA
Smad3 5′-R: AAAGACCTCCCCTCCGATGTAGTAG
Smad3 3′-F: ACACATTGGGAGAGGTGTGC
Smad3 3′-R: GCAAGGGTCCATTCAGGTGT
Smad3 shRNA-2-F: GCCTTACCACTATCAGAGAG
Smad3 shRNA-2-R: AACCTTACCTCATCAGAGAG
Smad3 shRNA-3-F: AACTCTCCAATGTCAACCG
Smad3 shRNA-3-R: GCTGTCCAATGTCAACCG
Real-time PCR gel analysis
Four percent agarose gel was made using agarose from Fisher Scientific and TAE buffer from Biosciences. 6× loading dye from Fermentas was added to the final amplification product and 15 ul of the mixture was added to the gel. QuickLoad 100 bp DNA Ladder from New England Biolabs was used as reference. Pictures of the gel were taken with a BioRad GelDoc/ChemiDoc Imager and Quantity One software. Pixel density was then analyzed with ImageJ (NIH) by subtracting the background pixel density and normalizing each gene of interest to GAPDH respectively.
Western blot analysis
Cells were lysed in RIPA buffer containing 1 mM PMSF, 1 mM sodium orthovanadate, PhosSTOP phosphatase inhibitor cocktail (Roche), and complete protease inhibitor cocktail (Roche). 30 ug of total protein extract in Laemmli buffer were resolved by SDS-PAGE on 10% precast gels (TGX, Bio-Rad) and were transferred to PVDF membranes (Millipore). Membranes were blocked for 1 h in 5% non-fat milk in PBST at room temperature. Primary antibodies against GAPDH (1:2000), actin (1:1000), pERK1/2 (1:1000), total ERK1/2 (1:1000), pSmad2/3 (1:1000), total Smad2/3 (1:1000), and OXTR (1:1000) were diluted in 5% non-fat milk in PBST. PVDF membranes were incubated in antibody solutions either overnight at 4 °C or 2 h at room temperature. Horseradish peroxidase-conjugated secondary antibodies were diluted 1:2000 in 1% BSA, and membranes were incubated for 1 h at room temperature. Blots were developed using ECL reagents (Advansta and Thermo Scientific) and were analyzed with Bio-Rad Gel Doc/Chemi Doc Imaging System and Quantity One software. Results of multiple assets were quantified by digitalizing the data and by normalizing pixel density of examined proteins by GAPDH pixel density with ImageJ (NIH).
Smad2/3: Cell Signaling Technology #3102
pSmad 2/3: Cell Signaling Technology #8828
ERK1/2: abcam ab184699
pERK1/2: Cell Signaling Technology #9101S
GAPDH: abcam ab9485
beta-actin: (ThermoFisher MA5-15739)
OXTR: proteintech 23045-1-AP
Secondary antibodies (all from Santa Cruz Biotechnology)
Bovine anti-rabbit IgG-HRP: sc-2370
Goat anti-rabbit IgG-HRP: sc-2004
Donkey anti-goat IgG-HRP: sc-2020
Bovine anti-goat IgG-HRP: sc-2350
Bovine anti-mouse IgG-HRP: sc-2371
Goat anti-mouse IgG-HRP: sc-2005
Goat anti-rat IgG-HRP: sc-2006.
Myoblasts were cultured in chamber slides (LabTek CC2 coated glass) and were transduced with lentivirus as mentioned above. For OXTR staining, the slides were fixed in cold 4% paraformaldehyde (PFA) on ice for 5 min and were washed three times in PBS. For Ki67/Pax7 staining, the slides were fixed in 70% cold ethanol at 4 °C overnight. The slides were then permeabilized with 0.25% Triton-100 in PBS for 5 min, and were blocked in staining buffer (1% calf serum in PBS) for 2 h. Primary antibodies were added into staining buffer and were incubated for at least 4 h at room temperature or overnight at 4 °C. Three PBS washes were performed before secondary antibody were added and incubated overnight at 4 °C. Slides were then washed three times with PBS before mounted for fluorescence imaging.
Invitrogen Alexa Fluor 488 Donkey anti-Rabbit IgG (H + L) Secondary Antibody
Invitrogen Alexa Fluor 488 Goat anti-Mouse IgG (H + L) Secondary Antibody
GEO dataset analysis
Patient studies for viral infections were identified from Gene Expression Omnibus datasets (GEO https://www.ncbi.nlm.nih.gov/geo/). Every dataset with information on OXTR, TGFβR1, or TGFβR2 were analyzed, regardless of cell types and types of viral infection. The fold changes of each gene measured in patients were normalized to non-infected control groups in a human study. For datasets with time-course experiments, time points starting and after 72 h were calculated. The accessions of GEO datasets used are as follows: GSE18816, GSE23031, GSE16593, GSE23031, GSE2067, GSE2067, GSE3292, GSE13597, GSE6802, GSE2815, GSE48466, GSE24533, GSE3397, GSE22589, GSE49954, GSE18816, GSE4785, GSE20755, GSE20948, GSE3397, GSE30719, and GSE27131.
A non-paired, two-tailed t test was performed on all of the respective data. Quantified data were presented as means (SD). P values of < 0.05 were considered statistically significant.
Control lentiviral vectors downregulates OXTR in mouse myogenic cells in vitro and mouse muscle stem cells in vivo
Lentiviral-mediated RNAi and transgenic gene expression have proven to be highly efficient and typically use empty vector backbones and non-target shRNA as controls. However, the influences of the vector controls on genes other than targeted were not well studied.
We transduced primary mouse myoblasts with Smad3-targeting shRNA and non-mammalian, non-target shRNA controls and compared the expression levels of Smad3 and OXTR with those of primary mouse myoblasts receiving a mock transduction of no viral vectors. As expected, the non-mammalian shRNA control served as a reliable reference for the downregulation of Smad3 mRNA by the Smad3 targeting shRNA (Fig. 1a). However, interestingly, the non-target shRNA vectors significantly downregulated OXTR expression by ~70%, as compared to the un-transduced primary mouse myoblasts (Fig. 1a). Furthermore, the downregulation of Smad3 by the Smad3-targeted shRNAs further diminished OXTR expression, as compared to non-target shRNAs (Fig. 1a, Additional file 1: Figure S1A).
To examine if a similar phenomenon occurs in vivo, old C57.B6 mice, in which TGF-β/pSmad3 pathway becomes elevated with age, were injected intramuscularly with non-target non-mammalian (GFP) shRNA viral particles, versus three different Smad3-targeting shRNAs as in . The transduced muscle satellite cells were isolated from the in vivo injected muscles 3 days later, as we have published . The in vivo transduced satellite cells showed the expected downregulation of Smad3 when muscle was injected with Smad3-targeting shRNA viral particles, as compared to the non-target shRNA vector (Fig. 1b). Very interestingly, the in vitro observed trend of OXTR downregulation by viral transduction was also detected in these in vivo studies (Fig. 1b, Additional file 1: Figure S1B). Similar to the in vitro results, Smad3 shRNA in vivo additionally downregulated OXTR, as compared to the non-target shRNA viral vector (Fig. 1b).
While non-mammalian control shRNA vectors are designed to not target any known, canonical mammalian sequences, the presence of shRNA will engage with and activate the RNA-induced silencing complex (RISC). Interested in whether the downregulation of OXTR comes from viral transduction or RISC activation, primary mouse myoblasts were transduced with an empty vector (encoding just the viral backbone) that does not contain a hairpin insert, or with the above non-mammalian (GFP) shRNA control particles. OXTR expression levels were compared with those of non-transduced primary myoblasts. Our data shows that both control vectors caused a significant (~70%) downregulation in OXTR expression even at MOI = 0.05, which is much lower than the titer that is typically recommended for viral transductions; and that control shRNA particles downregulated OXTR faster than the viral backbone particles (Additional file 1: Figure S2A, B).
These data establish that transduction with control viral vectors perturbs key cell signaling molecules, such as OXTR at low viral titers.
OXTR downregulation upon lentiviral transduction is evolutionarily conserved between mice and humans, and manifests upon different viral infections in human populations
Interested in exploring the conservation of the observed phenomena between mouse and human cells, we repeated the above in vitro experiments with primary human myoblasts. Transduction of these human muscle precursors with non-mammalian, non-target (GFP) shRNA controls showed significant downregulation in OXTR expression, suggesting that OXTR downregulation upon lentiviral transduction is evolutionarily conserved between mouse and human. And similar to the outcomes with mouse myoblasts, this attenuation became stronger when a Smad3-targeting shRNA vector was used (e.g., when human Smad3 expression was experimentally diminished (Fig. 2B)). An examination of the complementarity demonstrated that mouse and human Smad3 loci have very high homology in the area targeted by this specific shRNA (Additional file 1: Figure S3) consistently with our data that this specific Smad3 shRNA targets Smad3 in both species.
Extrapolating these data, we asked whether OXTR downregulation would generally be caused by viral infections in humans and analyzed published GEO databases from human studies that focused on different viral infections (HIV, SIV, influenza virus, etc.). Specifically, we data-mined previously published GEO databases that report the levels of OXTR, Smad3, TGFBR1, and TGFBR2, which were analyzed with statistical rigor in the presence of various viral infections (in various cell types), as compared to the non-infected control groups within each human study.
Interestingly, this analysis confirmed and extrapolated our conclusions by revealing that OXTR becomes significantly downregulated upon different unrelated viral infections in human populations (Fig. 3a). In addition to pooling the results from all databases, the data from each individual database were also examined: for each group, OXTR showed downregulation upon viral infections ranging from 47 to 82%, as compared with the level of this receptor in healthy individuals (Fig. 3b).
These results demonstrate an evolutionary conservation of OXTR attenuation by experimental and natural viral infections and indicate a significance of the uncovered phenomenon for human health.
Lentiviral transductions do not downregulate cell surface receptors broadly
To answer whether the downregulation of receptors by non-target shRNA vectors applies to only certain proteins exemplified by OXTR, or manifests more generally, we examined the levels of TGF-β receptors 1 and 2 (TGFBR1 and TGFBR2) in primary mouse and human myoblasts by real-time qRT-PCR and western blotting. Interestingly, in contrast to the OXTR, expression levels of TGFBR1 and 2 were not significantly altered upon control lentivirial transduction in either mouse or human myogenic cells, and TGFBR2 actually became elevated in human myoblasts (Fig. 4a, b).
Moreover, while OXTR was significantly downregulated (Fig. 3), the expression levels of TGFBR1, TGFBR2, and Smad3 genes were not diminished in people afflicted with viral infections, as per our data mining study with the above-described human GEO databases (Fig. 4c). As for Smad3 expression, some human studies showed an obvious upregulation in expression upon viral infection, while others showed no change from healthy controls (Fig. 4c).
The expression levels of OXTR were further examined at protein levels through western blotting and were found to be reduced after the control lentiviral transduction, in a time-dependent manner (Fig. 4d). This downregulation of OXTR was further verified by immunofluorescence (Additional file 1: Figure S4). Furthermore, a key signaling molecule down-stream of OXTR, pERK , became attenuated concordantly with OXTR and with the same timing (Fig. 4d). In contrast, the protein levels of TGFBR2 did not change significantly, in agreement with our qRT-PCR results.
These results establish that viral infection promotes selected downregulation of certain, but not all cell surface receptors, which may perturb signaling networks. More specifically, these data suggest a skewing from OXTR/pERK to TGF-β/pSmad3 signaling upon viral infections, which is reminiscent of the age-imposed changes that associate with broad degenerative pathologies . Interestingly, in human populations such conserved pathway skewing took place regardless of the specific viral infections or their distinct mechanisms of action or diseases.
Viral infections decrease myogenic activity in vitro
Since OXTR signaling plays a major role in maintaining myogenic activity and muscle health , we hypothesized that control viral transductions, which attenuate OXTR, would also diminish the myogenicity of even young muscle progenitor cells.
To examine this hypothesis, we compared the in vitro proliferation of myogenic (Pax7+) cells that were derived from young mice. Specifically, primary mouse myoblasts were transduced with either empty vector (encoding just the viral backbone) or with non-mammalian (GFP-encoding) shRNA particles and were compared with the mock-transduced cohort. All cells were fixed at 5 days after the transduction and were analyzed by co-immunodetection of Pax7+ and the proliferation marker Ki67 (Fig. 5a). The percentage of Ki67+/Pax7+ double positive, as well as Ki67+ and Pax7+ single positive myoblasts out of all Hoechst + cells were quantified and compared between the cohorts (Fig. 5b).
Interestingly and in complete agreement with our above observations, as compared to the un-transduced myoblasts, transduction with either the viral backbone or with the control shRNA particles dramatically and with high statistical significance reduced the percent of proliferating primary muscle progenitor cells from ~85% down to ~20%. Additionally, the numbers of cells expressing Pax7+ also significantly declined, proliferating or not (Fig. 5b), suggesting that viruses inhibit this key myogenic transcriptional factor . When measured separately from the myogenic marker Pax7, the cell proliferation, e.g., percent of Ki67+ cells, has also declined upon the control viral transduction (Fig. 5b), which is consistent with the inhibition of p21 by OXTR . Since so few cells were proliferating, it is unlikely that Pax7− cells overtook the culture during this experiment.
These data suggest that viral infections might play a negative role in skeletal muscle maintenance by suppressing myogenic proliferation and may lead to the loss of muscle stem cell identity through downregulation of Pax7 .
Numerous studies have been conducted with viral backbone and non-targeting shRNA vectors in mammalian cells and in vivo, and these have been assumed to control for changes in cell behavior and gene expression as compared to the experimental gene delivery. The data presented in our work shifts this paradigm by demonstrating that non-target non-mammalian shRNA encoding and even viral backbone-encoding viruses significantly downregulate a canonical cell-signaling pathway (OXTR-ERK) and perturb the TGF-β/OXT balance. Notably, our results show that cells and animals exposed to control viruses (in vitro and in vivo) will differ from un-transduced cohorts in a number of ways.
OXTR positively regulates muscle stem cell responses through the MAPK/ERK pathway, which among other broad effects attenuates CDK inhibitors and promote myogenic proliferation ; and in addition to its role in muscle health, oxytocin-OXTR is a determinant in trust and bonding, child bearing, mental well-being, and metabolic homeostasis, and it is also implicated in bone health and metabolism [21,22,23,24,25,26]. In this regard, our data suggest that in addition to known pathological manifestations, viral infection may play a less direct role in decreasing the health and regeneration of skeletal muscle, bone, brain-psychological well-being and inducing obesity by skewing key cell signaling networks that regulate homeostatic tissue maintenance.
While we show that viral infections and experimental transductions do not broadly attenuate all cell-surface receptors, it is likely that more than just OXTR becomes downregulated. For example, Pax7 rapidly declined in primary myoblasts transduced with the backbone or GFP shRNA control vectors. In this regard, it would be interesting to examine whether and how Pax7 and other myogenic cell-fate control genes  experience a virally induced attenuation (OXTR-dependent or by a parallel pathway).
Finally, our study establishes that TGF-β/pSmad3 and OXTR pathways actively interact and imply that sustained attenuation of the TGF-β/pSmad3 pathway (elevated with age) might be not entirely positive for tissue health and maintenance, as with time it will block the function of the hormone oxytocin.
Lentiviral control vectors downregulate the OXTR/ERK signaling, allowing the TGF-β/pSmad pathway to become prevalent. This phenomenon is evolutionary conserved between mice and humans; and it manifests in people who suffer from different viral infections. Skewing from OXTR to TGF-β/pSmad is predicted to diminish the maintenance and repair of skeletal muscle (and possibly, other tissues), and accordingly, control viral infections diminish myogenicity.
Mothers against decapentaplegic homolog 3
TGFβ receptor 1
TGFβ receptor 2
Transforming growth factor beta
STEWART S. Lentivirus-delivered stable gene silencing by RNAi in primary cells. RNA. 2003;9:493–501.
Artavanis-Tsakonas S. Notch signaling: cell fate control and signal integration in development. Science. 1999;284:770–6.
Brack A, Conboy I, Conboy M, Shen J, Rando T. A temporal switch from notch to Wnt signaling in muscle stem cells is necessary for normal adult myogenesis. Cell Stem Cell. 2008;2:50–9.
Canning C, Lee L, Luo S, Graham A, Jones C. Neural tube derived Wnt signals cooperate with FGF signaling in the formation and differentiation of the trigeminal placodes. Neural Dev. 2008;3:35.
Carlson M, Conboy M, Hsu M, Barchas L, Jeong J, Agrawal A, Mikels A, Agrawal S, Schaffer D, Conboy I. Relative roles of TGF-Î21 and Wnt in the systemic regulation and aging of satellite cell responses. Aging Cell. 2009;8:676–89.
Heine U. Role of transforming growth factor-beta in the development of the mouse embryo. J Cell Biol. 1987;105:2861–76.
Kook S, Son Y, Choi K, Lee H, Chung W, Hwang I, Lee J. Cyclic mechanical stress suppresses myogenic differentiation of adult bovine satellite cells through activation of extracellular signal-regulated kinase. Mol Cell Biochem. 2007;309:133–41.
Lewis J. Notch signalling and the control of cell fate choices in vertebrates. Semin Cell Dev Biol. 1998;9:583–9.
Pelton R. Immunohistochemical localization of TGF beta 1, TGF beta 2, and TGF beta 3 in the mouse embryo: expression patterns suggest multiple roles during embryonic development. J Cell Biol. 1991;115:1091–105.
Shi D, Bourdelas A, Umbhauer M, Boucaut J. Zygotic Wnt/β-catenin signaling preferentially regulates the expression of Myf5 gene in the mesoderm of xenopus. Dev Biol. 2002;245:124–35.
Wagner J, Schmidt C, Nikowits W, Christ B. Compartmentalization of the somite and myogenesis in chick embryos are influenced by Wnt expression. Dev Biol. 2000;228:86–94.
Volonte D, Liu Y, Galbiati F. The modulation of caveolin-1 expression controls satellite cell activation during muscle repair. FASEB J. 2005;19(2):237–9.
Brionne T, Tesseur I, Masliah E, Wyss-Coray T. Loss of TGF-β1 Leads to increased neuronal cell death and microgliosis in mouse brain. Neuron. 2003;40:1133–45.
Brack A, Conboy M, Roy S, Lee M, Kuo C, Keller C, Rando T. Increased Wnt signaling during aging alters muscle stem cell fate and increases fibrosis. Science. 2007;317:807–10.
Carlson M, Suetta C, Conboy M, Aagaard P, Mackey A, Kjaer M, Conboy I. Molecular aging and rejuvenation of human muscle stem cells. EMBO Molec Med. 2009;1:381–91.
Yousef H, Conboy M, Morgenthaler A, Schlesinger C, Bugaj L, Paliwal P, Greer C, Conboy I, Schaffer D. Systemic attenuation of the TGF-β pathway by a single drug simultaneously rejuvenates hippocampal neurogenesis and myogenesis in the same old mammal. Oncotarget. 2015;6:11959–78.
Conboy I, Conboy M, Rebo J. Systemic problems: a perspective on stem cell aging and rejuvenation. Aging. 2015;7:754–65.
Carlson ME, Hsu M, Comboy IM. Imbalance between pSmad3 and Notch induces CDK inhibitors in old muscle stem cells. Nature. 2008;454:528–32.
Carlson M, Silva H, Conboy I. Aging of signal transduction pathways, and pathology. Exp Cell Res. 2008;314:1951–61.
Elabd C, Cousin W, Upadhyayula P, Chen R, Chooljian M, Li J, Kung S, Jiang K, Conboy I. Oxytocin is an age-specific circulating hormone that is necessary for muscle maintenance and regeneration. Nat Comm. 2014;5. https://www.nature.com/articles/ncomms5082.
Beranger G, Pisani D, Castel J, Djedaini M, Battaglia S, Amiaud J, Boukhechba F, Ailhaud G, Michiels J, Heymann D, Luquet S, Amri E. Oxytocin reverses ovariectomy-induced osteopenia and body fat gain. Endocrinology. 2014;155:1340–52.
Camerino C. Low sympathetic tone and obese phenotype in oxytocin-deficient mice. Obesity. 2009;17:980–4.
Deblon N, Veyrat-Durebex C, Bourgoin L, Caillon A, Bussier A, Petrosino S, Piscitelli F, Legros J, Geenen V, Foti M, Wahli W, Di Marzo V, Rohner-Jeanrenaud F. Mechanisms of the anti-obesity effects of oxytocin in diet-induced obese rats. PLoS ONE. 2011;6:e25565.
Elabd C, Basillais A, Beaupied H, Breuil V, Wagner N, Scheideler M, Zaragosi L, Massiéra F, Lemichez E, Trajanoski Z, Carle G, Euller-Ziegler L, Ailhaud G, Benhamou C, Dani C, Amri E. Oxytocin controls differentiation of human mesenchymal stem cells and reverses osteoporosis. Stem Cells. 2008;26:2399–407.
Takayanagi Y, Kasahara Y, Onaka T, Takahashi N, Kawada T, Nishimori K. Oxytocin receptor-deficient mice developed late-onset obesity. NeuroReport. 2008;19:951–5.
Tamma R, Colaianni G, Zhu L, DiBenedetto A, Greco G, Montemurro G, Patano N, Strippoli M, Vergari R, Mancini L, Colucci S, Grano M, Faccio R, Liu X, Li J, Usmani S, Bachar M, Bab I, Nishimori K, Young L, Buettner C, Iqbal J, Sun L, Zaidi M, Zallone A. Oxytocin is an anabolic bone hormone. Proc Natl Acad Sci. 2009;106:7149–54.
Zou K, Lisio MD, Huntsman HD, Pincu Y, Mahmassani Z, Miller M, et al. Laminin-111 improves skeletal muscle stem cell quantity and function following eccentric exercise. Stem Cells Translat Med. 2014;3:1013–22.
Lee H, Macbeth A, Pagani J, Young W: Oxytocin: The Great facilitator of life. Prog Neurobiol. 2009;88(2):127–151.
Modahl C, Green L, Fein D, Morris M, Waterhouse L, Feinstein C, Levin H. Plasma oxytocin levels in autistic children. Biol Psychiatry. 1998;43:270–7.
Ruggeri B, Sarkans U, Schumann G, Persico A. Biomarkers in autism spectrum disorder: the old and the new. Psychopharmacology. 2013;231:1201–16.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 − ΔΔCT method. Methods. 2001;25:402–8.
Conboy MJ, Conboy IM. Preparation of adult muscle fiber-associated stem/precursor cells. Methods Mol Biol. 2010;621:149–63. doi:10.1007/978-1-60761-063-2_10.
Rudnicki MA, Jaenisch R. The MyoD family of transcription factors and skeletal myogenesis. BioEssays. 1995;17:203–9.
Berzsenyi M, Woollard D, McLean C, Preiss S, Perreau V, Beard M, Scott Bowden D, Cowie B, Li S, Mijch A, Roberts S. Down-regulation of intra-hepatic T-cell signaling associated with GB virus C in a HCV/HIV co-infected group with reduced liver disease. J Hepatol. 2011;55:536–44.
Blackham S, Baillie A, Al-Hababi F, Remlinger K, You S, Hamatake R, McGarvey M. Gene expression profiling indicates the roles of host oxidative stress, apoptosis, lipid metabolism, and intracellular transport genes in the replication of hepatitis C virus. J Virol. 2010;84:5404–14.
do Valle T, Billecocq A, Guillemot L, Alberts R, Gommet C, Geffers R, Calabrese K, Schughart K, Bouloy M, Montagutelli X, Panthier J. A new mouse model reveals a critical role for host innate immunity in resistance to Rift valley fever. J Immunol. 2010;185:6146–56.
Gerlach R, Camp J, Chu Y, Jonsson C. Early host responses of seasonal and pandemic influenza A viruses in primary well-differentiated human lung epithelial cells. PLoS ONE. 2013;8:e78912.
Gersten M, Alirezaei M, Marcondes M, Flynn C, Ravasi T, Ideker T, Fox H. An integrated systems analysis implicates EGR1 Downregulation in Simian immunodeficiency virus encephalitis-induced neural dysfunction. J Neurosci. 2009;29:12467–76.
Huang Y, Li Z, Hyseni X, Schmitt M, Devlin R, Karoly E, Soukup J. Identification of gene biomarkers for respiratory syncytial virus infection in a bronchial epithelial cell line. Genom Med. 2008;2:113–25.
Lee S, Chan R, Gardy J, Lo C, Sihoe A, Kang S, Cheung T, Guan Y, Chan M, Hancock R, Peiris M. Systems-level comparison of host responses induced by pandemic and seasonal influenza A H1N1 viruses in primary human type I-like alveolar epithelial cells in vitro. Respir Res. 2010;11:147.
Lund A, Lucero J, Herbert L, Liu Y, Naik J. Human immunodeficiency virus transgenic rats exhibit pulmonary hypertension. Am J Physiol Lung Cell Mol Physiol. 2011;301:L315–26.
Manel N, Hogstad B, Wang Y, Levy D, Unutmaz D, Littman D. A cryptic sensor for HIV-1 activates antiviral innate immunity in dendritic cells. Nature. 2010;467:214–7.
Ndolo T, George M, Nguyen H, Dandekar S. Expression of simian immunodeficiency virus Nef protein in CD4+ T cells leads to a molecular profile of viral persistence and immune evasion. Virology. 2006;353:374–87.
Nguyen H, Sankaran S, Dandekar S. Hepatitis C virus core protein induces expression of genes regulating immune evasion and anti-apoptosis in hepatocytes. Virology. 2006;354:58–68.
Shornick L, Wells A, Zhang Y, Patel A, Huang G, Takami K, Sosa M, Shukla N, Agapov E, Holtzman M. Corrections: airway epithelial versus immune cell Stat1 function for innate defense against respiratory viral infection. J Immunol. 2011;187:2834-2834.
Tyner J, Uchida O, Kajiwara N, Kim E, Patel A, O'Sullivan M, Walter M, Schwendener R, Cook D, Danoff T, Holtzman M. CCL5-CCR5 interaction provides antiapoptotic signals for macrophage survival during viral infection. Nat Med. 2005;11:1180–7.
Verhoeven D, George M, Hu W, Dang A, Smit-McBride Z, Reay E, Macal M, Fenton A, Sankaran-Walters S, Dandekar S. Enhanced innate antiviral gene expression, IFN-, and cytolytic responses are predictive of mucosal immune recovery during simian immunodeficiency virus infection. J Immunol. 2014;192:3308–18.
Zhao B, Zheng S, Gong L, Wang Y, Chen C, Jin W, Zhang D, Yuan X, Guo J, Duan Z, He Y. T Lymphocytes from chronic HCV-infected patients are primed for activation-induced apoptosis and express unique pro-apoptotic gene signature. PLoS ONE. 2013;8:e77008.
We would like to thank Yan Liu and Michael Conboy for help with experimental systems, helpful discussion of this work and editing of the manuscript.
This work was supported by the NIH R56, Rogers’ family, Calico/QB3, Packer Endowment, and SENS foundation awards to IMC.
Availability of data and materials
The datasets used and/or analyzed during the study are available from the corresponding author on reasonable requests.
YL performed the experiments shown in all figures, analyzed the data, and contributed to writing the manuscript. IMC interpreted the data and wrote the manuscript. Both authors read and approved the final manuscript.
The authors declare that they have no competing interest.
Consent for publication
Ethics approval and consent to participate
All procedures were performed in accordance with the administrative panel of the Office of Laboratory Animal Care, UC Berkeley. The protocol was approved by the UC Berkeley Animal Care and Use Committee. Old C57BL/6 mice (22–24 months) were obtained from NIH (Bethesda, MD, USA) and used in accordance with the approved protocol.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional comparative qRT PCR data. a. The in vitro studies with mouse primary myoblasts that were transduced with shRNA 1, 2, 3, and their mix. b. Comparative qRT PCR on un-transduced myoblasts vs. satellite sells, vs. control shRNA in vivo transduced satellite cells. N = 3, *,**, *** P < 0.05. Figure S2. Dose curve and timing of OXTR downregulation by control viruses. a. Non-target shRNA and empty vector particles downregulate OXTR in a dose-dependent manner. OXTR downregulation increases with viral dosage with both empty vector particles and non-target shRNA vectors. Myoblasts transduced with non-target shRNAs were collected at 72 h, and those transduced with empty vector particles at 120 h. OXTR downregulation is significant even at low viral dosage. N = 3 for each cohort. ***P < 0.05. b. Empty vector particles take longer to decrease OXTR expression compared with non-target shRNAs, but in both cases significant downregulation is detected by 120 h post transduction. N = 3 for each cohort and each time point. ***P < 0.05. Figure S3. High sequence homology between human and mouse Smad3 loci in the area targeted by shRNA are observed and are demonstrated by the alignment. Figure S4. Lentiviral transfection downregulates OXTR as determined by the immunofluorescence. Immunofluorescence for OXTR (green) was performed with un-transduced primary human myoblasts versus human myoblasts transduced with control lentiviral shRNA vectors. Hoechst (blue) labels all nuclei. OXTR-specific immunofluorescence was diminished in primary human myoblasts transfected with non-targeted shRNA vectors, as compared to the un-transduced cells. Figure S5. Gel confirmation of qRT PCR results. Amplification products from the real-time RT-PCR reactions were run on a gel, and the pixel intensities of each band were quantified by normalizing to GAPDH. *P ≤ 0.05; **P ≤ 0.01. (PDF 4347 kb)