|
|
||||||||
1Division of Allergy, Pulmonary, and Critical Care Medicine, Vanderbilt University, Nashville, Tennessee; 2Division of Pulmonary Sciences and Critical Care Medicine, University of Colorado Health Sciences Center, Denver, Colorado; 3Ottawa Health Research Institute, Ottawa, Ontario, Canada; and 4Center for Lung Biology, Department of Medicine, University of South Alabama, Mobile, Alabama
Submitted 1 April 2008 ; accepted in final form 18 August 2008
| ABSTRACT |
|---|
|
|
|---|
pulmonary circulation and disease; gene expression; genetically altered mice
The familial form is usually associated with mutations in the type II receptor for the bone morphogenetic protein pathway, BMPR2 (17). To determine the phenotype and in vivo molecular effect of BMPR2 mutation, we previously created a doxycycline-inducible transgenic mouse, which expressed an intracellular domain truncation of BMPR2, BMPR2delx4+, in smooth muscle on induction. These animals had increased right ventricular systolic pressures (RVSP) (36), increased cytokines (12, 33), loss of smooth muscle differentiation (33), and defects in vasoreactivity and vasoreactivity pathways (33, 42). However, pulmonary structural changes in this model were limited (36).
The mutation used in this original BMPR2delx4+ mouse was derived from family UK21 and consisted of a "T" base insertion at base 504, resulting in a premature stop 18 amino acids into the kinase domain (21). This resulted in loss of the entire intracellular domain and associated functions. The intracellular domain consists of a kinase domain at amino acids 200–500 responsible for phosphorylation of the type I BMP receptor leading to activation of SMAD transcription factors (4) and a long COOH-terminal tail from amino acids 500–1038. The kinase domain may also interact with other targets, including Rack1 (44). The BMPR2 tail has several poorly understood functions, including regulation of p38 and p42/44 MAPK (27, 40) and interaction with LIMK (10), c-Src (37), and Tctex (22). Some patient mutations are within this cytoplasmic tail and leave SMAD functions intact while causing disruptions to tail functions (21, 26). For example, an arginine to termination mutation at amino acid 899 (R899X) is found in family US33 (17).
We hypothesized that this mutation would result in a subset of the manifestations associated with the BMPR2delx4+ mutation, since SMAD activity remains intact with the R899X mutation. To test this hypothesis, we created transgenic mice that express the BMPR2 R899X mutation (BMPR2R899X) in smooth muscle when induced by doxycycline. After 9 wk of induction in adult mice, mice were hemodynamically phenotyped, fluorescent microangiography (FMA) was used to assess vascular pruning, and tissue sections were examined by immunohistochemistry. We found that expression of the BMPR2R899X mutation resulted in increased RVSP accompanied by substantial pulmonary vascular pruning and remodeling. Gene arrays were performed on transgenic mice with both high and low RVSP; we found that stress response, muscle organization and function, proliferation, apoptosis, and developmental pathways were all changed before RVSP increased. After elevation of RVSP, there were additional changes in these pathways as well as increased muscle structural genes and increased angiogenesis-related genes.
| METHODS |
|---|
|
|
|---|
Heart catheterization and FMA.
Mice are given tribromoethanol (500 mg/kg ip) to induce a surgical plane of anesthesia. The animals are then shaved to expose the surgical area. Mice are placed on a heated surgical table (872/1, 872/H; Harvard Apparatus, Holliston, MA) and secured with surgical tape. Systemic blood pressure and pulse is measured via a tail cuff and pulse transducer run through a PowerLab NIBP Controller (ADInstruments, Colorado Springs, CO). The surgical site is viewed using a Zeiss OPMI-1 surgical microscope. An incision of
1 in. in length is made extending from the animal's chin down to the right armpit. The thyroid gland is then blunt-dissected upward to expose the underlying tissue and the right jugular vein. The jugular vein is then separated from surrounding tissue using dissecting forceps until the body of the vessel is completely free from adherent tissues. The cranial end of the jugular is tied off completely, and a loose tie is then made at the caudal end of the exposed jugular using 4-0 braided silk suture. Four-inch microdissecting scissors are then used to make a small incision in the medial aspect of the right jugular vein. A Millar 1.4 French pressure-volume microtip catheter transducer (SPR-839; Millar Instruments, Houston, TX) connected to a PowerLab/8s (ADInstruments) is then inserted through the incision and gently threaded down into the right ventricle. Proper placement within the ventricle is determined through observation of the pressure-volume loop obtained from the catheter. The loose caudal suture is then tightened to secure the catheter in place. Once the catheter is properly placed, data are collected using Chart 5 (ADInstruments).
Once blood pressure and volume data are collected, the caudal suture is reloosened, and the catheter removed. One hundred units of heparin is then injected through the jugular incision to prevent clotting. The suture is then retightened to prevent bleeding. The animals are then removed from the surgical table to a dissecting area.
Plane of anesthesia is redetermined postsurgery, and an overdose of sodium pentobarbital is administered if the withdrawal reflex is returning. After complete loss of withdrawal reflex, an incision is made just below the xiphoid process, and the ventral portion of the rib cage is cut away being careful not to nick the lungs. A small cut is then made at the vertex of the left atrium. A needle is inserted into the right ventricle, and
10 ml of a flushing solution (10 mM sodium phosphate, 127 mM sodium chloride, 10 mM EDTA, 5 U/ml heparin) is flushed through the lungs and out of the left atrium. A 45°C solution of 10% (1:10 vol/vol) fluorescent microbeads (FluoSpheres F8811; Molecular Probes, Eugene, OR) in 1% low melting point (LMP) agarose in PBS is then slowly flushed via the right ventricle through the lungs until it runs out of the left atrium. Lungs are simultaneously inflated with a 0.8% solution of LMP agarose in PBS via the trachea. The chest is packed off with ice until the LMP agarose has congealed. The heart and lungs are removed en bloc and fixed 48 h in 4% paraformaldehyde. Sections for FMA are cut using a vibratome.
Luciferase reporter assay. A7r5 [CRL-1444, American Type Culture Collection (ATCC)] cells were plated and grown in 12-well plates in DMEM with 10% FBS and 1% penicillin-streptomycin. When the cells were 60–70% confluent, wells were washed, and the media replaced with 0.5% FBS DMEM overnight. The following day, the cells were transfected with 3 µg of pBRE, 1.2 µg of TK Renilla (for internal control), and 3 µg of either control vector (pGEM), wild-type Bmpr2, or BMPR2 with R899X mutation, using FuGENE 6 (Roche, Indianapolis, IN) as transfecting medium at a 3:1 FuGENE-to-DNA ratio.
After 24 h, cells received 0, 5, 15, or 50 ng/ml recombinant BMP4. After 48 h from transfection, the cells were lysed directly in the plates using 200 µl of 1x lysis buffer from Dual-Luciferase Assay Kit (Promega). Cells were rocked for 15 min subjected to one round of freeze-thaw before analysis. Twenty microliters of the lysate and 100 µl each of Luciferase Assay Reagent II and Stop & Glo Reagent were used for each assay. The assays were carried out using a Turner BioSystems Luminometer with dual injectors.
Western blot and immunohistochemistry. Western blots were performed as previously described (12) using primary antibodies at 1:200 dilution and secondary antibodies at 1:3,000 dilution. Primary antibodies included Id1 (sc-488; Santa Cruz Biotechnology, Santa Cruz, CA), phospho-Smad1 (cat. no. 9511; Cell Signaling Technology, Danvers, MA), and phospho-p38 (sc-7973, Santa Cruz Biotechnology).
Immunohistochemistry was performed as previously described (33) using the following primary antibodies and dilutions: Flag Ab-1 NeoMarkers (Fremont, CA) rabbit polyclonal, 1:250; actin A5441 Sigma (St. Louis, MO) mouse monoclonal, 1:1,000; actin Abcam (Cambridge, MA) ab5694 rabbit polyclonal, 1:1,000; von Willebrand factor (vWF) Dako A0082 rabbit polyclonal, 1:1,000; CD45 sc-25590 Santa Cruz Biotechnology rabbit polyclonal, 1:100; CD133 sc-30219 Santa Cruz Biotechnology rabbit polyclonal, 1:200; CD3 epsilon ab49943 Abcam rabbit polyclonal, 1:200; Mac-3 550292 BD Pharmingen (San Jose, CA) mouse monoclonal, 1:200; CD11b 550282 BD Pharmingen mouse monoclonal, 1:200; phospho-p42/44 9106 Cell Signaling Technology mouse monoclonal, 1:50.
Affymetrix arrays. Total RNA was isolated using an RNeasy Mini Kit (QIAGEN, Valencia, CA) from 30 mg of right lung from six individual mice. Samples were prepared for Affymetrix arrays using 2.5 µg of total RNA. First and second strand complimentary DNA was synthesized using standard techniques. Biotin-labeled antisense complimentary RNA was produced by an in vitro transcription reaction. Mouse Genome 430 2.0 microarrays (Affymetrix, Foster City, CA) were hybridized with 20 µg of cRNA. Target hybridization, washing, staining, and scanning probe arrays were done following an Affymetrix GeneChip Expression Analysis Manual. All array results have been submitted to the NCBI gene expression and hybridization array data repository (GEO; http://www.ncbi.nlm.nih.gov/geo/) as series GSE11018 [NCBI GEO] .
Array analysis. Affymetrix CEL files were loaded into dChip 2007 array analysis software. The dChip algorithm is capable of detecting significant differences at signal strengths lower than those usable in Microarray Suite (Affymetrix, Santa Clara, CA; Ref. 19). Overall signal strength from arrays was normalized to the median array, and expression levels were determined using the perfect match/mismatch (PM/MM) algorithm. Gene ontology was determined using the Classify Genes tool within dChip (45) with gene ontology files downloaded from the Gene Ontology Consortium (http://www.geneontology.org; Ref. 10a). Our definition for changed genes was a 1.2x change for each of four pairwise comparisons with a minimum absolute change of 100, which results in a false discovery rate of 27%.
| RESULTS |
|---|
|
|
|---|
-smooth muscle actin in both pulmonary artery and airway smooth muscle (Fig. 1B). Columnar epithelial cells and elastic lamina are false positive as shown by secondary antibody-only control (Supplemental Fig. 1A, available in the data supplement online at the AJP-Lung Cellular and Molecular Physiology web site). Supplemental Fig. 1B shows serial sections double-stained with either FLAG tag and actin or vWF (endothelial marker) and actin, showing that nonmuscularized vessels do not express FLAG tag, further supporting smooth muscle specificity. Although we cannot rule out expression below the threshold of detection in other tissue types, these data indicate that the transgene is substantially made in the correct cell type and is correctly expressed on the cell surface.
|
|
1.5x normalized to β-actin. We also found an increase in p42/44 phosphorylation by immunohistochemistry (Supplemental Fig. 2C), apparently in both endothelium and smooth muscle, implying a paracrine effect. Contrary to expectation, we did not see alteration in phosphorylation of cofilin, a LIMK target (Supplemental Fig. 2D) in whole lung homogenates. We cannot rule out a tissue-specific effect obscured by use of whole lung, as array data below suggest alterations in actin organization pathways. Nine weeks of BMPR2R899X expression results in elevated RVSP in some mice and vascular pruning in all mice. Nineteen 4-wk-old SM22-rtTA x TetO7-BMPR2R899X mice and 14 SM22-rtTA-only littermate controls were fed doxycycline for 9 wk, systemic pressures were determined by tail cuff, right ventricular pressures were assessed by closed-chest Millar catheter, FMA were performed on some, and the frozen and fixed tissue were collected. We found that 5 of 14 SM22-rtTA x TetO7-BMPR2R899X mice from which RVSP could be obtained had elevated RVSP (Fig. 3A); right ventricular muscularization was consistent with this (P < 0.02 by correlation z-test, data not shown). Counts of <100-µm muscularized vessels per field by immunohistochemistry showed an almost perfect correlation with RVSP (Fig. 3B). Thus there was strong support in both heart and pulmonary vessel muscularization for the pressures measured. Systemic pressures were not different between transgenic mice and controls and did not correlate with RVSP (data not shown).
|
|
|
Muscular lesions are surrounded by macrophages and T cells and include some CD133+ cells. The adventitial component of the vascular lesions contained many CD45+ cells, indicating that they are derived from circulating cells (data not shown). Further staining for T cells and macrophages showed that the vessels were surrounded closely by macrophages (Fig. 6A) with a large number of T cells making up the extended lesion (Fig. 6B). Some of the adventitial cells are also positive for progenitor marker CD34 (Supplemental Fig. 4). Lungs from SM22-rtTA-only mice have only a few alveolar macrophages and no apparent T cells (Fig. 6, A and B, right). Although the muscle making up the central lesion did not stain positive for CD45 (data not shown), there appeared to be a substantial number of cells that were positive for both actin and CD133 (Fig. 6C), suggesting that at least some of the muscle cells differentiated from circulating cells. Vessels that were filled primarily with CD133-positive cells were also CD45 positive (Fig. 6D); from our current study design, we cannot determine whether this is an earlier stage in the development of lesions, such as those in Fig. 6C, or a different type of lesion. In lungs from control mice, only a few isolated alveolar macrophages can be found, and both T cells and CD133 staining were essentially never found.
|
|
Comparison of controls to BMPR2R899X mice with normal RVSP allows examination of pathways dysregulated by mutation and compensation for the mutation without the complicating factor of additional changes caused by elevated pulmonary pressure. In this comparison, we found 181 unique, named genes with altered expression using criteria with a false discovery rate of 27%; our goal was sensitivity over specificity. When sorted into gene ontology groups, we found that the majority fell into four broad categories: developmental genes, stress response, muscle structure and function, and cell cycle (Fig. 8A and Supplemental Table 1).
|
Statistically overrepresented gene ontology groups are likely to have a lower number of falsely discovered genes. For instance,
4% of genes on the array are involved with actin dynamics; with 181 genes altered, one would expect about (0.04 x 181) seven to be actin dynamics-related if they were randomly distributed, with two of these spurious. Since we actually have 22 genes in this category, the category is highly overrepresented and thus likely to be a real change; in addition, we still only expect two of the genes in this category to be falsely discovered.
Comparing genes dysregulated in the BMPR2R899X model to those we have previously reported in our SM22-rtTA x TetO7-BMPR2delx4+ model, we found both substantial overlap as well as major differences, which allows us to make guesses as to domain function (Fig. 8C). The most striking differences are genes dysregulated in BMPR2delx4+ mice but not BMPR2R899X mice. BMPR2R899X mice do not show the broad loss of smooth muscle markers seen in the earlier model (33), suggesting that smooth muscle differentiation state is driven by SMAD activity. More surprisingly, given the apparently inflammatory phenotype seen in BMPR2R899X mice, they do not have the broad increase in cytokines, CXC chemokines, or complement activation seen in the BMPR2delx4+ mice (33), at least at this time point. Conversely, BMPR2R899X mice had several categories of genes changed not seen in BMPR2delx4+ mice, including cell cycle and Wnt pathway genes. Other pathways were dysregulated in both models and had substantial overlap in specific genes altered (marked in Supplemental Table 1).
BMPR2R899X mice with elevated RVSP show additional changes in angiogenesis, vasoreactivity, and injury response.
The BMPR2R899X mice with elevated RVSP had substantial additional changes to gene expression; using the same criteria as above, an additional 269 changed (Supplemental Table 2). These include many genes that are clearly suggestive of response to injury. In Fig. 9A, these include heme oxygenase-1 (Hmox1), an oxidative stress response (28), CD44, a cell adhesion marker characteristic of response to endothelial barrier injury (31), and myotrophin (Mtpn), a pressure and stretch sensor that triggers hypertrophy (30). Angiogenesis factors also increase with elevated RVSP (Fig. 9B), including canonical angiogenesis genes VegfA and angiopoietin-1 but also PDGF receptor-
(PdgfRA), required for fibroblast recruitment in angiogenesis (7). Finally, the advent of high right ventricular pressure results in increases in vasoreactivity-related genes (Fig. 9C) such as endothelin receptor-
(EdnRA), nitric oxide signal transducer guanylate cyclase (Gucy1a3), and thioredoxin interacting protein (Txnip), a peroxisome proliferator-activated receptor-
(PPAR
)-interacting protein that regulates oxidative stress (38) and contractile energy homeostasis (41). For all of these genes, expression in the normal RVSP group is also trending upward, although not yet significantly; we assume that this is because these also are starting to have difficulty regulating pressure. Broadly, additional genes differentially regulated in BMPR2R899X mice with elevated RVSP but not normal RVSP likely indicate genes that are regulated in response to pressure or stretch, not in response to BMPR2 dysfunction.
|
| DISCUSSION |
|---|
|
|
|---|
Gene arrays for SM22-rtTA x TetO7-BMPR2R899X mice assayed before RVSP increases show alterations in expression of genes involved in development, stress response, cell cycle, and muscle structure and function (Fig. 8). Although it is impossible to determine whether the change in any individual gene is etiologic or compensatory, changes in entire classes of genes implicate those pathways in disease development, thereby narrowing the focus for downstream effector pathways for mutant BMPR2. Of equal importance, they also exclude pathways; angiogenesis and conventional vasoreactivity pathways do not show substantial altered expression until after RVSP increases (Fig. 9), suggesting that they are a response to increased pulmonary vascular pressure, not the cause of it.
Integrating the data from this study with our (33, 36) past studies and recent data from the literature (15, 32), we hypothesize that a decrease in BMP signal results in activation of injury repair mechanisms with each downstream BMPR2 effector pathway, SMAD, MAPK, and actin organization, playing a different role. Decreased SMAD signaling through BMPR2 results in both a switch in smooth muscle from a contractile to a synthetic state and increased cytokine expression (12, 33). Increased p38MAPK phosphorylation (Fig. 2D) drives downstream stress responses, possibly including proliferation (Ref. 39; Fig. 8, A and B, and Supplemental). Direct BMPR2-mediated alterations in actin organization through LIMK1 and Tctex may result in multiple effects seen in PAH in patients, including muscle migration (18), matrix deposition (2), and inflammatory cell adhesion (8, 29). Different BMPR2 mutations may cause dysfunction in any or all of these elements, but one can easily imagine how overactive injury response through any of these pathways would predispose to PAH.
These studies were done at a single time point and so cause and effect relationships between these elements are essentially a matter for speculation. However, the differences between those animals that had normal RVSP and those that had high RVSP leads us to the following hypothesis for pathogenesis.
Loss of proper tail domain function leads to alterations in actin organization pathways through LIMK and Tctex. These either directly lead to recruitment of circulating cells through alteration in adhesion factors or they lead to pulmonary vascular stiffening or constriction. In the latter case, this is insufficient to directly increase RVSP, but since shear stress is exquisitely sensitive to vessel diameter (25), even subtle changes in vessel diameter can radically alter shear stress and thus endothelial injury. Some combination of direct effects of BMPR2 mutation, signaling from the circulating cells, and direct endothelial injury result in increased proliferation. This proliferation results in blocked small vessels in the pulmonary vasculature (Fig. 7). When pruning is sufficiently advanced, RVSP starts to increase, leading to conventional muscularization of all pulmonary vessels. The combination of increased pressure and either the Bmpr2 mutation or the presence of inflammatory cells (Fig. 6) leads to the larger lesions; although these happen in secondary forms of human PAH, they are never seen purely as the result of increased pressure in mice. The fact that some of the concentric muscular lesions are CD133-positive suggests that they are least partially populated by recruitment of circulating progenitors.
At the start of this study, our goal was to find out which subset of the phenotype found in the earlier SM22-rtTA x TetO7-BMPR2delx4+ mouse model was due to loss of the tail domain and, by deduction, the loss of the kinase domain. The earlier model had increased RVSP likely due to defects in vasoreactivity and loss of complete smooth muscle differentiation (33, 36, 42). Although it also had defects in cytokine signaling (12, 33), these had no obvious impact on phenotype. The BMPR2delx4+ mouse looked like it was best considered a model of asymptomatic mutation carriers, who may have some baseline defects in vasoreactivity (11) but are in need of a second hit to develop disease.
The much more penetrant phenotype seen in the SM22-rtTA x TetO7-BMPR2R899X mice, which at least nominally have a less severe mutation, was thus somewhat surprising. There are two potential explanations for this. First, it may just be a dose effect; the BMPR2delx4+ mice downregulated their own promoter and had approximately threefold lower expression of transgene on average than BMPR2R899X mice (data not shown), which may have limited their phenotype. The more interesting hypothesis is that loss of proper tail domain function has a stronger phenotype when kinase domain function is preserved; this is suggested by the recent data from the model in which expression of a small interfering RNA to BMPR2 causes defects in pulmonary smooth muscle but not PAH (20).
This model system has several limitations. Although inducible overexpression of a dominantly acting mutant overcomes limitations of BMPR2 knockouts or constitutive mutations, it is not representative of the stoichiometry found in any human patient, which may alter the relative importance of different BMPR2 downstream pathways. It is also possible that the native BMPR2 R899X mutation found in patients is subject to exon skipping or nonsense mediated decay. Furthermore, the doxycycline that we use to induce the transgene itself inhibits matrix metalloproteinases and through these angiogenesis (9); inhibition of angiogenesis exacerbates phenotype in the rat models of PAH (34, 35). On the other hand, doxycycline has been proposed as a treatment for pulmonary diseases such as lymphangioleiomyomatosis (24) and Marfan syndrome (6). We have attempted to control for potential effects by administering doxycycline to control mice as well, but the effect of matrix may be over- or underrepresented in this model. Furthermore, limitation of the transgene to smooth muscle removes potentially important effects in endothelium, circulating cells, or other tissues. This model is thus best understood as one mechanism by which BMPR2 mutation can directly lead to PAH rather than the only mechanism.
In summary, this is the first report of an animal model demonstrating extensive pulmonary vascular structural changes resulting from Bmpr2 mutation in vivo. The power of this model is that it will now be possible to dissect the molecular etiology of the disease in detail with correlation of small molecule inhibitors to alleviation of specific subsets of the molecular, hemodynamic, and histological phenotype. This model should be extremely useful to the research community both in examining early molecular and physical events in the development of PAH and as a platform to validate potential treatments.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. R. Stenmark, B. Meyrick, N. Galie, W. J. Mooi, and I. F. McMurtry Animal models of pulmonary arterial hypertension: the hope for etiological discovery and pharmacological cure Am J Physiol Lung Cell Mol Physiol, December 1, 2009; 297(6): L1013 - L1032. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Kugathasan, J. B. Ray, Y. Deng, E. Rezaei, D. J. Dumont, and D. J. Stewart The angiopietin-1-Tie2 pathway prevents rather than promotes pulmonary arterial hypertension in transgenic mice J. Exp. Med., September 28, 2009; 206(10): 2221 - 2234. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Huang, D. Wang, K. Ihida-Stansbury, P. L. Jones, and J. F. Martin Defective pulmonary vascular remodeling in Smad8 mutant mice Hum. Mol. Genet., August 1, 2009; 18(15): 2791 - 2801. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Majka and D. F. Alvarez Neovascular capacity of endothelial progenitor cells in the adult pulmonary circulation Am J Physiol Lung Cell Mol Physiol, June 1, 2009; 296(6): L868 - L869. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |