Objective: To report the antibody identification, blood management during pregnancy and the monitoring process of fetal hemolytic disease of fetus and newborn (HDFN) in a pregnant woman with a history of blood transfusion and pregnancy who developed anti-Jr . Methods: Saline tube technique and anti-human globulin technique were used for maternal blood typing, unexpected antibody screening and identification, as well as for determining antibody titer and IgG subclasses. PCR-SSP was employed for genotyping of 18 blood group systems. Next-generation sequencing (NGS) was utilized for gene sequencing of 38 blood group systems. Sanger sequencing was applied to verify rare blood group mutations detected by NGS and to investigate the corresponding rare blood group genes in family members. Blood preparation was achieved through anemia management in prenatal clinics and autologous blood collection during pregnancy. The newborn underwent the three primary tests for HDFN and plasma IgG subclass testing. Results: The pregnant woman's blood type was B, RhD positive, with a positive unexpected antibody screen, and the antibody identification pattern was consistent with a high-frequency antigen antibody. Gene sequencing revealed a homozygous ABCG2 c.376C>T mutation in the woman, resulting in the Jr(a-) phenotype, and anti-Jr antibody was present in her plasma. No compatible Jr(a-) blood was found among family members. The maternal anti-Jr IgG titer remained stable at 256 during pregnancy, with no detectable IgG1 or IgG3 subclasses against the Jr antigen. A total of 800 mL of autologous blood was collected in two stages during pregnancy. The newborn was B, RhD positive, Jr(a+), with a positive unexpected antibody screen (anti-Jr ). IgG subclass typing detected no IgG1 or IgG3. The direct antiglobulin test was positive, while the acid elution test was negative. Conclusion: The combination of serology and blood group genetic analysis provides a diagnostic basis for identifying antibodies to high-frequency antigens. Managing perinatal anemia and implementing staged autologous blood storage can secure blood supply for the perioperative period. IgG antibody subclass typing offers a reference for clinical assessment and prevention of HDFN.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTo clarify whether METTL14 mediates the core role of acupuncture at Neiguan (PC6) in promoting myelination and improving behavior in young autistic rats through gene intervention technology. MethodsThe ASD model was established by intraperitoneal injection of valproic acid (VPA) in pregnant rats. Male offspring were intracerebroventricularly injected with adenovirus-packaged METTL14 shRNA (sh-METTL14) or its control (sh-NC) on postnatal day 1, with a model group set as well. Subsequently, the juvenile rats were divided into model group, acupuncture group, acupuncture+sh-NC group, and acupuncture+sh-METTL14 group. The acupuncture group received acupuncture at Neiguan (PC6) from postnatal day 7, once daily for 21 consecutive days. Neurobehavioral changes were evaluated by behavioral tests; METTL14 knockdown efficiency and the expression of METTL14, METTL3, and PTEN were detected by quantitative real-time PCR (qRT-PCR) and Western blot (WB); PTEN m⁶A levels were measured by RNA immunoprecipitation-qPCR (RIP-qPCR); myelin ultrastructure, expression of myelin basic protein (MBP) and neurofascin 155 (NF155), and dendritic spine density were observed using transmission electron microscopy (TEM), enzyme-linked immunosorbent assay (ELISA), immunofluorescence, qRT-PCR, and primary neuron culture. ResultsBehaviorally, knockdown of METTL14 significantly counteracted the beneficial effects of acupuncture in improving self-grooming, open field exploration, three-chamber social interaction, and Morris water maze learning and memory (P<0.05, P<0.01). Compared with the acupuncture+sh-NC group, the acupuncture+sh-METTL14 group showed significantly decreased mRNA and protein expression of hippocampal METTL14 (P<0.01), and the upregulating effects of acupuncture on METTL3 and PTEN expression were reversed (P<0.01). Meanwhile, knockdown of METTL14 significantly inhibited the acupuncture-induced increase in PTEN m⁶A levels (P<0.01). Morphologically, knockdown of METTL14 attenuated the improvement of myelin structure by acupuncture, reversed the downregulation of MBP and upregulation of NF155 induced by acupuncture, and blocked the increase in dendritic spine density (P<0.05, P<0.01). ConclusionMETTL14 is a key molecule mediating the therapeutic effect of acupuncture at Neiguan. Acupuncture at Neiguan upregulates METTL14, thereby enhancing m⁶A methylation modification of PTEN mRNA to stabilize its expression, ultimately promoting myelin development and improving behavioral symptoms in ASD juvenile rats. This preliminarily reveals the modern biological connotation of “opening Xuanfu and dredging myelin”.

