Hypoxia inducible factor -1 regulates WSSV-induced glycolytic genes in the white shrimp Litopenaeus vannamei

Comprehensive review on glucose 6 phosphate dehydrogenase: A critical immunometabolic and redox switch in insects

Mounting an active immune response is energy intensive and demands the reallocation of nutrients to maintain the body's resistance and tolerance against infections. Central to this metabolic adaptation is Glucose-6-phosphate dehydrogenase (G6PDH), a housekeeping enzyme involve in pentose phosphate pathway (PPP). PPP play an essential role in generating ribose, which is critical for nicotinamide adenine dinucleotide phosphate (NADPH). It is vital for physiological and cellular processes such as generating nucleotides, fatty acids and reducing oxidative stress. The G6PDH is extremely conserved enzyme across species in PP shunt. The deficiency of enzymes leads to serious consequences on organism, particularly on adaptation and development. Acute deficiency can lead to impaired cell development, halted embryonic growth, reduce sensitivity to insulin, hypertension and increase inflammation. Historically, research focusing on G6PDH and PPP have primarily targeted diseases on mammalian. However, our review has investigated the unique functions of the G6PDH enzyme in insects and greatly improved mechanistic understanding of its operations. This review explore how G6PDH in insects plays a crucial role in managing the redox balance and immune related metabolism. This study aims to investigate the enzyme's role in different metabolic adaptations.

WSSV early protein WSSV004 enhances viral replication by suppressing LDH activity

White spot syndrome virus (WSSV) is known to upregulate glycolysis to supply biomolecules and energy for the virus's replication. At the viral genome replication stage, lactate dehydrogenase (LDH), a glycolytic enzyme, shows increased activity without any increase in expression. In the present study, yeast 2-hybrid screening was used to identify WSSV proteins that interacted with LvLDH isoform 1 and 2, and these included the WSSV early protein WSSV004. The interaction between WSSV004 and LvLDH1/2 was confirmed by co-immunoprecipitation. Immunofluorescence showed that WSSV004 co-localized with LvLDH1/2 in the cytoplasm. dsRNA silencing experiments showed that WSSV004 was crucial for WSSV replication. However, although WSSV004 silencing led to the suppression of total LvLDH gene expression during the viral late stage, there was nevertheless a significant increase in LvLDH activity at this time. We also used affinity purification-mass spectrometry to identify cellular proteins that interact with WSSV004, and found a total of 108 host proteins and 3 WSSV proteins with which it potentially interacts. Bioinformatics analysis revealed that WSSV004 and its interacting proteins might be responsible for various biological pathways during infection, including vesicular transport machinery and RNA-related functions. Collectively, our study suggests that WSSV004 serves as a multifunctional modulator to facilitate WSSV replication.

Impact of hypoxia on glucose metabolism and hypoxia signaling pathways in juvenile horseshoe crabs Tachypleus tridentatus

