Improving the activity and thermostability of GH2 β-glucuronidases via domain reassembly

Glycoside hydrolase family 2 (GH2) enzymes are generally composed of three domains: TIM-barrel domain (TIM), immunoglobulin-like β-sandwich domain (ISD), and sugar-binding domain (SBD). The combination of these three domains yields multiple structural combinations with different properties. Theoretically, the drawbacks of a given GH2 fold may be circumvented by efficiently reassembling the three domains. However, very few successful cases have been reported. In this study, we used six GH2 β-glucuronidases (GUSs) from bacteria, fungi, or humans as model enzymes and constructed a series of mutants by reassembling the domains from different GUSs. The mutants PGUS-At, GUS-PAA, and GUS-PAP, with reassembled domains from fungal GUSs, showed improved expression levels, activity, and thermostability, respectively. Specifically, compared to the parental enzyme, the mutant PGUS-At displayed 3.8 times higher expression, the mutant GUS-PAA displayed 1.0 time higher catalytic efficiency (k_cat /K_m ), and the mutant GUS-PAP displayed 7.5 times higher thermostability at 65°C. Furthermore, two-hybrid mutants, GUS-AEA and GUS-PEP, were constructed with the ISD from a bacterial GUS and SBD and TIM domain from fungal GUSs. GUS-AEA and GUS-PEP showed 30.4% and 23.0% higher thermostability than GUS-PAP, respectively. Finally, molecular dynamics simulations were conducted to uncover the molecular reasons for the increased thermostability of the mutant.

Switchable dual-band to broadband terahertz metamaterial absorber incorporating a VO2 phase transition

We design and demonstrate a thermally switchable terahertz metamaterial absorber consisting of an array of orthogonal coupled split-ring metal resonators involving a VO₂ phase transition. Numerical results indicate that the active metamaterial always absorbs the TE wave in dual-band regardless of insulating and metallic VO₂, while the insulator-to-metal phase transition enables a switchable effect between dual-band and broadband absorption of the TM wave with the resonant frequency tunability of 33%. Especially under the metallic VO₂ state, the absorption properties are polarization-dependent and exhibit a switching effect between dual-band and broadband absorption with the increase of the polarization angle. The tunable absorption mechanism can be explained by effective impedance theory and electric energy density distributions. The proposed dual-band to broadband metamaterial switching absorber may have broad applications in sensors, imaging and emitters.

Autopsy and statistical evidence of disturbed hemostasis progress in COVID-19: medical records from 407 patients

BACKGROUND

The progression of coagulation in COVID-19 patients with confirmed discharge status and the combination of autopsy with complete hemostasis parameters have not been well studied.

OBJECTIVE

To clarify the thrombotic phenomena and hemostasis state in COVID-19 patients based on epidemiological statistics combining autopsy and statistical analysis.

METHODS

Using autopsy results from 9 patients with COVID-19 pneumonia and the medical records of 407 patients, including 39 deceased patients whose discharge status was certain, time-sequential changes in 11 relevant indices within mild, severe and critical infection throughout hospitalization according to the Chinese National Health Commission (NHC) guidelines were evaluated. Statistical tools were applied to calculate the importance of 11 indices and the correlation between those indices and the severity of COVID-19.

RESULTS

At the beginning of hospitalization, platelet (PLT) counts were significantly reduced in critically ill patients compared with severely or mildly ill patients. Blood glucose (GLU), prothrombin time (PT), activated partial thromboplastin time (APTT), and D-dimer levels in critical patients were increased compared with mild and severe patients during the entire admission period. The International Society on Thrombosis and Haemostasis (ISTH) disseminated intravascular coagulation (DIC) score was also high in critical patients. In the relatively late stage of nonsurvivors, the temporal changes in PLT count, PT, and D-dimer levels were significantly different from those in survivors. A random forest model indicated that the most important feature was PT followed by D-dimer, indicating their positive associations with disease severity. Autopsy of deceased patients fulfilling diagnostic criteria for DIC revealed microthromboses in multiple organs.

CONCLUSIONS

Combining autopsy data, time-sequential changes and statistical methods to explore hemostasis-relevant indices among the different severities of the disease helps guide therapy and detect prognosis in COVID-19 infection.

Transformation pathways of chlorinated paraffins relevant for remediation: a mini-review

