VEGF 189 Expression Is Highly Related to Adaptation of the Plateau Pika ( Ochotona curzoniae ) Inhabiting High Altitudes

Abstract

Li, Hongge, Songchang Guo, Yongming Ren, Depeng Wang, Honghao Yu, Wenjing Li, Xinquan Zhao, and Zhijie Chang. VEGF₁₈₉ expression is highly related to adaptation of the plateau pika (Ochotona curzoniae) inhabiting high altitudes. High Alt Med Biol 14:395–404, 2013.— The plateau pika (Ochotona curzonia) has adapted to high-altitude hypoxia during evolution. Higher microvessel density in specific tissues and a blunted hypoxic pulmonary vasoconstriction response are the critical components of this adaptation. VEGF, vascular endothelial growth factor, has proved to be a key regulator of angiogenesis in response to tissue hypoxia and to play an important role in vascular vasodilation. However, the role of VEGF in adaptation to high-altitude hypoxia in the plateau pika remains unknown. In this study, we cloned cDNAs for VEGF₁₆₅ and VEGF₁₈₉ and examined their expression in pikas inhabiting altitudes of 3200 and 4750 m. Phylogenetic analysis reveals that pika VEGF₁₆₅ and VEGF₁₈₉ are evolutionarily conserved. Real-time PCR analysis demonstrates that VEGF₁₆₅ and VEGF₁₈₉ display tissue and altitude-specific expression patterns. Interestingly, we found that the levels of VEGF₁₈₉ mRNA are significantly higher than those of VEGF₁₆₅ in the brain and muscle tissues of the pika, which is different from what was previously observed in sea-level mammals. VEGF₁₈₉ mRNA levels in brain, muscle, and lung of the pika increased with increased habitat altitude, whereas VEGF₁₆₅ shows less change. Our study suggests an important role for VEGF₁₈₉ in adaptation to hypoxia by the plateau pika in the high-altitude environment.

Introduction

Hypoxia and cold climate are the two most important ecological challenges for plateau animals. Native high-altitude animals have developed unique cardiovascular, hematological, and respiratory systems that ensure a successful adaptation to the harsh environment. Higher microvessel density, increased red blood corpuscle count, and blunted hypoxic pulmonary vasoconstriction response, which help animals to adapt to high-altitude hypoxia, have been reported in many plateau mammals (Anand et al., 1988, Anand et al., 1986, Azwai et al., 2007, Benavides et al., 1989, Durmowicz et al., 1993, Ge et al., 1998, Jiang et al., 1991, Jiang et al., 1992, Qi et al., 2008, Wang et al., 2008, Wei et al., 2006, Ye et al., 1994).

The plateau pika (Ochotona cuzoniae) is a small, nonhibernating mammal that inhabits alpine meadows above 3000 m, and is a keystone species on the Qinghai–Tibetan plateau ecosystem (Fig. 1). Over the course of evolution, the plateau pika has acquired a remarkable tolerance to hypoxia and the cold plateau environment by improving the utilization of oxygen (Du et al., 1982), increasing the resting metabolic rate, and having a nonshivering thermogenesis (Li et al., 2001). A number of strategies that are used by the pika to adapt to the harsh environment have been elucidated. These strategies include: 1) an increase in erythrocyte count, reduction in the mean corpuscular volume (Ye et al., 1994), changes in hemoglobin (Hb) and 2,3-diphosphoglycerate concentrations (Ge et al., 1998), and an increase in the oxygen affinity to Hb (He et al., 1994); 2) adaptation to its environment in terms of hepatic function, metabolism (Li et al., 1986), and neuroendocrine function (Du et al., 1983, Wu et al., 2001); 3) blunted hypoxic pulmonary vasoconstriction and thin-walled pulmonary arterioles (Ge et al., 1998); and 4) higher capillary densities, mitochondrial density, and concentration of myoglobin (for heart) in heart and lung (Qi et al., 2008, Wang et al., 2008). All these strategies reveal that the pika's adaptation to high altitudes occurs at both the physiological and biochemical levels.

FIG. 1.

The plateau pika (Ochotona curzoniae), also known as the black-lipped pika, is a species of mammal in the Ochotonidae family commonly known as pika. They prefer to live in elevation of 3000 to 5000 m, mostly in the Tibetan Plateau.

The molecular basis of these adaptations in the pika have occurred because of evolutionary changes in several critical factors, including HIF-1α (Li et al., 2008, Zhao et al., 2004), leptin (Yang et al., 2006, 2008), and hemoglobin (Yingzhong et al., 2007). These factors have been reported to be expressed in a tissue-specific manner (Zhao et al., 2004) and with changes that are related to altitude (Li et al., 2008, Yang et al., 2006).