Marine hypoxia poses a significant challenge in the contemporary marine environment. The horseshoe crab, an ancient benthic marine organism, is confronted with the potential threat of species extinction due to hypoxia, making it an ideal candidate for studying hypoxia tolerance mechanisms. In this experiment, juvenile Tachypleus tridentatus were subjected to a 21-day trial at DO:2 mg/L (hypoxia) and DO:6 mg/L conditions. The experimental timeline included a 14-day exposure phase followed by a 7-day recovery period. Sampling occurred on days 0, 7, 14, and 21, where the period from day 14 to day 21 corresponds to seven days of recuperation. Several enzymatic activities of important proteins throughout this investigation were evaluated, such as succinate dehydrogenase (SDH), phosphofructokinase (PFK), hexokinase (HK), lactate dehydrogenase (LDH), and pyruvate kinase (PK). Concurrently, the relative expression of hexokinase-1 (HK), hypoxia-inducible factor 1-alpha inhibitor (FIH), and hypoxia-inducible factor 1-alpha (HIF-1α), pyruvate dehydrogenase phosphatase (PDH), succinate dehydrogenase assembly factor 4 (SDH), and Glucose-6-phosphatase (G6Pase) were also investigated. These analyses aimed to elucidate alterations in the hypoxia signaling pathway and respiratory energy metabolism. It is revealed that juvenile T. tridentatus initiated the HIF pathway under hypoxic conditions, resulting in an upregulation of HIF-1α and FIH-1 gene expression, which in turn, influenced a shift in metabolic patterns. Particularly, the activity of glycolysis-related enzymes was promoted significantly, including PK, HK, PKF, LDH, and the related HK gene. In contrast, enzymes linked to aerobic respiration, PDH, and SDH, as well as the related PDH and SDH genes, displayed down-regulation, signifying a transition from aerobic to anaerobic metabolism. Additionally, the activity of gluconeogenesis-related enzymes such as PK and G6Pase gene expression were significantly elevated, indicating the activation of gluconeogenesis and glycogenolysis pathways. Consequently, juvenile T. tridentatus demonstrated an adaptive response to hypoxic conditions, marked by changes in respiratory energy metabolism modes and the activation of hypoxia signaling pathways.

HIF-1-mediated regulation of LDH gene unravels key insights into MCDV-1 pathogenesis in mud crabs Scylla paramamosain

Hypoxia-inducible factors -1 (HIF-1) is a crucial transcription factor that regulates the expression of glycolytic genes. Our previous study proved that the Mud crab dicistrovirus-1 (MCDV-1) can induce aerobic glycolysis that favors viral replication in mud crab Scylla paramamosain. However, the role of HIF-1 on key glycolytic genes during the MCDV-1 infection has not been examined. In this study, the intricate interplay between HIF-1 and the key glycolysis enzyme, lactate dehydrogenase (LDH), was investigated after MCDV-1 infection. The expression of LDH was significant increased after MCDV-1 infection. Additionally, the expression of HIF-1α was upregulated following MCDV-1 infection, potentially attributed to the downregulation of prolyl hydroxylase domains 2 expression. Subsequent examination of the SpLDH promoter identified the presence of hypoxia response elements (HREs), serving as binding sites for HIF-1α. Intriguingly, experimental evidence demonstrated that SpHIF-1α actively promotes SpLDH transcription through these HREs. To further elucidate the functional significance of SpHIF-1α, targeted silencing was employed, resulting in a substantial reduction in SpLDH expression, activity, and lactate concentrations in MCDV-1-infected mud crabs. Notably, SpHIF-1α-silenced mud crabs exhibited higher survival rates and lower viral loads in hepatopancreas tissues following MCDV-1 infection. These results highlight the critical role of SpHIF-1α in MCDV-1 pathogenesis by regulating LDH gene dynamics, providing valuable insights into the molecular mechanisms underlying the virus-host interaction.

Respiratory syncytial virus co-opts hypoxia-inducible factor-1α-mediated glycolysis to favor the production of infectious virus

IMPORTANCE

Respiratory syncytial virus (RSV) is the leading etiological agent of lower respiratory tract illness. However, efficacious vaccines or antiviral drugs for treating RSV infections are currently not available. Indeed, RSV depends on host cells to provide energy needed to produce progeny virions. Glycolysis is a series of oxidative reactions used to metabolize glucose and provide energy to host cells. Therefore, glycolysis may be helpful for RSV infection. In this study, we show that RSV increases glycolysis by inducing the stabilization, transcription, translation, and activation of hypoxia-inducible factor (HIF)-1α in infected cells, which is important for the production of progeny RSV virions. This study contributes to understanding the molecular mechanism by which HIF-1α-mediated glycolysis controls RSV infection and reveals an effective target for the development of highly efficient anti-RSV drugs.

The pyruvate kinase of the whiteleg shrimp Litopenaeus vannamei: Gene structure and responses to short term hypoxia

