Evaluating eDNA and eRNA metabarcoding for aquatic biodiversity assessment: From bacteria to vertebrates.

The monitoring and management of aquatic ecosystems depend on precise estimates of biodiversity. Metabarcoding analyses of environmental nucleic acids (eNAs), including environmental DNA (eDNA) and environmental RNA (eRNA), have garnered attention for their cost-effective and non-invasive biomonitoring capabilities. However, the accuracy of biodiversity estimates obtained through eNAs can vary among different organismal groups. Here we evaluate the performance of eDNA and eRNA metabarcoding across nine organismal groups, ranging from bacteria to terrestrial vertebrates, in three cross-sections of the Yangtze River, China. We observe robust complementarity between eDNA and eRNA data. The relative detectability of eNAs was notably influenced by major taxonomic groups and organismal sizes, with eDNA providing more robust signals for larger organisms. Both eDNA and eRNA exhibited similar cross-sectional and longitudinal patterns. However, the detectability of larger organisms declined in eRNA metabarcoding, possibly due to differential RNA release and decay among different organismal groups or sizes. While underscoring the potential of eDNA and eRNA in large river biomonitoring, we emphasize the need for differential interpretation of eDNA versus eRNA data. This highlights the importance of careful method selection and interpretation in biomonitoring studies.

An In Vivo Assessment of a Novel Temperature Control Flexible Ureteroscope System for Monitoring and Controlling Intrarenal Temperature During Flexible Ureteroscopy.

OBJECTIVE: To evaluate the efficacy of a novel temperature control flexible ureteroscope system in the precise monitoring and control of intrarenal temperature during ureteroscopy. METHODS: We developed a novel temperature control flexible ureteroscope system (PT-Scope), including a temperature-monitoring ureteroscope and a irrigation-suction platform for temperature regulation. A porcine thermometry model was established to observe temperature changes under varying holmium laser powers (10, 20, 30W) and irrigation rates (0, 20, 50 ml/min), utilizing PCN thermometry and PT-Scope measurements, with subsequent evaluation of temperature variations at different distances from the laser fiber tip. A porcine kidney stone model was established while porcine were randomly assigned to two groups: In the temperature control group, PT-Scope was connected to the irrigation-suction platform with temperature regulation, while in the non-temperature control group without temperature regulation. Comparative analysis was performed to evaluate differences in intrarenal temperature between the two groups. RESULTS: Across various laser powers and irrigation rates, the temperature measurement capability of the PT-Scope was precise, demonstrating consistency with PCN temperature measurements. The temperature obtained from the PT-Scope reflect the temperature approximately 0.05 cm away from the fiber tip, whereas temperatures close to fiber tip were significantly higher. The peak temperature of the temperature control group vs non-temperature control group were 31.70±2.609℃and 44.37±3.318℃, respectively (p<0.01). The mean temperature of the temperature control group vs non-temperature control group was 27.40±2.107℃ vs 35.9±1.921℃(p <0.01). CONCLUSION: PT-Scope has demonstrated the capability to precisely monitor and control intrarenal temperature within a safe threshold.

Citizen science in action: Time-resolved immunofluorescence-based field detection of antibiotics with portable analytical kit.

Citizen scientist-based environmental monitoring and public education are becoming increasingly popular. However, current technologies for antibiotic-based novel contaminant identification are still restricted to laboratory sample collection and analysis due to detection methodologies and apparatus limitations. This study developed a time-resolved immunofluorescence-based simultaneous field-based assay for ciprofloxacin (CIP) and enrofloxacin (ENR) that matches test results to geographic locations. The assay helps the public understand the potential levels of antibiotic exposures in their environments and helps them take appropriate action to reduce risk. The assay was developed using smartphones and social software in addition to rapid testing. The method uses a portable, low-cost analytical kit with a smartphone app to build a field-based detection platform for the detection and analysis of ENR and CIP in water and aquatic products. The methodological evaluation was good, with detection limits of 0.4 ng/mL and 0.5 ng/g for ENR in water and fish, and quantification limits of 1.2 ng/mL and 1.4 ng/g, with recoveries of 89.0 %-101.0 % and 78.0 %-97.0 %. For CIP in water and fish, the limits of detection were 0.3 ng/mL and 0.4 ng/g, the limits of quantification were 0.9 ng/mL and 1.2 ng/g, and the recoveries were 75.0 %-91.0 % and 72.0 %-89.0 %, both with coefficients of variation <15 %. These limits were sufficient to prevent the two antibiotics from crossing over during simultaneous detection. The assay was validated using real samples to assess the effectiveness of the assay platform in field deployments, and the results were consistent with those obtained through liquid chromatography-tandem mass spectrometry (LC-MS) and enzyme-linked immunoassay (ELISA) techniques. In addition, the TRFIA assay process requires less time, uses more portable instruments, and is less complex than traditional methods. This study provides a new scientific, accurate, and rapid detection method for antibiotic detection by citizen scientists, helping scientists to obtain a wider range of data and providing more opportunities to solve scientific problems.