Vascular endothelial growth factor (VEGF) plays a key role in angiogenesis, the formation of new capillaries from pre-existing vasculature, during development and in response to tissue hypoxia. In addition, VEGF has been reported to cause vasodilatation in vitro in a dose-dependent manner (Ku et al., 1993) and to have effects on hematopoiesis and bone marrow-derived cells (Broxmeyer et al., 1995, Gerber et al., 2002, Hattori et al., 2001). Furthermore, VEGF expression is tightly regulated at the levels of transcription, mRNA stability (Levy et al., 1998) and translation (Stein et al., 1998) under hypoxia. Hypoxia inducible factor 1 (HIF-1), a key mediator of the hypoxic response, induces transcription of VEGF (Forsythe et al., 1996).

VEGF is expressed as several different isoforms due to alternative splicing of the primary transcript depending on the presence or absence of exons 6 and 7 (Yamazaki et al., 2006). These isoforms have differing affinity for heparin, therefore are either secreted (soluble, VEGF₁₂₁ and VEGF₁₆₅) or cell- or matrix-associated (VEGF₁₈₉, VEGF₂₀₆, and partially VEGF₁₆₅). The different isoforms also have differing binding affinity to the VEGF receptors (Gitay-Goren et al., 1992), which results in the diversity of their bioactivities (Yamazaki et al., 2006).

Since VEGF has a fundamental role in angiogenesis, the response to hypoxia and its effects on vasodilation, we hypothesize that VEGF may have a different expression pattern in the plateau pika. In this study, we cloned cDNAs for the VEGF₁₆₅ and VEGF₁₈₉ isoforms from pika, and determined their expression patterns in the plateau pika. Our results show that VEGF₁₆₅ and VEGF₁₈₉ display tissue and altitude-specific expression patterns in this animal.

Materials and Methods

Animal tissue preparation

Plateau pikas were captured near the Haibei Alpine Meadow Research Station, Chinese Academy of Sciences (altitude 3200 m) and Hoh Xil region near Kunlun Mountain (altitude 4750 m) where vegetation type is alpine meadow in Qinghai province, China. The annual mean air temperatures at these two sites are −1.7°C and −11.7°C, respectively, for the last decade. Zhou Le et al. found that genetic distance and geographical distance of plateau pika population have no significant correlation (Zhou et al., 2007). Yang et al. (2008) found that altitude has no significant effect on substitution rates of the gene-like leptin in pika. Here ten individuals of plateau pika from each site were used for mRNA analysis. Animals were killed by cervical dislocation and immediately dissected at collection. Heart, lung, liver, spleen, kidney, brain, and muscle tissues were rapidly removed and frozen in liquid nitrogen. All procedures involving the handling and care of animals were in accordance with China's Practice for the Care and Use of Laboratory Animals and were approved by the Chinese Zoological Society.

RNA and cDNA preparation

Total RNA was extracted and purified from heart, lung, liver, spleen, kidney, brain, and muscle tissues of the plateau pika using TRIZOL reagent (Invitrogen, Carlsbad, CA, USA). RNA samples were then treated with RNase-free DNase I (TaKaRa Biotechnology Co. Ltd., Dalian, China), and the concentration was determined with an Ultrospec 3000. RNA integrity was tested by electrophoresis through a formaldehyde-denaturing 1% agarose gel. Four micrograms of total RNA treated with DNase I were used for first-strand cDNA synthesis using the RevertAid™ H Minus First Strand cDNA Synthesis kit (MBI, Fermantase, Opelstr., Germany) with oligo(dT)₁₈ primers in a final volume of 20 μL. Once synthesized, cDNA samples were diluted 20-fold with nuclease-free water and used for standard PCR or real-time PCR reactions.

Cloning

Primers were designed based on the full-length coding sequence (CDS) of human, mouse, and Norway rat VEGF_165/164. Primers are shown in Figure 2. Using plateau pika brain and muscle cDNA as templates, amplification was carried out with an initial denaturation at 94°C for 5 min, followed by 30 cycles of 94°C for 40 sec, 52°C for 30 sec and 72°C for 1 min, and a final extension at 72°C for 10 min. Amplified DNA fragments were subcloned into the pGEM-T Easy Vector (Promega, Madison, WI, USA) and sequenced.

FIG. 2.

PCR primers used for RT-PCR and real-time RT-PCR. (A) Schematic of the plateau pika VEGF₁₆₅ and VEGF₁₈₉ open reading frames (ORFs). Exons are indicated with 1, 2, 3, 4, 5, 6, 7, and 8. The black, pink, and blue arrows indicate the positions of the PCR primers for cloning, RT-PCR, and real-time RT-PCR analysis, respectively. The red part is exon 6 of the pika VEGF₁₈₉, which is not found in pika VEGF₁₆₅. (B) Primers used in the clone, RT-PCR, and real-time RT-PCR analysis.

Sequence analysis

The CDS of plateau pika VEGF was translated to obtain the protein sequence using BioEdit software. Multiple sequence alignments were performed with the deduced VEGF_165/164 and VEGF_189/188 protein sequences of Ochotona curzoniae, Homo sapiens (VEGF₁₆₅, AAA35789; VEGF₁₈₉, CAC19513), Oryctolagus cuniculus (Watkins et al., 1999), Rattus norvegicus (VEGF₁₆₄, AAL07526; VEGF₁₈₈, AAL07528), and Mus musculus (VEGF₁₆₄, NP_033531; VEGF₁₈₈, NP_001020421) using CLUSTAL X 1.81 (Thompson et al., 1994). Phylogenetic trees were constructed from the amino acid sequences of VEGF_165/164 and VEGF_189/188 of plateau pika and other species using the neighbor-joining method with MEGA version 4.0 (Tamura et al., 2007).