The shrimp Litopenaeus vannamei is the main farmed crustaceans worldwide. This crustacean suffers environmental changes in oxygen availability that affect its energy metabolism. Pyruvate kinase (PK) catalyzes the last reaction of glycolysis and is key for the regulation of glycolysis and gluconeogenesis. There is ample knowledge about mammalian PK, but in crustaceans, the information is very scarce. In this study, we analyzed in silico the structures of the PK gene and protein. Also, the effects of hypoxia on gene expression, enzymatic activity, glucose, and lactate in hepatopancreas and muscle were analyzed. The PK gene is 15,103 bp and contains 11 exons and 10 introns, producing four mRNA variants by alternative splicing and named PK1, PK2, PK3 and PK4, and two proteins with longer C-terminus and two with a 12 bp insertion. The promoter contains putative binding sites for transcription factors (TF) that are typically involved in stress responses. The deduced amino acid sequences contain the classic domains, binding sites for allosteric effectors and potential reversible phosphorylation residues. Protein modeling indicates a homotetramer with highly conserved structure. The effect of hypoxia for 6 and 12 h showed tissue-specific patterns, with higher expression, enzyme activity and lactate in muscle, but higher glucose in hepatopancreas. Changes in response to hypoxia were detected at 12 h in expression with induction in muscle and reduction in hepatopancreas, while enzyme activity was maintained, and glucose and lactate decreased. These results show rapid changes in expression and metabolites, while enzyme activity was maintained to cope with short-term hypoxia.

WSSV exploits AMPK to activate mTORC2 signaling for proliferation by enhancing aerobic glycolysis

AMPK plays significant roles in the modulation of metabolic reprogramming and viral infection. However, the detailed mechanism by which AMPK affects viral infection is unclear. The present study aims to determine how AMPK influences white spot syndrome virus (WSSV) infection in shrimp (Marsupenaeus japonicus). Here, we find that AMPK expression and phosphorylation are significantly upregulated in WSSV-infected shrimp. WSSV replication decreases remarkably after knockdown of Ampkα and the shrimp survival rate of AMPK-inhibitor injection shrimp increases significantly, suggesting that AMPK is beneficial for WSSV proliferation. Mechanistically, WSSV infection increases intracellular Ca2+ level, and activates CaMKK, which result in AMPK phosphorylation and partial nuclear translocation. AMPK directly activates mTORC2-AKT signaling pathway to phosphorylate key enzymes of glycolysis in the cytosol and promotes expression of Hif1α to mediate transcription of key glycolytic enzyme genes, both of which lead to increased glycolysis to provide energy for WSSV proliferation. Our findings reveal a novel mechanism by which WSSV exploits the host CaMKK-AMPK-mTORC2 pathway for its proliferation, and suggest that AMPK might be a target for WSSV control in shrimp aquaculture.

Integrated analysis of intestinal microbiota and metabolomic reveals that decapod iridescent virus 1 (DIV1) infection induces secondary bacterial infection and metabolic reprogramming in Marsupenaeus japonicus

In recent years, with global warming and increasing marine pollution, some novel marine viruses have become widespread in the aquaculture industry, causing huge losses to the aquaculture industry. Decapod iridescent virus 1 (DIV1) is one of the newly discovered marine viruses that has been reported to be detected in a variety of farmed crustacean and wild populations. Several previous studies have found that DIV1 can induce Warburg effect-related gene expression. In this study, the effects of DIV1 infection on intestinal health of shrimp were further explored from the aspects of histological, enzymatic activities, microorganisms and metabolites using Marsupenaeus japonicus as the object of study. The results showed that obvious injury in the intestinal mucosa was observed after DIV1 infection, the oxidative and antioxidant capacity of the shrimp intestine was unbalanced, the activity of lysozyme was decreased, and the activities of digestive enzymes were disordered, and secondary bacterial infection was caused. Furthermore, the increased abundance of harmful bacteria, such as Photobacterium and Vibrio, may synergized with DIV1 to promote the Warburg effect and induce metabolic reprogramming, thereby providing material and energy for DIV1 replication. This study is the first to report the changes of intestinal microbiota and metabolites of M. japonicus under DIV1 infection, demonstrating that DIV1 can induce secondary bacterial infection and metabolic reprogramming. Several bacteria and metabolites highly associated with DIV1 infection were screened, which may be leveraged for diagnosis of pathogenic infections or incorporated as exogenous metabolites to enhance immune response.