The initiation of adaptive immune responses relies on the precise recognition and interpretation of antigenic information. In this process, the specific binding of T cell receptors (TCRs) to peptide-major histocompatibility complex (pMHC) molecules represents one of the key molecular events in the initiation of adaptive immune responses. Accordingly, the structural features of TCR-pMHC complexes provide a fundamental basis for dissecting antigen recognition mechanisms and support rational vaccine design, therapeutic target discovery in TCR-based immunotherapy, and TCR identification and optimization. However, experimental determination of TCR-pMHC structures remains costly, time-consuming, and limited in coverage, making computational approaches essential for rapidly obtaining reliable structural information. Computational methods for predicting the structures of TCR-pMHC complexes have advanced rapidly in recent years, driven by progress in deep learning-based modeling frameworks and the increasing availability of structural and sequence resources. Despite these developments, most existing tools do not adequately distinguish the key structural and biophysical differences between MHC class I (MHC-I) and MHC class II (MHC-II) complexes during model construction. As a consequence, their predictive performance differs substantially between class I and class II complexes. In general, structural predictions for class I complexes outperform those for class II complexes. This discrepancy may be related to several fundamental differences between the two systems, including the architecture of the peptide-binding groove, the distribution of peptide lengths, and the properties of peptide flanking residues (PFRs). Compared with MHC-I molecules, MHC-II molecules usually bind longer antigenic peptides, which typically range from 13 to 25 amino acids in length. PFRs at both termini of these peptides participate in regulating the overall conformation of TCR-pMHC class II complexes and exert a pronounced effect on the geometric and physicochemical characteristics of the TCR-pMHC binding interface. Furthermore, within the TCR recognition interface, the complementarity-determining regions (CDRs) consist of segments that differ markedly in conformational behavior. They commonly include regions that are relatively rigid and structurally stable, together with highly flexible segments exhibiting substantial conformational plasticity. These rigidity-flexibility features constitute an essential structural basis enabling TCRs to recognize diverse peptide-MHC ligands and to accommodate conformational heterogeneity at the interface. However, many current modeling tools, in an effort to enforce global conformational stability or reduce structural noise, tend to over-constrain intrinsically flexible regions. Such oversimplification may lead to inappropriate rigidification of flexible CDR loops, resulting in local structural distortions, compromised interface geometry, or even complete modeling failure for specific complexes. Against this background, the review approaches the field from the perspective of computational differences between MHC-I and MHC-II complexes. We first systematically organize and summarize available resources related to TCRs and pMHCs, including structural datasets, sequence databases, prediction tools, and benchmarking studies. We then focus on five representative tools capable of predicting both class I and class II complexes—AlphaFold2, AlphaFold3, TCRmodel2, tFold-TCR, and TCR-pHLA_ModellerS. After excluding structures present in the training sets of these tools, we constructed a benchmark dataset comprising 25 class I and 10 class II TCR-pMHC complexes in the bound state and conducted a systematic evaluation using this dataset. We first employ widely used general evaluation metrics, including All-Atom Root Mean Square Deviation (All-Atom RMSD), Backbone RMSD, Template Modeling score (TM-score), and DockQ, to assess the global conformational accuracy and interface modeling quality of class I and class II complexes. For class II complexes, we propose for the first time a peptide flanking residue deviation index, including the PFRs-Deviation Index (PFRs-DI), N-PFR-Deviation Index (N-PFR-DI), and C-PFR-Deviation Index (C-PFR-DI), to quantitatively characterize conformational deviations in PFRs. In addition, we propose the CDR conformational consistency index (CCC) designed to qualitatively evaluate the ability of prediction tools to capture TCR CDR conformational flexibility. These metrics collectively assess a tool’s ability to model both overall conformation and critical functional regions, thereby addressing the limitations of existing evaluation criteria that overemphasize global structure while inadequately capturing modeling quality in key functional areas. This establishes a unified analytical framework for MHC-I and MHC-II complexes to guide data resource selection, modeling strategy formulation, and evaluation system development. The framework further advances computational modeling and provides crucial support for multi-scale analysis of TCR-pMHC recognition mechanisms and their biological functions.