RT-PCR analysis

Expression of total VEGF, VEGF_165, and VEGF₁₈₉ mRNAs were determined in pika heart, lung, liver, spleen, kidney, brain, and muscle by a one-step RT-PCR method. Aliquots of total RNA were reverse-transcribed at 50°C for 30 min. The primers (sense, 5′TTGCTGCTCTACCTCCAC3′; antisense, 5′ATGTCCACCAAGGTCTCG3′) for the amplification of total VEGF were designed in the common coding regions (nucleotides 1–422) of all VEGF isoforms. One pair of primers was designed to determine both VEGF₁₆₅ and VEGF₁₈₉. PCR products from this pair of primers have different sizes and can be separated on gel electrophoresis. These primers are shown in Figure1.

Real-time PCR analysis

PCR amplification was performed by using SYBR GreenER^™ qPCR Supermix (Invitrogen) in a final reaction volume of 20 μL. Each reaction contained 0.8 μL of diluted cDNA samples or plasmid along with 200 nM primers. For standard curves, dilutions of plasmid-DNA constructs containing the amplicons of each gene were used. The PCR plate was incubated at 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 sec and 60°C for 1 min. Absolute gene quantification was carried out using the ABI PRISM 7000 Sequence Detection System (Applied Biosystems, Foster City, CA, USA). Samples were isolated from individual animals, and each sample was tested in triplicate. Normalization was done by dividing the average expression measurement of each gene by that of β-actin in each sample.

In order to specifically amplify pika VEGF₁₆₅ or VEGF₁₈₉, reverse primers specific to each sequence were designed, which also do not allow genomic DNA to be amplified. All primers are shown in Figure 2.

Data analysis

Statistical analyses were performed using SPSS 13.0 software. A student's t-test was used to evaluate the significance of differences between groups of experiments. All values were expressed as mean±S.D. and P<0.05 was considered to be statistically significant.

Results

Cloning and characterization of VEGF cDNAs in the plateau pika

We cloned two VEGF isoforms from the brain and muscle of plateau pika with full-length CDS of 576 bp and 648 bp. These two isoforms of pika VEGF were predicted to contain 191 and 215 amino acids, both with an apparent signal peptide (26 amino acids). The nucleotide sequences and the predicted amino acid sequences are presented in Figure 3 and were deposited in GenBank (cDNA: EU262733 and EU262734; protein: ACA23169 and ACA23170, respectively). The mature secreted proteins are predicted to be 165 and 189 amino acids in length; therefore we named the two isoforms pVEGF₁₆₅ and pVEGF₁₈₉, respectively. The molecular weights (MW) of the two isoforms were calculated to be 19.43 KDa and 22.32 KDa and have isoelectric points (pI) of 7.93 and 9.27. The pika VEGFs were predicted to have one N-glycosylation site, one N-myristoylation site, one casein kinase II phosphorylation site, five protein kinase C phosphorylation sites, and one cAMP- and cGMP-dependent protein kinase phosphorylation site. In addition, there is one bipartite nuclear localization signal profile in VEGF₁₈₉.

FIG. 3.

Plateau pika VEGF cDNA and deduced amino acid sequence. The sequence represents the combined sequence of two cDNA derived from pika VEGF₁₆₅ and VEGF₁₈₉. Predicted motifs and the signal peptide sequence are shaded with a gray background with abbreviated names under them. The blue font indicates a bipartite nuclear localization signal profile (NLS_BP) in VEGF₁₈₉. The part of the cDNA (72bp) and amino acid sequence (24aa) missing in VEGF₁₆₅ is boxed. Oligonucleotide primers used to generate the clones are underlined. The abbreviation for the motifs are: CK-2, casein kinase II phosphorylation site; N-glycosylation, N-glycosylation site; N-myristoylation, N-myristoylation site; PKC, protein kinase C phosphorylation site; CAMP, cAMP- and cGMP-dependent protein kinase phosphorylation site.

Analysis of the amino acid sequence reveals that the pika VEGF₁₆₅ shares 94, 96, 90, 90, 95, 96, 75, and 76% of identity to orthologs in human, rabbit, mouse, Norway rat, bovine, dog, chicken, and African clawed frog. The amino acid sequence identity of pika VEGF₁₈₉ to orthologs in human, rabbit, mouse, Norway rat, dog, and chicken is also high (94, 96, 90, 90, 95, and 76%, respectively). Cysteine residues, which can form disulfide bonds, at positions 26, 51, 57, 60, 61, 68, 102, and 104 of the receptor-binding domain; and 141, 144, 159, 161, 163, 170, 182, and 184 in the heparin-binding domain of pika VEGF₁₆₅ and VEGF₁₈₉ are highly conserved among species (Fig. 4).