White Spot Syndrome Virus Triggers a Glycolytic Pathway in Shrimp Immune Cells (Hemocytes) to Benefit Its Replication

White spot syndrome virus (WSSV) is the causative agent of a shrimp disease that inflicts in huge economic losses in shrimp-farming industry. WSSV triggers aerobic glycolysis in shrimp immune cells (hemocytes), but how this virus regulates glycolytic enzymes or pathway is yet to be characterized. Therefore, mRNA levels and activity of four important glycolytic enzymes, Hexokinase (HK), Phosphofructokinase (PFK), Pyruvate kinase (PK), and Lactate dehydrogenase (LDH), were measured in WSSV-infected shrimp hemocytes. Gene expression of HK and PFK, but not LDH or PK, was increased at the viral genome replication stage (12 hpi); furthermore, activity of these enzymes, except HK, was concurrently increased. However, there was no increased enzyme activity at the viral late stage (24 hpi). In vivo dsRNA silencing and glycolysis disruption by 2-DG further confirmed the role of glycolysis in virus replication. Based on tracing studies using stable isotope labeled glucose, glycolysis was activated at the viral genome replication stage, but not at the viral late stage. This study demonstrated that WSSV enhanced glycolysis by activating glycolytic enzyme at the viral genome replication stage, providing energy and biomolecules for virus replication.

The Possible Role of Glucose-6-Phosphate Dehydrogenase in the SARS-CoV-2 Infection

Glucose-6-phosphate dehydrogenase (G6PD) is the second rate-limiting enzyme of the pentose phosphate pathway. This enzyme is present in the cytoplasm of all mammalian cells, and its activity is essential for an adequate functioning of the antioxidant system and for the response of innate immunity. It is responsible for the production of nicotinamide adenine dinucleotide phosphate (NADPH), the first redox equivalent, in the pentose phosphate pathway. Viral infections such as SARS-CoV-2 may induce the Warburg effect with an increase in anaerobic glycolysis and production of lactate. This condition ensures the success of viral replication and production of the virion. Therefore, the activity of G6PD may be increased in COVID-19 patients raising the level of the NADPH, which is needed for the enzymatic and non-enzymatic antioxidant systems that counteract the oxidative stress caused by the cytokine storm. G6PD deficiency affects approximately 350–400 million people worldwide; therefore, it is one of the most prevalent diseases related to enzymatic deficiency worldwide. In G6PD-deficient patients exposed to SARS-CoV-2, the amount of NADPH is reduced, increasing the susceptibility for viral infection. There is loss of the redox homeostasis in them, resulting in severe pneumonia and fatal outcomes.

The potential role of eyestalk in the immunity of Litopenaeus vannamei to Vibrio parahaemolyticus infection II. From the perspective of long non-coding RNA