ObjectiveTo analyze the expression level and clinical significance of Y‑box binding protein‑1（YB1） in tissues and cells of oral squamous cell carcinoma （OSCC）. MethodsBioinformatics analysis， immunohistochemistry， Western blot and qRT-PCR were used to detect the expression of YB1 in OSCC tissues， and the prognosis and survival analysis were performed. The effects of YB1 on the proliferation， migration and invasion of OSCC cells were evaluated by CCK-8， scratch assay and Transwell assay， and the expression changes of phosphoinositide 3-kinase/protein kinase B （PI3K/AKT） signaling pathway and epithelial-mesenchymal transition （EMT） related markers after YB1 knockdown were detected by Western blot. ResultsBioinformation analysis showed that the expression level of YB1 in OSCC samples was higher than that in normal samples （P0.001）， and its expression was correlated with gender （P=0.022） and tumor differentiation degree （P0.001）. Kaplan-Meier plotter survival analysis showed that the survival rate of YB1 was low （Log-rank P=0.012）. Immunohistochemistry， Western blot and qRT-PCR experiments verified that the expression level of YB1 in OSCC tissues was higher than that of normal tissues. The results of cell experiments showed that knockdown of YB1 inhibited the malignant biological behaviors such as proliferation， migration and invasion of OSCC cells. Western blot results showed that the knockdown of YB1 could reduce the protein expression levels of p-PI3K and p-AKT， and inhibit EMT. ConclusionYB1 is highly expressed in tissues and cells in OSCC， and reveals the important role of the YB1/PI3K/AKT axis in the process of EMT of OSCC， which is expected to become a new tumor marker for early diagnosis and treatment and prognosis evaluation of OSCC.

Background: and Purpose Vertebrobasilar fusiform aneurysms (VBFAs) carry substantial morbidity and mortality, but optimal management for unruptured VBFAs remains unclear. We compared the safety and efficacy of conservative management (CM), stent-assisted coiling (SAC), and flow diverters (FDs) in patients with unruptured VBFAs, focusing on long-term prognosis. Methods: This study included data from a nationwide Chinese cohort of patients with vertebrobasilar dissecting aneurysms. Inverse probability of treatment weighting (IPTW) balanced confounders across groups. The primary outcome was poor prognosis (modified Rankin Scale score >2). Secondary outcomes included aneurysm rupture, ischemic stroke, compression symptoms, and VBFA-related deaths. Logistic regression estimated odds ratios (ORs) and 95% confidence intervals (CIs). Subgroup and sensitivity analyses were performed. Results: Among 1,115 patients with unruptured VBFAs, 838 (median age, 54 years; 655 men) were included. After IPTW, baseline characteristics were balanced. Median follow-up was 54 months. FD was associated with a lower risk of poor prognosis than CM (OR, 0.48 [95% CI, 0.30 to 0.77]; p=0.002), with no difference between CM and SAC. FD also reduced aneurysm rupture (OR, 0.20 [95% CI, 0.07 to 0.60]; p=0.004) and compression symptoms (OR, 0.30 [95% CI, 0.13 to 0.68]; p=0.004) versus CM. Time-to-event analyses further revealed significant differences in vertebral artery lesions and Type I–II VBFAs, whereas no significant differences were observed in basilar or vertebrobasilar junction lesions or in Type III–IV VBFAs. Conclusions Compared with CM, FD was associated with improved long-term outcomes in unruptured VBFAs, particularly in vertebral artery lesions and Type I–II VBFAs, although residual confounding cannot be excluded.

ObjectiveHepatocellular carcinoma (HCC) represents 90% of all primary liver cancers. The main risk factors associated with HCC include viral hepatitis (B and/or C), alcohol abuse, and metabolic dysfunction-associated steatotic liver disease (MASLD), which progressively advance to liver fibrosis, cirrhosis, and ultimately evolve into HCC. Surgical resection represents the most effective treatment for HCC, while recent advances in immunotherapy, including immune checkpoint inhibitors and adoptive cell therapies, have provided improved treatment prospects for patients with unresectable HCC. However, the complex metabolic heterogeneity of HCC limits the therapeutic efficacy. Metabolic intermediates acyl-CoA not only provide energy and substrates for numerous biochemical reactions but also serve as donors for protein lysine acylation, a major class of post-translational modification (PTM). Therefore, a deeper understanding of the molecular mechanisms underlying protein lysine acylation and hepatocarcinogenesis is urgently needed. MethodsThe levels of protein lysine acylation and silence information regulator 5 (SIRT5) expression levels in clinical HCC samples were analyzed by Western blot. Quantitative malonylome and succinylome of HCC samples were analyzed by antibody-based affinity enrichment coupled with tandem mass spectrometry. The proliferation of HCC cells was analyzed with Cell Counting Kit-8 (CCK-8) assays, the apoptosis was quantified by Annexin V-FITC/propidium iodide (PI) staining coupled with flow cytometry, and the ability of cells to migrate was assayed by Transwell assays. The enzymatic activity of glutathione S-transferase Mu 1 (GSTM1) was quantified. Transgenic mice with hepatic overexpression of SIRT5 were constructed using CRISPR-Cas9, and primary hepatocarcinogenesis was induced by administration of diethylnitrosamine. ResultsWestern blot analysis indicated that the expression level of SIRT5 was elevated in clinical samples from HCC patients, and the levels of lysine malonylation, glutarylation, and succinylation were significantly reduced in HCC tissues. Knockout of SIRT5 in MHCC-97H and MHCC-97L hepatoma cells suppressed cell proliferation, and increased the percentage of apoptotic cells significantly. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the differentially malonylome and succinylome of HCC samples revealed significant enrichment in two major classes of biological processes: core energy metabolism (e.g., glycolysis/gluconeogenesis, tricarboxylic acid metabolic process, fatty acid beta oxidation) and detoxification and oxidative stress response (e.g., response to toxic substance, chemical carcinogenesis, reactive oxygen species (ROS)). SIRT5 removes malonylation from lysine residues in GSTM1 and restores its detoxification activity, which is crucial for the survival of hepatocytes under stressed conditions. More importantly, in vivo experiment indicated that hepatic-specific overexpression of SIRT5 in mice accelerated diethylnitrosamine-induced liver fibrosis and hepatocarcinogenesis, indicating the critical role of SIRT5 in HCC progression. ConclusionThis study highlights the previously unrecognized SIRT5-GSTM1 axis as a key regulator in hepatocarcinogenesis, and suggests a potential target for the treatment of patients with HCC.