Prevalence and impact of viral myocarditis in patients with severe fever with thrombocytopenia syndrome.

To explore the association and impact between viral myocarditis and mortality in patients with severe fever with thrombocytopenia syndrome. A dynamic analysis was conducted between fatal group and nonfatal group regarding the daily epidemiology data, clinical symptoms, and electrocardiogram (ECG), echocardiogram, and laboratory findings. Outcomes of patients with and without viral myocarditis were compared. The association between viral myocarditis and mortality was analyzed. Among 183 severe fever with thrombocytopenia syndrome patients, 32 were in the fatal group and 151 in the nonfatal group; there were 26 (81.25%) with viral myocarditis in the fatal group, 66 (43.70%) with viral myocarditis in the nonfatal group (p < 0.001), 79.35% of patients had abnormal ECG results. The abnormal rate of ECG in the fatal group was 100%, and in the nonfatal group was 74.83%. Univariate analysis found that the number of risk factors gradually increased on Day 7 of the disease course and reached the peak on Day 10. Combined with the dynamic analysis of the disease course, alanine aminotransferase, aspartate aminotransferase, creatine kinase, creatine kinase fraction, lactate dehydrogenase, hydroxybutyrate dehydrogenase, neutrophil count, serum creatinine, Na, Ca, carbon dioxide combining power, amylase, lipase, activated partial thromboplastin time and thrombin time had statistically significant impact on prognosis. The incidence of fever with thrombocytopenia syndrome combined with viral myocarditis is high, especially in the fatal group of patients. Viral myocarditis is closely related to prognosis and is an early risk factor. The time point for changes in myocarditis is Day 7 of the course of the disease.

Construction of an early differentiation diagnosis model for patients with severe fever with thrombocytopenia syndrome and hemorrhagic fever with renal syndrome.

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging infectious disease with a high mortality rate. Differentiating between SFTS and hemorrhagic fever with renal syndrome (HFRS) is difficult and inefficient. Retrospective analysis of the medical records of individuals with SFTS and HFRS was performed. Clinical and laboratory data were compared, and a diagnostic model was developed based on multivariate logistic regression analyzes. Receiver operating characteristic curve analysis was used to evaluate the diagnostic model. Among the 189 patients, 113 with SFTS and 76 with HFRS were enrolled. Univariate analysis revealed that more than 20 variables were significantly associated with SFTS. Multivariate logistic regression analysis revealed that gender, especially female gender (odds ratio [OR]: 4.299; 95% confidence interval [CI]: 1.163-15.887; p = 0.029), age ≥65 years (OR: 16.386; 95% CI: 3.043-88.245; p = 0.001), neurological symptoms (OR: 12.312; 95% CI: 1.638-92.530; p = 0.015), leukopenia (<4.0 × 109/L) (OR: 17.355; 95% CI: 3.920-76.839; p < 0.001), and normal Cr (OR: 97.678; 95% CI: 15.483-616.226; p < 0.001) were significantly associated with SFTS but not with HFRS. The area under the curve of the differential diagnostic model was 0.960 (95% CI: 0.936-0.984), which was significantly better than that of each single factor. In addition, the model exhibited very excellent sensitivity and specificity (92.9% and 85.5%, respectively). In cases where HFRS and SFTS are endemic, a diagnostic model based on five parameters, such as gender, age ≥65 years, neurological symptoms, leukopenia and normal Cr, will facilitate the differential diagnosis of SFTS and HFRS in medical institutions, especially in primary care settings.

Enhancing amphibian biomonitoring through eDNA metabarcoding.