Long non-coding RNAs (lncRNAs) have been linked to immunological modulation. Unfortunately, little is known about the processes of immune control in shrimp. In crustaceans such as Litopenaeus vannamei, a prominent aquaculture species, the X-organ-sinus gland complex (XO-SG) in the eyestalk is an essential neuroendocrine regulatory organ. Eyestalk ablation is commonly employed in aquaculture to accelerate ovarian maturation in shrimp. It does, however, have a negative impact on the shrimps' immunocompetence and causes death. As a result, we used RNA-seq to profile the transcriptomes of L. vannamei hemocytes infected with Vibrio parahaemolyticus after the eyestalk ablation. Following strict transcript screening procedures, 2307 lncRNAs were identified from L. vannamei hemocytes in this study. Pearson correlation analysis was finally used to uncover 535 DElncRNAs and 1566 DEmRNA targets. According to the Venn diagram analysis, 326 non-eyestalk regulatory lncRNAs (NElncRNAs) with a target of 1014 non-eyestalk regulatory genes (NEmRNAs), 47 eyestalk negative regulatory lncRNAs (ENRlncRNAs) with a target of 95 eyestalk negative regulatory genes (ENRmRNAs), and 162 eyestalk positive regulatory lncRNAs (EPRlncRNAs) with a target of 457 eyestalk positive regulatory genes (EPRmRNAs) were screened. The bioinformatics analysis revealed that lncRNAs were associated with Axon regeneration, Rap1 signaling pathway, Thyroid hormone signaling pathway, TGF-beta signaling pathway, and PI3K-Akt signaling pathway, implying that lncRNAs may play a role in the regulation of the neuroendocrine-immune (NEI) system. Furthermore, several lncRNAs targeting HSP70, YWHAZ, FER2, HIF1α, and Notch were discovered and verified by qRT-PCR. These findings showed that regulation of lncRNAs in hemocytes which were controlled by the eyestalk might be one of the impact variables in controlling the differential expression of mRNAs associated with immune response in L. vannamei infected with V. parahaemolyticus.

Characterization of phosphofructokinase (PFK) from mud crab Scylla paramamosain and its role in mud crab dicistrovirus-1 proliferation

Phosphofructokinase (PFK), the key enzyme of glycolysis, can catalyze the irreversible transphosphorylation of fructose-6-phosphate forming fructose-1, 6-biphosphate. In the present study, a PFK gene from the mud crab Scylla paramamosain, named SpPFK, was cloned and characterized. The full length of SpPFK contained a 5'untranslated region (UTR) of 249 bp, an open reading frame of 2,859 bp, and a 3'UTR of 1,248 bp. The mRNA of SpPFK was highly expressed in the gill, followed by the hemocytes and muscle. The expression of SpPFK was significantly up-regulated after mud crab dicistrovirus-1 (MCDV-1) infection. Knocking down SpPFK in vivo by RNA interference significantly reduced the expression of lactate dehydrogenase after MCDV-1 infection. Furthermore, silencing of SpPFK in vivo increased the survival rate of mud crabs and decreased the MCDV-1 copies in the gill and hepatopancreas after MCDV-1 infection. All these results suggested that SpPFK could play an important role in the process of MCDV-1 proliferation in mud crab.

Correlation between acute brain injury and brain metabonomics in dichlorvos-poisoned broilers

Dichlorvos (DDVP) is an insecticide with neurotoxicity that is widely used in agricultural production and life. However, the effects of acute DDVP poisoning on brain tissue remain underinvestigated. The purpose of this study was to evaluate the differences within 15 min-6 h in plasma biochemical indexes, brain histology and metabolites among three groups of commercial broilers orally administered different dosages of DDVP one time: (1) high-dose group (11.3 mg/kg), (2) low-dose group (2.48 mg/kg) and (3) control group (0 mg/kg). The results of biochemical indexes showed that acute DDVP poisoning could cause hyperglycemia and oxidative stress in poisoned broilers. Histological examination showed that DDVP could induce brain edema, abnormal expression of glial fibrillary acidic protein (GFAP) and neuronal mitochondrial damage in broilers. Whole-brain metabolism showed that DDVP could significantly change the secretion of neurotransmitters, energy metabolism, amino acid metabolism and nucleotide metabolism. Correlation analysis showed that metabolites such as hypoxanthine, acetylcarnitine and glucose 6-phosphate were significantly correlated with blood glucose, biomarkers of oxidative stress and brain injury pathology. The results of this study provide new insights into the molecular mechanism of brain tissue responses to acute DDVP exposure in broilers and deliver important information for clinical research on neurodegenerative diseases caused by acute DDVP poisoning.