In the past decades, the environmental presence and ecological risks of chlorinated paraffins (CPs), an emerging class of organic halogen compounds, have been receiving increasing attention worldwide. Short-chain CPs (SCCPs) and medium-chain CPs (MCCPs) constitute the important CPs of considerable concern. In this review article, the state-of-the-art research status on the environmental transformation of CPs, including thermal decomposition, photolytic and photocatalytic degradation, biological metabolism, and atmospheric transformation, was summarized and integrated in detail. The degradation efficiency and transformation products of CPs in these environmental processes were evaluated, in which dechlorination was considered as the major reaction pathway. Notably, waste incineration of CPs has been demonstrated to generate a variety of persistent chlorinated aromatic hydrocarbons such as polychlorinated biphenyls and polychlorinated naphthalenes, which have more significant environmental impacts. Additionally, photodegradation and photocatalysis are suggested as the feasible techniques for efficient removal of SCCPs from water matrices. Overall, the current transformation studies of CPs could facilitate the comprehensive understanding of their environmental behaviors and fate as well as the development of promising remediation strategies for pollution control.

Cellular and molecular landscape of mammalian sinoatrial node revealed by single-cell RNA sequencing

Bioelectrical impulses intrinsically generated within the sinoatrial node (SAN) trigger the contraction of the heart in mammals. Though discovered over a century ago, the molecular and cellular features of the SAN that underpin its critical function in the heart are uncharted territory. Here, we identify four distinct transcriptional clusters by single-cell RNA sequencing in the mouse SAN. Functional analysis of differentially expressed genes identifies a core cell cluster enriched in the electrogenic genes. The similar cellular features are also observed in the SAN from both rabbit and cynomolgus monkey. Notably, Vsnl1, a core cell cluster marker in mouse, is abundantly expressed in SAN, but is barely detectable in atrium or ventricle, suggesting that Vsnl1 is a potential SAN marker. Importantly, deficiency of Vsnl1 not only reduces the beating rate of human induced pluripotent stem cell - derived cardiomyocytes (hiPSC-CMs) but also the heart rate of mice. Furthermore, weighted gene co-expression network analysis (WGCNA) unveiled the core gene regulation network governing the function of the SAN in mice. Overall, these findings reveal the whole transcriptome profiling of the SAN at single-cell resolution, representing an advance toward understanding of both the biology and the pathology of SAN.

The spontaneous bioelectrical activity of pacemaker cells in sinoatrial node (SAN) triggers the heartbeats. Here, the authors perform single-cell RNA sequencing in the mouse SAN and identify molecular and cellular features of the SAN conserved in rabbit and cynomolgus monkey, identifying a new potential SAN marker.

Pulmonary thromboembolism after distal ulna and radius fractures surgery: A case report and a literature review

BACKGROUND

Pulmonary thromboembolism (PTE) is a serious postoperative complication that can occur after a fracture. Generally, PTE is caused by the falling off of lower extremity deep vein thrombosis (LEDVT) after lower limb fracture surgery. LEDVT and PTE after upper extremity fracture surgery are very rare. PTE is one of the most common clinical causes of sudden death. Venous thromboembolism includes PTE and DVT. We experienced one case of LEDVT and PTE after distal ulna and radius fracture surgery. The purpose of our report is to raise awareness for orthopedic surgeons that PTE can occur after distal ulna and radius fracture surgery, and patients with high risk factors should be considered for prevention and treatment of thrombosis in a timely manner.

CASE SUMMARY

We report a 51-year-old Chinese male who had severe fractures of the left distal ulna, radius and little finger after a motorcycle accident. The patient underwent external fixation, open reduction and internal fixation. On the third post-operative day, computed tomographic pulmonary angiography showed PTE. Doppler ultrasonography showed thrombus formation in the bilateral posterior tibial veins. After a period of anticoagulation therapy, on the 25^th d after the PTE, computed tomographic pulmonary angiography showed that thrombus in both sides of the pulmonary artery disappeared. Furthermore, about 4 mo after the PTE, thrombosis in the deep veins of the lower limbs disappeared. About 1 year after the surgery, X-rays showed good fracture healing, and the function of the wrist joint recovered well.

CONCLUSION

Though rare, PTE can occur after distal ulna and radius fracture surgery and patients with high risk factors should be considered for prevention and treatment of thrombosis in a timely manner.

EMG Signal Filtering Based on Variational Mode Decomposition and Sub-Band Thresholding