Surveying biodiversity has taken a quantum leap with environmental DNA (eDNA) metabarcoding, an immensely powerful approach lauded for its efficiency, sensitivity, and non-invasiveness. This approach emerges as a game-changer for the elusive realm of endangered and rare species-think nocturnal, environmentally elusive amphibians. Here, we have established a framework for constructing a reliable metabarcoding pipeline for amphibians, covering primer design, performance evaluation, laboratory validation, and field validation processes. The Am250 primer, located on the mitochondrial 16S gene, was optimal for the eDNA monitoring of amphibians, which demonstrated higher taxonomic resolution, smaller species amplification bias, and more extraordinary detection ability compared to the other primers tested. Am250 primer exhibit an 83.8% species amplification rate and 75.4% accurate species identification rate for Chinese amphibians in the in silico PCR and successfully amplified all tested species of the standard samples in the in vitro assay. Furthermore, the field-based mesocosm experiment showed that DNA can still be detected by metabarcoding even days to weeks after organisms have been removed from the mesocosm. Moreover, field mesocosm findings indicate that eDNA metabarcoding primers exhibit different read abundances, which can affect the relative biomass of species. Thus, appropriate primers should be screened and evaluated by three experimental approaches: in silico PCR simulation, target DNA amplification, and mesocosm eDNA validation. The selection of a single primer set or multiple primers' combination should be based on the monitoring groups to improve the species detection rate and the credibility of results.

Pancreatic Injury Is Associated with Poor Prognosis in Severe Fever with Thrombocytopenia Syndrome.

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging infectious disease. Previous studies have primarily focused on the epidemiological and clinical characteristics of patients with SFTS, whereas pancreatic injury has received little attention. This study investigated the effects of pancreatic injury on the prognosis of patients with SFTS. A total of 156 patients diagnosed with SFTS between April 2016 and April 2022 were included in the analysis. Multivariate logistic regression analysis showed that pancreatic injury (odds ratio [OR] = 3.754, 95% confidence interval [CI]: 1.361-79.036, P = 0.024) and neurological symptoms (OR = 18.648, 95% CI: 4.921-70.668, P < 0.001) were independent risk factors for mortality. The receiver operating characteristic curve indicated that serum pancreatic enzymes were predictive of progression to death in patients with SFTS. The area under the curve (AUC) for amylase was 0.711, with an optimal cutoff value of 95.5 U/L, sensitivity of 96.4%, and specificity of 35.9%. Lipase had an AUC of 0.754, an optimal cutoff value of 354.75 U/L, sensitivity of 75%, and specificity of 67.2%. Thus, pancreatic injury was associated with a poor prognosis of SFTS and can be used as an important reference for SFTS determination and prognostic assessment.

Clinicopathological characteristics and its association with digestive system tumors of 1111 patients with Schistosomiasis japonica.

Schistosomiasis japonicum can cause different degrees of organ damage and complex human immune pathological reactions, which often invade the intestine and liver. The purpose of this study was to explore the pathological types and pathological changes of Schistosomiasis and their correlation with some digestive system tumors. Hematoxylin eosin staining was performed on the diseased tissues of 1111 Schistosomiasis cases. We counted the deposition sites of Schistosoma eggs, analyzed the pathological characteristics, and compared the clinicopathological characteristics of Schistosomiasis associated digestive system tumors and non-Schistosomiasis digestive system tumors. We found that Schistosoma japonicum can cause multi organ and multi system damage, with 469 cases of inflammation, 47 cases of adenoma, and 519 cases of adenocarcinoma. Other types include cysts, stromal tumors, malignant lymphomas, and neuroendocrine tumors. Schistosomiasis associated tumors, including gastric cancer, liver cancer, colon cancer and rectal cancer, were compared with non-Schistosomiasis tumors. There were significant differences in age, gender and tumor differentiation between the two groups. Our study shows Schistosomiasis is a systemic disease, causing multiple organ and system damage in the human body. Its clinicopathological types are diverse, and there may be a pathological change process of "Inflammation-adenoma-carcinoma". Schistosomiasis associated digestive system tumors differ from non-Schistosomiasis tumors in some clinicopathological features.

T Cell Activating Thermostable Self-Assembly Nanoscaffold Tailored for Cellular Immunity Antigen Delivery.

