Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia
© Zhao et al.; licensee BioMed Central Ltd. 2013
Received: 22 April 2013
Accepted: 16 September 2013
Published: 18 September 2013
Skip to main content
© Zhao et al.; licensee BioMed Central Ltd. 2013
Received: 22 April 2013
Accepted: 16 September 2013
Published: 18 September 2013
Responses to hypoxia have been investigated in many species; however, comparative studies between conspecific geographical populations at different altitudes are rare, especially for invertebrates. The migratory locust, Locusta migratoria, is widely distributed around the world, including on the high-altitude Tibetan Plateau (TP) and the low-altitude North China Plain (NP). TP locusts have inhabited Tibetan Plateau for over 34,000 years and thus probably have evolved superior capacity to cope with hypoxia.
Here we compared the hypoxic responses of TP and NP locusts from morphological, behavioral, and physiological perspectives. We found that TP locusts were more tolerant of extreme hypoxia than NP locusts. To evaluate why TP locusts respond to extreme hypoxia differently from NP locusts, we subjected them to extreme hypoxia and compared their transcriptional responses. We found that the aerobic metabolism was less affected in TP locusts than in NP locusts. RNAi disruption of PDHE1β, an entry gene from glycolysis to TCA cycle, increased the ratio of stupor in TP locusts and decreased the ATP content of TP locusts in hypoxia, confirming that aerobic metabolism is critical for TP locusts to maintain activity in hypoxia.
Our results indicate that TP and NP locusts have undergone divergence in hypoxia tolerance. These findings also indicate that insects can adapt to hypoxic pressure by modulating basic metabolic processes.
Animals with wide altitudinal distribution are subjected to different ambient oxygen levels, representing a promising system for disclosing the mechanism of hypoxia adaptation. These animals have been effectively used to elucidate the molecular and physiological foundations of hypoxia adaptation. For example, the bar-headed geese, Anser indicus, that can fly at 9 km elevation during migration, possesses a nonsynonymous mutation of COX3 gene (Trp116→Arg) leading to a higher COX affinity and COX activity as compared to low-altitude geese . The deer mouse, Peromyscus maniculatus, is continuously distributed from sea-level elevation to high altitude above 4,300 m. DNA sequences of two α-globin genes revealed that Asp64→Gly substitution appeared in high-altitude populations of P. maniculatus at unusually high frequency and could substantially increase oxygen-binding affinity of α-globin . Several hypoxia-responsive genes, including EPAS1, EGLN1 and PPARA, were under positive selection in Tibetan natives compared to lowlanders [3–5]. Although this approach holds much promise, currently it is mainly applied to vertebrates rather than invertebrates.
The migratory locust, Locusta migratoria, is the most widely distributed locust species in the world. Despite of their wide geographic distribution, the migratory locust around the world is taxonomically the same species . A recent research revealed that the locust population of L. migratoria in the world could be divided into two distinct lineages, the Northern lineage and the Southern lineage, both of which were originated from Africa and then separated 895,000 years ago . TP locusts, previously identified as L. migratoria tibetensis, belong to the Southern lineage, whereas NP locusts belong to the Northern lineage. Multilocus microsatellite genotyping analysis revealed that TP locusts are genetically distinct from NP locusts and they have experienced adaptive differentiation coupled to Quaternary glaciations events . A Bayesian phylogram based on locust mitochondrial sequences also supported that TP locusts had no frequent genetic exchange with NP locusts . TP locusts likely expanded into Tibet 34,000 and 40,000 years ago and have since inhabited the Tibetan Plateau . Although the partial oxygen pressure (pO2) in Tibet (~13 kPa) is above the critical oxygen pressure of locusts (~3 kPa) [9, 10], the Tibetan pO2 still exerts selective pressures on locusts . Long-term habitation to the Tibetan Plateau probably has adapted TP locusts to the chronic low-oxygen atmosphere, conferring on them superior tolerance of hypoxia.
The evolution is ultimately attributed to genetic mutations which may occur in the coding region or the regulatory region of a gene. If mutations occur in the coding region, they may modify the molecular function of the gene; if in the regulatory region, they may alter its gene expression levels. The genetic variations in the coding regions are successfully evaluated in Tibetans and lowlanders by using human SNP arrays [3, 4] and human exon capture arrays . However, such arrays are not available for locusts. Instead, we have developed the locust microarray platform, which has been proved to be an effective approach to study locust gene expression profiles [12–14]. Comparison of gene expression patterns using locust microarrays is an alternate approach to understand the natural selection and adaptive evolution of TP locusts [15, 16].
In this study, we confirmed the superior capacity of TP locusts to tolerate hypoxia through a series of phenotypic assays. Then we examined gene expression profiles of field and laboratory locusts using locust microarrays and found that TP locusts utilized the aerobic metabolism in a better way than NP locusts, which was further confirmed by RNAi experiments. Our data demonstrate that TP and NP locusts have diverged in hypoxia tolerance by different modulation of basic metabolism.
To understand how field TP locusts respond to the high-altitude environments, we identified the significantly enhanced pathways in field-collected TP locusts (Additional file 2). Most of the enhanced pathways involve amino acid metabolism such as glycine, serine and threonine metabolism. Of particular interest is the selenoamino acid metabolism, which may help reduce the oxidative stress  and decrease the radiation damage to amino acid and proteins . Citrate cycle (TCA cycle) was also up-regulated, indicating that the aerobic pathway was possibly enhanced to produce sufficient energy molecules in field TP locusts. Interestingly, olfactory transduction was also strengthened, which is possibly related to food foraging and feeding in the field TP locusts.
To investigate whether the gene PDHE1β undergoes genetic variations between TP and NP locusts, we sequenced the coding region of PDHE1β in 15 individuals of TP and NP locusts (Additional file 5). Multiple sequence alignment revealed no significant genetic variation in the coding region of PDHE1β (Additional file 8).
The wide distribution of L. migratoria makes it possible that locusts from different habitats may have different adaptations to specific environments. TP locusts likely expanded into Tibet 34,000 and 40,000 years ago  and have since inhabited the Tibetan Plateau. Long-term habitation to the Tibetan Plateau probably has adapted TP locusts to the chronic low-oxygen atmosphere, conferring on them superior tolerance of hypoxia.
We examined four phenotypic traits in TP and NP locusts (Figure 2), all of which demonstrate that TP locusts are better prepared for hypoxia tolerance. Smaller body size in TP locusts is possibly a response to the chronic hypoxia as noted in fruit flies [25, 26]. Using hypoxic stupor - a common response to extreme hypoxia in insects  - as criteria, we found that TP locusts were more tolerant of extreme hypoxia. Recovery time from stupor, once used as a criteria to screen for hypoxia-sensitive fruit flies , revealed that NP locusts were sensitive to hypoxic injury, possibly related to longer polysynaptic transmission in central nervous system . The respiration rates under a series of hypoxia suggested that hypoxia suppressed aerobic respiration drastically when pO2 was near the critical value as reported in American grasshopper, Schistocerca americana; however, TP locusts maintained higher respiration rates than NP locusts, indicating TP locusts were able to utilize oxygen more efficiently under extreme hypoxia. Taken together, these results indicate that the two geographical populations have different thresholds for hypoxia tolerance.
The gene expression profiles of the field-collected NP and TP locusts showed that multiple factors contributed to their transcriptional differences (Figure 3). One major difference between the field habitats of NP and TP locusts is the oxygen content, which decreases by 39.2% in Tibet compared to in North China Plain. The HIF target genes and antioxidant genes were differentially expressed between field NP and TP locusts, possibly due to the different oxygen content. The antioxidant genes were suppressed in TP locusts, which is possibly due to a decreased need to avoid oxygen toxicity . Another important difference is the UV radiation, which is more severe in Tibet. UV radiation can induce a variety of mutagenic and cytotoxic DNA lesions, which are harmful to the genetic stability . The up-regulation of genes repairing UV-damaged DNA reflects the different UV radiation in the field. In addition, other factors such as temperature and nutrition can also affect the lipid metabolism and cell growth [31, 32]. To identify the genes that contributed to the enhanced tolerance of hypoxia in TP locusts, we took measures to minimize the effects of factors other than oxygen.
TP locusts showed better tolerance of extreme hypoxia than NP locusts (Figure 2) and the capacity of TP locusts to modulate basic metabolism contributed to this superiority (Figures 4 and 5). Altering the metabolic pathways to produce more ATP is a common strategy to cope with hypoxia for living organisms. For example, metabolomic and computational analyses of fruitfly flight muscles under 0.5% oxygen levels showed that conversion of pyruvate from lactate to alanine and acetate could convey hypoxia tolerance by improving ATP-producing efficiency per glucose . Flux-balance analyses showed that flies adapted to 4% oxygen levels produced more ATP per glucose by lowering pyruvate carboxylase flux and preferring the usage of Complex I to Complex II [26, 34].
The metabolic processes of TCA and ETC. were less suppressed in TP locusts (Figures 4 and 5). These results indicate that TP locusts could better utilize aerobic metabolism than NP locusts to produce ATP under extreme hypoxia, and is consistent with the fact that TP locusts had a higher respiration rate under extreme hypoxia (Figure 2D). It appears that higher efficiency to utilize limited oxygen to produce more ATP helps TP locusts conquer the extreme hypoxia. As observed in mammals, the COX subunit was switched from COX4-1 to COX4-2 to maintain optimal efficiency of mitochondrial respiration under 1% oxygen levels . We confirmed the importance of aerobic metabolism in the tolerance of extreme hypoxia for TP locusts through RNA interference (Figure 6). Our results suggest that disrupting PDHE1β, a subunit of E1 component in pyruvate dehydrogenase complex (PDH), dramatically suppressed energy charge and increased stupor in hypoxia for TP locusts, indicating aerobic metabolism was critical for TP locusts to resist extreme hypoxia. PDH catalyzes an irreversible reaction that links glycolysis in the cytoplasm and TCA cycle in the mitochondrion, controlling the entry from anaerobic metabolism to aerobic metabolism . Disrupting PDHE1β would prevent transfer of substrates through flip-flop action , and thus prohibit pyruvate from entry into the TCA cycle, in which NADH and FADH2 are produced for following production of ATP. Missense mutations in PDHE1α gene, which weaken the transfer of substrates as we did using RNAi, lead to PDH deficiency and are associated with lactic acidosis and central nervous system dysfunction [36, 38]. This suggests that decreased expression of PDHE1β in NP locusts under hypoxia probably contributes to their vulnerability to extreme hypoxia and unaffected expression of PDHE1β help TP locusts conquer extreme hypoxia.
TP locusts adapted to altitude hypoxia and evolved superior capacity to produce ATP efficiently through aerobic metabolism under hypoxia. The locust strategy is well consistent with that taken by vertebrates. For example, genetic comparison between Tibetan natives and lowlanders indicated that one SNP in EPAS1 gene showed significant difference in frequency among the two populations, which is associated with erythrocyte abundance and hemoglobin concentration that is proportional to oxygen transport [4, 5]. Genomic data of Tibetan antelope showed signals of adaptive evolution and expansion of gene families associated with energy metabolism and oxygen transmission . Population analysis of deer mice from different altitudes revealed that Asp64→Gly substitution in α-globin gene is the dominant genotype in the deer mice on the high-altitude localities, which substantially increases oxygen-binding affinity of α-globin . These data show that hypoxia adaptation has genetic foundations and that selection forces favor genetic variations that help to generate more energy molecules under hypoxia.
We compared the coding sequences of PDHE1β between TP and NP locusts and found few genetic variations (Additional file 8); possibly the major difference of this gene lies in the promoter region rather than the coding region. To find the genetic evidence underlying the different hypoxia tolerance of TP and NP locusts, it is necessary to compare their genetic variations on a genome-wide scale, as was done in human using exon capture arrays  and SNP arrays [3, 4]. Our ongoing locust genome project will help decipher the genomic evidence of hypoxia adaptation in TP locusts in the future. Apart from genetic variations, an increasing number of studies show that epigenetic regulation also plays a crucial role in the cellular response to hypoxia [40–42]. To fully understand the different tolerance of hypoxia between TP and NP locusts, studies in terms of epigenetic regulation of hypoxic responses should be conducted in the future.
We confirmed the enhanced capacity of TP locusts to tolerate hypoxia, compared their gene expression profiles in field and lab populations, and validated the critical role of aerobic metabolism in TP locusts’ hypoxia tolerance. This study demonstrates that TP and NP locusts have undergone divergence in hypoxia tolerance due to their long-term inhabiting different altitude; it also indicates that insects can adapt to hypoxic pressure by modulating basic metabolic processes. We speculate that genetic variations in basic metabolic processes contribute to their different tolerance of hypoxia, which is worth examining in the future.
Locusts were collected in the field in 2005 and then raised under experimental conditions as stock cultures. TP locusts were collected in Naidong County, Duilongdeqing County, and Mozhugongka County, whose altitudes were around 4,000 m and pO2 was around 13 kPa (Figure 1). NP locusts were collected in Wudi County Shandong Province, Huanghua Hebei Province, and Tianjin City, whose altitudes were around 26 m and pO2 was around 21 kPa (Figure 1). The locusts were taken back to Beijing and raised gregariously as described before . TP locusts have a phenomenon of compulsory diapause and their eggs require about three months in 4°C to terminate diapause, while NP locusts reproduce continuously through the year. We hatched in batches the TP locust eggs belonging to the same generation and there were at most two generations in a year for TP locusts. All the experiments were done with male locusts 10 days after eclosion (Figure 2A).
Four assays were done to determine the phenotypic differences between NP and TP locusts. (1) Femur length. The left or right femurs were randomly chosen for both male and female locusts and measured using a vernier caliper (2) Stupor ratio. Locusts were subjected to hypoxia in a glass tank placed in a hypoxia chamber; compressed air was diluted with nitrogen to decrease the oxygen content to the required levels. The oxygen level was monitored by an electrode in a hypoxia chamber during the whole process. Five partial oxygen pressures (1.2 kPa, 1.6 kPa, 2.0 kPa, 2.4 kPa and 2.8 kPa) were tested, and six groups of locusts with ten individuals in each group were used at each oxygen level. When the oxygen level was very low, locusts lost coordination, collapsed on the bottom of the glass tank and fell into a motionless state of stupor. Locusts were kept at the test oxygen level for 2 hours and during this period we recorded the number of locusts in stupor every 5 minutes. The number of locusts in stupor fluctuated during the initial minutes of the treatment, but became very stable after 1 hour. We reported the number of locusts in stupor at the end of two-hour treatment. The stupor ratios were calculated at each oxygen levels for TP and NP locusts, and compared using student’s t-test after arcsine-square-root transformation. (3) Recovery time. We placed male adult locusts in syringes and then flushed the syringes for five minutes with nitrogen gas so as to clear the atmosphere (Additional file 9). After that, the syringes were sealed tightly and incubated in 30°C for 20 min. Then, the locusts were returned to normoxia, and the time taken by the locusts to become conscious (i.e. stand upright) was recorded as recovery time. Thirty individuals were tested respectively for TP and NP locusts. (4) CO2 production. CO2 production rates were measured as previously described with minor modifications  (Additional files 9 and 10). Five oxygen levels, i.e. 1.6 kPa, 3.2 kPa, 6.5 kPa, 13 kPa and 21 kPa, were selected for comparison. Male locusts subjected to the specified oxygen levels were incubated in a thermostat container for 30 min in 30°C. The locusts were weighed immediately after their CO2 production was measured.
Thoracic muscles of field populations were sampled in the field and preserved immediately in RNAlater (Ambion) to prevent RNA degradation. For laboratory populations, after locusts were subjected to 1.6 kPa pO2 for 2 hours, thoracic muscles were immediately sampled and frozen in the liquid nitrogen. Total RNA was extracted following the protocol of the RNeasy Mini Kit (Qiagen). Then, RNA was reversely transcribed into cDNA, labeled with Cy3 and Cy5, and finally hybridized to locust microarrays as previously described [11–14]. Microarrays were scanned by using a GenePix 4000B microarray scanner and image analysis was performed by using GenePix Pro 6.0 (Axon Instruments). Six independent hybridizations with biological replicates were conducted with a direct-comparison strategy. DEGs were identified with the limma package and the cutoff for DEGs was set as FDR < 0.05 and fold change > 1.5 . DEGs were mapped to gene ontology (GO) terms by using WEGO  and to pathways by using KOBAS . All microarray data were MIAME compliant. Both raw and processed data were deposited in the NCBI Gene Expression Omnibus [GEO: GSE43327].
The standard curve method was used to measure the relative RNA expression level as previously described [11, 46]. RT-qPCR amplifications were conducted with an MX3000P spectrofluorometric thermal cycler (Stratagene) and a RealMasterMix (SYBR Green) kit (Tiangen). Melting curve analysis was performed to confirm the specificity of amplification. Primer sequences were presented in Additional file 5.
The enzyme activities of pyruvate dehydrogenase complex (PDH) and NADH dehydrogenase (complex I) in locust flight muscle was measured as previously described . Aconitase activity in locust flight muscle submitochondrial particles was measured by monitoring the increasing rate of cis-aconitase at 240 nm . Briefly, the buffer (50 mM Tris, pH 7.4, 0.6 mM MnCl2) was preheated in 30°C for 3 min and then isocitrate (20 mM) was added to the buffer. The mitochondrial protein (final concentration 50 μg · ml-1) was then added to the mixture and the absorbance was monitored for 2 min at 30°C to calculate the increasing absorbance rate of cis-aconitate.
Double-stranded RNA (dsRNA) of green fluorescent protein (GFP) and PDHE1β was prepared using Promega T7 RiboMAX™ Express RNAi System following the manufacturer’s instruction (Additional file 5) [12, 13]. 18 μg (6 μg/μl) of double-stranded GFP (dsGFP) and dsPDHE1β was injected into male adult locusts at the second ventral segment of the abdomen. Three days later the expression levels of PDHE1α and β were both examined by RT-qPCR to evaluate the interfering efficiency and specificity.
The concentrations of AMP, ADP and ATP were determined using previously established method . Briefly, the samples (i.e., muscle homogenate and isolated mitochondria) were first fixed on ice with a 0.5 volume of 2.3 M perchloric acid for 30 min. The insoluble material was removed via centrifugation, and the supernate was neutralized with 2.5 M KHCO3. Afterward, the supernate was filtered through a 0.45 μm membrane (Agilent, USA) and analyzed via high-performance liquid chromatography (HPLC) at 254 nm using an Agilent 1100 HPLC system (Agilent, USA). Protein concentration was determined via BCA protein assay (Pierce, USA) and used to normalize enzymatic activities and concentrations of AMP, ADP and ATP.
North China population
Differentially expressed gene
Partial oxygen pressure.
This research was supported by grants from the National Basic Research Program of China (No. 2012CB114102) and National Natural Science Foundation of China (No. 31210103915 and 30830022). This study was also supported by the Research Network of Computational Biology ‘RNCB’ and the Supercomputing Center at Beijing Institutes of Life Science, Chinese Academy of Sciences.
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.