Surface electromyography (EMG) signals are inevitably contaminated by various noise components, including powerline interference (PLI), baseline wandering (BW), and white Gaussian noise (WGN). These noises directly degrade the efficiency of EMG processing and affect the accuracy and robustness of further applications. Currently, most of the EMG filters only target one category of noise. Here, we propose a novel filter to remove all three types of noise. The noisy EMG signal is first decomposed into an ensemble of band-limited modes using variational mode decomposition (VMD). Each category of noise is located within specific modes and is separately removed in sub-bands. In particular, WGN is suppressed by soft thresholding with a noise level-dependent threshold. The denoising performance was assessed from simulated and experimental signals using three performance metrics: the root mean square error ([Formula: see text]), the improvement in signal-to-noise ratio ([Formula: see text]), and the percentage reduction in the correlation coefficient ( η). Other methods, including traditional infinite impulse response (IIR) filters, empirical mode decomposition (EMD) method, and ensemble empirical mode decomposition (EEMD) method, were examined for comparison. The proposed method achieved the best performance to remove BW or WGN. It also effectively reduced PLI noise when the signal-to-noise ratio (SNR) was low. The SNR was improved by 18.6, 19.2, and 8.0 dB for EMG signals corrupted with PLI, BW, and WGN at -6 dB SNR, respectively. The experimental results illustrated that noise was completely removed from resting states, and obvious spikes were distinguished from action states. For two of the ten subjects, the improved SNR reached 20 dB. This study explores the special characteristics of VMD and demonstrates the feasibility of using the VMD-based filter to denoise EMG signals. The proposed filter is efficient at removing three categories of noise and can be used for any application that requires EMG signal filtering at the preprocessing stage, such as gesture recognition and EMG decomposition.

Cadmium exposure alters expression of protective enzymes and protein processing genes in venom glands of the wolf spider Pardosa pseudoannulata

Cadmium (Cd) pollution is currently the most serious type of heavy metal pollution throughout the world. Previous studies have shown that Cd elevates the mortality of paddy field spiders, but the lethal mechanism remains to be explored profoundly. In the present study, we measured the activities of protective enzymes (acetylcholinesterase, glutathione peroxidase, phenol oxidase) and a heavy metal chelating protein (metallothionein) in the pond wolf spider Pardosa pseudoannulata after Cd exposure. The results indicated that Cd initially increased the enzyme activities and protein concentration of the spider after 10- and 20-day exposure before inhibiting them at 30-day exposure. Further analysis showed that the enzyme activities in the cephalothorax were inhibited to some extent. Since the cephalothorax region contains important venom glands, we performed transcriptome sequencing (RNA-seq) analysis of the venom glands collected from the spiders after long-term Cd exposure. RNA-seq yielded a total of 2826 differentially expressed genes (DEGs), and most of the DEGs were annotated into the process of protein synthesis, processing and degradation. Furthermore, a mass of genes involved in protein recognition and endoplasmic reticulum (ER) -associated protein degradation were down-regulated. The reduction of protease activities supports the view that protein synthesis and degradation in organelles and cytoplasm were dramatically inhibited. Collectively, our outcomes illustrate that Cd poses adverse effects on the expression of protective enzymes and protein, which potentially down-regulates the immune function in the venom glands of the spiders via the alteration of protein processing and degradation in the ER.

Facile construction of a molecularly imprinted polymer-based electrochemical sensor for the detection of milk amyloid A

A molecularly imprinted electrochemical sensor for the detection of serum amyloid A (MAA) in milk was established for early diagnosis of subclinical mastitis in dairy cows. The electrochemical sensor was initially constructed using a nanocomposite material (reduced graphene oxide/gold nanoparticles, AuNPs@rGO) to modify the working electrode. The template protein, MAA, was then immobilized using pyrrole as the functional monomer to carry out the electropolymerization. Finally, the template protein was removed to form a molecular imprint film with the capability to qualitatively and quantitatively signaling of MAA. Cyclic voltammetry (CV), differential pulse voltammetry (DPV), and scanning electron microscopy (SEM) were used to characterize the modification process of the molecularly imprinted electrochemical sensors. Under optimized conditions, the sensor shows two well-behaved linear relationships in the MAA concentration range 0.01 to 200 ng/mL. A lower detection limit was estimated to be 5 pg/mL (S/N = 3). Other parameters including the selectivity, reproducibility (RSD 3.2%), and recovery rate (96.1-103%) are all satisfactory. Compared with the traditional methods, detection of MAA to determine the subclinical mastitis of dairy cows can efficiently be diagnosed and hence prevent an outbreak of dairy cow mastitis. The electrochemical sensor can detect MAA more rapidly, sensitively, and inexpensively than the ELISA-based MAA detection. These advantages indicate that the method is promising for early diagnosis of dairy cows.