Antigen delivery based on non-virus-like particle self-associating protein nanoscffolds, such as Aquifex aeolicus lumazine synthase (AaLS), is limited due to the immunotoxicity and/or premature clearance of antigen-scaffold complex resulted from triggering unregulated innate immune responses. Here, using rational immunoinformatics prediction and computational modeling, we screen the T epitope peptides from thermophilic nanoproteins with the same spatial structure as hyperthermophilic icosahedral AaLS, and reassemble them into a novel thermostable self-assembling nanoscaffold RPT that can specifically activate T cell-mediated immunity. Tumor model antigen ovalbumin T epitopes and the severe acute respiratory syndrome coronavirus 2 receptor-binding domain are loaded onto the scaffold surface through the SpyCather/SpyTag system to construct nanovaccines. Compared to AaLS, RPT -constructed nanovaccines elicit more potent cytotoxic T cell and CD4+ T helper 1 (Th1)-biased immune responses, and generate less anti-scaffold antibody. Moreover, RPT significantly upregulate the expression of transcription factors and cytokines related to the differentiation of type-1 conventional dendritic cells, promoting the cross-presentation of antigens to CD8+ T cells and Th1 polarization of CD4+ T cells. RPT confers antigens with increased stability against heating, freeze-thawing, and lyophilization with almost no antigenicity loss. This novel nanoscaffold offers a simple, safe, and robust strategy for boosting T-cell immunity-dependent vaccine development.

Unsupervised biological integrity assessment by eDNA biomonitoring of multi-trophic aquatic taxa.

The biological integrity of global freshwater ecosystems is threatened by ever-increasing environmental stressors due to increased human activities, such as land-use change, eutrophication, toxic pollutants, overfishing, and exploitation. Traditional ecological assessments of lake or riverine ecosystems often require human supervision of a pre-selected reference area, using the current state of the reference area as the expected state. However, selecting an appropriate reference area has become increasingly difficult with the expansion of human activities. Here, an unsupervised biological integrity assessment framework based on environmental DNA metabarcoding without a prior reference area is proposed. Taxon richness, species dominance, co-occurrence network density, and phylogenetic distance were used to assess the aquatic communities in the Taihu Lake basin. Multi-gene metabarcoding revealed comprehensive biodiversity at multiple trophic levels including algae, protists, zooplankton, and fish. Fish sequences were mainly derived from 12S, zooplankton mainly from mitochondrial cytochrome C oxidase subunit I, and algae and protists mainly from 18S. There were significant differences in community composition among lakes, rivers, and reservoirs but no significant differences in the four fundamental biological indicators. The algal and zooplankton integrities were positively correlated with protist and fish integrities, respectively. Additionally, the algal integrity of lakes was found to be significantly lower than that of rivers. The unsupervised assessment framework proposed in this study allows different ecosystems, including the same ecosystem in different seasons, to adopt the same indicators and assessment methods, which is more convenient for environmental management and decision-making.

Destabilizing Effects of Environmental Stressors on Aquatic Communities and Interaction Networks across a Major River Basin.

Human-driven environmental stressors are increasingly threatening species survival and diversity of river systems worldwide. However, it remains unclear how the stressors affect the stability changes across aquatic multiple communities. Here, we used environmental DNA (eDNA) data sets from a human-dominated river in China over 3 years and analyzed the stability changes in multiple communities under persistent anthropogenic stressors, including land use and pollutants. First, we found that persistent stressors significantly reduced multifaceted species diversity (e.g., species richness, Shannon's diversity, and Simpson's diversity) and species stability but increased species synchrony across multiple communities. Second, the structures of interaction networks inferred from an empirical meta-food web were significantly changed under persistent stressors, for example, resulting in decreased network modularity and negative/positive cohesion. Third, piecewise structural equation modeling proved that the persistent stress-induced decline in the stability of multiple communities mainly depended upon diversity-mediated pathways rather than the direct effects of stress per se; specifically, the increase of species synchrony and the decline of interaction network modularity were the main biotic drivers of stability variation. Overall, our study highlights the destabilizing effects of persistent stressors on multiple communities as well as the mechanistic dependencies, mainly through reducing species diversity, increasing species synchrony, and changing interaction networks.

Citizen science meets eDNA: A new boom in research exploring urban wetland biodiversity.

•eDNA citizen science provides a comprehensive picture of the biodiversity.•Non-native species reduced the local fish diversity in urban wetlands.•Expanding water area can improve wetland biodiversity.

Asymmetric Assembly and Self-Adjuvanted Antigen Delivery Platform for Improved Antigen Uptake and Antitumor Effect.