FIG. 4.

Multiple alignment of the amino acid sequence of plateau pika VEGF_165/189 with those from human (Homo sapiens), rabbit (Oryctolagus cuniculus), Norway rat (Rattus norvegicus), and mouse (Mus musculus). Identical amino acids are shaded with a gray background. The receptor-binding domain and heparin-binding domain are boxed. Conserved cystein residues are indicated by the asterisks.

Phylogenetic trees were constructed based on the VEGF₁₆₅ and VEGF₁₈₉ sequences of plateau pika and the other vertebrates, respectively. The results indicate that both VEGF₁₆₅ and VEGF₁₈₉ from pika are most closely related to those from the rabbit (Fig. 5A and B). These phylogenetic trees are in accordance with those obtained from hemoglobin (Yingzhong et al., 2007), HIF-1α, and Cytb (not published).

FIG. 5.

Phylogenetic relationships of the plateau pika and other species VEGF_165/164 (A) and VEGF_{189 /188} (B).

Pika and human VEGF₁₆₅ and VEGF₁₈₉ show some differences, although they are highly conserved. Differences occur in the receptor-binding domain for both VEGF₁₆₅ and VEGF₁₈₉ (Glu6Gly, Asp8Gly, His9Gln, Lys10Asn, Pro11His, Arg22Gln, Arg27His, Ala58Gly, Ser65Gly, and Phe74Ser in mature protein) as well as a change in Ala139Val, next to the heparin-binding domain of VEGF₁₈₉ (Fig. 6).

FIG. 6.

Comparison of the amino acid sequences of VEGF₁₆₅ and VEGF₁₈₉ between plateau pika and human. Differences between the sequences are shown in the white region. A receptor-binding domain and a heparin-binding domain are boxed.

Expression of total VEGF in plateau pika

Total VEGF mRNA was expressed in a tissue-specific mode in the plateau pika. The mRNA levels of VEGF were the highest in lungs, and the levels were higher in hearts, livers, kidneys, and muscles than in spleens and brains. This expression pattern in the pika is similar with that in mouse living at sea level (data not shown).

Expression patterns of plateau pika VEGF₁₆₅ and VEGF₁₈₉

To study the expression patterns of pika VEGF₁₆₅ and VEGF₁₈₉, RT-PCR was performed using mRNA from different tissues. Results demonstrate that VEGF₁₆₅ and VEGF₁₈₉ were widely expressed in the pika (heart, lung, liver, spleen, kidney, brain, and muscle), but the relative amounts of VEGF₁₆₅ and VEGF₁₈₉ mRNAs varied among tissues (Fig. 7A). To determine the relative levels of VEGF₁₆₅ and VEGF₁₈₉ mRNA accurately, a real-time RT-PCR analysis was performed. The results of this analysis showed that the expression of the two isoforms was quite distinct in the different tissues of the pika. VEGF₁₈₉ was most highly expressed in the lung, but also highly abundant in the heart. The expression of VEGF₁₆₅ was not as high as VEGF₁₈₉ in the lung and heart and showed similar levels in the heart, lung, liver, and kidney. The relative expression levels of VEGF₁₆₅ and VEGF₁₈₉ mRNA showed some interesting patterns: VEGF₁₈₉ was the dominant isoform in the heart, lung, and muscle tissues, while both were roughly equally abundant in liver and kidney tissues and both were expressed weakly in the spleen and brain tissues (Fig. 7B). Despite low expression, the level of VEGF₁₈₉ is higher than that of VEGF₁₆₅ in the brain.

FIG. 7.

Expression of total VEGF,VEGF₁₆₅ and VEGF₁₈₉ mRNA in different plateau pika tissues. (A) RT-PCR analysis of total VEGF, VEGF₁₆₅, and VEGF₁₈₉ expression in the pika heart, lung, liver, spleen, kidney, brain, and muscle. (B) Transcript levels of VEGF₁₆₅ and VEGF₁₈₉ assessed by real time RT-PCR. Expression levels were normalized to β-actin. All results are expressed as means±S.D. (n=10).

Expression of VEGFs in plateau pika in relation to habitat altitude

Since VEGF plays an important role in the tissue response to hypoxia, we questioned whether the expression of VEGF₁₆₅ and VEGF₁₈₉ would be different in the pikas inhabiting different altitudes. To resolve this question, we determined the mRNA levels of VEGF₁₆₅ and VEGF₁₈₉ in the brain, muscle, and lung tissues of pikas inhabiting altitudes of 3200 and 4750 m. Then VEGF₁₆₅ and VEGF₁₈₉ mRNA levels were normalized by β-actin since β-actin, as a housekeeping gene, has no significant difference between these two groups (Li et al., 2009).The data show that VEGF₁₆₅ and VEGF₁₈₉ mRNA levels were significantly increased in the brains of the animals inhabiting the altitude of 4700 m, compared with those at 3200 m (p<0.01 and p<0.01, respectively). The levels of VEGF₁₈₉ appear to be increased more than those of VEGF₁₆₅ (Fig. 8A). In the muscle (Fig. 8B) and lung (Fig. 8C), VEGF₁₈₉ levels were also significantly increased in the animals inhabiting an altitude of 4750 m compared to those at 3200 m (p<0.01 and p<0.001, respectively), while VEGF₁₆₅ shows no significant difference. These data indicated that VEGF₁₈₉ mRNA levels are increased in the animals that inhabit higher altitudes.