[Corrigendum] Nicorandil protects mesenchymal stem cells against hypoxia and serum deprivation-induced apoptosis

After the publication of the above article, the authors have realized that Figs. 2 and 4 in their paper were published with incorrect images; regarding Fig. 2, the data featured in Fig. 2A (for the H/SD + Nico 1000 µM panel) were repeated with those featured in Fig. 1C (the 6 h H/SD panel), and the data shown for Bcl-2 in Fig. 4C were selected incorrectly. These errors arose inadvertently as a consequence of misassembling the figures. The revised versions of Figs. 2 and 4, featuring the corrected data panels for the above‑mentioned experiments, are shown on the next page. Note that the revised data shown for these Figures do not affect the overall conclusions reported in the paper. The authors express their gratitude to the Editor of International Journal of Molecular Medicine for allowing them the opportunity to publish this corrigendum, and apologize to the readership for any inconvenience caused. [the original article was published in International Journal of Molecular Medicine 36: 415-423, 2015; DOI: 10.3892/ijmm.2015.2229].

Pulmonary oligometastases treated by stereotactic body radiation therapy (SBRT): a single institution's experience

Background

To investigate the effect of stereotactic body radiation therapy (SBRT) on pulmonary oligometastases and to analyze the clinical factors and dose parameters affecting local recurrence-free survival (LRFS) and overall survival (OS).

Methods

This study retrospectively enrolled a total of 84 patients (148 lesions) treated in our department from May 2015 to November 2018. Pulmonary oligometastases was defined as up to 5 metastatic lesions in the lung and with both the primary tumor and any extra-thoracic metastases being controlled. Patients receiving a BED₁₀ (biological effective dose, α/β =10) of SBRT ≥75 Gy and a dose/fraction ≥4 Gy were enrolled. The patient group consisted of 52 men (61.9%) and 32 women (38.1%), with a median age 56 years (range, 29-80 years). Median tumor diameter was 1.71cm (range, 1.2-5.0 cm). The BED₁₀ was 75-119 Gy in 4-15 fractions. Univariate and multivariate Cox regression analyses were performed on factors predicting the outcomes.

Results

All patients completed the treatment as planned, and the median follow-up time was 20.3 months. The median OS for the entire group was 34.3 months, with an actuarial 1-, 2-, 3- and 5-year OS of 74.7%, 59.4%, 49.7%, and 36.8%, respectively. Among the 148 lesions in the whole group, 19 (12.8%) lesions had local recurrence (LR). The median LRFS time for all patients was 56.9 months. The LRFS rate was 93.6%, 83.5%, 81.4%, and 76.6% at 1, 2, 3, and 5 years, respectively. No patient developed acute grade 3 or 4 toxicity. On univariate analysis, age ≥63 years old, primary site of colorectal cancer, BED₁₀ <85.2 Gy, pathological type of adenocarcinoma, planning target volume (PTV) min BED₁₀ <76.6 Gy, and gross tumor volume (GTV) ≥8.8 cc, were significantly associated with poorer LRFS. Multivariate analysis showed that age ≥63 years old, primary site of colorectal cancer, and PTV min BED₁₀ <76.6 Gy were significant risk factors affecting LRFS.

Conclusions

SBRT is feasible for pulmonary oligometastasis with favorable local control and minimal toxicity. Multiple dose parameters, instead of a prescription dose only, in combination with clinical parameters, should be considered for optimal local control.

Design of easy-manufacturing superdirective antenna: a theoretical study

An easy-fabricating superdirective antenna with needlelike radiation in the microwave band is proposed for the first time to the best of our knowledge. The electromagnetic field of the superdirective antenna excited by an electric line source is calculated. The parameters of the superdirective antenna are obtained by optimization algorithms. A two-step optimization method is carried out to reduce the difficulty of fabrication. The mechanism of superdirectivity in the microwave band is explained briefly by traditional multilayered dielectric resonances. A machinable 10-layered superdirective antenna, using this method, is proposed. Unlike traditional superdirective antennas, the proposed antenna, after two-step optimization in this paper, has all-positive integral permittivities, which reduce the fabrication difficulty significantly.

Transcriptome sequencing reveals the effects of cadmium toxicity on the cold tolerance of the wolf spider Pirata subpiraticus