The development of effective tumor vaccines is an important direction in the field of cancer prevention/immunotherapy. Efficient antigen delivery is essential for inducing effective antitumor responses for tumor vaccines. Lumazine synthase (BLS) from Brucella spp. is a decameric protein with delivery and adjuvant properties, but its application in tumor vaccines is limited. Here, we developed an antigen delivery platform by combining a BLS asymmetric assembly and the Plug-and-Display system of SpyCatcher/SpyTag. An asymmetric assembly system consisting of BLSke and BLSdr was developed to equally assemble two molecules. Then, the MHC-I-restricted ovalbumin peptide (OVA(257-264) SIINFEKL) was conjugated with BLSke, and a cell-penetrating peptide (CPP) KALA was conjugated with BLSdr using the SpyCatcher/SpyTag system. KALA modification enhanced internalization of OVA peptides by DCs as well as promoted the maturation of DCs and the cross-presentation of SIINFEKL. Moreover, the immunotherapy of a KALA-modified vaccine suppressed tumor growth and enhanced CD8+ T cell responses in E.G7-OVA tumor-bearing mice. In the prophylactic model, KALA-modified vaccination showed the most significant protective effect and significantly prolonged the survival period of tumor challenged mice. In conclusion, the asymmetric assembly platform equally assembles two proteins or peptides, avoiding their spatial or functional interference. This asymmetric assembly and Plug-and-Display technology provide a universal platform for rapid development of personalized tumor vaccines.

Multiple regression analysis of risk factors related to radiation pneumonitis.

BACKGROUND: Radiation pneumonitis (RP) is a severe complication of thoracic radiotherapy that may lead to dyspnea and lung fibrosis, and negatively affects patients' quality of life. AIM: To carry out multiple regression analysis on the influencing factors of radiation pneumonitis. METHODS: Records of 234 patients receiving chest radiotherapy in Huzhou Central Hospital (Huzhou, Zhejiang Province, China) from January 2018 to February 2021, and the patients were divided into either a study group or a control group based on the presence of radiation pneumonitis or not. Among them, 93 patients with radiation pneumonitis were included in the study group and 141 without radiation pneumonitis were included in the control group. General characteristics, and radiation and imaging examination data of the two groups were collected and compared. Due to the statistical significance observed, multiple regression analysis was performed on age, tumor type, chemotherapy history, forced vital capacity (FVC), forced expiratory volume in the first second (FEV1), carbon monoxide diffusion volume (DLCO), FEV1/FVC ratio, planned target area (PTV), mean lung dose (MLD), total number of radiation fields, percentage of lung tissue in total lung volume (vdose), probability of normal tissue complications (NTCP), and other factors. RESULTS: The proportions of patients aged ≥ 60 years and those with the diagnosis of lung cancer and a history of chemotherapy in the study group were higher than those in the control group (P < 0.05); FEV1, DLCO, and FEV1/FVC ratio in the study group were lower than those in the control group (P < 0.05), while PTV, MLD, total field number, vdose, and NTCP were higher than in the control group (P < 0.05). Logistic regression analysis showed that age, lung cancer diagnosis, chemotherapy history, FEV1, FEV1/FVC ratio, PTV, MLD, total number of radiation fields, vdose, and NTCP were risk factors for radiation pneumonitis. CONCLUSION: We have identified patient age, type of lung cancer, history of chemotherapy, lung function, and radiotherapy parameters as risk factors for radiation pneumonitis. Comprehensive evaluation and examination should be carried out before radiotherapy to effectively prevent radiation pneumonitis.

Using Defect Control To Break the Stability-Activity Trade-Off in Enzyme Immobilization via Competitive Coordination.

Immobilization of enzymes within metal-organic frameworks is a powerful strategy to enhance the long-term usability of labile enzymes. However, the thus-confined enzymes suffer from the trade-off between enhanced stability and reduced activity because of the contradiction between the high crystallinity and the low accessibility. Here, by taking laccase and zeolitic imidazolate framework-8 (ZIF-8) as prototypes, we disclosed an observation that the stability-activity trade-off could be solved by controlling the defects via competitive coordination. Owing to the presence of competitive coordination between laccase and the ligand precursor of ZIF-8, there existed a three-stage process in the de novo encapsulation: nucleation-crystallization-recrystallization. Our results show that the biocomposites collected before the occurrence of recrystallization possessed both increased activity and enhanced stability. The findings here shed new light on the control of defects through the subtle use of competitive coordination, which is of great significance for the engineering application of biomacromolecules.

The rubber tree RALF peptide hormone and its receptor protein kinase FER implicates in rubber production.

