Genome-wide analysis of transcriptional changes in the thoracic muscle of the migratory locust, Locusta migratoria, exposed to hypobaric hypoxia

Integrated transcriptome and metabolome analysis reveals the molecular responses of Pardosa pseudoannulata to hypoxic environments

Terrestrial organisms are likely to face hypoxic stress during natural disasters such as floods or landslides, which can lead to inevitable hypoxic conditions for those commonly residing within soil. Pardosa pseudoannulata often inhabits soil crevices and has been extensively studied, yet research on its response to hypoxic stress remains unclear. Therefore, we investigated the adaptive strategies of Pardosa pseudoannulata under hypoxic stress using metabolomics and transcriptomics approaches. The results indicated that under hypoxic stress, metabolites related to energy and antioxidants such as ATP, D-glucose 6-phosphate, flavin adenine dinucleotide (FAD), and reduced L-glutathione were significantly differentially expressed. Pathways such as the citric acid (TCA) cycle and oxidative phosphorylation were significantly enriched. Transcriptome analysis and related assessments also revealed a significant enrichment of pathways associated with energy metabolism, suggesting that Pardosa pseudoannulata primarily copes with hypoxic environments by modulating energy metabolism and antioxidant-related substances.

In silico promoter analysis and functional validation identify CmZFH, the co-regulator of hypoxia-responsive genes CmScylla and CmLPCAT

Oxygen (O₂) plays an essential role in aerobic organisms including terrestrial insects. Under hypoxic stress, the cowpea bruchid (Callosobruchus maculatus) ceases feeding and growth. However, larvae, particularly 4th instar larvae exhibit very high tolerance to hypoxia and can recover normal growth once brought to normoxia. To better understand the molecular mechanism that enables insects to cope with low O₂ stress, we performed RNA-seq to distinguish hypoxia-responsive genes in midguts and subsequently identified potential common cis-elements in promoters of hypoxia-induced and -repressed genes, respectively. Selected elements were subjected to gel-shift and transient transfection assays to confirm their cis-regulatory function. Of these putative common cis-elements, AREB6 appeared to regulate the expression of CmLPCAT and CmScylla, two hypoxia-induced genes. CmZFH, the putative AREB6-binding protein, was hypoxia-inducible. Transient expression of CmZFH in Drosophila S2 cells activated CmLPCAT and CmScylla, and their induction was likely through interaction of CmZFH with AREB6. Binding to AREB6 was further confirmed by bacterially expressed CmZFH recombinant protein. Deletion analyses indicated that the N-terminal zinc-finger cluster of CmZFH was the key AREB6-binding domain. Through in silico and experimental exploration, we discovered novel transcriptional regulatory components associated with gene expression dynamics under hypoxia that facilitated insect survival.

Comparative transcriptome analysis reveals molecular strategies of ghost moth Thitarodes armoricanus in response to hypoxia and anoxia

Hypoxia or anoxia greatly impact the survival of many animal species. The ghost moth Thitarodes armoricanus is distributed in the Tibetan Plateau at an average elevation of approximate 4 km above sea level and has probably evolved a superior capacity to tolerate low oxygen levels. In this study, transcriptome analysis using high-throughput RNA-seq revealed common and different adaptation strategies of T. armoricanus in response to hypoxia (11% O₂) or anoxia. T. armoricanus adopted three common strategies for adaptation to hypoxia or anoxia: Up-regulated signal transduction pathways essential for cellular survival, strengthened cell and organelle structure and activity, and activated immune system. Under hypoxia, T. armoricanus might develop a strategy to adapt to hypoxia by suppressing TCA, oxidative phosphorylation pathways, and hypoxanthine catabolism. T. armoricanus larvae kept active under hypoxia but became coma under anoxia, probably relating to up-regulated or suppressed dopamine synthesis pathway. Furthermore, the HIF system seemed not to be essential for regulating the hypoxic and anoxic responses of this insect in Tibetan Plateau. This study provides a global view of gene expression profiles and suggests common and different adaptation strategies of T. armoricanus under hypoxic and anoxic conditions. The results are helpful for understanding the mechanism responsible for the low oxygen level tolerance of this insect species.

Genetic variation in PTPN1 contributes to metabolic adaptation to high-altitude hypoxia in Tibetan migratory locusts

Animal and human highlanders have evolved distinct traits to enhance tissue oxygen delivery and utilization. Unlike vertebrates, insects use their tracheal system for efficient oxygen delivery. However, the genetic basis of insect adaptation to high-altitude hypoxia remains unexplored. Here, we report a potential mechanism of metabolic adaptation of migratory locusts in the Tibetan Plateau, through whole-genome resequencing and functional investigation. A genome-wide scan revealed that the positively selected genes in Tibetan locusts are predominantly involved in carbon and energy metabolism. We observed a notable signal of natural selection in the gene PTPN1, which encodes PTP1B, an inhibitor of insulin signaling pathway. We show that a PTPN1 coding mutation regulates the metabolism of Tibetan locusts by mediating insulin signaling activity in response to hypoxia. Overall, our findings provide evidence for the high-altitude hypoxia adaptation of insects at the genomic level and explore a potential regulatory mechanism underlying the evolved metabolic homeostasis.

