Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia
BMC Genomics volume 14, Article number: 631 (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 (p O2) in Tibet (~13 kPa) is above the critical oxygen pressure of locusts (~3 kPa) [9, 10], the Tibetan p O2 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.
Phenotypic differences between NP and TP locusts
The migratory locust Locusta migratoria is widely distributed in China, both on the low-altitude North China Plain (NP) and on the high-altitude Tibetan Plateau (TP) (Figure 1A). The natural habitats of NP and TP locusts show great differences, one of which is the partial oxygen pressure (p O2) which decreases around 40% in Tibet compared to that in North China (Figure 1B). To test whether NP and TP locusts have evolved different sensitivity to hypoxia, we compared their phenotypic traits from morphological, behavioral and physiological perspectives. First, we compared their body size by using femur length as an index. The average femur length of male and female NP locusts was 21.59 mm and 25.21 mm; 13.7% and 17% larger than that of male and female TP locusts, respectively (Figure 2A). Second, to directly compare the capacity of NP and TP locusts to tolerate hypoxia, we measured the ratio of locusts exhibiting stupor by subjecting them to five different levels of extreme hypoxia. When p O2 decreased from 2.8 kPa to 1.2 kPa, the stupor ratio of NP locusts rose from 52% to 96% while that of TP locusts increased from 6% to 84%. Under all the tested oxygen levels, the stupor ratios of TP locusts were smaller than those of NP locusts, with the biggest difference at 1.6 kPa p O2 (Figure 2B). Third, to evaluate the effects of hypoxic injury on TP and NP locusts, we compared the time required to recover from stupor. It took TP locusts 4.2 min to recover from stupor, significantly less than 6.2 min required by NP locusts (Figure 2C). Finally, to evaluate the effect of hypoxia on aerobic metabolism, we compared metabolic rates under a series of p O2 (1.6 kPa, 3.2 kPa, 6.5 kPa, 13 kPa and 21 kPa) by measuring their CO2 production rates. Locusts did not show significant differences in metabolic rates when the p O2 was equal to or above 3.2 kPa. When p O2 reached 1.6 kPa, CO2 production rates of both populations plummeted drastically; however, the TP locusts still showed higher CO2 production rates than the NP locusts (Figure 2D). Taken together, these phenotypic differences suggest that TP locusts have evolved a better capacity to cope with hypoxia than NP locusts.
Gene expression profiles of field NP and TP locusts
To explore how NP and TP locusts diverged at the transcriptional level after adapting to different habitats, we compared the gene expression of thoracic muscles between the field populations. We identified 324 differentially expressed genes (DEGs) up-regulated in TP locusts and 555 DEGs up-regulated in NP locusts (Additional file 1). Of the 879 genes, only 24 were among the 235 previously identified target genes of hypoxia-inducible factor . Based on GO, KEGG and NCBI annotations, DEGs were further classified into 14 functional categories (Figure 3, Additional file 1), which include a variety of basic functions, e.g. cellular structure, lipid metabolism, antioxidant processes, DNA repair and gene expression. In the field-caught TP locusts, structural genes involved in the assembly of cytoskeleton, muscles and cuticle were suppressed (Figure 3; Category 1, Additional file 1), such as tubulin (LM04131, LM00737, LM03863, Additional file 1), tropomyosin (LM00621, Additional file 1) and endocuticle structural glycoprotein (LM00591, LM00412, Additional file 1). Genes related to fatty acid synthesis and β oxidation were also suppressed in TP locusts (Figure 3; Category 3, Additional file 1), such as acetyl-coenzyme A synthetase (LM06760 and LM06776, Additional file 1) which activates acetate so that it can be used for lipid synthesis  and delta(3,5)-delta(2,4)-dienoyl-CoA isomerase (LM01097, Additional file 1) which catalyzes β oxidation of unsaturated fatty acids . Besides, the genes in the antioxidant systems were suppressed in TP locusts as well (Figure 3; Category 7, Additional file 1), e.g. superoxide dismutase (LM03935, Additional file 1), peroxiredoxin (LM00341, Additional file 1), and glutaredoxin (LM00278, Additional file 1). Genes involved in genetic information processing outnumbered any other categories (Figure 3; Category 11, Additional file 1). These genes were involved in DNA repair, transcription, translation, and protein degradation. Of particular interest, two genes involved in repairing UV-damaged DNA were up-regulated in TP locusts, namely ubiquitin-conjugating enzyme E2 (LM02929, Additional file 1) which is required for postreplication repair of UV-damaged DNA  and UV excision repair protein RAD23 homolog B (LM02985, Additional file 1) which is involved in global genome nucleotide excision repair by acting as component of the XPC complex . In contrast, a gene involved in recombination repair was up-regulated in NP locusts (LM03194, Additional file 1), which promotes homologous recombination at sites of DNA damage . The gene expression profiles showed that multiple factors contributed to the transcriptional differences between the field NP and TP locusts.
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.
Gene expression profiles of laboratory NP and TP locusts under extreme hypoxia
NP and TP locusts responded differently to extreme hypoxia (Figure 2B-D). To understand what contributed to the superior tolerance of extreme hypoxia for TP locusts, we compared the transcriptional responses of their thoracic muscles to extreme hypoxia (1.6 kPa p O2). We identified 247 hypoxia-responsive genes from NP locusts and only 111 from TP locusts (Figure 4A, Additional files 3 and 4). The expression levels of 10 genes, including both DEGs and non-DEGs, were confirmed by real-time quantitative PCR (RT-qPCR) (Additional files 5 and 6). The genes were classified into 14 gene ontology (GO) categories. Metabolic processes had the largest number of DEGs and the biggest difference in DEGs number between the NP and TP locusts (Figure 4B). Further investigation of the DEGs in metabolic processes showed that these genes were mainly involved in glycolysis, the citrate cycle, and the electron transfer chain (Figure 4C, D). Genes involved in anaerobic metabolism were up-regulated in both NP and TP locusts (Figure 4C, D), including lactate dehydrogenase (LDH), glucose transporter (GluT), hexokinase type 2 (HK2), phosphoglycerate kinase (PGK), pyruvate kinase (PK), and 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 4 (PFK). In contrast, genes involved in aerobic processes were down-regulated only in NP locusts, including complex I subunit B13, pyruvate dehydrogenase (PDH) E1 component subunit beta and aconitase (Figure 4C).
Role of aerobic metabolism in the tolerance of extreme hypoxia
Based on the facts that TP locusts showed higher aerobic respiration rates than the NP locusts under extreme hypoxia (Figure 2D) and the aerobic metabolic processes were not as severely suppressed in TP locusts as in NP locusts (Figure 4C-D), we speculate that the aerobic metabolism plays an important role in the adaptation of TP locusts to extreme hypoxia. We first examined the gene expression levels of three DEGs under extreme hypoxia: PDH E1 component subunit beta, aconitase, and complex I subunit B13. Consistent with microarray results, the three genes were all suppressed in NP locusts and not significantly affected in TP locusts (Figure 5A). We then measured the enzymatic activities of PDH complex, aconitase, and complex I. We found in NP locusts that enzymatic activities of PDH complex and aconitase were reduced under extreme hypoxia while that of complex I were not significantly altered (Figure 5B). In TP locusts, the activity of aconitase was also reduced, but it remained higher than that in NP locusts; similarly, the activity of PDH complex was also higher in TP locusts than in NP locusts under extreme hypoxia (Figure 5B).
To further confirm our speculation, we blocked the aerobic metabolic branch by disrupting PDHE1β using RNAi , an entry gene from glycolysis to TCA cycle (Figure 4D). The expression level of PDHE1β was drastically reduced by 20-fold while expression of PDHE1α, partner of PDHE1β, was not significantly affected in both TP and NP locusts (Figure 6A). A behavioral test under hypoxia showed that the stupor ratio of TP locusts with PDHE1β disrupted was 33.3% higher than that in control locusts (Figure 6B). To further evaluate the effect of RNAi, we measured the contents of AMP, ADP, and ATP, and calculated the energy charge (EC), which is an indicator of cellular energy status. The contents of these molecules did not change significantly in both TP and NP locusts under normoxia after RNAi (Additional file 7A-C). However, when locusts were placed under extreme hypoxia (1.6 kPa p O2), the contents of AMP, ADP, and ATP changed significantly in TP locusts after RNAi compared to control (dsGFP). When PDHE1β was disrupted in TP, AMP was 1.7 nmol/mg higher than control, ADP 2.6 nmol/mg lower, and ATP 1.9 nmol/mg higher (Additional file 7D-F). Under normoxia, RNAi did not alter EC significantly after disrupting PDHE1β (Figure 6C). However, under hypoxia, EC values of NP and TP locusts decreased in both control and treated locusts; TP locusts with PDHE1β disrupted had lower EC values than the control (Figure 6C). Taken together, gene expression, enzymatic activities, and RNAi data support that aerobic metabolism plays an important role in the adaptation of TP locusts to extreme hypoxia.
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 p O2 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 p O2 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 p O2 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.
Microarray experiments and data analysis
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 p O2 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].
Real-time quantitative PCR (RT-qPCR)
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.
Measurement of AMP, ADP, and ATP using HPLC
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.
Scott GR, Schulte PM, Egginton S, Scott ALM, Richards JG, Milsom WK: Molecular evolution of cytochrome c oxidase underlies high-altitude adaptation in the bar-headed goose. Mol Biol Evol. 2011, 28: 351-363. 10.1093/molbev/msq205.
Storz JF, Sabatino SJ, Hoffmann FG, Gering EJ, Moriyama H, Ferrand N, Monteiro B, Nachman MW: The molecular basis of high-altitude adaptation in deer mice. PLoS Genet. 2007, 3: e45-10.1371/journal.pgen.0030045.
Peng Y, Yang Z, Zhang H, Cui C, Qi X, Luo X, Tao X, Wu T, Ouzhuluobu , Basang , Ciwangsangbu , Danzengduojie , Chen H, Shi H, Su B: Genetic variations in Tibetan populations and high altitude adaptation at the Himalayas. Mol Biol Evol. 2011, 28: 1075-1081. 10.1093/molbev/msq290.
Simonson TS, Yang Y, Huff CD, Yun H, Qin G, Witherspoon DJ, Bai Z, Lorenzo FR, Xing J, Jorde LB, Prchal JT, Ge R: Genetic evidence for high-altitude adaptation in Tibet. Science. 2010, 329: 72-75. 10.1126/science.1189406.
Yi X, Liang Y, Huerta-Sanchez E, Jin X, Cuo ZXP, Pool JE, Xu X, Jiang H, Vinckenbosch N, Korneliussen TS, Zheng H, Liu T, He W, Li K, Luo R, Nie X, Wu H, Zhao M, Cao H, Zou J, Shan Y, Li S, Yang Q, Asan , Ni P, Tian G, Xu J, Liu X, Jiang T, Wu R, Zhou G, Tang M, Qin J, Wang T, Feng S, Li G, Huasang , Luosang J, Wang W, Chen F, Wang Y, Zheng X, Li Z, Bianba Z, Yang G, Wang X, Tang S, Gao G, Chen Y, Luo Z, Gusang L, Cao Z, Zhang Q, Ouyang W, Ren X, Liang H, Zheng H, Huang Y, Li J, Bolund L, Kristiansen K, Li Y, Zhang Y, Zhang X, Li R, Li S, Yang H, Nielsen R, Wang J, Wang J: Sequencing of 50 human exomes reveals adaptation to high altitude. Science. 2010, 329: 75-78. 10.1126/science.1190371.
Ma C, Yang P, Jiang F, Chapuis M-P, Shali Y, Sword GA, Kang LE: Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust. Mol Ecol. 2012, 21: 4344-4358. 10.1111/j.1365-294X.2012.05684.x.
Chen YL: A new subspecies of Locusta migratoria, Tibetan migratory locust (Locusta migratoria tibetensis subsp.n.). Acta Ent Sin. 1963, 12: 463-475.
Zhang D-X, Yan L-N, Ji Y-J, Hewitt G, Huang Z-S: Unexpected relationships of substructured populations in Chinese Locusta migratoria. BMC Evol Biol. 2009, 9: 144-10.1186/1471-2148-9-144.
Arieli R, Lehrer C: Recording of locust breathing frequency by barometric method exemplified by hypoxic exposure. J Insect Physiol. 1988, 34: 325-10.1016/0022-1910(88)90143-6.
Greenlee KJ, Harrison JF: Development of respiratory function in the American locust Schistocerca americana I. Across-instar effects. J Exp Biol. 2004, 207: 497-508. 10.1242/jeb.00767.
Zhao DJ, Zhang ZY, Harrison J, Kang L: Genome-wide analysis of transcriptional changes in the thoracic muscle of the migratory locust, Locusta migratoria, exposed to hypobaric hypoxia. J Insect Physiol. 2012, 58: 1424-1431. 10.1016/j.jinsphys.2012.08.006.
Ma Z, Guo W, Guo X, Wang X, Kang L: Modulation of behavioral phase changes of the migratory locust by the catecholamine metabolic pathway. Proc Natl Acad Sci U S A. 2011, 108: 3882-3887. 10.1073/pnas.1015098108.
Guo W, Wang X, Ma Z, Xue L, Han J, Yu D, Kang L: CSP and Takeout genes modulate the switch between attraction and repulsion during behavioral phase change in the migratory locust. PLoS Genet. 2011, 7: e1001291-10.1371/journal.pgen.1001291.
Wang H, Ma Z, Cui F, Wang X, Guo W, Lin Z, Yang P, Kang L: Parental phase status affects the cold hardiness of progeny eggs in locusts. Funct Ecol. 2012, 26: 379-389. 10.1111/j.1365-2435.2011.01927.x.
Gilad Y, Oshlack A, Rifkin SA: Natural selection on gene expression. Trends Genet. 2006, 22: 456-461. 10.1016/j.tig.2006.06.002.
Ranz JM, Machado CA: Uncovering evolutionary patterns of gene expression using microarrays. Trends Ecol Evol. 2006, 21: 29-37. 10.1016/j.tree.2005.09.002.
Luong A, Hannah VC, Brown MS, Goldstein JL: Molecular characterization of human acetyl-CoA synthetase, an enzyme regulated by sterol regulatory element-binding proteins. J Biol Chem. 2000, 275: 26458-26466. 10.1074/jbc.M004160200.
Gurvitz A, Wabnegger L, Yagi AI, Binder M, Hartig A, Ruis H, Hamilton B, Dawes IW, Hiltunen JK, Rottensteiner H: Function of human mitochondrial 2,4-dienoyl-CoA reductase and rat monofunctional Delta3-Delta2-enoyl-CoA isomerase in beta-oxidation of unsaturated fatty acids. Biochem J. 1999, 344: 903-914. 10.1042/0264-6021:3440903.
Koken M, Reynolds P, Bootsma D, Hoeijmakers J, Prakash S, Prakash L: Dhr6, a Drosophila homolog of the yeast DNA-repair gene RAD6. Proc Natl Acad Sci U S A. 1991, 88: 3832-3836. 10.1073/pnas.88.9.3832.
Sugasawa K, Ng JM, Masutani C, Maekawa T, Uchida A, van der Spek PJ, Eker AP, Rademakers S, Visser C, Aboussekhra A, Wood RD, Hanaoka F, Bootsma D, Hoeijmakers JH: Two human homologs of Rad23 are functionally interchangeable in complex formation and stimulation of XPC repair activity. Mol Cell Biol. 1997, 17: 6924-6931.
Sander M, Lowenhaupt K, Rich A: Drosophila Rrp1 protein: an apurinic endonuclease with homologous recombination activities. Proc Natl Acad Sci U S A. 1991, 88: 6780-6784. 10.1073/pnas.88.15.6780.
Higuchi A, Takahashi K, Hirashima M, Kawakita T, Tsubota K: Selenoprotein P controls oxidative stress in Cornea. PLoS ONE. 2010, 5: e9911-10.1371/journal.pone.0009911.
Shimazu F, Tappel AL: Selenoamino acids: decrease of radiation damage to amino acids and proteins. Science. 1964, 143: 369-371. 10.1126/science.143.3604.369.
He Z-B, Cao Y-Q, Yin Y-P, Wang Z-K, Chen B, Peng G-X, Xia Y-X: Role of hunchback in segment patterning of Locusta migratoria manilensis revealed by parental RNAi. Dev Growth Differ. 2006, 48: 439-445. 10.1111/j.1440-169X.2006.00881.x.
Klok CJ, Harrison JF: Atmospheric hypoxia limits selection for large body size in insects. PLoS ONE. 2009, 4: e3876-10.1371/journal.pone.0003876.
Zhou D, Xue J, Chen J, Morcillo P, Lambert JD, White KP, Haddad GG: Experimental selection for Drosophila survival in extremely low O2 environment. PLoS ONE. 2007, 2: e490-10.1371/journal.pone.0000490.
Krishnan SN, Sun Y-A, Mohsenin A, Wyman RJ, Haddad GG: Behavioral and electrophysiologic responses of Drosophila melanogaster to prolonged periods of anoxia. J Insect Physiol. 1997, 43: 203-210. 10.1016/S0022-1910(96)00084-4.
Haddad GG, Sun Y, Wyman RJ, Xu T: Genetic basis of tolerance to O2 deprivation in Drosophila melanogaster. Proc Natl Acad Sci U S A. 1997, 20: 10809-10812.
Contreras HL, Bradley TJ: Metabolic rate controls respiratory pattern in insects. J Exp Biol. 2009, 212: 424-428. 10.1242/jeb.024091.
Rastogi RP, Richa , Kumar A, Tyagi MB, Sinha RP: Molecular mechanisms of ultraviolet radiation-induced DNA damage and repair. J Nucleic Acids. 2010, 2010: 592980-
Mandal SB, Sen PC, Chakrabarti P, Sen K: Effect of respiratory deficiency and temperature on the mitochondrial lipid metabolism of Aspergillus niger. Can J Microbiol. 1978, 24: 586-592. 10.1139/m78-095.
Francis D, Barlow PW: Temperature and the cell cycle. Symp Soc Exp Biol. 1988, 42: 181-201.
Feala JD, Coquin L, McCulloch AD, Paternostro G: Flexibility in energy metabolism supports hypoxia tolerance in Drosophila flight muscle: metabolomic and computational systems analysis. Mol Syst Biol. 2007, 3: 99-
Feala J, Coquin L, Zhou D, Haddad G, Paternostro G, McCulloch A: Metabolism as means for hypoxia adaptation: metabolic profiling and flux balance analysis. BMC Syst Biol. 2009, 3: 91-10.1186/1752-0509-3-91.
Fukuda R, Zhang H, Kim J-w, Shimoda L, Dang CV, Semenza GL: HIF-1 regulates cytochrome oxidase subunits to optimize efficiency of respiration in hypoxic cells. Cell. 2007, 129: 111-122. 10.1016/j.cell.2007.01.047.
Brown RM, Head RA, Boubriak II, Leonard JV, Thomas NH, Brown GK: Mutations in the gene for the E1b subunit: a novel cause of pyruvate dehydrogenase deficiency. Hum Genet. 2004, 115: 123-127.
Ciszak EM, Korotchkina LG, Dominiak PM, Sidhu S, Patel MS: Structural basis for flip-flop action of thiamin pyrophosphate-dependent enzymes revealed by human pyruvate dehydrogenase. J Biol Chem. 2003, 278: 21240-21246. 10.1074/jbc.M300339200.
Wexler ID, Kerr DS, Ho L, Lusk MM, Pepin RA, Javed AA, Mole JE, Jesse BW, Thekkumkara TJ, Pons G: Heterogeneous expression of protein and mRNA in pyruvate dehydrogenase deficiency. Proc Natl Acad Sci U S A. 1988, 85: 7336-7340. 10.1073/pnas.85.19.7336.
Ge R-L, Cai Q, Shen Y-Y, San A, Ma L, Zhang Y, Yi X, Chen Y, Yang L, Huang Y, He R, Hui Y, Hao M, Li Y, Wang B, Ou X, Xu J, Zhang Y, Wu K, Geng C, Zhou W, Zhou T, Irwin DM, Yang Y, Ying L, Bao H, Kim J, Larkin DM, Ma J, Lewin HA: Draft genome sequence of the Tibetan antelope. Nat Commun. 2013, 4: 1858-
Nguyen MP, Lee S, Lee YM: Epigenetic regulation of hypoxia inducible factor in diseases and therapeutics. Arch Pharm Res. 2013, 36: 252-263. 10.1007/s12272-013-0058-x.
Tsai YP, Wu KJ: Epigenetic regulation of hypoxia-responsive gene expression: Focusing on chromatin and DNA modifications. Int J Cancer. 2013, DOI: 10.1002/ijc.28190
Watson JA, Watson CJ, McCann A, Baugh J: Epigenetics, the epicenter of the hypoxic response. Epigenetics. 2010, 5: 293-296. 10.4161/epi.5.4.11684.
Kang L, Chen X, Zhou Y, Liu B, Zheng W, Li R, Wang J, Yu J: The analysis of large-scale gene expression correlated to the phase changes of the migratory locust. Proc Natl Acad Sci U S A. 2004, 101: 17611-17615. 10.1073/pnas.0407753101.
Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34: W293-W297. 10.1093/nar/gkl031.
Wu J, Mao X, Cai T, Luo J, Wei L: KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006, 34: W720-W724. 10.1093/nar/gkl167.
Zhao D-J, Guo K, Kang L: Identification of condition-specific reference genes from microarray data for locusts exposed to hypobaric hypoxia. FEBS Open Bio. 2012, 2: 235-240.
Zhang Z-Y, Chen B, Zhao D-J, Kang L: Functional modulation of mitochondrial cytochrome c oxidase underlies adaptation to high-altitude hypoxia in a Tibetan migratory locust. P Roy Soc Lond B Bio. 2013, 280: 20122758-10.1098/rspb.2012.2758.
Krebs HA, Holzach O: The conversion of citrate into cis-aconitate and isocitrate in the presence of aconitase. Biochem J. 1952, 52: 527-528.
Schweinsberg PD, Loo TL: Simultaneous analysis of ATP, ADP, AMP, and other purines in human erythrocytes by high-performance liquid chromatography. J Chromatogr. 1980, 181: 103-107.
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.
The authors declare that they have no competing interests.
DZ and LK designed the study. DZ, AC and JH designed and performed the respiration experiments. DZ performed the phenotypic assays, microarray experiments, QPCR and RNAi. DZ and ZZ performed measurement of enzymatic activities and HPLC experiments. DZ drafted the manuscript. LK, JH, AC and ZZ revised the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Differentially expressed genes between the field populations of locusts (TP/NP). The genes were classified into 14 categories according their functions. The possible targets of hypoxia-inducible factor were indicated with asterisk (*). (XLS 130 KB)
Additional file 2:Significantly enhanced pathways in field-collected TP locusts. These pathways were identified using KO-Based Annotation System (KOBAS) with all the genes on the locust microarrays as background. (XLS 24 KB)
Additional file 3:Hypoxia-responsive genes in NP locusts. NP locusts were subjected to 1.6 kPa p O2 for two hours and then thoracic muscles were sampled. The experiments were done with six biological replicates. (XLS 50 KB)
Additional file 4:Hypoxia-responsive genes in TP locusts. TP locusts were subjected to 1.6 kPa p O2 for two hours and then thoracic muscles were sampled. The experiments were done with six biological replicates. (XLS 36 KB)
Additional file 6:RT-qPCR validation of microarray results. Ten genes, including non-differentially expressed genes, were randomly selected for validating the gene expression levels. The fold changes of six biological replicates were averaged and consistent with the microarray results. (EPS 650 KB)
Additional file 7:Concentration of AMP, ADP and ATP in locust flight muscle. The concentrations of AMP, ADP and ATP show no difference between TP and NP locusts under normoxia (A-C); however, they showed difference between TP and NP locusts under extreme hypoxia (D-F) (*p < 0.05, t-test, n = 6). (EPS 565 KB)
Additional file 8:Multiple sequence alignment of PDHE1β. The full coding sequence from 15 individuals of TP and NP locusts were sequenced respectively. The nucleotide sequences were translated into peptides using ORF Finder and were aligned with ClustalX. Positions with the same amino acid in all the individuals were labeled with asterisk (*). (PDF 8 KB)
Additional file 9:Diagram of respirometry setup. Gas stream is denoted in thick lines with arrows indicating flow direction. DR, drierite, to remove water vapor; AS, ascarite (Sigma-Aldrich) to scrub CO2; FM, flowmeter, to control the flow rate to 120 ml/min; LI-COR, an equipment to quantify CO2, connected to UI relaying data to a computer; pump, to pull the air through the respirometry at a constant speed of 100 ml/min; S1, syringe No.1, its plunger is pulled out and then it is plugged with cotton ball; S2, syringe No.2, a locust individual is positioned inside S2 and the air in S2 is evacuated by the gas from the air tank to treat the locust in the specific oxygen; S3, syringe No.3, the respired gas was injected into the respirometry system to quantify the produced CO2 using LI-COR. (EPS 819 KB)
Authors’ original submitted files for images
About this article
Cite this article
Zhao, D., Zhang, Z., Cease, A. et al. Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia. BMC Genomics 14, 631 (2013). https://doi.org/10.1186/1471-2164-14-631