Objective To investigate the epidemiological characteristics of malaria and provide a basis for developing improved prevention and control measures. Methods Data were obtained from the Chinese Disease Prevention and Control Information System. Malaria surveillance data for Nanning City from January 1,2023, to December 31,2024, were exported from the Infectious Disease Reporting Information Management Subsystem. The characteristics of the two groups of malaria cases were compared. Results A total of 103 imported malaria cases were reported in Nanning City in 2024, representing a 38.32% decrease compared with the same period of the previous year. No statistically significant difference were observed between cases reported in 2023 and 2024 in terms of average age, gender ratio, proportion of parasite species, and monthly reporting distribution;however, statistically significant differences were found in the proportion of reporting areas and current residence areas(χ2= 13.572 and 10.355, respectively; P = 0.001 and 0.035, respectively). The proportion of cases reported in Shanglin County and the proportion of cases residing in Shanglin County were both lower than those during the same period of the previous year. Conclusions The high aggregation of imported malaria cases in Nanning City has decreased. Medical institutions in areas other than Shanglin County should strengthen their vigilance against malaria.

Objective To gain insights into the biological characteristics of different populations of Culex tritaeniorhynchus within and around the Yellow Sea wetland from the perspective of the occurrence of resistance, we investigated the levels of resistance to deltamethrin and kdr gene mutation in the wetland and its peripheral areas. Methods Specimens were collected from Cx. tritaeniorhynchus populations at two monitoring sites in the Rare Bird National Nature Reserve and Tiaozi Ni Wetland Scenic Area, and also from two populations in Yancheng City and the Liuhe District of Nanjing, and the resistance of these mosquitoes to deltamethrin was determined using the CDC biotest bottle method. For each concentration of deltamethrin assessed, a random subset of exposed specimens was selected for amplification of the kdr gene fragment, followed by Sanger sequencing to identify and analyze resistance-associated mutations. Results The LC50 levels of deltamethrin among mosquitoes from the four populations in Luhe, Yancheng, the Rare Bird National Nature Reserve and the Tiaozi Ni Wetland Scenic Area were 2.048 5, 7.798 2, 3.473 3, and 17.695 5 mg/mL, respectively, with corresponding concentrations of deltamethrin ranging from 0.005 to 5.000,0.050 to 50.000,0.050 to 25.000 and 0.050 to 50.000 mg/mL, respectively. Furthermore, the ranges of the KT50 values were 11.76-107.43, 67.05-216.30,29.77-107.43 and 28.40-329.51 min; the 1-h knockdown rates were 34.58%-99.15%, 9.52%-43.80%, 55.09%-73.01%, and 10.09%-68.07%; and the 24-h mortality rates were 12.15%-67.52%,9.52%-79.56%,13.17%-82.21%, and 11.01%-78.99%, respectively. With respect to kdr gene mutation, we assayed a total of 63,70,59, and 57 mosquitoes for the four populations, for which we detected L1014F mutation frequencies of 14.29%, 35.00%, 20.34%, and 31.58%, respectively, with a majority of these mutations being heterozygous for resistance. In addition, five adult mosquitoes were identified has having synonymous mutations at site 1011[i. e. , AAT(asparagine)mutation to AAC(asparagine)]. Conclusions Our findings revealed the clear resistance of Cx. tritaeniorhynchus to deltamethrin in the Yancheng region of the Yellow Sea wetland, and the resistance phenotype and kdr frequency of Cx. tritaeniorhynchus in the wetland environment were comparable to those of Cx. tritaeniorhynchus in the wetland environment, thereby indicating that the resistance of different populations of Cx. tritaeniorhynchus was homogeneous under the pressure of different insecticide selection within and around the wetland. However, the underlying mechanisms need to be further studied.