FIG. 8.

Comparison of VEGF₁₆₅ and VEGF₁₈₉ mRNA levels in plateau pikas inhabiting different altitudes. mRNA expression levels of VEGF₁₆₅ and VEGF₁₈₉ in brain (A), muscle (B), and lung (C) of the pika collected at two different altitudes (3200 and 4750 m) are shown. Expression levels were normalized to β-actin. All the results are expressed as means±S.D. (n=10). Asterisks indicate significant differences compared with the altitude of 3200 m (**p<0.01, ***p<0.001).

Discussion

To adapt to a high-altitude environment, plateau animals have developed unique characteristics. The plateau pika, as a typical plateau mammal, has been extensively studied to understand this adaptation. HIF-1α (Li et al., 2008, Zhao et al., 2004) and leptin (Yang et al., 2006) expression, and leptin (Yang et al., 2008) and Hb (Yingzhong et al., 2007) amino acid substitutions have been linked to adaptation in the pika. Here we present for the first time the isolation and characterization of VEGF, another gene that is sensitive to the environment, from the plateau animal. The novel expression pattern of VEGF₁₆₅ and VEGF₁₈₉ may be involved in adaptation of the pika since phylogenetic analysis reveals the conservation of pika VEGF₁₆₅ and VEGF₁₈₉.

Several studies have shown that VEGF_165/164 and VEGF_189/188 isoforms mRNA display tissue-specific expression. VEGF₁₈₈ is the most abundant isoform in lung and heart, whereas VEGF₁₆₄ is more highly expressed in murine muscle, brain, liver, spleen, and kidney (Bacic et al., 1995, Ng et al., 2001). Watkins et al. (1999) reported that VEGF₁₈₉ is most highly expressed in rabbit lung, whereas VEGF₁₆₅ is most abundant in liver, spleen, and kidney. The predominance of VEGF₁₈₈ in heart and lung distinguishes these two tissues from all others examined. This suggests that VEGF_189/188 may play a specific or exclusive role in cardiopulmonary function. It has been demonstrated that VEGF₁₈₈ produced by the rat pulmonary epithelium mediates the assembly and stabilization of the highly organized vessel network that surrounds the alveoli (Ng et al., 2001).

In our study, RT-PCR and real time RT-PCR-based expression analysis revealed that VEGF₁₈₉ mRNA expression was significantly higher than VEGF₁₆₅ in the plateau pika brain and muscle. This is a novel finding because earlier studies had shown that VEGF₁₆₅ mRNA expression is higher than VEGF₁₈₉ in the brain and muscle of mice and rat, animals found at sea level (Bacic et al., 1995, Ng et al., 2001). The recent study (Zheng et al., 2011) found that the microvessel density is significantly increased in the pika brain compared with sea-level rat brain (602.1±66.98 vs. 426.75±58.25 microvessels per mm²). Therefore, higher-level VEGF₁₈₉ may be required for more microvessels in pika brain to adapt to hypoxia since VEGF₁₈₉ was reported to have the capability to induce the most-dense, small, sprouting microvessels (Yuan et al., 2011). However, there is a lower microvessel density but higher Mb (myoglobin) in the pika skeletal muscle compared to sea-level rat muscle (microvessel density: 603.70±115.88 vs. 1009.63±138.48 microvessels per mm²; Mb: 0.91±0.09 vs. −0.08±0.49 logC_Mb/logC_GAPDH). Yan, et al.(2009) found the muscle-specific Mb overexpressing transgenic mice had a lower vascular density, although it, at VEGF level, has similar response to ischemia injury with wild-type mouse. VEGF overexpression has been reported to increase the Mb expression rather than angiographic score of collateral arteries in the mouse muscle (van Weel et al., 2004). Therefore, the higher-level VEGF₁₈₉ in the pika muscle may be involved in its adaptation to high altitude hypoxia through increasing Mb expression since Mb has higher affinity to O₂.

Changes in the ratio of VEGF isoforms have been demonstrated to be related to responses to systemic hypoxia (Hofstaetter et al., 2004) or glucose starvation (Zhang et al., 2002). Systemic hypoxia has been shown to differentially affect expression of VEGF isoform mRNA (Marti et al., 1998, Watkins et al., 1999). Birot et al. (2004) have reported that VEGF₁₈₈ levels are transiently elevated in the rat heart in response to systemic hypoxia, whereas no change was found in VEGF₁₆₄ mRNA levels. Transient elevation of VEGF has been demonstrated to occur at the time of initial adaptation in the mouse brain in response to chronic hypoxia (Kuo et al., 1999). Hypoxia is well known to be a major feature of the plateau environment. In this study, we found that VEGF₁₈₉ mRNA levels remain higher in the brain and muscle of native pika inhabiting higher altitudes. This may account for the higher abundance of VEGF₁₈₉ compared with VEGF₁₆₅ in the pika brain and muscle. It should be noted that the increase is permanent, not transient. This suggests that pika may need to retain relatively higher levels of VEGF₁₈₉ than VEGF₁₆₅ expression in the brain and muscle in order to adapt to the higher-altitude environments, as changes in VEGF₁₆₅ expression were less than those for VEGF₁₈₉ levels between pikas inhabiting different altitudes.