Vertebrate adaptation to high-altitude life has been extensively investigated, while invertebrates are less well-studied. Here, the authors find signals of adaptive evolution in genomes of migratory locusts from the Tibetan Plateau, and implicate a PTPN1 coding mutation in their hypoxia response.

Inhibition of mitochondrial respiration under hypoxia and increased antioxidant activity after reoxygenation of Tribolium castaneum

Regulating the air in low-oxygen environments protects hermetically stored grains from storage pests damage. However, pests that can tolerate hypoxic stress pose a huge challenge in terms of grain storage. We used various biological approaches to determine the fundamental mechanisms of Tribolium castaneum to cope with hypoxia. Our results indicated that limiting the available oxygen to T. castaneum increased glycolysis and inhibited the Krebs cycle, and that accumulated pyruvic acid was preferentially converted to lactic acid via anaerobic metabolism. Mitochondrial aerobic respiration was markedly suppressed for beetles under hypoxia, which also might have led to mitochondrial autophagy. The enzymatic activity of citrate synthase decreased in insects under hypoxia but recovered within 12 h, which suggested that the beetles recovered from the hypoxia. Moreover, hypoxia-reperfusion resulted in severe oxidative damage to insects, and antioxidant levels increased to defend against the high level of reactive oxygen species. In conclusion, our findings show that mitochondria were the main target in T. castaneum in response to low oxygen. The beetles under hypoxia inhibited mitochondrial respiration and increased antioxidant activity after reoxygenation. Our research advances the field of pest control and makes it possible to develop more efficient strategies for hermetic storage.

Adaptive evidence of mitochondrial genomes in Dolycoris baccarum (Hemiptera: Pentatomidae) to divergent altitude environments

Given mitochondrion is the 'energy and oxygen usage factories', adaptive signatures of mitochondrial genes have been extensively investigated in vertebrates from different altitudes, but few studies focus on insects. Here, we sequenced the complete mitochondrial genome (mitogenome) of Dolycoris. baccarum living in the Tibetan Plateau (DBHC, ∼3200 m above sea level (asl)) and conducted a detailed comparative analysis with another D. baccarum mitogenome (DBQY) from relatively low altitude (∼1300 m asl). All the 37 mitochondrial genes were highly conserved and under purifying selection, except for two mitochondrial protein-coding genes (MPCGs) (atp6 and nad5) that showed positively selected signatures. We therefore further examined non-synonymous substitutions in atp6 and nad5, by sequencing more individuals from three populations with different altitudes. We found that these non-synonymous substitutions were polymorphic in these populations, likely due to relaxed selection constraints in different altitudes. Purifying selection in all mitochondrial genes may be due to their functional importance for the precision of ATP production usually. Length difference in mitochondrial control regions between DBHC and DBQY was also conversed at the population level, indicating that sequence size adjustments in control region may be associated with adaptation to divergent altitudes. Quantitatively real-time PCR analysis for 12 MPCGs showed that gene expression patterns had a significant change between the two populations, suggesting that expression levels of MPCGs could be modulated by divergent environmental pressures (e.g. oxygen content and ambient temperature). These results provided an important guide for further uncovering genetic mechanisms of ecological adaptation in insects.

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

The impact of metagenomic interplay on the mosquito redox homeostasis

Mosquitoes are exposed to oxidative challenges throughout their life cycle. The primary challenge comes from a blood meal. The blood digestion turns the midgut into an oxidative environment, which imposes pressure not only on mosquito fecundity and other physiological traits but also on the microbiota in the midgut. During evolution, mosquitoes have developed numerous oxidative defense mechanisms to maintain redox homeostasis in the midgut. In addition to antioxidants, SOD, catalase, and glutathione system, sufficient supply of the reducing agent, NADPH, is vital for a successful defense against oxidative stress. Increasing evidence indicates that in response to oxidative stress, cells reconfigure metabolic pathways to increase the generation of NADPH through NADP-reducing networks including the pentose phosphate pathway and others. The microbial homeostasis is critical for the functional contributions to various host phenotypes. The symbiotic microbiota is regulated largely by the Duox-ROS pathway in Drosophila. In mosquitoes, Duox-ROS pathway, heme-mediated signaling, antimicrobial peptide production and C-type lectins work in concert to maintain the dynamic microbial community in the midgut. Microbial mechanisms against oxidative stress in this context are not well understood. Emerging evidence that microbial metabolites trigger host oxidative response warrants further study on the metagenomic interplay in an oxidative environment like mosquito gut ecosystem. Besides the classical Drosophila model, hematophagous insects like mosquitoes provide an alternative model system to study redox homeostasis in a symbiotic metagenomic context.

