Skip to main content

Gene expression analysis of a murine model with pulmonary vascular remodeling compared to end-stage IPAH lungs



Idiopathic pulmonary arterial hypertension (IPAH) continues to be one of the most serious intractable diseases that might start with activation of several triggers representing the genetic susceptibility of a patient. To elucidate what essentially contributes to the onset and progression of IPAH, we investigated factors playing an important role in IPAH by searching discrepant or controversial expression patterns between our murine model and those previously published for human IPAH. We employed the mouse model, which induced muscularization of pulmonary artery leading to hypertension by repeated intratracheal injection of Stachybotrys chartarum, a member of nonpathogenic and ubiquitous fungus in our envelopment.


Microarray assays with ontology and pathway analyses were performed with the lungs of mice. A comparison was made of the expression patterns of biological pathways between our model and those published for IPAH.


Some pathways in our model showed the same expression patterns in IPAH, which included bone morphogenetic protein (BMP) signaling with down-regulation of BMP receptor type 2, activin-like kinase type 1, and endoglin. On the other hand, both Wnt/planar cell polarity (PCP) signaling and its downstream Rho/ROCK signaling were found alone to be activated in IPAH and not in our model.


Activation of Wnt/PCP signaling, in upstream positions of the pathway, found alone in lungs from end stage IPAH may play essential roles in the pathogenesis of the disease.


Pulmonary hypertension is a hemodynamic state characterized by elevation of the mean pulmonary arterial pressure leading to right ventricular (RV) failure and premature death. Pulmonary arterial hypertension (PAH) affects the small muscular arteries and arterioles in the lung and is histologically characterized by endothelial and smooth muscle cell proliferation, medial thickening, and thrombosis in situ. Idiopathic pulmonary arterial hypertension (IPAH), one of 6 subcategories proposed by Dana Point Classification [1], accounts for approximately half of PAH cases [2] and up to 40% of patients with no family history carries mutations in the bone morphogenetic protein receptor type 2 (BMPR2) gene [1]. 7% of patients with IPAH has a family history [3], and about 70% of these have long been recognized and are usually due to mutations in BMPR2, or much less commonly, 2 other members of the transforming growth factor superfamily, activin-like kinase type 1 (ALK1) and endoglin (ENG) [1]. While BMPR2 mutation strongly predisposes to IPAH, only 20% of mutation carriers develop a clinical disease [4]. This finding suggests that the development of IPAH first requires a genetic susceptibility, followed by one or several secondary triggering factors such as modifier genes and some sort of stimulus [5]. However, the pathogenesis of IPAH remains unclear.

To elucidate the pathophysiology of IPAH, we conducted genome-wide analysis of RNA expression profiles in lungs obtained from the murine model, which showed a favorable reproducibility to remodel pulmonary arteries induced by inoculation of Stachybotrys chartarum, an ubiquitous fungus in the surrounding environment [6]. This was followed in the study by exploration for factors playing a significant role in the onset of IPAH by searching discrepant or controversial expression patterns between in the model and those in IPAH previously published.


Detection of S. genes in the human lung tissue of IPAH patients and controls by nested polymerase chain reaction (PCR)