The two sites in the present study have an extremely cold and hypoxic plateau climate. The most obvious differences between these two sites are altitude and ambient temperature. The altitudes of the two sites are 3200 m and 4750 m, with average temperatures of −1.7°C and −11.7°C, respectively. As altitude increases and temperature decreases, the partial pressure of oxygen decreases. Thus the partial pressure of O₂ (Po₂) at higher altitude is lower than the Po₂ at lower altitude. In our study, we demonstrate that under in vivo conditions, in ambient atmosphere, VEGF₁₈₉ mRNA is expressed at higher levels in lungs in pika inhabiting a higher altitude compared with those at a lower altitude. This trend is similar to the changes in HIF-1α (the hypoxia inducible subunit of HIF-1) protein expression in pika lung with altitudes (Li et al., 2008). It is well known that hypoxia-induced expression of VEGF is under the control of HIF-1 in other species (Damert et al., 1997), therefore it can be deduced that the higher constitutive level of VEGF₁₈₉ mRNA in the lung may be supported by the higher constitutive levels of HIF-1α protein in the pika inhabiting higher altitudes. However, Rissanen et al. (2006) found that temperature is also involved in regulation of HIF-1α in crucian carp. Here, low temperature may also play an important role in altitude-related VEGF regulation (Rissanen et al., 2006), since cold is reported to be involved in angiogenesis through upregulating VEGF expression by HIF in mouse adipose tissue (Xue et al., 2009). Therefore, the two prime ecological factors, hypoxia and cold, both may play an important role in the adaption of pika to harsh altitude climates through VEGF.

In summary, we have cloned VEGF₁₆₅ and VEGF₁₈₉ isoforms from the plateau pika, and found that they were expressed in a tissue-specific manner. Specifically, VEGE₁₈₉ has a higher expression than VEGF₁₆₅ in pika brain and muscle, which may be involved in pika brain's and muscle's adaptation to hypoxia through different mechanisms. Further, VEGE₁₈₉ levels in pika brain, muscle, and lung are correlated to the altitudes of their habitat. All these changes could allow the animal to adapt to the hypoxia challenges of its natural habitat.

Footnotes

Acknowledgments

We thank Dr. David M. Irwin (University of Toronto) for reading and revising our manuscript, and Dr. Tong Sun and Dr. Yan-Jie Li for their help in the real-time PCR experiment

Author Disclosure Statement

The authors declare no conflicting financial interests.

This work was greatly supported by a grant from The Knowledge Innovation project of the Chinese Academy of Sciences (0254051211).

References

Anand

, Heath

, Williams

, Deen

, Ferrari

, Bergel

, and Harris

. (1988). The pulmonary circulation of some domestic animals at high altitude. Intl J Biometeorol, 32:56–64.

Anand

, Harris

, Ferrari

, Pearce

, and Harris

. (1986). Pulmonary haemodynamics of the yak, cattle, and cross breeds at high altitude. Thorax, 41:696–700.

Azwai

, Abdouslam

, Al-Bassam

, Dawek

AMA

, and Al-Izzi

SAL

. (2007). Morphological characteristics of blood cells in clinically normal adult llamas (lama glama). Veterinarski Arhiv, 77:69–79.

Bacic

, Edwards

, and Merrill

. (1995). Differential expression of vascular endothelial growth factor (vascular permeability factor) forms in rat tissues. Growth Factors, 12:11–15.

Benavides

, Perez

, Espinoza

, Cabello

, Riquelme

, Parer

, and Llanos

. (1989). Cardiorespiratory functions in the fetal llama. Respir Physiol, 75:327–334.

Birot

, Peinnequin

, Simler

, van Cuyck-Gandre

, Hamel

, and Bigard

. (2004). Vascular endothelial growth factor expression in heart of rats exposed to hypobaric hypoxia: Differential response between mrna and protein. J Cell Physiol, 200:107–115.

Broxmeyer

, Cooper

, Li

, Lu

, Song

, Kwon

, Warren

, and Donner

. (1995). Myeloid progenitor cell regulatory effects of vascular endothelial cell growth factor. Int J Hematol, 62:203–215.

Damert

, Ikeda

, and Risau

. (1997). Activator-protein-1 binding potentiates the hypoxia-induciblefactor-1-mediated hypoxia-induced transcriptional activation of vascular-endothelial growth factor expression in c6 glioma cells. Biochem J, 327:419–423.

, and Li