RAPID ALKALINIZATION FACTORs (RALFs), which are secreted peptides serving as extracellular signals transduced to the inside of the cell, interact with the receptor-like kinase FERONIA (FER) and participates in various biological pathways. Here, we identified 23 RALF and 2 FER genes in Hevea brasiliensis (para rubber tree), and characterized their expression patterns in different tissues, across the process of leaf development, and in response to the rubber yield-stimulating treatments of tapping and ethylene. Four Hevea latex (the cytoplasm of rubber-producing laticifers)-abundant RALF isoforms, HbRALF19, HbRALF3, HbRALF22, and HbRALF16 were listed with descending expression levels. Of the four HbRALFs, expressions of HbRALF3 were markedly regulated in an opposite way by the treatments of tapping (depression) and ethylene (stimulation). All of the four latex-abundant RALFs specifically interacted with the extracellular domain of HbFER1. Transgenic Arabidopsis plants overexpressing these HbRALFs displayed phenotypes similar to those reported for AtRALFs, such as shorter roots, smaller plant architecture, and delayed flowering. The application of HbRALF3 and HbRALF19 recombinant proteins significantly reduced the pH of Hevea latex, an important factor regulating latex metabolism. An in vitro rubber biosynthesis assay in a mixture of latex cytosol (C-serum) revealed a positive role of HbFER1 in rubber biosynthesis. Taken together, these data provide evidence for the participation of the HbRALF-FER module in rubber production.

A specific reverse complement sequence for distinguishing Brucella canis from other Brucella species.

Canine brucellosis is primarily caused by Brucella canis, but other Brucella species can also cause the disease. Identifying sequences specific to B. canis and establishing PCR assays that can distinguish between B. canis and other Brucella species is essential to determine the etiology of canine brucellosis and the source of infection and to achieve effective control. We analyzed the gaps and SNPs of genomes I and II from B. canis strain RM6/66 and B. melitensis strain 16M using the Mauve genome alignment software, and the specificity of each of these differential regions was analyzed by BLAST. A 132 bp specific sequence was found between the DK60_915 (glycosyl hydrolase 108 family protein) and DK60_917 (aldose 1-epimerase) loci in B. canis chromosome 1. Further comparative analysis revealed that this is a reverse complement sequence between B. canis and other Brucella species. Then, three primers were designed based on the sequence that could detect B. canis with a 310 bp amplification product or other Brucella species with a 413 bp product. The PCR based on these primers had reasonable specificity and a sensitivity of 100 copies of Brucella DNA. The detection results for the blood samples of the aborted dogs showed a favorable accordance with the Bruce-ladder multiplex PCR assay. In conclusion, we found a specific reverse complement sequence between B. canis and other Brucella and developed a PCR method that allows a more comprehensive identification of the pathogen involved in canine brucellosis. These findings provide an effective means for preventing and controlling brucellosis.

Trehalose 6-Phosphate/SnRK1 Signaling Participates in Harvesting-Stimulated Rubber Production in the Hevea Tree.

Trehalose 6-phosphate (T6P), the intermediate of trehalose biosynthesis and a signaling molecule, affects crop yield via targeting sucrose allocation and utilization. As there have been no reports of T6P signaling affecting secondary metabolism in a crop plant, the rubber tree Hevea brasiliensis serves as an ideal model in this regard. Sucrose metabolism critically influences the productivity of natural rubber, a secondary metabolite of industrial importance. Here, we report on the characterization of the T6P synthase (TPS) gene family and the T6P/SNF1-related protein kinase1 (T6P/SnRK1) signaling components in Hevea laticifers under tapping (rubber harvesting), an agronomic manipulation that itself stimulates rubber production. A total of fourteen TPS genes were identified, among which a class II TPS gene, HbTPS5, seemed to have evolved with a function specialized in laticifers. T6P and trehalose increased when the trees were tapped, this being consistent with the observed enhanced activities of TPS and T6P phosphatase (TPP) and expression of an active TPS-encoding gene, HbTPS1. On the other hand, SnRK1 activities decreased, suggesting the inhibition of elevated T6P on SnRK1. Expression profiles of the SnRK1 marker genes coincided with elevated T6P and depressed SnRK1. Interestingly, HbTPS5 expression decreased significantly with the onset of tapping, suggesting a regulatory function in the T6P pathway associated with latex production in laticifers. In brief, transcriptional, enzymatic, and metabolic evidence supports the participation of T6P/SnRK1 signaling in rubber formation, thus providing a possible avenue to increasing the yield of a valuable secondary metabolite by targeting T6P in specific cells.