Upon getting the informed consent from patient’s family, a part of lung tissue obtained at autopsy from 9 children (under 18 years of age) with IPAH and 9 age-matched controls for detection of S. DNA in the tissue of human lung. It had been confirmed that none of them had inflammatory changes or other kind of infection in lungs by histopathological examination before DNA preparation. This protocol was approved by ethical committee of Toho University School of Medicine (# 20029). Using DNA isolated from the lung tissue, nested PCR and gel electrophoresis were performed in the usual manner. The following primers were used for identification of S. gene from autopsied lungs [7, 8]: a primary primer of a two-step nested assay, forward primer S-Chart-tri5-F2-5-TACACCAGGGAGGAGCGTGTT-3 and reverse primer S-Chart-tri5-R2-5-GCCGACAATGGTTCGAAGGGA-3 produced a product of 138 bp. For the nested PCR primers, forward primer S-Chart-tri5-F2n-5-GAAAATCTCCAGTATGCCT-3 and reverse primer S-Chart-tri5-R2n-5-AGCCTCCAGTCTCTTGGGG-3 produced a product of 96 bp. For β-actin, the control for DNA quality, a two-step nested assay was also used with the following primers: forward primer B3-5-ACACAACTGTGTTCACTAGC-3 and reverse primer B4-5-CAACTTCATCCACGTTCACC-3 produced a product of 110 bp. For the nested PCR primers, forward primer b-globin-nestF34-52-5-AACCTCAAACAGACACCAT-3 and reverse primer b-globin-nestR103-85-5-TTGCCCCACAGGGCAGTAA-3 produced a product of 70 bp. The PCR amplification was performed as follows: initial denaturation at 94°C for 2 min, followed by 35 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, with an extension at 72°C for 30 s, and a final extension at 72°C for 7 min. The second-round PCR reactions were performed in a manner identical to that applied for the first strand PCR, except for using different sets of primers. The PCR products were analyzed by electrophoresis on an agarose gel stained with ethidium bromide upon preparation.

Fungal preparation and intratracheal injection

S. chartarum (IFM 53637), which produces trichothecene mycotoxins, was isolated from house dust in Japan, and has been stored in the culture collection of the Medical Mycology Research Center, Chiba University. The fungus was grown on potato dextrose agar (PDA) slants for 3 weeks at 25°C. Spores were collected in RPMI1640 medium (Sigma, St. Louis, MO, USA) and the concentration was adjusted to 4 × 105 spores / ml. Spore concentrations and appearance of the suspension were evaluated under light microscopy before use. Six-week-old male ddY mice (Tokyo Laboratory Animals Science, Tokyo, Japan) were employed in this study. Mice were lightly anaesthetized with an intraperitoneal injection of ketamine (65 mg/kg BW) and xylazine (13 mg/kg BW). Their mean weight was 27.4 ± 1.21 g. The mice were placed in a supine position and a 24 G intravascular catheter (Insyte-W; Becton-Dickinson, Sandy, UT, USA) was then inserted intratracheally. The spore suspension (25 μl / mouse) containing 1 × 104 spores was injected through the catheter into the trachea of each mouse 12 times at 4–5 day intervals for 8 weeks (n = 3) as described previously [6]. Control mice (n = 3) were injected with the same volume of RPMI-1640 medium rather than the spore suspension. All mice were cared for in accordance with the rules and regulations set out by the Prime Minister’s Office of Japan. Animal protocols were approved by the Special Committee on Animal Welfare of Chiba University. (DOU: 21–65).

Histopathology and morphometric analysis of pulmonary arteries

Mice were sacrificed using by an overdose of diethyl ether inhalation 7 days after the last injection. Lungs were removed and fixed with a 10% formaldehyde solution, embedded in paraffin, cut into 3 μm-thick sections, and stained with hematoxylin and eosin for histopathological examination. Elastic fiber was stained with Elastica-van Gieson staining (Muto pure chemicals, Tokyo, Japan). Morphometric analyses were performed to determine the luminal stenosis of the pulmonary arteries. At least 200 pulmonary arteries per lung section from each mouse were chosen at random and examined. Cross-sections of arteries observed in the section were used to measure the distance between the external elastic lamina, internal elastic lamina, and intravascular lumen. All images were analyzed using IMAGE J 1.36b software (National Institutes of Health, Bethesda, MD, USA). The stricture rate was therefore calculated. Morphometry measurements were performed according to the techniques described in Dail and Hammar’s Pulmonary Pathology [9]. The thickness of media was calculated by subtracting the distance between the internal elastic lamina from that of the external lamina, and the thickness of the intima was calculated by subtracting the distance between the intravascular lumen from that of the internal elastic lamina. The distance between the external elastic lamina of the artery was defined as the diameter. Arteries were divided into three groups according to diameter: 50 < μm, 50–100 μm, and &100 μm. Data are given as mean ± SD. Statistical analyses were performed using Mann-Whitney's U test. Differences were considered significant at P < 0.05.

RNA isolation and quality identification

RNA was isolated from the whole lung homogenates for both microarray analysis and Real-time (RT) Quantitative PCR with the RNeasy Lipid Tissue Mini Kit (Qiagen, Alameda, CA, USA) according to the manufacturer’s instructions and stored at −80°C. Total RNA quality was assessed and confirmed using the Agilent Bioanalyzer 2100 (Agilent Technologies, Palo Alto, CA, USA) for visualization of the 28S and 18S rRNA bands. RNA concentration and purity were also assessed and confirmed using the UV spectrophotometer NanoDrop™ND-1000 (NanoDrop Technologies, Wilmington, DE, USA), which calculates 260/280 ratios.

RNA preparation for microarray analysis

cDNA preparation and microarray analysis were conducted at Bio Matrix Research (Chiba, Japan) using the Affymetrix system (Santa Clara, CA, USA). Isolated total RNA (100 μg) was converted into double-stranded cDNA using 3IVT Express kit (Affymetrix, Santa Clara, CA, USA), which was purified using a GeneChip Sample Cleanup Module (Affymetrix, Santa Clara, CA, USA). In vitro transcription reactions were performed using a GeneChip IVT Labeling Kit, which includes T7 RNA polymerase and biotin-labeled ribonucleotides. Biotin-labeled cRNA was purified using a GeneChip Sample Cleanup Module. The concentration of cRNA was calculated from light absorbance at 260 nm using a UV spectrophotometer. cRNA (15 mg) was then fragmented at 94°C in the presence of a fragmentation buffer (Affymetrix, Santa Clara, CA, USA). The labeled cRNA was purified, fragmented, and spiked with in vitro transcription controls.

Microarray analysis

Mouse Genome 430 2.0 microarrays (Affymetrix, Santa Clara, CA, USA) were hybridized with 12.5 μg of cRNA. The array was incubated for 16 hr at 45°C, and automatically washed and stained with the GeneChip Hybridization, Wash and Stain Kit (Affymetrix, Santa Clara, CA, USA) on an Affymetrix GeneChip Fluidics station. The arrays were analyzed using the GeneChip Scanner 3000. All preparations were run on quality-controlled chips and had 3/5 signal ratios of less than 3. The expression value of the transcript was computed using Affymetrix® GeneChip® Command Console® Software (AGCC) with the MAS5 algorithm [10], in which the probabilities of the values of each transcript were indicated as the “Flag” Present (p 0 to <0.04), Marginal (p 0.04 to <0.06), and Absent (p 0.06 to <0.5). Further analysis was performed with probes that had a present call in all analyzed samples. For analysis, the data were normalized using GeneSpring® GX 10.0 (Agilent Technologies, Palo Alto, CA, USA) data-mining software, per-chip normalization to the 50th percentile of the measurements for the array, and per-gene by normalizing to the median measurement for the gene across all the arrays in the data set. In addition, fold changes were calculated by this software for each gene between the experimental groups and controls. Statistically significant differences were investigated by means of unpaired t-tests. Gene expression differences with p < 0.05 and at least a ±1.3-fold change were considered statistically significant. We employed the Biological Networks Gene Ontology (GO) tool BINGO [11] ( to find statistically over- or under-represented GO categories in biologic data as the tool for GO analysis of the stored genes. The analysis was done using the ‘hyper geometric test’, and all GO biological process terms that were significant with P < 0.05 (after correcting for multiple term testing using Benjamin and Hochberg false discovery rate corrections) were selected as over-represented and under-represented. Furthermore, the open access and curated pathway database REACTOME [12] ( was used to determine which events (reactions and/or pathways) were statistically overrepresented in a set of genes.

RT Quantitative PCR for evaluating the microarray results

cDNA preparation and RT-PCR were performed at Bio Matrix Research (Chiba, Japan). 2 μg of isolated RNA was used for cDNA synthesis (40 μl) using a Superscript III First Strand Synthesis System (Invitrogen, Carlsbad, CA, USA). To evaluate the concentration and purity of cDNA, 260/280 ratios were calculated using the UV spectrophotometer NanoDrop™ND-1000 (NanoDrop Technologies, Wilmington, DE, USA). PCR was performed in a 15 μl reaction mixture containing 1 μl of sample cDNA, 0.75 μl of TaqMan® Gene Expression Assays, 7.50 μl of TaqMan Universal PCR Master Mix, and 5.75 μl of RNase / DNase free water (Applied Biosystems, Foster City, CA, USA). 15 μl of PCR reaction mix was transferred into a 384-well reaction plate. The primer sequences are given in Table 1. Gene expression was measured on a 7900HT Fast RT-PCR system (Applied Biosystems, Foster City, CA, USA) with cycle conditions of 50°C/2 min, 95°C/10 min, 95°C/15 sec, and 60°C/1 min (steps 3–4 were repeated 40 times). Assay results were collected and analyzed using SDS 2.2 software (Applied Biosystems, Foster City, CA, USA). Each value is the mean of three biological replicates. Data are given as mean ± SD. Statistical analyses were performed by non-paired t-tests. Differences were considered significant at P < 0.05.

Table 1 Primers used in RT-PCR

Literature search concerning gene expression pattern in IPAH

We used PubMed ( to search for previous studies published since 2000 that analyzed biological molecules with altered expression using lung samples isolated from IPAH, compared with normal control or secondary PAH. From microarray studies, we referred to open access data stored in Gene Expression Omnibus: GEO ( and identified the genes that were differentially expressed following the methods described in each study. GO and Pathway analysis were performed on the genes in the same way as we done on our data.

Comparing expression patterns of molecules between IPAH and experimental model

The biological molecules reported in various studies were collected and divided into groups according to the gene ontology biological process and pathways using the Gene Ontology Annotation Database GOA [13] (, REACTOME [12], and KEGG [14] ( The expression pattern of each group in IPAH was compared with our PAH models.


Detection of S. genes in human lung tissue from patients with IPAH and controls by nested PCR

S. chartarum DNA was detected in 6 of 9 lung samples among both two groups, children with IPAH and age-matched controls (Figure 1). There was no difference in the frequency (approximately 70%) of the detection in children, with or without of IPAH.

Figure 1

Detection of Stachybotrys chartarum gene in the human lung tissue of IPAH patients and controls. Agarose gel electrophoresis showed the amplification products after nested polymerase chain reaction for S. chartarum DNA using the standard assay. P1- P2: Positive control (96 bp and 70 bp); S1-S9: Children with IPAH; C1-C9: Age-matched control; N: Negative control (No DNA). S. chartarum and β-globin bands are 96 and 70 bp, respectively. 96 bp–specific band was detectable in P1-2, S1-2, S4, S6, S8-9, C1, C3-4, C6, C8-9.

Pathological findings in experimental PAH

Diffuse symmetric thickening of intima and media in the pulmonary artery was shown in the experimental group (Figure 2). The thickened intima and media were accompanied by proliferation of myointimal and smooth muscle cells, respectively. None of arteries showed alterations corresponding to necrosis, thrombosis, and plexiform lesions. The changes developed in arteries of small and medium-sized were mostly uniform. The thickened intima and media were statistically significant regardless of the size of vessels (Figure 3).

Figure 2

Pulmonary vascular remodeling in Stachybotrys chartarum -exposed mice. Peripheral pulmonary arteries in experimental group (A, C, and E) and control group (B, D, and F). A and B: hematoxylin and eosin double stain, C, D, E, and F: Elastic Van Gieson stain. (A) Pulmonary vascular thickening was present over the whole lung. (C) The arterial wall shows symmetrical thickening of the intima and media. (E) This vascular remodeling almost totally obliterates the lumens. (B, D, and F) Structure of the lung from control group is unaltered. Scale bars: 500 μm (A, B), 250 μm (C, D), and 50 μm (E, F).

Figure 3

Morphometric analysis evaluating luminal stenosis of pulmonary arteries. Significant thickening of both intima and media were observed in the exposure group when compared with the control regardless of the size of vessels. Data are presented as means with standard error of the means. Statistical analyses were performed using Mann-Whitney's U test. A bar with ** is significantly different (P < 0.01).

In this model, no venous canals were altered. None were found of components from injected fungus (spores and hyphae) and changes likely induced by the injection, such as perivascular cuffing and intraalveolar inflammatory exudates.

Gene expression in the PAH model mouse

Upon normalizing the expression values for the samples, the scatter plot of log intensity values was obtained as shown in Figure 4. 337 and 503 genes were found to exhibit up- and down-regulation from the lungs of mice with inoculation of the fungus in comparison to those from the control group. The most markedly up- and down- regulated genes are shown in Table 2. Down regulation of BMPR2, ALK1, and ENG was found in our PAH model. Down regulation of SMAD family member 6 expression was also observed, while the expressions of other genes involved in bone morphogenetic protein (BMP) signaling were unchanged. All the microarray data are MIAME compliant and the complete microarray data were deposited in GEO (accession number GSE23178). 696 biological process terms were detected by GO analysis in up-regulated genes. GO terms related to the immune system and cytokines accounted for about 80% of total terms. Among the remaining terms, estrogen receptor signaling pathway and serotonin transport/secretion were included. A statistically significant biological process was not found for down-regulated genes. Pathway analysis revealed that reactions and pathways related to the immune system, Janus kinase/signal transducers and activators of the transcription (JAK/STAT) pathway, and hemostasis etc. were detected (Table 3) in up-regulated genes. Additionally, vascular endothelial growth factor (VEGF), platelet-derived growth factor (PDGF), apoptosis, BMP signaling, etc. were detected in down-regulated genes.

Figure 4

Scatter plot of microarray dataset. Gene expression signal of control group (Y-axis) plotted against experimental group (X-axis). The expression data are the average of three independent gene chips for each group.

Table 2 The top 10 most up-regulated and down-regulated genes
Table 3 Events (reactions and/or pathways) associated with up-regulated and down-regulated genes

Validation of microarray results by RT quantitative PCR analysis

RT-PCR was used to validate 14 selected genes that were induced or suppressed by the exposure. The correlation of fold changes in gene expression between the arrays and PCR is shown in Table 4. The results demonstrate completely the same gene expression pattern between both methods. The alterations of gene expressions were statistically significant in BMPR2, ENG, Vascular endothelial growth factor A, Platelet-derived growth factor alpha polypeptide, matrix metallopeptidase (MMP) 19, MMP12, eosinophil-associated ribonuclease A family member 11, and chemokine ligand 9.

Table 4 Quantitative PCR Validation of Microarray Expression data

Altered expression of biological molecules and genes in lungs from IPAH

The biological molecules reported in previous studies [1528] were listed in Additional file 1. Microarray data of patients with IPAH and normal controls were referable from 3 individual studies [15, 16, 27]. The numbers of IPAH investigated in each study were 2 [16], 7 [27], and 18 [15], the mean age and its standard deviation were 44 ± 10 [16], 29 ± 16 [27], and 44 ± 18 [15], respectively. We analyzed the microarray data in each study and tried to extract the genes that showed common expression patterns through these three studies. However, there were few genes, and none of significant GO or pathways was identified among the previous reports regarding to IPAH.

Comparing expression patterns of molecules between IPAH and experimental model

Events (reactions and/or pathways) and the expression patterns in IPAH extracted from the previous reports are listed with comparison to those resulted from our experimental model [1528] (Table 5).

Table 5 Events (reactions and/or pathways) and the expression patterns in IPAH and in our PAH model


Since PAH is a progressive disease of unknown cause involving pulmonary arterial remodeling, characterized by relentless deterioration and death, intense investigations have been conducted in a variety of animal models [29] to know pathophysiology. The most commonly used were rats exposed to either hypoxia or monocrotaline, and newer models were introduced that involved modification of these approaches using rodents including transgenic mice [29, 30].There were at least three genomewide studies conducting rat models among them, but little have been discussed with comparison to those in the human disease with pathway and GO analyses [3133]. We have therefore aimed to elucidate a part of pathophysiology of PAH accompanied by pulmonary arterial remodeling with comparison in gene expression pattern between those previously known in end-stage IPAH and our murine model, of which muscularization in media and intima of pulmonary arteries was induced by inoculation of nonpathogenic fungus [6, 34]. It was found that a large frequency of S. chartarum gene in the lung of both children with IPAH and age-matched controls in autopsy cases, whereas the prior histological examination had revealed no inflammatory changes with an association to fungal infection. The result suggests that the airway of human generally exposed by the ubiquitous fungus. Accordingly, unknown intrinsic factors may play a significant role in the onset of IPAH.

It was explored in the study to search the pathways which show common expression pattern among IPAH patients using microarray data previously published as well as up-loaded at the open access sites [15, 16, 27], but none were elucidated. A part of the result may reflect clinical differences in subjects such as age, gender, stage of the disease, and etc. On the other hand, the present murine model showed similar expression pattern of genes which has been generally accepted as some biological molecules associated with the pathogenesis of IPAH, such as BMP signaling [1, 35, 36]: BMPR2, ALK1, and ENG. Although some of previous genome-wide analyses [3133] conducting rat models reported alterations of signaling pathway (Table 6), the present model revealed more favorable similarity in number of pathways expressing same manner to IPAH. However, the result from our study might be affected by recruited cellular components and/or reformed extracellular matrix, because RNA was extracted from the whole lung of the mice. Although RNA extraction was done one week after the latest inoculation of the fungus when no inflammatory change had been previously confirmed [6], in vitro work evaluating cells from the lungs of the mice would be required to look at functional alterations. In addition, whereas multiple time point’s evaluation could highlight the difference in pathophysiology of models of between our mice and rats previously published, we have evaluated once at the point when the most significant muscularization had been confirmed [6], since a priority in the study would have been set in comparing gene expression patterns between end-stage IPAH and our model. It has been generally accepted that the altered expression of BMP signaling was one of the important molecular reactions when pulmonary vascular remodeling developed as the response to some sort of stimulus [2, 3]. This can be supported by the facts that only a few cases of BMPR2 mutation carriers developed clinical disease [4] and BMPR2-knockdown mice did not develop pulmonary artery medial hypertrophy, spontaneously [29]. Since remodeling of pulmonary artery in our model may be a sequel to inoculation of nonpathogenic fungus, down regulation of BMPR2 signaling should be simply understood as a consequence of muscularization which might be induced by alteration of signaling pathways at the upper stream. Besides BMP signaling, four more pathways known as identical expression patterns in IPAH were found in our PAH model. It has been suggested that inflammation might participate in the onset and propagation of pulmonary vascular remodeling in PAH via the JAK/STAT pathway [37, 38], because elevated levels of inflammatory cytokines could trigger inflammation that is characteristic of PAH of both connective tissue disease-associated [39] and Virus-associated [40]. While PAH is up to threefold more prevalent in women than men [41], the increased expression of molecules associated with the estrogen signaling pathway are reported in both genders of patients with IPAH [15]. In addition, the evidence implicating serotonin has been discussed with correlation to the anorexigenic drugs aminorex and fenfluramine [42]. As with alteration of BMP signaling, the altered expression of these pathways would be simply associated with the consequence of pulmonary vascular remodeling developed as the response to some sort of stimulus. It might be better to understand as a sequel to vascular remodeling because that coagulation activity could be activated by various stimuli such as a cytokine and endothelial dysfunction [43].

Table 6 Summarized events (reactions and/or pathways) and the expression patterns in previously rat models

On the other hand, it emerged that some discrepant gene expression patterns between those previously known in IPAH and our model. Four pathways were identified as those altered alone in IPAH which comprised up-regulations of the Wnt/planar cell polarity (PCP) signaling pathway [27], the Ras homolog (Rho)/Rho- Associated Coiled-Coil-Forming Protein Kinase (ROCK) pathway [26, 27], and the hypoxia response pathway [23, 24], and down-regulations of the transforming growth factor beta (TGFB) signaling pathway [16]. Both PDGF signaling pathway [25] and VEGF signaling pathway [15, 26] were found as those up-regulated in IPAH and down-regulated in this model. Interestingly, not detected was a pathway showing the reverse pattern.

Among these various pathways, we would focus on the Wnt/PCP pathway as the essential pathway of pathogenesis in IPAH because this pathway is presently known to be in the upper levels of hierarchical pathways regulating other related pathways [1214]. Laumanns et al. reported that microarray analysis of lung tissue from patients with IPAH demonstrated the contribution of this pathway to the pathogenesis of IPAH [27]. It has been reported Wnt family of signaling proteins is essential for organ development in general, and lung morphogenesis in particular [44]. Especially, the PCP pathway signals through activation of the Rho/ROCK signaling pathway are implicated in cytoskeletal organization and epithelial cell polarity. In addition, this pathway has been shown to be required for normal lung development, and reports are beginning to emerge of links between PCP pathways and lung disease [45]. Besides IPAH, PCP gene expression changes were observed within isolated pulmonary vasculature in patients with pulmonary fibrosis. Studies of patients with IPAH have shown significant up-regulation of PCP signaling [27] and down-regulation of TGFB signaling [16]. Investigations have also shown inhibition of canonical Wnt signaling by PCP ligands [45], crosstalk between canonical Wnt signaling and TGFB signaling [46], and recruitment of both canonical and non-canonical Wnt pathways are required in BMP2 mediated angiogenesis in human pulmonary artery endothelial cells [47]. A part of results from the present study, no change of pathways around this system can support that altered interaction between the canonical Wnt pathway and the TGFB signaling pathway play an essential role for human disease and further clarification is expected of the role of the PCP pathway in the pathogenesis of IPAH. On the other hand, some previous reports conducted human disease reported the Rho/ROCK signaling pathway is also activated, but it is known that many other stimuli can activate this system [48]. An activation of this pathway may conduct to smooth muscle cell proliferation in pulmonary artery [49], we wish to understand this event as secondary episode to alteration in PCP signaling pathways.

Although, the details in the mechanism is unclear, various factors (e.g., Hypoxia, inflammation, and shear stress) are known as the inducer of growth factors such as VEGF and PDGF [50, 51]. Accordingly, their activation may not be essential because they usually play at the lower level of hierarchical pathways, which are activated by many sorts of stimuli.


Discrepancy in gene expression pattern between this model and the human disease previously reported suggests that activation of Rho/ROCK signaling via Wnt/PCP signaling plays an essential role in pathogenesis of IPAH.


  1. 1.

    Simonneau G, Robbins IM, Beghetti M, Channick RN, Delcroix M, Denton CP, Elliott CG, Gaine SP, Gladwin MT, Jing ZC, Krowka MJ, Langleben D, Nakanishi N, Souza R: Updated clinical classification of pulmonary hypertension. J Am Coll Cardiol. 2009, 54: S43-S54. 10.1016/j.jacc.2009.04.012.

    PubMed  Article  Google Scholar 

  2. 2.

    Peacock AJ, Murphy NF, McMurray JJ, Caballero L, Stewart S: An epidemiological study of pulmonary arterial hypertension. Eur Respir J. 2007, 30: 104-109. 10.1183/09031936.00092306.

    PubMed  CAS  Article  Google Scholar 

  3. 3.

    Sankelo M, Flanagan JA, Machado R, Harrison R, Rudarakanchana N, Morrell N, Dixon M, Halme M, Puolijoki H, Kere J, Elomaa O, Kupari M, Räisänen-Sokolowski A, Trembath RC, Laitinen T: BMPR2 mutations have short lifetime expectancy in primary pulmonary hypertension. Hum Mutat. 2005, 26: 119-124. 10.1002/humu.20200.

    PubMed  CAS  Article  Google Scholar 

  4. 4.

    Newman JH, Trembath RC, Morse JA, Grunig E, Loyd JE, Adnot S, Coccolo F, Ventura C, Phillips JA, Knowles JA, Janssen B, Eickelberg O, Eddahibi S, Herve P, Nichols WC, Elliott G: Genetic basis of pulmonary arterial hypertension: current understanding and future directions. J Am Coll Cardiol. 2004, 43: 33S-39S. 10.1016/j.jacc.2004.02.028.

    PubMed  CAS  Article  Google Scholar 

  5. 5.

    Firth AL, Mandel J, Yuan JX: Idiopathic pulmonary arterial hypertension. Dis Model Mech. 2010, 3: 268-273. 10.1242/dmm.003616.

    PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Ochiai E, Kamei K, Watanabe A, Nagayoshi M, Tada Y, Nagaoka T, Sato K, Sato A, Shibuya K: Inhalation of Stachybotrys chartarum causes pulmonary arterial hypertension in mice. Int J Exp Pathol. 2008, 89: 201-208. 10.1111/j.1365-2613.2008.00585.x.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  7. 7.

    Black JA, Dean TR, Foarde K, Menetrez M: Detection of Stachybotrys chartarum using rRNA, tri5, and beta-tubulin primers and determining their relative copy number by real-time PCR. Mycol Res. 2008, 112: 845-851. 10.1016/j.mycres.2008.01.006.

    PubMed  CAS  Article  Google Scholar 

  8. 8.

    Black JA, Foarde KK: Comparison of four different methods for extraction of Stachybotrys chartarum spore DNA and verification by real-time PCR. J Microbiol Methods. 2007, 70: 75-81. 10.1016/j.mimet.2007.03.017.

    PubMed  CAS  Article  Google Scholar 

  9. 9.

    Groshong SD, Tomashefski JF, Cool CD: Pulmonary Hypertensive vascular Disease. Dail and Hammar’s pulmonary pathology. Volume 1. Edited by: Tomashefski JF Jr. 2008, New York: Springer, 1032-1087. 3

    Google Scholar 

  10. 10.

    Pepper SD, Saunders EK, Edwards LE, Wilson CL, Miller CJ: The utility of MAS5 expression summary and detection call algorithms. BMC Bioinforma. 2007, 8: 273-10.1186/1471-2105-8-273.

    Article  Google Scholar 

  11. 11.

    Maere S, Heymans K, Kuiper M: BiNGO: a Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21: 3448-3449. 10.1093/bioinformatics/bti551.

    PubMed  CAS  Article  Google Scholar 

  12. 12.

    Croft D, O'Kelly G, Wu G, Haw R, Gillespie M, Matthews L, Caudy M, Garapati P, Gopinath G, Jassal B, Jupe S, Kalatskaya I, Mahajan S, May B, Ndegwa N, Schmidt E, Shamovsky V, Yung C, Birney E, Hermjakob H, D'Eustachio P, Stein L: Reactome: a database of reactions, pathways and biological processes. Nucleic Acids Res. 2011, 39: D691-D697. 10.1093/nar/gkq1018.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  13. 13.

    Barrell D, Dimmer E, Huntley RP, Binns D, O'Donovan C, Apweiler R: The GOA database in 2009-an integrated gene ontology annotation resource. Nucleic Acids Res. 2009, 37: D396-D403. 10.1093/nar/gkn803.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  14. 14.

    Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M: KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010, 38: D355-D360. 10.1093/nar/gkp896.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  15. 15.

    Rajkumar R, Konishi K, Richards TJ, Ishizawar DC, Wiechert AC, Kaminski N, Ahmad F: Genomewide RNA expression profiling in lung identifies distinct signatures in idiopathic pulmonary arterial hypertension and secondary pulmonary hypertension. Am J Physiol Heart Circ Physiol. 2010, 298: H1235-H1248. 10.1152/ajpheart.00254.2009.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  16. 16.

    Fantozzi I, Huang W, Zhang J, Zhang S, Platoshyn O, Remillard CV, Thistlethwaite PA, Yuan JX: Divergent effects of BMP-2 on gene expression in pulmonary artery smooth muscle cells from normal subjects and patients with idiopathic pulmonary arterial hypertension. Exp Lung Res. 2005, 31: 783-806. 10.1080/01902140500461026.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  17. 17.

    Sanchez O, Marcos E, Perros F, Fadel E, Tu L, Humbert M, Dartevelle P, Simonneau G, Adnot S, Eddahibi S: Role of endothelium-derived CC Chemokine ligand 2 in idiopathic pulmonary arterial hypertension. Am J Respir Crit Care Med. 2007, 176: 1041-1047. 10.1164/rccm.200610-1559OC.

    PubMed  CAS  Article  Google Scholar 

  18. 18.

    Masri FA, Xu W, Comhair SA, Asosingh K, Koo M, Vasanji A, Drazba J, Anand-Apte B, Erzurum SC: Hyperproliferative apoptosis-resistant endothelial cells in idiopathic pulmonary arterial hypertension. Am J Physiol Lung Cell Mol Physiol. 2007, 293: L548-L554. 10.1152/ajplung.00428.2006.

    PubMed  CAS  Article  Google Scholar 

  19. 19.

    Hecker M, Zaslona Z, Kwapiszewska G, Niess G, Zakrzewicz A, Hergenreider E, Wilhelm J, Marsh LM, Sedding D, Klepetko W, Lohmeyer J, Dimmeler S, Seeger W, Weissmann N, Schermuly RT, Kneidinger N, Eickelberg O, Morty RE: Dysregulation of the IL-13 receptor system a novel pathomechanism in pulmonary arterial hypertension. Am J Respir Crit Care Med. 2010, 182: 805-818. 10.1164/rccm.200909-1367OC.

    PubMed  CAS  Article  Google Scholar 

  20. 20.

    Courboulin A, Paulin R, Giguère NJ, Saksouk N, Perreault T, Meloche J, Paquet ER, Biardel S, Provencher S, Côté J, Simard MJ, Bonnet S: Role for miR-204 in human pulmonary arterial hypertension. J. Exp. Med. 2011, 208: 535-548. 10.1084/jem.20101812.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  21. 21.

    Dewachter L, Adnot S, Fadel E, Humbert M, Maitre B, Barlier-Mur AM, Simonneau G, Hamon M, Naeije R, Eddahibi S: Angiopoietin/Tie2 pathway influences smooth muscle hyperplasia in idiopathic pulmonary hypertension. Am J Respir Crit Care Med. 2006, 174: 1025-1033. 10.1164/rccm.200602-304OC.

    PubMed  CAS  Article  Google Scholar 

  22. 22.

    Yu Y, Keller SH, Remillard CV, Safrina O, Nicholson A, Zhang SL, Jiang W, Vangala N, Landsberg JW, Wang JY, Thistlethwaite PA, Channick RN, Robbins IM, Loyd JE, Ghofrani HA, Grimminger F, Schermuly RT, Cahalan MD, Rubin LJ, Yuan JX: A functional single-nucleotide polymorphism in the TRPC6 gene promoter associated with idiopathic pulmonary arterial hypertension. Circulation. 2009, 119: 2313-2322. 10.1161/CIRCULATIONAHA.108.782458.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  23. 23.

    Lepetit H, Eddahibi S, Fadel E, Frisdal E, Munaut C, Noel A, Humbert M, Adnot S, D'Ortho MP, Lafuma C: Smooth muscle cell matrix metalloproteinases in idiopathic pulmonary arterial hypertension. Eur Respir J. 2005, 25: 834-842. 10.1183/09031936.05.00072504.

    PubMed  CAS  Article  Google Scholar 

  24. 24.

    Eddahibi S, Humbert M, Fadel E, Raffestin B, Darmon M, Capron F, Simonneau G, Dartevelle P, Hamon M, Adnot S: Serotonin transporter overexpression is responsible for pulmonary artery smooth muscle hyperplasia in primary pulmonary hypertension. J Clin Invest. 2001, 108: 1141-1150.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  25. 25.

    Perros F, Montani D, Dorfmüller P, Durand-Gasselin I, Tcherakian C, Le Pavec J, Mazmanian M, Fadel E, Mussot S, Mercier O, Hervé P, Emilie D, Eddahibi S, Simonneau G, Souza R, Humbert M: Platelet-derived growth factor expression and function in idiopathic pulmonary arterial hypertension. Am J Respir Crit Care Med. 2008, 178: 81-88. 10.1164/rccm.200707-1037OC.

    PubMed  CAS  Article  Google Scholar 

  26. 26.

    Do e Z, Fukumoto Y, Takaki A, Tawara S, Ohashi J, Nakano M, Tada T, Saji K, Sugimura K, Fujita H, Hoshikawa Y, Nawata J, Kondo T, Shimokawa H: Evidence for Rho-kinase activation in patients with pulmonary arterial hypertension. Circ J. 2009, 73: 1731-1739. 10.1253/circj.CJ-09-0135.

    PubMed  Article  Google Scholar 

  27. 27.

    Laumanns IP, Fink L, Wilhelm J, Wolff JC, Mitnacht-Kraus R, Graef-Hoechst S, Stein MM, Bohle RM, Klepetko W, Hoda MA, Schermuly RT, Grimminger F, Seeger W, Voswinckel R: The noncanonical WNT pathway is operative in idiopathic pulmonary arterial hypertension. Am J Respir Cell Mol Biol. 2009, 40: 683-691. 10.1165/rcmb.2008-0153OC.

    PubMed  CAS  Article  Google Scholar 

  28. 28.

    Golpon HA, Geraci MW, Moore MD, Miller HL, Miller GJ, Tuder RM, Voelkel NF: HOX genes in human lung: altered expression in primary pulmonary hypertension and emphysema. Am J Pathol. 2001, 158: 955-966. 10.1016/S0002-9440(10)64042-4.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  29. 29.

    Stenmark KR, Meyrick B, Galie N, Mooi WJ, McMurtry IF: Animal models of pulmonary arterial hypertension: the hope for etiological discovery and pharmacological cure. Am J Physiol Lung Cell Mol Physiol. 2009, 297: L1013-L1032. 10.1152/ajplung.00217.2009.

    PubMed  CAS  Article  Google Scholar 

  30. 30.

    Rabinovitch M: Cellular and molecular pathobiology of pulmonary hypertension conference summary. Chest. 2005, 128: 642S-646S. 10.1378/chest.128.6_suppl.642S.

    PubMed  Article  Google Scholar 

  31. 31.

    Moreno-Vinasco L, Gomberg-Maitland M, Maitland ML, Desai AA, Singleton PA, Sammani S, Sam L, Liu Y, Husain AN, Lang RM, Ratain MJ, Lussier YA, Garcia JG: Genomic assessment of a multikinase inhibitor, sorafenib, in a rodent model of pulmonary hypertension. Physiol Genomics. 2008, 33: 278-291. 10.1152/physiolgenomics.00169.2007.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  32. 32.

    Kwapiszewska G, Wilhelm J, Wolff S, Laumanns I, Koenig IR, Ziegler A, Seeger W, Bohle RM, Weissmann N, Fink L: Expression profiling of laser-microdissected intrapulmonary arteries in hypoxia-induced pulmonary hypertension. Respir Res. 2005, 6: 109-10.1186/1465-9921-6-109.

    PubMed  PubMed Central  Article  Google Scholar 

  33. 33.

    Vaszar LT, Nishimura T, Storey JD, Zhao G, Qiu D, Faul JL, Pearl RG, Kao PN: Longitudinal transcriptional analysis of developing neointimal vascular occlusion and pulmonary hypertension in rats. Physiol Genomics. 2004, 17: 150-156. 10.1152/physiolgenomics.00198.2003.

    PubMed  CAS  Article  Google Scholar 

  34. 34.

    Pestka JJ, Yike I, Dearborn DG, Ward MD, Harkema JR: Stachybotrys chartarum, trichothecene mycotoxins, and damp building-related illness: new insights into a public health enigma. Toxicol Sci. 2008, 104: 4-26. 10.1093/toxsci/kfm284.

    PubMed  CAS  Article  Google Scholar 

  35. 35.

    Chida A, Shintani M, Nakayama T, Furutani Y, Hayama E, Inai K, Saji T, Nonoyama S, Nakanishi T: Missense mutations of the BMPR1B (ALK6) gene in childhood idiopathic pulmonary arterial hypertension. Circ J. 2012, 76: 1501-1508. 10.1253/circj.CJ-11-1281.

    PubMed  CAS  Article  Google Scholar 

  36. 36.

    Shintani M, Yagi H, Nakayama T, Saji T, Matsuoka R: A new nonsense mutation of SMAD8 associated with pulmonary arterial hypertension. J Med Genet. 2009, 46: 331-337. 10.1136/jmg.2008.062703.

    PubMed  CAS  Article  Google Scholar 

  37. 37.

    Price LC, Wort SJ, Perros F, Dorfmüller P, Huertas A, Montani D, Cohen-Kaminsky S, Humbert M: Inflammation in pulmonary arterial hypertension. Chest. 2012, 141: 210-221. 10.1378/chest.11-0793.

    PubMed  CAS  Article  Google Scholar 

  38. 38.

    Pfeifer AC, Timmer J, Klingmüller U: Systems biology of JAK/STAT signalling. Essays Biochem. 2008, 45: 109-120. 10.1042/BSE0450109.

    PubMed  CAS  Article  Google Scholar 

  39. 39.

    Nicolls MR, Taraseviciene-Stewart L, Rai PR, Badesch DB, Voelkel NF: Autoimmunity and pulmonary hypertension: a perspective. Eur Respir J. 2005, 26: 1110-1118. 10.1183/09031936.05.00045705.

    PubMed  CAS  Article  Google Scholar 

  40. 40.

    Cicalini S, Chinello P, Petrosillo N: HIV infection and pulmonary arterial hypertension. Expert Rev Respir Med. 2011, 5: 257-266. 10.1586/ers.11.10.

    PubMed  Article  Google Scholar 

  41. 41.

    Thenappan T, Shah SJ, Rich S, Gomberg-Maitland M: A USA-based registry for pulmonary arterial hypertension: 1982–2006. Eur Respir J. 2007, 30: 1103-1110. 10.1183/09031936.00042107.

    PubMed  CAS  Article  Google Scholar 

  42. 42.

    Maclean MR, Dempsie Y: The serotonin hypothesis of pulmonary hypertension revisited. Adv Exp Med Biol. 2010, 661: 309-322. 10.1007/978-1-60761-500-2_20.

    PubMed  CAS  Article  Google Scholar 

  43. 43.

    Johnson SR, Granton JT, Mehta S: Thrombotic arteriopathy and anticoagulation in pulmonary hypertension. Chest. 2006, 130: 545-552. 10.1378/chest.130.2.545.

    PubMed  CAS  Article  Google Scholar 

  44. 44.

    Pongracz JE, Stockley RA: Wnt signalling in lung development and diseases. Respir Res. 2006, 7: 15-10.1186/1465-9921-7-15.

    PubMed  PubMed Central  Article  Google Scholar 

  45. 45.

    Yates LL, Dean CH: Planar polarity: a new player in both lung development and disease. Organogenesis. 2011, 7: 209-216. 10.4161/org.7.3.18462.

    PubMed  PubMed Central  Article  Google Scholar 

  46. 46.

    Königshoff M, Eickelberg O: WNT signaling in lung disease: a failure or a regeneration signal?. Am J Respir Cell Mol Biol. 2010, 42: 21-31. 10.1165/rcmb.2008-0485TR.

    PubMed  Article  Google Scholar 

  47. 47.

    de Jesus Perez VA, Alastalo TP, Wu JC, Axelrod JD, Cooke JP, Amieva M, Rabinovitch M: Bone morphogenetic protein 2 induces pulmonary angiogenesis via Wnt-beta-catenin and Wnt-RhoA-Rac1 pathways. J Cell Biol. 2009, 184: 83-99. 10.1083/jcb.200806049.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  48. 48.

    Buchsbaum RJ: Rho activation at a glance. J Cell Sci. 2007, 120: 1149-1152. 10.1242/jcs.03428.

    PubMed  CAS  Article  Google Scholar 

  49. 49.

    Fukumoto Y, Matoba T, Ito A, Tanaka H, Kishi T, Hayashidani S, Abe K, Takeshita A, Shimokawa H: Acute vasodilator effects of a Rho-kinase inhibitor, fasudil, in patients with severe pulmonary hypertension. Heart. 2005, 91: 391-392. 10.1136/hrt.2003.029470.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

  50. 50.

    Norrby K: Mast cells and angiogenesis. APMIS. 2002, 110: 355-371. 10.1034/j.1600-0463.2002.100501.x.

    PubMed  CAS  Article  Google Scholar 

  51. 51.

    Chen L, Endler A, Shibasaki F: Hypoxia and angiogenesis: regulation of hypoxia-inducible factors via novel binding factors. Exp Mol Med. 2009, 41: 849-857. 10.3858/emm.2009.41.12.103.

    PubMed  CAS  PubMed Central  Article  Google Scholar 

Download references


This work was supported by the Health Science Research Grants for Research on Emerging and Re-emerging Infectious Diseases (H22-Shinkou-Ippan-8 and H23-Shinkou-Ippan-018) from the Ministry of Health, Labor and Welfare of Japan, a grant from the Strategic Basis on Research Grounds for Non-governmental Schools at Heisei 20th from the Ministry of Education, Culture, Sports, Science and Technology of Japan to KS, a grant from the Strategic Research Foundation Grant-aided Project for Private schools at Heisei 23rd from Ministry of Education, Culture, Sports, Science and Technology of Japan to KS, 2011–2015, a Toho University project grant #23-19, 21, and 28, Yokohama Foundation for Advancement of Medical Science grant to YO and MS, Dr Yanase’s grant of Toho University Medical School to YO, and Grant-in-Aid for Young Scientists B from the Ministry of Education, Culture, Sports, Science and Technology of Japan (KAKENHI, #24790364) to YO.

Author information



Corresponding author

Correspondence to Megumi Wakayama.

Additional information

Competing interests

Dr. Shibuya reports receiving research grants from Pfizer Inc., Janssen Pharmaceutical K.K., Dainippon Sumitomo Pharma Co., Astellas Pharmaceutical Company, Taiho Pharmaceutical, and Pola pharma. All authors declare that they have no conflict of interest.

Authors’ contributions

KS designed the study, participated in the data collections and the interpretation of the results, and drafted the manuscript as the first author. YO carried out the experiments and the histopathological evaluation and revised the manuscript as an equal contributor of the present study. EO, HK, and MS carried out the experiments and participated in the data collections. HN and TI carried out the histopathological evaluation and performed the statistical analysis. DS, NT, and TN carried out the histopathological evaluation and prepared figures. MW participated in the interpretation of the results, revised the manuscript, and gave final approval to the manuscript as a corresponding author. TS conceptualized the study and advised the first author on pulmonary hypertension as a clinical doctor. KK conceptualized the study and helped to draft the manuscript. KS conceptualized and designed the study, carried out histopathological and statistical evaluation, and revised the manuscript as a last author. All authors contributed to conceptualizing and writing this study. Furthermore, all authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Shimodaira, K., Okubo, Y., Ochiai, E. et al. Gene expression analysis of a murine model with pulmonary vascular remodeling compared to end-stage IPAH lungs. Respir Res 13, 103 (2012).

Download citation


  • Pulmonary Vascular Remodeling
  • Stachybotrys chartarum
  • BMP signaling
  • BMPR2
  • PCP pathway