HIF-1 is involved in the regulation of expression of metallothionein and apoptosis incidence in different oxygen conditions in the white shrimp Litopenaeus vannamei

The white shrimp Litopenaeus vannamei is exposed to hypoxic conditions in natural habitats and in shrimp farms. Hypoxia can retard growth, development and affect survival in shrimp. The hypoxia-inducible factor 1 (HIF-1) regulates many genes involved in glucose metabolism, antioxidant proteins, including metallothionein (MT) and apoptosis. In previous studies we found that the L. vannamei MT gene expression changed during hypoxia, and MT silencing altered cell apoptosis; in this study we investigated whether the silencing of HIF-1 affected MT expression and apoptosis. Double-stranded RNA (dsRNA) was used to silence HIF-1α and HIF-1β under normoxia, hypoxia, and hypoxia plus reoxygenation. Expression of HIF-1α, HIF-1β and MT, and apoptosis in hemocytes or caspase-3 expression in gills, were measured at 0, 3, 24 and 48 h of hypoxia and hypoxia followed by 1 h of reoxygenation. The results showed that hemocytes HIF-1α expression was induced during hypoxia and reoxygenation at 3 h, while HIF-1β decreased at 24 and 48 h. In normoxia, HIF-1 silencing in hemocytes increased apoptosis at 3 h and decreased at 48 h; while in gills, caspase-3 increased at 3, 24 and 48 h. In hypoxia, HIF-1 silencing decreased apoptosis in hemocytes at 3 h, but caspase-3 increased in gills. During reoxygenation, apoptosis in hemocytes and caspase-3 in gills increased. During normoxia in hemocytes, silencing of HIF-1 decreased MT expression, but in gills, MT increased. During hypoxia and reoxygenation, silencing induced MT in hemocytes and gills. These results indicate HIF-1 differential participation in MT expression regulation and apoptosis during different oxygen conditions.

Integrative Transcriptomic Analysis Reveals the Immune Mechanism for a CyHV-3-Resistant Common Carp Strain

Anti-disease breeding is becoming the most promising solution to cyprinid herpesvirus-3 (CyHV-3) infection, the major threat to common carp aquaculture. Virus challenging studies suggested that a breeding strain of common carp developed resistance to CyHV-3 infection. This study illustrates the immune mechanisms involved in both sensitivity and anti-virus ability for CyHV3 infection in fish. An integrative analysis of the protein-coding genes and long non-coding RNAs (lncRNAs) using transcriptomic data was performed. Tissues from the head kidney of common carp were extracted at days 0 (the healthy control) and 7 after CyHV-3 infection (the survivors) and used to analyze the transcriptome through both Illumina and PacBio sequencing. Following analysis of the GO terms and KEGG pathways involved, the immune-related terms and pathways were merged. To dig out details on the immune aspect, the DEGs were filtered using the current common carp immune gene library. Immune gene categories and their corresponding genes in different comparison groups were revealed. Also, the immunological Gene Ontology terms for lncRNA modulation were retained. The weighted gene co-expression network analysis was used to reveal the regulation of immune genes by lncRNA. The results demonstrated that the breeding carp strain develops a marked resistance to CyHV-3 infection through a specific innate immune mechanism. The featured biological processes were autophagy, phagocytosis, cytotoxicity, and virus blockage by lectins and MUC3. Moreover, the immune-suppressive signals, such as suppression of IL21R on STAT3, PI3K mediated inhibition of inflammation by dopamine upon infection, as well as the inhibition of NLRC3 on STING during a steady state. Possible susceptible factors for CyHV-3, such as ITGB1, TLR18, and CCL4, were also revealed from the non-breeding strain. The results of this study also suggested that Nramp and PAI regulated by LncRNA could facilitate virus infection and proliferation for infected cells respectively, while T cell leukemia homeobox 3 (TLX3), as well as galectin 3 function by lncRNA, may play a role in the resistance mechanism. Therefore, immune factors that are immunogenetically insensitive or susceptible to CyHV-3 infection have been revealed.