Characterization and analysis of a de novo transcriptome from the pygmy grasshopper Tetrix japonica

The pygmy grasshopper Tetrix japonica is a common insect distributed throughout the world, and it has the potential for use in studies of body colour polymorphism, genomics and the biology of Tetrigoidea (Insecta: Orthoptera). However, limited biological information is available for this insect. Here, we conducted a de novo transcriptome study of adult and larval T. japonica to provide a better understanding of its gene expression and develop genomic resources for future work. We sequenced and explored the characteristics of the de novo transcriptome of T. japonica using Illumina HiSeq 2000 platform. A total of 107 608 206 paired-end clean reads were assembled into 61 141 unigenes using the trinity software; the mean unigene size was 771 bp, and the N50 length was 1238 bp. A total of 29 225 unigenes were functionally annotated to the NCBI nonredundant protein sequences (Nr), NCBI nonredundant nucleotide sequences (Nt), a manually annotated and reviewed protein sequence database (Swiss-Prot), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. A large number of putative genes that are potentially involved in pigment pathways, juvenile hormone (JH) metabolism and signalling pathways were identified in the T. japonica transcriptome. Additionally, 165 769 and 156 796 putative single nucleotide polymorphisms occurred in the adult and larvae transcriptomes, respectively, and a total of 3162 simple sequence repeats were detected in this assembly. This comprehensive transcriptomic data for T. japonica will provide a usable resource for gene predictions, signalling pathway investigations and molecular marker development for this species and other pygmy grasshoppers.

De novo transcriptome characterization of the ghost moth, Thitarodes pui, and elevation-based differences in the gene expression of its larvae

Thitarodes pui larvae are the hosts of a medicinal fungus, Ophiocordyceps sinensis, and are naturally distributed at an altitude of 4100-4650 m on Segrila Mountain of the Tibetan Plateau. Here, we conducted transcriptome profiling of T. pui larvae based on the Illumina high-throughput sequencing platform. Subsequently, we explored elevation-based differences of T. pui by comparing gene expression profiles between larvae at high-altitude (natural conditions) and larvae exposed to short-term (2months) low-altitude conditions. A total of 105,935,208 clean reads were assembled into 70,048 unigenes with a mean length of 639 bp. All unigenes were searched against public databases, with 51.26% unigenes being successfully annotated in the NR, SWISS-PROT, EuKaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genome (KEGG) databases. A total of 11,846 unigenes were functionally classified into 239 KEGG pathways. Metabolism was the most represented pathway, with 4271 unigenes (36.05%). Using the transcriptome data as a reference, 629 (311 up-regulated/318 down-regulated) genes were differentially expressed by low-altitude larvae when compared with those of high-altitude larvae. The most significantly differentially expressed genes were annotated in the processes of carbohydrate metabolism, lipid metabolism, and respiration. This report provides valuable information about the T. pui transcriptome for future genomic studies, including how gene expression is altered in larvae reared at different elevations.

Acute high-altitude hypoxic brain injury: Identification of ten differential proteins

Hypobaric hypoxia can cause severe brain damage and mitochondrial dysfunction, and is involved in hypoxic brain injury. However, little is currently known about the mechanisms responsible for mitochondrial dysfunction in hypobaric hypoxic brain damage. In this study, a rat model of hypobaric hypoxic brain injury was established to investigate the molecular mechanisms associated with mitochondrial dysfunction. As revealed by two-dimensional electrophoresis analysis, 16, 21, and 36 differential protein spots in cerebral mitochondria were observed at 6, 12, and 24 hours post-hypobaric hypoxia, respectively. Furthermore, ten protein spots selected from each hypobaric hypoxia subgroup were similarly regulated and were identified by mass spectrometry. These detected proteins included dihydropyrimidinase-related protein 2, creatine kinase B-type, isovaleryl-CoA dehydrogenase, elongation factor Ts, ATP synthase beta-subunit, 3-mercaptopyruvate sulfurtransferase, electron transfer flavoprotein alpha-subunit, Chain A of 2-enoyl-CoA hydratase, NADH dehydrogenase iron-sulfur protein 8 and tropomyosin beta chain. These ten proteins are all involved in the electron transport chain and the function of ATP synthase. Our findings indicate that hypobaric hypoxia can induce the differential expression of several cerebral mitochondrial proteins, which are involved in the regulation of mitochondrial energy production.

Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia

Background

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.

Results

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.

Conclusions

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.