As the predominant predator of pests in rice fields, spiders have been exposed to cadmium (Cd) pollution for a long time. The livability of spiders during the overwintering period is closely related to population growth in spring, but the effects of Cd on spider's survival of cold hardness and the underlining mechanism remain unclear. In the present study, we found that some growth parameters (body length, width, mass and livability) in the wolf spider Pirata subpiraticus were altered distinctively under Cd stress. To investigate the effects of Cd toxicity on the spider at molecular levels, RNA-sequencing was performed on the spiderlings undergoing ambient temperature alterations. Transcriptome data showed that a total of 807 differentially expressed genes (DEGs) were yielded in the comparison. The obtained DEGs were mainly linked with metabolism-related process, including oxidoreductase activity and lipid transport, and 25 DEGs were associated with the reported cryoprotectants, including glycerol, arginine, cysteine, heat shock protein, glucose and mannose. Growth factors (insulin growth factor, platelet-derived growth factor and transforming growth factor) and cytochrome P450 encoding genes were dramatically expressed in the spider. Furthermore, transcriptional factors (TFs) family were characterized according to the transcriptomic profile, and ZBTB TFs were represented the most distinctive alterations in the characterized genes. Collectively, our study illustrated that Cd poses disadvantageous effects on the growth of P. subpiraticus at cold ambient temperature, and the spiders are capable of responding to the adverse Cd stress by expressing the genes involved in the metabolism of energy substances, cryoprotectants and immune-related components.

Galactosylation of Monosaccharide Derivatives of Glycyrrhetinic Acid by UDP-Glycosyltransferase GmSGT2 from Glycine max

Glycyrrhetinic acid (GA), a pentacyclic triterpenoid aglycone, is the major functional component in licorice which mainly exists in the form of functional glycosides in licorice. The introduction of a sugar moiety to the C-3 OH of GA to yield glycosylated derivatives has been reported, but the late-stage glycosylation of GA-3-O-sugar to form rare GA glycosides with more complexed glycosyl decoration has been rarely reported. In this study, a unique UDP-galactosyltransferase GmSGT2 from Glycine max was found to transfer a galactose to the C2 position of the sugar moiety of GA-3-O-monoglucuronide (GAMG) and GA-3-O-monoglucose. In addition to UDP-galactose, GmSGT2 also recognizes UDP-glucose, UDP-xylose, and UDP-arabinose with relative activities of 32.1-89.2%. Based on a test of 12 typical natural products, GmSGT2 showed high specificity toward the pentacyclic triterpenoid skeleton as the sugar acceptor. Molecular docking was performed to elucidate the substrate recognition mechanism of GmSGT2 toward GAMG.

Correction to: Roles of microRNA-34a targeting SIRT1 in mesenchymal stem cells

An amendment to this paper has been published and can be accessed via the original article.

Chiral Organophosphorous Pesticide Fosthiazate: Absolute Configuration, Stereoselective Bioactivity, Toxicity, and Degradation in Vegetables

Fosthiazate is a widely used chiral organophosphorous nematicide with four stereoisomers. The present study systemically assessed the stereoselectivity of fosthiazate for the first time, including absolute configuration confirmation, stereoselective bioactivity toward nematode and aphid, toxicity to honeybees, and stereoselective degradation in cucumber and pepper under field conditions. The absolute configurations of the four stereoisomers that eluted on the Superchiral IG-3 column were confirmed as (1S,3R)-(-)-fosthiazate, (1S,3S)-(-)-fosthiazate, (1R,3S)-(+)-fosthiazate, and (1R,3R)-(+)-fosthiazate. In comparison to the other two stereoisomers, (1S,3R)-fosthiazate and (1S,3S)-fosthiazate possess more than 100 times bioactivity and 10 times toxicity toward the target and non-target organisms, respectively. The molecular docking found that (1S,3R)-fosthiazate and (1S,3S)-fosthiazate had shorter binding distances and lower energies with acetylcholinesterase (AChE), which illuminated the mechanism of the experimental results. In addition, both of the high-bioactive stereoisomers had faster degradation rates in cucumber and pepper. On the basis of the results of bioactivity, toxicity, and degradation behavior, the stereoisomer mixture with (1S,3R)-fosthiazate and (1S,3S)-fosthiazate will be a better option than racemic fosthiazate to increase the bioactivity and reduce application rates.