How do abiotic environmental conditions influence shrimp susceptibility to disease? A critical analysis focussed on White Spot Disease

White Spot Syndrome Virus (WSSV) causes White Spot Disease (WSD) and is historically the most devastating disease in the shrimp industry. Global losses from this disease have previously exceeded $3 bn annually, having a major impact on a global industry worth US$19 bn per annum. Shrimp are cultured predominantly in enclosed ponds that are subject to considerable fluctuations in abiotic conditions and WSD outbreaks are increasingly linked to periods of extreme weather, which may cause major fluctuations in pond culture conditions. Combined with the intensity of production in these systems, the resulting suboptimal physicochemical conditions have a major bearing on the susceptibility of shrimp to infection and disease. Current knowledge indicates that pond temperature and salinity are major factors determining outbreak severity. WSSV appears to be most virulent in water temperatures between 25 and 28 °C and salinities far removed from the isoosmotic point of shrimp. Elevated temperatures (>30 °C) may protect against WSD, depending on the stage of infection, however the mechanisms mediating this effect have not been well established. Other factors relating to water quality that may play key roles in determining outbreak severity include dissolved oxygen concentration, nitrogenous compound concentration, partial pressure of carbon dioxide and pH, but data on their impacts on WSSV susceptibility in cultured shrimps is scarce. This illustrates a major research gap in our understanding of the influence of environmental conditions on disease. For example, it is not clear whether temperature manipulations can be used effectively to prevent or mitigate WSD in cultured shrimp. Therefore, developing our understanding of the impact of environmental conditions on shrimp susceptibility to WSSV may provide insight for WSD mitigation when, even after decades of research, there is no effective practical prophylaxis or treatment.

Evolutionary adaptation analysis of immune defense and hypoxia tolerance in two closely related Marsupenaeus species based on comparative transcriptomics

Kuruma shrimp, a major farmed shrimp species in the world, includes two cryptic or sibling species, Form I (Marsupenaeus japonicus) and Form II (Marsupenaeus pulchricaudatus). Due to the lack of genomic resources, little is known about the molecular mechanisms associated with immune defense and hypoxia tolerance. Here, we sequenced the transcriptomes of two closely related Marsupenaeus species and compared genomic divergence. This study obtained 77049 and 84561 unigenes with N50 values of 1281bp and 1244bp for M. japonicus and M. pulchricaudatus, respectively, and 5036 pairs of putative orthologs were identified between two Marsupenaeus species. Estimation of Ka/Ks ratios indicated that 165 orthologous genes may be under positive selection (Ka/Ks > 0.5), including 49 pairs with a Ka/Ks ratio >1. According to the peak of synonymous rates, the divergence time between M. japonicus and M. pulchricaudatus was about 0.26-0.69 Mya. These positively selected orthologous genes related to the immune process mainly comprised single VWC domain protein, legumain, ras-related C3 botulinum, caspase, C-type lectin and were enriched in functions related to immune (Toll-like receptor and PI3K-Akt signaling) and hypoxia signaling (HIF-1 signaling and VEGF signaling). In this study, dozens of caspase-like unigenes were screened from two Marsupenaeus transcriptomes. Among these, the PjCaspase orthologous gene was subjected to positive selection (Ka/Ks = 1.22), which had different secondary and three-dimensional structure prediction. Based on the single copy caspase gene, eight populations of Marsupenaeus species were divided into two phylogeographic lineages from the East and South China. We characterized the transcriptomes of the two Marsupenaeus species and obtained several key orthologs associated with immune defense and hypoxia tolerance, which provides new insights into the immunity and genetic divergence of the two varieties. Moreover, this study will facilitate further comparative genomic studies of the two varieties.