. (1982). Effect of simulated hypoxia acclimation on organism, organ and hematology in ochotona curzoniae and rats. Acta Theriolog Sinica, 2:35–42.

10.

, Li

, and Chen

. (1983). The changing of corticosterone level in native ochotona curzoniae. Acta Theriolog Sinica, 3: 47–52.

11.

Durmowicz

, Hofmeister

, Kadyraliev

, Aldashev

, and Stenmark

. (1993). Functional and structural adaptation of the yak pulmonary circulation to residence at high altitude. J Appl Physiol, 74:2276–2285.

12.

Forsythe

, Jiang

, Iyer

, Agani

, Leung

, Koos

, and Semenza

. (1996). Activation of vascular endothelial growth factor gene transcription by hypoxia-inducible factor 1. Mol Cell Biol, 16:4604–4613.

13.

, Kubo

, Kobayashi

, Sekiguchi

, and Honda

. (1998). Blunted hypoxic pulmonary vasoconstrictive response in the rodent Ochotona curzoniae (pika) at high altitude. Am J Physiol, 274:H1792–1799.

14.

Gerber

, Malik

, Solar

, et al. (2002). VEGF regulates haematopoietic stem cell survival by an internal autocrine loop mechanism. Nature, 417:954–958.

15.

Gitay-Goren

, Soker

, Vlodavsky

, and Neufeld

. (1992). The binding of vascular endothelial growth factor to its receptors is dependent on cell surface-associated heparin-like molecules. J Biol Chem, 267:6093–6098.

16.

Hattori

, Dias

, Heissig

, et al. (2001). Vascular endothelial growth factor and angiopoietin-1 stimulate postnatal hematopoiesis by recruitment of vasculogenic and hematopoietic stem cells. J Exp Med, 193:1005–1014.

17.

, Xu

, Meng

, Li

, and Wang

. (1994). Comparative analysis in transport and intake of oxygen between pikas (Ochotona curzoniae) and rats. J Prevent Med Chinese People's Liberation, 12:431–435.

18.

Hofstaetter

, Saad

, Samuel

, Wunderlich

, Choi

, and Glimcher

. (2004). Differential expression of VEGF isoforms and receptors in knee joint menisci under systemic hypoxia. Biochem Biophys Res Commun, 324:667–672.

19.

Jiang

, Gama

, and He

. (1991). A comparison on several hematologic values of yaks on tibet plateau at different altitudes. Acta Veterin Zootech Sinica, 22:20–26.

20.

Jiang

, He

, Gama

, and Bianba

. (1992). Comparison on several hematologic value of goat in tibet plateau at different altitude. Southwest China J Agricult Sci, 5:79–83.

21.

, Zaleski

, Liu

, and Brock

. (1993). Vascular endothelial growth factor induces EDRF-dependent relaxation in coronary arteries. Am J Physiol, 265:H586–592.

22.

Kuo

, Benhayon

, Przybylski

, Martin

, and LaManna

. (1999). Prolonged hypoxia increases vascular endothelial growth factor mrna and protein in adult mouse brain. J Appl Physiol, 86:260–264.

23.

Levy

, Chung

, Furneaux

, and Levy

. (1998). Hypoxic stabilization of vascular endothelial growth factor mrna by the rna-binding protein hur. J Biol Chem, 273:6417–6423.

24.

, Ren

, Guo

, et al. (2008). The protein level of hypoxia-inducible factor-1alpha is increased in the plateau pika (Ochotona curzoniae) inhabiting high altitudes. J Exp Zool, 311A:1–8.

25.

, Ren

, Guo

, et al. (2009). The protein level of hypoxia-inducible factor-1alpha is increased in the plateau pika (Ochotona curzoniae) inhabiting high altitudes. J Exp Zool A, Ecolog Genetics Physiol, 311:134–141.

26.

, Sun

, Huang

, et al. (2001). Cold adaptive thermogenesis in small mammals from different geographical zones of china. Comp Biochem Physiol A Mol Integr Physiol, 129:949–961.

27.

, You

, Chen

, and Du

. (1986). Effect of chronic hypoxia on the liver of Ochotona curzoniae and rat. Acta Theriolog Sinica, 6:261–266.

28.

Marti

, and Risau

. (1998). Systemic hypoxia changes the organ-specific distribution of vascular endothelial growth factor and its receptors. Proc Natl Acad Sci USA, 95:15809–15814.

29.

, Rohan

, Sunday

, Demello

, and D'Amore

. (2001). Differential expression of VEGF isoforms in mouse during development and in the adult. Dev Dyn, 220:112–121.

30.

, Wang

, Zhu

, Rao

, Wei

, and Wei

. (2008). Hypoxic adaptation of the hearts of plateau zokor (Myospalax baileyi) and plateau pika (Ochotona curzoniae). Acta Physiolog Sinica, 60:348–354.

31.

Rissanen

, Tranberg

, Sollid

, Nilsson

, and Nikinmaa

. (2006). Temperature regulates hypoxia-inducible factor-1 (hif-1) in a poikilothermic vertebrate, crucian carp (Carassius carassius). J Exp Biol, 209:994–1003.