Regulating Strategies for Producing Carbohydrate Active Enzymes by Filamentous Fungal Cell Factories

Filamentous fungi are important eukaryotic organisms crucial in substrate degradation and carbon cycle on the earth and have been harnessed as cell factories for the production of proteins and other high value-added products in recent decades. As cell factories, filamentous fungi play a crucial role in industrial protein production as both native hosts and heterologous hosts. In this review, the regulation strategies of carbohydrate active enzyme expression at both transcription level and protein level are introduced, and the transcription regulations are highlighted with induction mechanism, signaling pathway, and promoter and transcription factor regulation. Afterward, the regulation strategies in protein level including suitable posttranslational modification, protein secretion enhancement, and protease reduction are also presented. Finally, the challenges and perspectives in this field are discussed. In this way, a comprehensive knowledge regarding carbohydrate active enzyme production regulation at both transcriptional and protein levels is provided with the particular goal of aiding in the practical application of filamentous fungi for industrial protein production.

Network Pharmacology and Bioactive Equivalence Assessment Integrated Strategy Driven Q-markers Discovery for Da-Cheng-Qi Decoction to Attenuate Intestinal Obstruction

BACKGROUND

Intestinal obstruction (IO) is a kind of acute abdomen with high morbidity and mortality. Patients suffer from poor quality of life and tremendous financial pressure. Da-Cheng-Qi decoction (DCQD), a classical purgation prescription, has clinically been proven to be an effective treatment for IO.

PURPOSE

Network pharmacology integrated with bioactive equivalence assessment was used to discover the quality marker (Q-marker) of DCQD against IO.

METHODS

As there is hardly any targets recorded in database, thus the collection of IO targets was conducted by searching those of alternative diseases which have similar pathological symptoms with IO. In order to improve the reliability of the obtained targets, IO metabolomics data was introduced. Active compounds combination (ACC) was focused as potential Q-markers via component-target network analysis and function query from the identified components corresponding to the common targets. Bioequivalence between ACC and DCQD was assessed from the aspects of intestine motility (somatostatin secretion), inflammation (IL-6 secretion) and injury (wound healing assay) in vitro and was further validated in ileus rat model. PPI network analysis of core targets followed by gene pedigree classification and experimental validation confirmed the potential intervention pathway.

RESULTS

A combination of 11 ingredients, including emodin, physcion, aloe-emodin, rhein, chrysophanol, gallic acid, magnolol, honokiol, naringenin, tangeretin, and nobiletin was finally confirmed bioequivalence with DQCD to some extent and could serve as Q-markers for DCQD to attenuate IO. PI3K/AKT was verified as a possible affected pathway that DCQD exerted the effectiveness against IO.

CONCLUSION

For the disease with few recorded targets, searching those of alternative diseases which have similar pathological symptoms could be a feasible and effective approach. The proposed network pharmacology integrated bioactive equivalence evaluation paradigm is efficient to discover Q-marker of herbal formulae.

Mul-tiomics analysis of cadmium stress on the ovarian function of the wolf spider Pardosa pseudoannulata

Cadmium (Cd) pollution is widespread in paddy filed soil in China. In this study, the toxicity of Cd with regard to the female reproductive system of paddy spider Pardosa pseudoannulata was investigated by means of multi-omics analyses (transcriptome, proteome, and miRNAs). Decreased activities of detoxifying enzymes including peroxidase (POD), Glutathione S-transferases (GST), and superoxide dismutase were detected in the ovary of P. pseudoannulata. Of these, GST and POD were consistently down-regulated at the transcriptional and translational levels. Vitellogenin content and fecundity of the spider were also reduced by Cd burden. Five vitellogenin encodes genes were down-regulated while only vitellogenin-6 protein was up-regulated. But protein lipovitellin-1, the main composition of vitellin, was down-regulated. In addition, the correlation between the mitogen-activated protein kinase (MAPK) signaling pathway and Cd stress was identified. A down-regulated gene that encoding connector of kinase to AP-1 in the MAPK signaling pathway was regulated by the up-regulated miRNA (miRNA id: miRNA dan-miR- 318>der-miR-318>dgr-miR-318>dme-miR-318-3p > dmo-miR-318>dpe-miR-318>dps-miR-318>dse-miR-318>dsi-miR-318>dvi-miR-318>dwi-miR-318>dya-miR-318). In conclusion, Cd stress possesses distinct female reproductive toxicity on P. pseudoannulata through impairing antioxidant system and synthesis of vitellin.