32.

Stein

, Itin

, Einat

, Skaliter

, Grossman

, and Keshet

. (1998). Translation of vascular endothelial growth factor mrna by internal ribosome entry: Implications for translation under hypoxia. Mol Cell Biol, 18:3112–3119.

33.

Tamura

, Dudley

, Nei

, and Kumar

. (2007). Mega4: Molecular evolutionary genetics analysis (mega) software version 4.0. Mol Biol Evol, 24:1596–1599.

34.

van Weel

, Deckers

, Grimbergen

, et al. (2004). Vascular endothelial growth factor overexpression in ischemic skeletal muscle enhances myoglobin expression in vivo. Circ Res, 95:58–66.

35.

Wang

, Wei

, Qi

, and Zhu

. (2008). Characteristics of pulmonary acinus structure in the plateau zokor (Myospalax baileyi) and plateau pika (Ochotona curzoniae). Acta Zoolog Sinica, 54:531–539.

36.

Watkins

, D'Angio

, Ryan

, Patel

, and Maniscalco

. (1999). Differential expression of VEGF mRNA splice variants in newborn and adult hyperoxic lung injury. Am J Physiol, 276:L858–867.

37.

Wei

, Wei

, Zhang

, and Yu

. (2006). Blood-gas properties of plateau zokor (Myospalax baileyi). Comp Biochem Physiol A Mol Integr Physiol, 145:372–375.

38.

, and Du

. (2001). Study of activity of HPA in rat and Ochotona curzoniae under hypoxia. Acta Theriolog Sinica, 21:195–198.

39.

Xue

, Petrovic

, Cao

, et al. (2009). Hypoxia-independent angiogenesis in adipose tissues during cold acclimation. Cell Metabol, 9:99–109.

40.

Yamazaki

, and Morita

. (2006). Molecular and functional diversity of vascular endothelial growth factors. Mol. Divers, 10:515–527.

41.

Yang

, Zhao

, Guo

, Li

. Qi

, Wang

, and Cao

. (2006). Leptin cDNA cloning and its mRNA expression in plateau pikas (Ochotona curzoniae) from different altitudes on Qinghai-Tibet plateau. Biochem Biophys Res Commun, 345:1405–1413.

42.

Yang

, Wang

, Zhao

, et al. (2008). Natural selection and adaptive evolution of leptin in the ochotona family driven by the cold environmental stress. PloS One, 3:e1472.

43.

Yang

, Guan

, Liu

, and Li

. (2009). Myoglobin over-expression attenuates angiogenic response in hindlimb ischemia in mice. Chin Med J, 122:1056–1060.

44.

, Cao

, and Bai

. (1994). Blood indices of plateau pika and relationship with hypoxia adaptation. Animals Sci Sinica, 2:115–120.

45.

Yingzhong

, Yue

, Guoen

, et al. (2007). Molecular cloning and characterization of hemoglobin alpha and beta chains from plateau pika (Ochotona curzoniae) living at high altitude. Gene, 403:118–124.

46.

Yuan

, Lin

, Chou

, et al. (2011). Functional and structural characteristics of tumor angiogenesis in lung cancers overexpressing different VEGF isoforms assessed by DCE- and SSCE-MRI. PloS One, 6:e16062.

47.

Zhang

, Conejo-Garcia

, Yang

, et al. (2002). Different effects of glucose starvation on expression and stability of VEGF mRNA isoforms in murine ovarian cancer cells. Biochem Biophys Res Commun, 292:860–868.

48.

Zhao

, Ning

, Zhu

, et al. (2004). Cloning of hypoxia-inducible factor 1alpha cdna from a high hypoxia tolerant mammal-plateau pika (Ochotona curzoniae). Biochem Biophys Res Commun, 316:565–572.

49.

Zheng

, Zhu

, Wang

, Wei

, and Wei

. (2011). [gene coding and mRNA expression of vascular endothelial growth factor as well as microvessel density in brain of plateau zokor: Comparison with other rodents]. Sheng li xue bao: [Acta Physiolog Sinica], 63:155–163.

50.

Zhou

, Yang

, Yu

, et al. (2007). Genetic diversity and differentiation in four natural populations of plateau pika (Ochotona curzoniae). Acta Theriolog Sinica, 27:8.

VEGF 189 Expression Is Highly Related to Adaptation of the Plateau Pika ( Ochotona curzoniae ) Inhabiting High Altitudes

Abstract

Abstract

Introduction

Materials and Methods

Animal tissue preparation

RNA and cDNA preparation

Cloning

Sequence analysis

RT-PCR analysis

Real-time PCR analysis

Data analysis

Results

Cloning and characterization of VEGF cDNAs in the plateau pika

Expression of total VEGF in plateau pika

Expression patterns of plateau pika VEGF165 and VEGF189

Expression of VEGFs in plateau pika in relation to habitat altitude

Discussion

Footnotes

Acknowledgments

Author Disclosure Statement

References

Expression patterns of plateau pika VEGF₁₆₅ and VEGF₁₈₉