Genome-wide association analysis of eggshell color of an F2 generation population reveals candidate genes in chickens

Eggshell color is an important visual characteristic that affects consumer preferences for eggs. Eggshell color, which has moderate to high heritability, can be effectively enhanced through molecular marker selection. Various studies have been conducted on eggshell color at specific time points. However, few longitudinal data are available on eggshell color. Therefore, the objective of this study was to investigate eggshell color using the Commission International de L'Eclairage L*a*b* system with multiple measurements at different ages (age at the first egg and at 32, 36, 40, 44, 48, 52, 56, 60, 66, and 72 weeks) within the same individuals from an F2 resource population produced by crossing White Leghorn and Dongxiang Blue chicken. Using an Affymetrix 600 single nucleotide polymorphism (SNP) array, we estimated the genetic parameters of the eggshell color trait, performed genome-wide association studies (GWASs), and screened for the potential candidate genes. The results showed that pink-shelled eggs displayed a significant negative correlation between L* values and both a* and b* values. Genetic heritability based on SNPs showed that the heritability of L*, a*, and b* values ranged from 0.32 to 0.82 for pink-shelled eggs, indicating a moderate to high level of genetic control. The genetic correlations at each time point were mostly above 0.5. The major-effect regions affecting the pink eggshell color were identified in the 10.3-13.0 Mb interval on Gallus gallus chromosome 20, and candidate genes were selected, including SLC35C2, PCIF1, and SLC12A5. Minor effect polygenic regions were identified on chromosomes 1, 6, 9, 12, and 15, revealing 11 candidate genes, including MTMR3 and SLC35E4. Members of the solute carrier family play an important role in influencing eggshell color. Overall, our findings provide valuable insights into the phenotypic and genetic aspects underlying the variation in eggshell color. Using GWAS analysis, we identified multiple quantitative trait loci (QTLs) for pink eggshell color, including a major QTL on chromosome 20. Genetic variants associated with eggshell color may be used in genomic breeding programs.

Effects of manganese glycine on eggshell quality, eggshell ultrastructure, and elemental deposition in aged laying hens

Poor eggshell quality of eggs laid by aged laying hens is the major problem affecting the length of the rearing period in the laying hen industry. Trace elements are required and play vital roles in the eggshell quality of laying hens. Appropriate dose of organic microelements is environmentally friendly and sufficient to satisfy the needs of hens because of their greater bioavailability and lower excretion than inorganic forms. The aim of this experiment was to investigate the effects of manganese (Mn) glycine (MG) on eggshell quality, elemental deposition, and eggshell ultrastructure in aged laying hens. A total of 720 Hy-Line Brown hens 70 weeks old were assigned equally to four groups with six replicates of 30 birds each. The hens were fed basal diets (without Mn supplementation) supplemented with 120 mg/kg of Mn from manganese sulfate monohydrate (MSM), or 40, 80, or 120 mg/kg Mn from MG for 12 weeks. Dietary supplementation with 80 mg/kg Mn from MG resulted in the greatest eggshell strength after 6 weeks of treatment (P = 0.047), and in greater eggshell strength than observed in the MSM control after 12 weeks of treatment (P = 0.025). After 12 weeks of treatment, the eggs of hens in the MG groups showed lower mammillary layer thickness in the blunt end, equator, and acute end than observed in the MSM control group (P < 0.001). With the exception of the blunt ends of eggs from hens in the 120 mg/kg MG group, the eggs of hens in the MG groups, compared with the MSM control group, exhibited a lower mammillary layer ratio, and greater palisade layer ratio and effective layer ratio in the blunt end, equator, and acute end (P < 0.001). Dietary supplementation with 80 mg/kg Mn from MG, compared with the MSM control and 40 and 120 mg/kg MG, resulted in the greatest palisade layer thickness and effective layer thickness, and the lowest mammillary layer thickness in the equator (P < 0.001, P = 0.001, P < 0.001, respectively). Furthermore, supplementation with 80 mg/kg Mn from MG exhibited the greatest ratio of the palisade layer and effective layer, and the lowest mammillary layer ratio in the blunt end and equator (all P < 0.001). The Mn content of eggshells in hens-fed diets supplemented with 80 and 120 mg/kg Mn from MG was greater than that in the MSM control and 40 mg/kg MG groups (P = 0.035). Dietary supplementation with 80 or 120 mg/kg Mn from MG resulted in greater tibia Mn content than observed in the 40 mg/kg MG group (P = 0.019), and greater yolk Mn content than observed in the 40 mg/kg MG and MSM control groups (P = 0.018). In conclusion, dietary supplementation with 80 mg/kg Mn from MG, compared with the MSM control (120 mg/kg Mn), may increase the deposition efficiency of Mn, alter eggshell elemental composition, improve eggshell ultrastructure, and enhance eggshell strength in aged laying hens.

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at 3.67±0.05±0.15  PeV. Below the knee, the spectral index is found to be -2.7413±0.0004±0.0050, while above the knee, it is -3.128±0.005±0.027, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -0.1200±0.0003±0.0341. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components.

Effects of exogenous energy on synthesis of steroid hormones and expression characteristics of the CREB/StAR signaling pathway in theca cells of laying hen

Energy and the cAMP-response element binding protein (CREB)/steroidogenic acute regulatory protein (StAR) signaling pathway play important roles in steroid hormone production and follicular development in hens. This present study aimed to investigate the effects of exogenous energy on the synthesis of steroid hormones and the expression characteristics of the CREB/StAR signaling pathway in theca cells of laying hen. The primary theca cells of small yellow follicles were randomly divided into 6 treatments and cultured in medium with glucose concentrations of 1, 1.5, 3, 4.5, 6, and 7.5 mg/mL for 48 h. It was found that growth was robust and cell outlines were clear when cells were cultured with 1, 1.5, 3, and 4.5 mg/mL glucose, but cell viability was diminished and cell density decreased after exposure to glucose at 6 and 7.5 mg/mL for 48 h. Cell viability showed an increasing and then decreasing quadratic response to increasing glucose concentration in culture (r2 = 0.688, P < 0.001). The cell viability of theca cells cultured with 4.5 mg/mL glucose was greater than those cultured with 1, 1.5, 6, and 7.5 mg/mL glucose (P < 0.05). The concentration of estradiol in the medium containing 3 mg/mL glucose was higher than in medium containing 1, 1.5, and 6 mg/mL glucose (P < 0.05). There was an increasing and then decreasing quadratic correlation between progesterone concentrations and glucose concentrations (r2 = 0.522, P = 0.002). The concentration of progesterone in medium with 4.5 mg/mL glucose was higher than in medium with 1 and 7.5 mg/mL glucose (P < 0.05). There was an increasing and then decreasing quadratic correlation between the relative expression of CREB1 (r2 = 0.752, P < 0.001), StAR (r2 = 0.456, P = 0.002), CYP1B1 (r2 = 0.568, P < 0.001), and 3β-HSD (r2 = 0.319, P = 0.018) in theca cells of laying hens and glucose concentrations after treatment with different glucose concentrations for 48 h. After treatment with 4.5 mg/mL glucose, the expression of StAR, CYP1B1, and 3β-HSD genes were increased compared to treatment with 1, 1.5, 3, 6, and 7.5 mg/mL glucose (P < 0.001). There was an increasing and then decreasing quadratic correlation between glucose concentrations and protein expression of CREB1 (r2 = 0.819, P < 0.001), StAR (r2 = 0.844, P < 0.001), 3β-HSD (r2 = 0.801, P < 0.001), and CYP11A1 (r2 = 0.800, P < 0.001) in theca cells of laying hens. The protein expression of CREB1, StAR, and 3β-HSD in theca cells cultured with 4.5 mg/mL glucose was higher than in other groups (P < 0.001). The results indicate that the appropriate glucose concentration (4.5 mg/mL) can improve the synthesis of steroid hormones in theca cells of laying hens through the upregulation of key genes and proteins in the CREB/StAR signaling pathway.

Mineral composition controls the stabilization of microbially derived carbon and nitrogen in soils: Insights from an isotope tracing model

Evidence is emerging that microbial products and residues (necromass) contribute greatly to stable soil organic matter (SOM), which calls for the necessity of separating the microbial necromass from other SOM pools in models. However, the understanding on how microbial necromass stabilizes in soil, especially the mineral protection mechanisms, is still lacking. Here, we incubated 13 C- and 15 N-labelled microbial necromass in a series of artificial soils varying in clay minerals and metal oxides. We found the mineralization, adsorption and desorption rate constants of necromass nitrogen were higher than those of necromass carbon. The accumulation rates of necromass carbon and nitrogen in mineral-associated SOM were positively correlated with the specific surface area of clay minerals. Our results provide direct evidence for the protection role of mineral in microbial necromass stabilization and provide a platform for simulating microbial necromass separately in SOM models.

Relationships between Phyllosphere Bacterial Communities and Leaf Functional Traits in a Temperate Forest

As a vital component of biodiversity, phyllosphere bacteria in forest canopy play a critical role in maintaining plant health and influencing the global biogeochemical cycle. There is limited research on the community structure of phyllosphere bacteria in natural forests, which creates a gap in our understanding of whether and/or how phyllosphere bacteria are connected to leaf traits of their host. In this study, we investigated the bacterial diversity and composition of the canopy leaves of six dominant tree species in deciduous broad-leaved forests in northeastern China, using high-throughput sequencing. We then compare the differences in phyllosphere bacterial community structure and functional genes of dominant tree species. Fourteen key leaf functional traits of their host trees were also measured according to standard protocols to investigate the relationships between bacterial community composition and leaf functional traits. Our result suggested that tree species with closer evolutionary distances had similar phyllosphere microbial alpha diversity. The dominant phyla of phyllosphere bacteria were Proteobacteria, Actinobacteria, and Firmicutes. For these six tree species, the functional genes of phyllosphere bacteria were mainly involved in amino acid metabolism and carbohydrate metabolism processes. The redundancy and envfit analysis results showed that the functional traits relating to plant nutrient acquisition and resistance to diseases and pests (such as leaf area, isotope carbon content, and copper content) were the main factors influencing the community structure of phyllosphere bacteria. This study highlights the key role of plant interspecific genetic relationships and plant attributes in shaping phyllosphere bacterial diversity.

Diversity and biogeography of plant phyllosphere bacteria are governed by latitude-dependent mechanisms

Predicting and managing the structure and function of plant microbiomes requires quantitative understanding of community assembly and predictive models of spatial distributions at broad geographic scales. Here, we quantified the relative contribution of abiotic and biotic factors to the assembly of phyllosphere bacterial communities, and developed spatial distribution models for keystone bacterial taxa along a latitudinal gradient, by analyzing 16S rRNA gene sequences from 1453 leaf samples taken from 329 plant species in China. We demonstrated a latitudinal gradient in phyllosphere bacterial diversity and community composition, which was mostly explained by climate and host plant factors. We found that host-related factors were increasingly important in explaining bacterial assembly at higher latitudes while nonhost factors including abiotic environments, spatial proximity and plant neighbors were more important at lower latitudes. We further showed that local plant-bacteria associations were interconnected by hub bacteria taxa to form metacommunity-level networks, and the spatial distribution of these hub taxa was controlled by hosts and spatial factors with varying importance across latitudes. For the first time, we documented a latitude-dependent importance in the driving factors of phyllosphere bacteria assembly and distribution, serving as a baseline for predicting future changes in plant phyllosphere microbiomes under global change and human activities.

Measurement of Ultra-High-Energy Diffuse Gamma-Ray Emission of the Galactic Plane from 10 TeV to 1 PeV with LHAASO-KM2A

The diffuse Galactic γ-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this Letter, we report the measurements of diffuse γ rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner (15°10  TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of -2.99±0.04, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of ∼3 than the prediction. A similar spectrum with an index of -2.99±0.07 is found in the outer Galaxy region, and the absolute flux for 10≲E≲60  TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations.

Effects of energy restriction during growing phase on the productive performance of Hyline Brown laying hens aged 6 to 72 wk

The aim of this study was to assess the effects of energy-restricted feeding during growing phase on the productive performance of Hyline Brown laying hens aged 6 to 72 wk. A total of 720 six-week-old layer chicks were allocated equally to 3 groups with 6 replicates of 40 pullets each, and were fed 1 of 3 diets that were nutritionally similar except for the apparent metabolizable energy corrected for nitrogen (AMEn) content. At the age of 6 to 17 wk, the pullets in the control group were given diet with 2,850 kcal/kg AMEn, and were fed ad libitum. The levels of AMEn in diet of pullets in the experimental groups were 90% (2,565 [2,850 × 90%] kcal/kg) and 80% (2,280 [2,850 × 80%] kcal/kg) of that in control group, and the daily amount of feed was restricted to the absolute quantity of the diet consumed by pullets in control group. At the age of 18 to 72 wk, all the hens were fed with the same diets ad libitum. As energy restriction increased in the growing phase, body weight (BW) dropped at the ages of 12 and 15 to 23 wk (at 23 wk: P = 0.001; at other ages: P < 0.001), but it showed no significant difference at 24 wk (P = 0.071). At 20 wk, restricting energy induced a delay in the development of sexual organs, including the ovary stroma, oviduct, and small yellow follicle (P < 0.05), as well as a delay in sexual maturity (P < 0.05). Consequently, the laying rate in the first and second periods dropped linearly (P = 0.046, 0.030, and 0.038, P < 0.001, respectively). The coefficient of variation (CV) in the BW at 19, 20, and 21 wk (P = 0.040, 0.023, and 0.042, respectively), the CV of age at first egg (P < 0.001), and CV of individual egg number at age 18 to 72 wk (P < 0.001) decreased linearly. There was a linear increase in the laying rate of hens in the later periods (at age 32-72 wk, P < 0.05), as well as in the average total egg number per hen and average laying rate at the age of 18 to 72 wk (P = 0.006). The average egg mass also showed a linear increase with increasing levels of energy restriction (P < 0.001). In summary, although appropriate energy restriction during growing phase delayed sexual maturity and sexual organ development in early-laying Hyline Brown pullets, it improved uniformity of BW, age at first egg laying, and individual egg number, and increased egg number per hen, laying rate, average egg mass, and number of settable eggs from 18 to 72 wk of age.

Ontogeny influences tree growth response to soil fertility and neighbourhood crowding in an old-growth temperate forest

BACKGROUND AND AIMS

Abiotic and biotic factors simultaneously affect tree growth and thus shape community structure and dynamics. In particular, trees of different size classes show different growth responses to soil nutrients and neighbourhood crowding, but our understanding of how species' joint responses to these factors vary between size classes remains limited in multi-storied temperate forests. Here, we investigated size class differences in tree growth response to soil gradients and neighbourhood crowding in an old-growth temperate forest.

METHODS

We combined growth data over 15 years from 38 902 individuals of 42 tree species with trait data in a 25-ha temperate forest plot in northeast China. We built hierarchical Bayesian models of tree growth to examine the effects of soil gradients and neighbourhood crowding between size classes and canopy types.

KEY RESULTS

We found that soil and neighbours mainly acted separately in shaping tree growth in small and large trees. Soil total nitrogen and phosphorus increased tree growth in small trees, in particular of understorey species, but not in large trees. Neighbours reduced tree growth in both tree size classes, with stronger effects on large than small trees, and on canopy than understorey species. Furthermore, small trees with higher specific leaf area grew faster in fertile soils, and small trees with less seed mass grew faster in crowded environments. Large trees with higher specific leaf area, specific root length and less seed mass grew faster in crowded environments, while these traits had limited influence on tree growth response to soil gradients.

CONCLUSIONS

Our study highlights the importance of size class in modulating the response of tree growth to soil and neighbours, and the differential role of species canopy types and functional traits in capturing these effects in large vs. small trees.

Spatial Distribution of Pinus koraiensis Trees and Community-Level Spatial Associations in Broad-Leaved Korean Pine Mixed Forests in Northeastern China

Investigating the spatial distributions and associations of tree populations provides better insights into the dynamics and processes that shape the forest community. Korean pine (Pinus koraiensis) is one of the most important tree species in broad-leaved Korean pine mixed forests (BKMFs), and little is known about the spatial point patterns of and associations between Korean pine and community-level woody species groups such as coniferous and deciduous trees in different developmental stages. This study investigated the spatial patterns of Korean pine (KP) trees and then analyzed how the spatial associations between KP trees and other tree species at the community level vary in different BKMFs. Extensive data collected from five relatively large sample plots, covering a substantial area within the natural distribution range of KP in northeastern China, were utilized. Uni- and bivariate pair correlation functions and mark correlation functions were applied to analyze spatial distribution patterns and spatial associations. The DBH (diameter at breast height) histogram of KP trees in northeastern China revealed that the regeneration process was very poor in the Changbai Mountain (CBS) plot, while the other four plots exhibited moderate or expanding population structures. KP trees were significantly aggregated at scales up to 10 m under the HPP null model, and the aggregation scales decreased with the increase in size classes. Positive or negative spatial associations were observed among different life stages of KP trees in different plots. The life history stages of the coniferous tree group showed positive spatial associations with KP saplings and juvenile trees at small scales, and spatial independence or negative correlations with larger KP trees at greater scales. All broad-leaved tree groups (canopy, middle, and understory layers) exhibited only slightly positive associations with KP trees at small scales, and dominant negative associations were observed at most scales. Our results demonstrate that mature KP trees have strong importance in the spatial patterns of KP populations, and site heterogeneity, limited seed dispersal, and interspecific competition characterize the spatial patterns of KP trees and community-level spatial associations with respect to KP trees, which can serve as a theoretical basis for the management and restoration of BKMFs in northeastern China.

[Clinical and neuroelectrophysiological characteristics of primary peripheral nerve hyperexcitability syndrome]

Objective: To investigate the clinical and neuroelectrophysiological characteristics of patients with primary peripheral nerve hyperexcitability syndrome (PNHS). Methods: The clinical data of 20 patients who were diagnosed with PNHS in Beijing Tiantan Hospital from April 2016 to January 2023 were retrospectively collected. All patients underwent neuroelectrophysiological examinations. Clinical and electrophysiological characteristics were compared between the antibody positive and antibody negative groups, according to serum and cerebrospinal fluid anti-contactin-associated protein-like 2 (CASPR2) and/or anti-leucine-rich glioma-inactivated protein 1 (LGI-1) antibodies. Results: There were 12 males and 8 females, with a mean age of (44.0±17.2) years and the disease course of [M (Q₁, Q₃)] 2.3 (1.1, 11.5) months. Motor symptoms included fasciculations, myokymia, muscle pain, cramps, and stiffness. These symptoms were commonly seen in the lower limbs (17 patients), followed by upper limbs (11 patients), face (11 patients) and trunk (9 patients). Nineteen (19/20) patients had sensory abnormalities and/or autonomic dysfunction, 13 patients had central nervous system involvement, and 5 patients had concomitant lung cancer or thymic lesions. The characteristic spontaneous potentials on needle electromyography (EMG) were myokymia potential (19 patients), fasciculation potential (12 patients), spastic potential (3 patients), neuromyotonic potential (1 patients), etc, which were commonly seen in the lower limb muscles, especially the gastrocnemius muscle(12 patients). After-discharge potential was found in 8 patients, and 7 were in the tibial nerve. Seven patients had positive serum anti-CASPR2 antibodies, and 3 of them had concomitant anti-LGI1 antibodies. And 1 patient had positive serum anti-LGI1 antibody alone. Compared with patients in the antibody negative group (n=12), the patients who had anti-VGKC complex antibodies (n=8) had a shorter course of disease [M (Q₁, Q₃): 1.8 (1, 2) months vs 9.5 (3.3, 20.3) months, P=0.012], higher incidence of after-discharge potential (6/8 vs 2/12, P=0.019). The immunotherapy regimen (multi-dru, single-drug, no immunotherapy: 6, 2, 0 patients) in antibody-positive patients was different from the antibody-negative group (3, 6, 3 patients, U=21.00, P=0.023). Conclusions: The symptoms of motor nerve hyperexcitation, characteristic EMG spontaneous potentials and after-discharge potentials in PNHS patients are most commonly seen in the lower limbs. Attention should be paid to concomitant sensory and autonomic nerve hyperexcitation. PNHS patients with positive serum anti-CASPR2 antibodies may require immunotherapy with multiple drugs.

[Epidemiological characteristics and risk factors of sepsis development and death in patients with extremely severe burns]

Objective: To explore the epidemiological characteristics and risk factors of sepsis development and death in patients with extremely severe burns. Methods: A retrospective case series study was conducted. From January 2017 to December 2021, 135 patients with extremely severe burns who met the inclusion criteria were admitted to the Department of Burn and Wound Repair of the Second Affiliated Hospital of Zhejiang University School of Medicine, including 100 males and 35 females, aged 18-84 years. The incidence and diagnosis time of sepsis, the rate of positive microbial culture of blood samples (hereinafter referred to as positive blood culture), and the mortality rate of all patients, as well as the incidence of sepsis and the pathogen of infection in patients with positive blood culture were recorded (statistically analyzed with chi-square test or Fisher's exact probability test). According to the occurrence of sepsis, all patients were divided into sepsis group (58 cases) and non-sepsis group (77 cases), and the gender, age, body mass index, history of hypertension, history of diabetes, combination of inhalation injury, burn site, burn type, total burn area, and combined injury of patients were compared between the two groups. According to the outcome, all patients were divided into death group (37 cases) and survival group (98 cases), and the aforementioned data grouped according to sepsis as well as the stability of shock period and the combination of sepsis of patients were compared between the two groups. The aforementioned data between two groups were statistically analyzed with univariate analysis of independent sample t test, Wilcoxon rank-sum test, Mann-Whitney U test, chi-square test, or Fisher's exact probability test. Factors with P<0.1 were selected for multivariate logistic regression analysis to screen independent risk factors of sepsis and death in patients with extremely severe burns. Results: Among all patients, the incidence of sepsis was 42.96% (58/135), the diagnosis time of sepsis was 14 (7, 24) d after injury, the positive blood culture rate was 62.22% (84/135), and the mortality rate was 27.41% (37/135). The incidence of sepsis of patients with positive blood culture was 69.05% (58/84). The top 5 pathogenic bacteria in the detection rate of septic patients with positive blood culture were Acinetobacter baumannii, Klebsiella pneumoniae, Pseudomonas aeruginosa, Staphylococcus aureus, and Enterobacter cloacae, ranking from high to low, and the proportion of Acinetobacter baumannii infected was significantly higher than that of non-septic patients with positive blood culture (χ2=7.49, P<0.05). Compared with those in non-sepsis group, the proportion of combination of inhalation injury, the proportion of perineal burns, and the total burn area of patients in sepsis group increased significantly (with χ2 values of 11.08 and 17.47, respectively, Z=5.68, P<0.05), while the other indicators did not change significantly (P>0.05). Multivariate logistic regression analysis showed that combination of inhalation injury, total burn area ≥80% total body surface area (TBSA), and perineal burns were independent risk factors for patients with extremely severe burns developing sepsis (with odds ratios of 3.15, 7.24, and 3.24, respectively, with 95% confidence intervals of 1.07 to 9.29, 1.79 to 29.34, and 1.21 to 8.68, respectively, P<0.05). Compared with those in survival group, the proportion of combination of inhalation injury, the proportion of perineal burns, and the proportion of combination of sepsis (with χ2 values of 6.55, 11.64, and 22.26, respectively, P values all <0.05), total burn area (Z=4.25, P<0.05), and proportion of instability of shock period (P<0.05) of patients in death group all increased significantly, while the other indicators did not change significantly (P>0.05). Multivariate logistic regression analysis showed that the instability of shock period and combination of sepsis were independent risk factors for death of patients with extremely severe burns (with odds ratios of 4.87 and 3.45, respectively, with 95% confidence intervals of 1.21 to 19.57 and 1.28 to 9.33, respectively, P<0.05). Conclusions: Patients with extremely severe burns have a high incidence of sepsis and a high mortality rate. The peak period of sepsis onset is 2 weeks after injury, with Acinetobacter baumannii as the most prominent infectious pathogen. Combination of inhalation injury, total burn area ≥80% TBSA, and perineal burns are independent risk factors for extremely severe burn patients complicated with sepsis, and combination of sepsis and instability of shock period are independent risk factors for death of patients with extremely severe burns.

A tera-electron volt afterglow from a narrow jet in an extremely bright gamma-ray burst

Some gamma-ray bursts (GRBs) have a tera-electron volt (TeV) afterglow, but the early onset of this has not been observed. We report observations with the Large High Altitude Air Shower Observatory of the bright GRB 221009A, which serendipitously occurred within the instrument field of view. More than 64,000 photons >0.2 TeV were detected within the first 3000 seconds. The TeV flux began several minutes after the GRB trigger, then rose to a peak about 10 seconds later. This was followed by a decay phase, which became more rapid ~650 seconds after the peak. We interpret the emission using a model of a relativistic jet with half-opening angle ~0.8°. This is consistent with the core of a structured jet and could explain the high isotropic energy of this GRB.

[Effects of three-dimensional bioprinting antibacterial hydrogel on full-thickness skin defect wounds in rats]

Objective: To explore the effects of three-dimensional (3D) bioprinting gelatin methacrylamide (GelMA) hydrogel loaded with nano silver on full-thickness skin defect wounds in rats. Methods: The experimental research method was adopted. The morphology, particle diameter, and distribution of silver nanoparticles in nano silver solution with different mass concentrations and the pore structure of silver-containing GelMA hydrogel with different final mass fractions of GelMA were observed by scanning electron microscope and the pore size was calculated. On treatment day 1, 3, 7, and 14, the concentration of nano silver released from the hydrogel containing GelMA with final mass fraction of 15% and nano silver with final mass concentration of 10 mg/L was detected by mass spectrometer. At 24 h of culture, the diameters of inhibition zone of GelMA hydrogel containing final mass concentration of 0 (no nano silver), 25, 50, and 100 mg/L nano silver against Staphylococcus aureus and Escherichia coli were detected. Fibroblasts (Fbs) and adipose stem cells (ASCs) were isolated respectively by enzymatic digestion using the discarded prepuce after circumcision from a 5-year-old healthy boy who was treated in the Department of Urology of the Second Affiliated Hospital of Zhejiang University School of Medicine in July 2020, and the discarded fat tissue after liposuction from a 23-year-old healthy woman who was treated in the Department of Plastic Surgery of the Hospital in July 2020. The Fbs were divided into blank control group (culture medium only), 2 mg/L nano sliver group, 5 mg/L nano sliver group, 10 mg/L nano sliver group, 25 mg/L nano sliver group, and 50 mg/L nano sliver group, which were added with the corresponding final mass concentrations of nano sliver solution, respectively. At 48 h of culture, the Fb proliferation viability was detected by cell counting kit 8 method. The Fbs were divided into 0 mg/L silver-containing GelMA hydrogel group, 10 mg/L silver-containing GelMA hydrogel group, 50 mg/L silver-containing GelMA hydrogel group, and 100 mg/L silver-containing GelMA hydrogel group and then were correspondingly treated. On culture day 1, 3, and 7, the Fb proliferation viability was detected as before. The ASCs were mixed into GelMA hydrogel and divided into 3D bioprinting group and non-printing group. On culture day 1, 3, and 7, the ASC proliferation viability was detected as before and cell growth was observed by live/dead cell fluorescence staining. The sample numbers in the above experiments were all 3. Four full-thickness skin defect wounds were produced on the back of 18 male Sprague-Dawley rats aged 4 to 6 weeks. The wounds were divided into hydrogel alone group, hydrogel/nano sliver group, hydrogel scaffold/nano sliver group, and hydrogel scaffold/nano sliver/ASC group, and transplanted with the corresponding scaffolds, respectively. On post injury day (PID) 4, 7, 14, and 21, the wound healing was observed and the wound healing rate was calculated (n=6). On PID 7 and 14, histopathological changes of wounds were observed by hematoxylin eosin staining (n=6). On PID 21, collagen deposition of wounds was observed by Masson staining (n=3). Data were statistically analyzed with one-way analysis of variance, analysis of variance for repeated measurement, Bonferroni correction, and independent sample t test. Results: The sliver nano particles in nano silver solution with different mass concentrations were all round, in scattered distribution and uniform in size. The silver-containing GelMA hydrogels with different final mass fractions of GelMA all showed pore structures of different sizes and interconnections. The pore size of silver-containing GelMA hydrogel with 10% final mass fraction was significantly larger than that of silver-containing GelMA hydrogels with 15% and 20% final mass fractions (with P values both below 0.05). On treatment day 1, 3, and 7, the concentration of nano silver released from silver-containing GelMA hydrogel in vitro showed a relatively flat trend. On treatment day 14, the concentration of released nano silver in vitro increased rapidly. At 24 h of culture, the diameters of inhibition zone of GelMA hydrogel containing 0, 25, 50, and 100 mg/L nano silver against Staphylococcus aureus and Escherichia coli were 0, 0, 0.7, and 2.1 mm and 0, 1.4, 3.2, and 3.3 mm, respectively. At 48 h of culture, the proliferation activity of Fbs in 2 mg/L nano silver group and 5 mg/L nano silver group was both significantly higher than that in blank control group (P<0.05), and the proliferation activity of Fbs in 10 mg/L nano silver group, 25 mg/L nano silver group, and 50 mg/L nano silver group was all significantly lower than that in blank control group (P<0.05). Compared with the that of Fbs in 0 mg/L silver-containing GelMA hydrogel group, the proliferation activity of Fbs in 50 mg/L silver-containing GelMA hydrogel group and 100 mg/L silver-containing GelMA hydrogel group was all significantly decreased on culture day 1 (P<0.05); the proliferation activity of Fbs in 50 mg/L silver-containing GelMA hydrogel group was significantly increased (P<0.05), while the proliferation activity of Fbs in 100 mg/L silver-containing GelMA hydrogel group was significantly decreased on culture day 3 (P<0.05); the proliferation activity of Fbs in 100 mg/L silver-containing GelMA hydrogel group was significantly decreased on culture day 7 (P<0.05). The proliferation activity of ASCs in 3D bioprinting group show no statistically significant differences to that in non-printing group on culture day 1 (P>0.05). The proliferation activity of ASCs in 3D bioprinting group was significantly higher than that in non-printing group on culture day 3 and 7 (with t values of 21.50 and 12.95, respectively, P<0.05). On culture day 1, the number of dead ASCs in 3D bioprinting group was slightly more than that in non-printing group. On culture day 3 and 5, the majority of ASCs in 3D bioprinting group and non-printing group were living cells. On PID 4, the wounds of rats in hydrogel alone group and hydrogel/nano sliver group had more exudation, and the wounds of rats in hydrogel scaffold/nano sliver group and hydrogel scaffold/nano sliver/ASC group were dry without obvious signs of infection. On PID 7, there was still a small amount of exudation on the wounds of rats in hydrogel alone group and hydrogel/nano sliver group, while the wounds of rats in hydrogel scaffold/nano sliver group and hydrogel scaffold/nano sliver/ASC group were dry and scabbed. On PID 14, the hydrogels on the wound surface of rats in the four groups all fell off. On PID 21, a small area of wounds remained unhealed in hydrogel alone group. On PID 4 and 7, the wound healing rates of rats in hydrogel scaffold/nano sliver/ASC group were significantly higher than those of the other three groups (P<0.05). On PID 14, the wound healing rate of rats in hydrogel scaffold/nano sliver/ASC group was significantly higher than the wound healing rates in hydrogel alone group and hydrogel/nano sliver group (all P<0.05). On PID 21, the wound healing rate of rats in hydrogel alone group was significantly lower than that in hydrogel scaffold/nano sliver/ASC group (P<0.05). On PID 7, the hydrogels on the wound surface of rats in the four groups remained in place; on PID 14, the hydrogel in hydrogel alone group was separated from the wounds of rats, while some hydrogels still existed in the new tissue of the wounds of rats in the other three groups. On PID 21, the collagen arrangement in the wounds of rats in hydrogel alone group was out of order, while the collagen arrangement in the wounds of rats in hydrogel/nano sliver group, and hydrogel scaffold/nano sliver/ASC group was relatively orderly. Conclusions: Silver-containing GelMA hydrogel has good biocompatibility and antibacterial properties. Its three-dimensional bioprinted double-layer structure can better integrate with new formed tissue in the full-thickness skin defect wounds in rats and promote wound healing.

[Yin Chan Quan Shu, the Obstetrics and Gynecology Monograph by Wang Kentang]

Yin Chan Quan Shu (Obstetrics and gynecology monograph) is a monograph on obstetrics and gynecology compiled by Wang Kentang in the Ming Dynasty. It had four volumes and was published in the thirtieth year of Wanli (1602) in the Ming Dynasty after it was edited by Zhang Shoukong and others. It was found that Yin Chan Quan Shu has four versions remaining. They were the version printed by Shu Lin Qiao Shan Tang in the Ming Dynasty, held in the National Library of China and the Cabinet Library of Japanese Official Documents Library; the version revised according to the version of Shu Lin Qiao Shan Tang, held in the Library of Capital Medical University, Tianjin Medical College, Shanghai Branch of the Chinese Medical Association, the Library of Guangzhou University of Chinese Medicine and the Cabinet Library of the National Archives of Japan; the version based on the version of Shu Lin Qiao Shan Tang in the Ming Dynasty, transcribed in the fourth year of Wen Hua (1807), collected in the Cabinet Library of the National Archives of Japan; the version transcribed according to the revised version in the Ming Dynasty, collected in the Shanghai Branch of the Chinese Medical Association. It was found that there was no evidence to support the existence of the so-called "version of Kangxi in the Qing Dynasty". This means almost all versions remaining came from the versions published in the Ming Dynasty. The references of Yin Chan Quan Shu came from Pulse Classic (Mai Jing), Chan Bao, Fu Ren Da Quan Liang Fang and other works with the supplement and development by Wang Kentang.Yin Chan Quan Shu was the main sources and foundation of the Criteria of Syndrome Identification and Treatment in Gynecology (Nv Ke Zheng Zhi Zhun Sheng) by Wang Kentang.

Constraints on Heavy Decaying Dark Matter from 570 Days of LHAASO Observations

The kilometer square array (KM2A) of the large high altitude air shower observatory (LHAASO) aims at surveying the northern γ-ray sky at energies above 10 TeV with unprecedented sensitivity. γ-ray observations have long been one of the most powerful tools for dark matter searches, as, e.g., high-energy γ rays could be produced by the decays of heavy dark matter particles. In this Letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340 days of data from 1/2-KM2A and 230 days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after γ/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest γ-ray constraints on the lifetime of heavy dark matter particles with mass between 10^{5} and 10^{9}  GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

Seasonal Influence of Biodiversity on Soil Respiration in a Temperate Forest

Soil respiration in forests contributes to significant carbon dioxide emissions from terrestrial ecosystems but it varies both spatially and seasonally. Both abiotic and biotic factors influence soil respiration but their relative contribution to spatial and seasonal variability remains poorly understood, which leads to uncertainty in models of global C cycling and predictions of future climate change. Here, we hypothesize that tree diversity, soil diversity, and soil properties contribute to local-scale variability of soil respiration but their relative importance changes in different seasons. To test our hypothesis, we conducted seasonal soil respiration measurements along a local-scale environmental gradient in a temperate forest in Northeast China, analyzed spatial variability of soil respiration and tested the relationships between soil respiration and a variety of abiotic and biotic factors including topography, soil chemical properties, and plant and soil diversity. We found that soil respiration varied substantially across the study site, with spatial coefficients of variation (CV) of 29.1%, 27.3% and 30.8% in spring, summer, and autumn, respectively. Soil respiration was consistently lower at high soil water content, but the influence of other factors was seasonal. In spring, soil respiration increased with tree diversity and biomass but decreased with soil fungal diversity. In summer, soil respiration increased with soil temperature, whereas in autumn, soil respiration increased with tree diversity but decreased with increasing soil nutrient content. However, soil nutrient content indirectly enhanced soil respiration via its effect on tree diversity across seasons, and forest stand structure indirectly enhanced soil respiration via tree diversity in spring. Our results highlight that substantial differences in soil respiration at local scales was jointly explained by soil properties (soil water content and soil nutrients), tree diversity, and soil fungal diversity but the relative importance of these drivers varied seasonally in our temperate forest.

El Niño-Southern Oscillation affects the species-level temporal variation in seed and leaf fall in a mixed temperate forest

El Niño-Southern Oscillation (ENSO), the variation between anomalously cold (La Niña) and warm conditions (El Niño), is one of the most prominent large-scale climate patterns with worldwide effects. Elevated seed and leaf fall has been found at the positive phase of ENSO (El Niño) in tropical forests. However, how seed and leaf fall respond to ENSO at species level is understudied, especially in temperate forests. In this study, we monitored seed and leaf fall at the species-level at 150 points across a 25-ha temperate forest in northeastern China over a span of 12 years. Using time series and wavelet analyses, we assessed three hypotheses: 1) temperate tree species' seed and leaf fall are strongly, but differently, correlated with ENSO and, 2) community synchrony in seed and leaf occurred both at seasonal and ENSO scales; finally, 3) local climatic modulated the effects of ENSO on seed and leaf fall. We found that ENSO was significantly correlated with seed and leaf fall of all species, although correlation strength varied across species (r = 0.206-0.658). Specifically, ENSO indices (ENSO12 or ENSO34) accounted for the most variation in seed and leaf fall of Acer pseudo-sieboldianum (40 % and 34 %, respectively) and ranged 4 %-31 % in all other species. Leaf fall was synchronous with ENSO cycles with a period of 2-7 years, but community synchrony of seed fall was only detected at seasonal scales. ENSO influenced seed fall of Fraxinus mandshurica and Tilla amurensis by mediating rainfall and relative humidity, respectively, highlighting the interactive effects of local climate and ENSO. Our findings highlight the potential effects of ENSO on ecosystems outside of tropical regions and improve our ability to predict regeneration dynamics and nutrient cycling of temperate forests under the context of global change.

[Imaging and cerebrospinal fluid features of two cases with lymphomatosis cerebri]

This article reported the clinical experience of diagnosis and treatment for two patients with lymphomatosis cerebri. Case 1 was female and aged 53 years old, while case 2 was male and aged 69 years old. Progressive cognitive impairment was the main clinical manifestation in both patients. Brain magnetic resonance imaging (MRI) suggested leukoencephalopathy with patchy or mass enhancement. Cerebral blood flow was reduced on perfusion imaging in one patient. Brain biopsy confirmed diffuse large B-cell lymphoma in both cases. The concentration of interleukin-10 in cerebrospinal fluid (CSF) of two patients was significantly increased, however, the result of CSF flow cytology was negative. The current study suggests that interleukin-10 in CSF is an important biological indicator for the diagnosis of lymphomatosis cerebri, but CSF flow cytometry may not be helpful. Moreover, cerebral hypoperfusion can be present in patients with lymphomatosis cerebri.

[Mid-term efficacy of China Net Childhood Lymphoma-mature B-cell lymphoma 2017 regimen in the treatment of pediatric Burkitt lymphoma]

Objective: To analyze the clinical characteristics of children with Burkitt lymphoma (BL) and to summarize the mid-term efficacy of China Net Childhood Lymphoma-mature B-cell lymphoma 2017 (CNCL-B-NHL-2017) regimen. Methods: Clinical features of 436 BL patients who were ≤18 years old and treated with the CNCL-B-NHL-2017 regimen from May 2017 to April 2021 were analyzed retrospectively. Clinical characteristics of patients at disease onset were analyzed and the therapeutic effects of patients with different clinical stages and risk groups were compared. Survival analysis was performed by Kaplan-Meier method, and Cox regression was used to identify the prognostic factors. Results: Among 436 patients, there were 368 (84.4%) males and 68 (15.6%) females, the age of disease onset was 6.0 (4.0, 9.0) years old. According to the St. Jude staging system, there were 4 patients (0.9%) with stage Ⅰ, 30 patients (6.9%) with stage Ⅱ, 217 patients (49.8%) with stage Ⅲ, and 185 patients (42.4%) with stage Ⅳ. All patients were stratified into following risk groups: group A (n=1, 0.2%), group B1 (n=46, 10.6%), group B2 (n=19, 4.4%), group C1 (n=285, 65.4%), group C2 (n=85, 19.5%). Sixty-three patients (14.4%) were treated with chemotherapy only and 373 patients (85.6%) were treated with chemotherapy combined with rituximab. Twenty-one patients (4.8%) suffered from progressive disease, 3 patients (0.7%) relapsed, and 13 patients (3.0%) died of treatment-related complications. The follow-up time of all patients was 24.0 (13.0, 35.0) months, the 2-year event free survival (EFS) rate of all patients was (90.9±1.4) %. The 2-year EFS rates of group A, B1, B2, C1 and C2 were 100.0%, 100.0%, (94.7±5.1) %, (90.7±1.7) % and (85.9±4.0) %, respectively. The 2-year EFS rates was higher in group A, B1, and B2 than those in group C1 (χ²=4.16, P=0.041) and group C2 (χ²=7.21, P=0.007). The 2-year EFS rates of the patients treated with chemotherapy alone and those treated with chemotherapy combined with rituximab were (79.3±5.1)% and (92.9±1.4)% (χ²=14.23, P<0.001) respectively. Multivariate analysis showed that stage Ⅳ (including leukemia stage), serum lactate dehydrogenase (LDH)>4-fold normal value, and with residual tumor in the mid-term evaluation were risk factors for poor prognosis (HR=1.38,1.23,8.52,95%CI 1.05-1.82,1.05-1.43,3.96-18.30). Conclusions: The CNCL-B-NHL-2017 regimen show significant effect in the treatment of pediatric BL. The combination of rituximab improve the efficacy further.

[Research advances on the treatment of hydrofluoric acid burns]

Hydrofluoric acid is a highly dangerous and toxic inorganic acid, which is widely used in industrial fields and daily life. The risk of hydrofluoric acid burns is related to hydrofluoric acid mass fraction, duration of exposure to hydrofluoric acid, burn area, burn depth, and burn site, etc. Hydrofluoric acid has strong toxicity and tissue penetration ability. A small area of hydrofluoric acid burns can cause death in a short time. Therefore, improving the understanding of the mechanism of hydrofluoric acid burns and learning how to treat hydrofluoric acid burns in different sites can further improve the cure rate of hydrofluoric acid burns.

[Research advances on biomaterials for the delivery of growth factors to regulate wound repair]

Wound repair is a highly coordinated and mutually regulated complex process involving various kinds of cells, extracellular matrices and cytokines. A variety of growth factors play an important regulatory role in wound healing, and it is critical to achieve effective delivery and sustained function of growth factors. In recent years, the application of biomaterials in tissue engineering has shown great potential, and the effective delivery of growth factors by biomaterials has attracted increasing attention. Based on this, this paper introduces the mechanism of related growth factors in the process of wound healing, focusing on the recent progress of biomaterial delivery of growth factors to accelerate wound healing, in order to provide new enlightenment for clinical wound treatment.

Sensitivity of Parity-Violating Electron Scattering to a Dark Photon

We explore the sensitivity of the parity-violating electron scattering (PVES) asymmetry in both elastic and deep-inelastic scattering to the properties of a dark photon. Given advances in experimental capabilities in recent years, there are interesting regions of parameter space where PVES offers the chance to discover new physics in the near future. There are also cases where the existence of a dark photon could significantly alter our understanding of the structure of atomic nuclei and neutron stars as well as parton distribution functions.

[A case of chronic refractory wound in the abdomen caused by residual foreign body]

A 59-year-old male patient with local sinus tract formation due to residual foreign body was admitted to the Second Affiliated Hospital of Zhejiang University College of Medicine on December 17, 2018. The examination showed that the residual foreign body was the component of a sticky cloth implanted when the patient underwent appendectomy 27 years ago. Hypertrophic scar developed at the right-lower abdominal incision for appendectomy 23 years ago and the secondary infection after cicatrectomy resulted in non-healing of the wound. The chronic refractory wound healed completely after surgical treatment in our hospital after this admission. The postoperative pathological examination revealed local inflammatory granuloma. This case suggests that chronic refractory wound is likely to form when secondary infection occurs following the surgical procedure near the implant, and aggressive surgery is an effective way to solve this problem.

[Epidemiological characteristics and treatment outcome analysis of 229 patients with hydrofluoric acid burns in hands]

Objective: To explore the epidemiological characteristics and treatment outcomes of patients with hydrofluoric acid burns in hands. Methods: A retrospective observational study was conducted. The medical records of 229 patients with hydrofluoric acid burns in hands who were admitted to Zhejiang Quhua Hospital from January 2008 to December 2020 and met the inclusion criteria were collected. The following statistical data of patients were collected, including gender, age, type of affiliated enterprise, hydrofluoric acid mass fraction, injury site, total burn area, prehospital time, length of hospital stay, length of wound healing, whether hypocalcemia and hypomagnesemia occurred or not on admission, whether surgery intervention was performed or not, and whether scar sequelae occurred or not. Single factor and multivariate logistic regression analysis were used to screen out the risk factors impacting surgery intervention and scar sequelae of all the patients and patients whose hydrofluoric acid mass fraction was known. Single factor and multivariate linear regression analysis were used to screen out the risk factors impacting the length of wound healing of all the patients and patients whose hydrofluoric acid mass fraction was known. Results: The 229 patients included 206 males and 23 females, with the majority aged 30 to 50 years (139 patients). The type of affiliated enterprise of majority patients was non-fluorine chemical enterprise. The hydrofluoric acid mass fraction was known in only 91 patients, mainly medium. The majority injury site was in the middle and end of finger. The total burn area was below or equal to 1% total body surface area. The prehospital time was 19 (9, 29) h. The length of hospital stay was 2 (1, 7) d. The length of wound healing was 12 (8, 18) d. The proportions of hypocalcemia and hypomagnesemia were 0.9% (2/229) and 1.3% (3/229) on admission, respectively. Thirty-six patients had surgeries and 83 patients had scar sequelae. In 229 patients, single factor logistic regression analysis showed that both type of affiliated enterprise and prehospital time were the factors impacting surgery intervention (with odds ratio values of 7.86 and 51.35, respectively, 95% confidence intervals of 1.83-33.76 and 11.89-221.78, respectively, P<0.01) and scar sequelae of patients (with odds ratio values of 3.62 and 27.40, respectively, 95% confidence intervals of 1.76-7.43 and 13.25-56.68, respectively, P<0.01); multivariate logistic regression analysis showed that prehospital time was the independent risks factor impacting surgery intervention and scar sequelae of patients (with odds ratio values of 43.00 and 24.55, respectively, 95% confidence intervals of 9.89-187.03 and 11.78-51.16, respectively, P<0.01); single factor linear regression analysis showed that both type of affiliated enterprise and prehospital time were the factors impacting the length of wound healing of patients (with β values of 6.16 and 12.83, respectively, 95% confidence intervals of 3.38-8.93 and 10.72-14.93, respectively, P<0.01); multivariate linear regression analysis showed that both type of affiliated enterprise and prehospital time were the independent risk factors impacting the length of wound healing of patients (with β values of 2.81 and 12.16, respectively, 95% confidence intervals of 0.50-5.13 and 10.00-14.31, respectively, P<0.05 or P<0.01). In 91 patients whose hydrofluoric acid mass fraction was known, single factor logistic regression analysis showed that type of affiliated enterprise, hydrofluoric acid mass fraction (low and high), and prehospital time were all the factors impacting surgery intervention of patients (with odds ratio values of 9.10, 11.25, 10.69, and 0.04, respectively, 95% confidence intervals of 1.15-72.25, 1.39-90.93, 1.32-86.59, and 0.01-0.19, respectively, P<0.05 or P<0.01), type of affiliated enterprise, hydrofluoric acid mass fraction, and prehospital time were all the factors impacting scar sequelae of patients (with odds ratio values of 0.32, 0.21, and 36.80, respectively, 95% confidence intervals of 0.11-0.92, 0.06-0.73, and 11.03-122.79, respectively, P<0.05 or P<0.01); multivariate logistic regression analysis showed that both hydrofluoric acid mass fraction and prehospital time were the independent risk factors impacting surgery intervention of patients (with odds ratio values of 11.51 and 0.04, respectively, 95% confidence intervals of 1.22-108.26 and 0.01-0.25, respectively, P<0.05 or P<0.01), prehospital time was the independent risk factor impacting scar sequelae of patients (odds ratio=37.71, with 95% confidence interval of 9.97-142.69, P<0.01); single factor linear regression analysis showed that type of affiliated enterprise, hydrofluoric acid mass fraction (low and high), and prehospital time were all the factors impacting the length of wound healing of patients (with β values of 7.12, -5.63, -9.74, and 13.50, respectively, 95% confidence intervals of 2.43-11.81, -10.59--0.68, -14.78--4.70, and 10.14-16.86, respectively, P<0.05 or P<0.01); multivariate linear regression analysis showed that both hydrofluoric acid mass fraction and prehospital time were the independent risk factors impacting the length of wound healing of patients (with β values of -5.84 and 0.09, respectively, 95% confidence intervals of -10.59--1.08 and 0.05-0.12, respectively, P<0.05 or P<0.01). Conclusions: The majority of patients with hydrofluoric acid burns in hands are young and middle-aged males. Type of affiliated enterprise, hydrofluoric acid mass fraction and prehospital time are the risk factors that affect the treatment outcomes of patients with hydrofluoric acid burns in hands.

Exploring Lorentz Invariance Violation from Ultrahigh-Energy γ Rays Observed by LHAASO

Recently, the LHAASO Collaboration published the detection of 12 ultrahigh-energy γ-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV γ rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy γ rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of the LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^{5} times the Planck scale M_{Pl} and that the second-order LIV scale is >10^{-3}M_{Pl}. Both limits improve by at least one order of magnitude the previous results.

[Comparison of alignment and position of prosthesis between portable accelerometer-based navigation device and conventional instrumentation in total knee arthroplasty with valgus deformity]

Objective: To explore the differences of alignment and position of prosthesis between portable accelerometer-based navigation device (PAD) and conventional instrumentation (CI) in total knee arthroplasty (TKA) with valgus deformity. Methods: Patients with knee osteoarthritis and valgus deformity who underwent primary TKA in Peking University Third Hospital from January 2017 to December 2020 were enrolled in this retrospective study and were divided into PAD group and CI group according to the surgical instruments. Five male patients and 44 female patients were included with a mean age of (67.2±7.0) years. The differences in preoperative general data, preoperative and postoperative alignment between the two groups were studied. Results: A total of 49 patients (25 patients in the PAD group and 24 in the CI group) were enrolled in this study. There were no statistically significant differences in gender, age, height, weight, body mass index, surgical side, preoperative hip-knee-ankle (HKA) angle, preoperative HKA angle deviation, Keblish classification and Ranawat classification between the two groups (all P>0.05). There was no significant difference in the accuracy of postoperative HKA angle (2.0°±1.4° vs 3.0°±2.2°, P=0.082), coronal femoral component angle (CFCA) (1.5°±1.2° vs 2.1°±1.6°, P=0.144) and coronal tibial component angle (CTCA) (1.2°±0.8° vs 1.3°±1.0°, P=0.695) between the two groups; but the standard deviation of the above-mentioned three indices in PAD group were all smaller than those in CI group. The rate of outliers of postoperative HKA angle of the PAD group was smaller than that in the CI group (P<0.05), but there was no significant difference in the rate of outliers of CFCA and CTCA between the two groups (both P>0.05). Conclusion: TKA assisted by PAD can provide good alignment and prosthesis position in patients with valgus deformity, and it is superior to TKA with CI in terms of precision and rate of outliers of postoperative overall alignment of lower extremity.

[Effect of the acetabular cup positioning and leg length restoration after total hip arthroplasty using robotic-assisted surgery system]

Objective: To assess the improvement of the robotic-assisted total hip arthroplasty (rTHA) regarding the acetabular cup positioning and the leg length restoration. Methods: Clinical data of 246 patients undergoing primary unilateral THA from September 2019 to June 2021 in Peking University Third Hospital were retrospectively reviewed, including patients treated with rTHA (n=113) and conventional THA (cTHA) (n=133). Thirty-seven male patients along with 76 females were enrolled into the rTHA group with a mean age of (62±9) years. In comparison, the average age of cTHA group was (60±12) years with 59 males and 74 females. The postoperative inclination and anteversion of the acetabular cup and the rate of acetabular cup within the Lewinnek and Callanan "safe zone" were documented and analyzed. For the patient with normal or surgically restored contralateral hip, the leg length discrepancy was also measured and compared between the two groups. Results: The overall mean postoperative inclination and anteversion was 39.0°±5.5° and 14.1°±5.1°, respectively. No statistically significant difference was found between the rTHA and cTHA group regarding postoperative inclination (39.3°±2.9° vs 38.7°±7.0°, P=0.383) and anteversion (13.4°±4.3° vs 14.7°±5.6°, P=0.054). In rTHA group, 97.3%(110/113) of the cups were implanted within the Lewinnek"safe zone"(75.9%(101/133) in cTHA group) and 94.7%(107/113) were within the Callanan"safe zone"(66.2%(88/133) in cTHA group), and those were both higher in rTHA group (both P<0.01). There was no significant difference in postoperative leg length discrepancy between the two groups (P=0.445). Meanwhile, 29(23.0%) cases of cTHA group had leg length discrepancy more than 5 mm, and it was 9.7%(9/93) in rTHA group. Conclusion: It indicated the benefit of rTHA in acetabular cup positioning as well as the leg length restoration.

[Comparison of alignment between computer assisted surgery and personal specific instrumentation in total knee arthroplasty]

Objective: To explore the differences of alignment and operative time between computer assisted surgery (CAS) and personal specific instrumentation (PSI) in total knee arthroplasty (TKA). Methods: Patients with knee osteoarthritis who underwent unilateral primary TKA in Peking University Third Hospital from June 2019 to April 2021 were enrolled in this prospective study and were randomly divided into CAS group and PSI group with random number table. The differences in preoperative general data, operative time, preoperative and postoperative alignment between the two groups were studied. Results: A total of 80 patients were enrolled in this study, including 40 patients in the CAS group with an average age of (68±6) years (male/female, 6/34 cases) and 40 patients in the PSI group with an average age of (69±6) years (male/female, 5/35 cases). There was no significant differences in gender, age, height, weight, body mass index, surgical side, preoperative hip-knee-ankle (HKA) angle, and preoperative HKA angle deviation between the two groups (all P>0.05). The accuracy of sagittal femoral component angle (sFCA) (1.7°±1.2° vs 2.5°±1.4°, P<0.05) and sagittal tibial component angle (sTCA) (1.2°±0.9° vs 2.3°±1.8°, P<0.01) of the CAS group were both smaller than those in the PSI group, but there was no significant difference in the accuracy of postoperative HKA angle, mechanical lateral distal femoral angle (mLDFA) and mechanical medial proximal tibia angle (mMPTA) between the two groups (all P>0.05). The precision of postoperative HKA angle and sTCA of the CAS group were both smaller than those in the PSI group (2.0°±1.0° vs 2.6°±1.7°, 1.2°±0.9° vs 2.3°±1.8°, both P<0.01), but there was no significant differences in the precision of mLDFA, mMPTA and sFCA between the two groups (all P>0.05). The rate of outliers of postoperative HKA angle (10.0% vs 27.5%, P<0.05), sFCA(5.0% vs 20.0%, P<0.05) and sTCA (2.5% vs 22.5%, P<0.01) of the CAS group were all smaller than those in the PSI group, but there was no significant difference in the rate of outliers of mLDFA and mMPTA between the two groups (both P>0.05). CAS group had significantly longer surgical time than the PSI group ((81±12) min vs (52±8) min, P<0.01). Conclusions: CAS is better than PSI in overall coronal alignment and sagittal position of the prosthesis after TKA. However, CAS requires a longer operation time.

[The role of first-aid network construction in the early treatment of patients with critically severe hydrofluoric acid burns]

Objective: To explore the role of first-aid network construction in the early treatment of patients with critically severe hydrofluoric acid burns. Methods: Twenty-seven fluorine chemical enterprises distributed in Zhejiang province, Jiangxi Province, Fujian Province, and Inner Mongolia Autonomous Region and 22 hospitals with burn/plastic department or professional burn treatment group in Zhejiang province, including Zhejiang Quhua Hospital, and 5 hospitals outside Zhejiang province were involved in the first-aid network construction as member units. As the main unit, Zhejiang Quhua Hospital was responsible for the daily maintenance and technical guidance of the first-aid network. Zhejiang Quhua Hospital was assigned as the designated emergency hospital for 20 fluorine chemical enterprises, a near emergency hospital to the other 7 fluorine chemical enterprises was assigned as the designated hospital for them. Medical records of 56 patients (all males) with critically severe hydrofluoric acid burns who admitted to 5 first-aid network hospitals from January 2006 to June 2021, meeting the inclusion criteria, were involved in the retrospective cohort study. Based on whether the enterprise belonging to the first-aid network construction or not, the patients were divided into first-aid network group (27 cases, aged (41±9) years) and non first-aid network group (29 cases, aged (42±10) years). After the patients in the first-aid network group were injured, the enterprises and hospitals linked up immediately. The hospital where the patient was treated mobilize the treatment force, equipment, materials, and drugs in advance by the first-aid network, thereby realizing seamless joint between pre-hospital first-aid and in-hospital treatment. The hospital started the first-aid process and temporarily mobilized the rescue forces, equipment, materials, and drug after patients in non first-aid network group arrived at the department of emergency of the hospital. The time from injury to medical service, the first detection time of serum calcium, the time staying in department of emergency, the duration of hypocalcemia and hypomagnesemia, and the treatment outcome of patients in the two groups were recorded. Data were statistically analyzed with chi-square test, Fisher's exact probability test, independent-sample t test, and Wilcoxon rank-sum test. Results: The time from injury to medical service, the first detection time of serum calcium, and the time staying in department of emergency of patients in first-aid network group were 40.0 (30.0, 55.0), 23.0 (17.5, 37.5), and 42.0 (37.0, 53.0) min, which were significantly shorter than 180.0 (120.0, 240.0), 31.0 (22.5, 47.5), 61.0 (52.0, 65.5) min in non first-aid network group (Z=-6.17, -1.98, -4.15, P<0.05 or P<0.01). The duration of hypocalcemia and hypo- magnesemia of patients in first-aid network group were 1.2 (1.1, 1.6) and 1.9 (1.7, 2.1) h, which were significantly shorter than 4.6 (3.1, 6.2) and 3.2 (2.5, 4.6) h in non first-aid network group (Z=-5.80, -4.81, P<0.01). Three patients (11.1%) in first-aid network group died, among whom 2 patients died at 40 min after injury and 1 patient died 9.0 h after injury. Four patients (13.8%) died in non first-aid network group at 3.0, 3.0, 4.5, and 7.0 h after injury, respectively. The mortality rates of patients in the two groups were similar (P>0.05). Conclusions: Critically severe hydrofluoric acid burn is an extremely urgent situation encountered in clinical practice. The construction of a first-aid network creates condition for on-site treatment of patients and improves the first-aid efficiency, thereby gaining time to save lives.

[Analysing pathogenic bacterial spectrum and drug resistance of bloodstream infection in patients with allogeneic hematopoietic stem cell transplantation]

Objective: To elucidate the clinical characteristics of bloodstream infection in patients with allogeneic hematopoietic stem cell transplantation (allo-HSCT) in our hospital and improves the survival of transplant patients with bloodstream infection. Methods: Two hundred and ten patients with allo-HSCT from the Department of Hematology were retrospectively analyzed between October 2014 and September 2019. Pathogen distribution, drug resistance, risk factors, and outcomes were investigated in 49 allo-HSCT patients with bloodstream infections. Results: Forty-nine of 210 patients with allo-HSCT had bloodstream infection, and 59 pathogenic microorganisms were identified, mainly Gram-negative bacteria (67.8%) , of which E. coli had the highest incidence (23.7%) , CRO accounted for 42.5%, and Grampositive bacteria accounted for 23.7% (without vancomycin or linezolid-resistant strain) . Additionally, fungi accounted for 8.5%. Univariate analysis suggested that the risk factors of bloodstream infection were gender, pretransplant disease status, and conditioning regimen. In contrast, multivariate analysis showed that bloodstream infection was mainly related to conditioning regimens. Further grouping results showed that 77.6% of patients with neutropenia had bloodstream infections, and 22.4% of patients with non-neutropenia had bloodstream infections; 81.0% of patients with active infections before transplantation had bloodstream infections, while bloodstream infection occurred in 16.9% of patients without active infection. Survival analysis showed that long-term survival of patients with bloodstream infection is shorter than that of patients without bloodstream infection and long-term survival of patients with CRO infection is shorter than that of patients without CRO infection. The survival of patients with neutropenia longer than 14 d is shorter than that of patients with neutropenia shorter than 14 d. Furthermore, there is no correlation between whether there is an active infection before transplantation and whether they are in a neutropenic state at the time of infection and survival. Conclusion: Our results suggest that effective prevention of bloodstream infections from drug-resistant bacteria, particularly CRO, shortening the duration of neutropenia, eradication of potential infections before transplantation, and patient-adaptive conditioning could reduce transplant-related mortality and improve prognosis.

[Quotations from The List of Medical Doctors (Ming Yi Lu) in The Medical Theory (Yi Shuo) and The Mengqiu of Famous Medical Doctors in the Previous Dynasties (Li Dai Ming Yi Meng Qiu)]

The Medical Theory (Yi Shuo) and The Mengqiu of Famous Medical Doctors in the Past Dynasties (Li Dai Ming Yi Meng Qiu) are two medical books that originated in the Song Dynasty. The two books quoted from The List of Medical Doctors (Ming Yi Lu) many times while the author and the issue time of The List of Medical Doctors (Ming Yi Lu) has not been clearly known up to know. It was found that there were 25 quotations in The Medical Theory (Yi Shuo) and 31 quotations in The Mengqiu of Famous Medical Doctors in the Past Dynasties (Li Dai Ming Yi Meng Qiu) from The List of Medical Doctors (Ming Yi Lu). The contents quoted in the two books are also found in The Mysterious Medical List (Shen Mi Ming Yi Lu) written by Dang Yongnian in the Song Dynasty. In this sense, it can be inferred that The List of Medical Doctors (Ming Yi Lu) mentioned above should be The Mysterious Medical List (Shen Mi Ming Yi Lu).

Tree growth response to soil nutrients and neighborhood crowding varies between mycorrhizal types in an old-growth temperate forest

Forest dynamics are shaped by both abiotic and biotic factors. Trees associating with different types of mycorrhizal fungi differ in nutrient use and dominate in contrasting environments, but it remains unclear whether they exhibit differential growth responses to local abiotic and biotic gradients where they co-occur. We used 9-year tree census data in a 25-ha old-growth temperate forest in Northeast China to examine differences in tree growth response to soil nutrients and neighborhood crowding between tree species associating with arbuscular mycorrhizal (AM), ectomycorrhizal (EM), and dual-mycorrhizal (AEM) fungi. In addition, we tested the role of individual-level vs species-level leaf traits in capturing differences in tree growth response to soil nutrients and neighborhood crowding across mycorrhizal types. Across 25 species, soil nutrients decreased AM tree growth, while neighborhood crowding reduced both AM and EM tree growth, and neither soil nor neighbors impacted AEM tree growth. Across mycorrhizal types, individual-level traits were stronger predictors of tree growth than species-level traits. However, most traits poorly mediated tree growth response to soil nutrients and neighborhood crowding. Our findings indicate that mycorrhizal types strongly shape differences in tree growth response to local soil and crowding gradients, and suggest that including plant-mycorrhizae associations in future work offers great potential to improve our understanding of forest dynamics.

[A retrospective study of hemophilic pseudotumor in maxillofacial region]

Objective: To provide clinical references for the diagnosis and treatment of hemophilic pseudotumor (HPT) in maxillofacial region. Methods: Fourteen cases of HPT in maxillofacial region from the Department of Stomatology, Beijing Children's Hospital from Jan 2009 to Jan 2019 were collected. Two cases were lost for follow-up and 12 patient,all boys, were finally followed up and included in the study. The patients aged from 13 months to 10 years old. The medical history, clinic manefestitions and the features of the radiology examination were recorded. The patients were treated by using replacement treatment first. If the conservative treatment was not effective, the patients then received operation combined with pereoperation replacement thearapy. The patients were followed up for 13 months to 10 years.There were 11 cases of hemophilia A, and 1 case of hemophilia B. Two cases were severe type, the others (10/12) were mild and moderate types. Only 1 case was diagnosed as hemophilia initially. Nine cases (9/12) were misdiagnosed as malignant tumors, 1 case was misdiagnosed as osteomyelitis and 1 case was misdiagnosed as hemangioma. Only 3 cases had identified history of trauma before. Results: All cases were treated with replacement therapy first, among which 10 cases were effective, 8 cases were cured by conservative therapy, 1 case had residual soft tissue fistula after conservative treatment and 1 case recurrented after conservative treatment for 8 months. Two patients with poor efficacy to the replacement treatment were performed operations and finally were cured. Conclusions: The misdiagnosis rate of HPT in maxillofacial region was high. The conservative factor replacement therapy could achieve good results in most children and could be used as the preferred treatment. If the conservative treatment was not effective, the surgical treatment was also a safe option.

Peta-electron volt gamma-ray emission from the Crab Nebula

The Crab Nebula is a bright source of gamma rays powered by the Crab Pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind. We report the detection of gamma rays from this source with energies from 5 × 10^-4 to 1.1 peta-electron volts with a spectrum showing gradual steepening over three energy decades. The ultrahigh-energy photons imply the presence of a peta-electron volt electron accelerator (a pevatron) in the nebula, with an acceleration rate exceeding 15% of the theoretical limit. We constrain the pevatron's size between 0.025 and 0.1 parsecs and the magnetic field to ≈110 microgauss. The production rate of peta-electron volt electrons, 2.5 × 10³⁶ ergs per second, constitutes 0.5% of the pulsar spin-down luminosity, although we cannot exclude a contribution of peta-electron volt protons to the production of the highest-energy gamma rays.

Consequences of spatial patterns for coexistence in species-rich plant communities

Ecology cannot yet fully explain why so many tree species coexist in natural communities such as tropical forests. A major difficulty is linking individual-level processes to community dynamics. We propose a combination of tree spatial data, spatial statistics and dynamical theory to reveal the relationship between spatial patterns and population-level interaction coefficients and their consequences for multispecies dynamics and coexistence. Here we show that the emerging population-level interaction coefficients have, for a broad range of circumstances, a simpler structure than their individual-level counterparts, which allows for an analytical treatment of equilibrium and stability conditions. Mechanisms such as animal seed dispersal, which result in clustering of recruits that is decoupled from parent locations, lead to a rare-species advantage and coexistence of otherwise neutral competitors. Linking spatial statistics with theories of community dynamics offers new avenues for explaining species coexistence and calls for rethinking community ecology through a spatial lens.

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622+3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High-Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of 8.2σ for E>25  TeV assuming a Gaussian template. The best-fit location is (right ascension, declination) =(95.47°±0.11°,37.92°±0.09°), and the extension is 0.40°±0.07°. The energy spectrum can be described by a power-law spectrum with an index of -2.92±0.17_{stat}±0.02_{sys}. No clear extended multiwavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of ∼160  TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

Polyphosphate (PolyP) for alveolar cleft repair: study protocol for a pilot randomized controlled trial

Objective

Bone grafting is an important surgical procedure to restore missing bone in patients with alveolar cleft lip/palate, aiming to stabilize either sides of the maxillary segments by inducing new bone formation, and in bilateral cleft cases also to stabilize the pre-maxilla. Polyphosphate (PolyP), a physiological polymer composed of orthophosphate units linked together with high-energy phosphate bonds, is a naturally existing compound in platelets which, when complexed with calcium as Ca-polyP microparticles (Ca-polyP MPs), was proven to have osteoinductive properties in preclinical studies.

Aim

To evaluate the feasibility, safety, and osteoinductivity of Ca-polyP MPs as a bone-inducing graft material in humans.

Methods

This prospective non-blinded first-in-man clinical pilot study shall consist of 8 alveolar cleft patients of 13 years or older to evaluate the feasibility and safety of Ca-PolyP MPs as a bone-inducing graft material. Patients will receive Ca-polyP graft material only or Ca-polyP in combination with biphasic calcium phosphate (BCP) as a bone substitute carrier. During the trial, the participants will be investigated closely for safety parameters using radiographic imaging, regular blood tests, and physical examinations. After 6 months, a hollow drill will be used to prepare the implantation site to obtain a biopsy. The radiographic imaging will be used for clinical evaluation; the biopsy will be processed for histological/histomorphometric evaluation of bone formation.

Discussion

This is the first-in-man study evaluating the safety and feasibility of the polyP as well as the potential regenerative capacity of polyP using an alveolar cleft model.

Trial registration

Indonesian Trial Registry INA-EW74C1N. Registered on 12 June 2020

Temperature sensitivity of SOM decomposition is linked with a K-selected microbial community

Temperature sensitivity (Q₁₀ ) of soil organic matter (SOM) decomposition is a crucial parameter to predict the fate of soil carbon (C) under global warming. Nonetheless, the response pattern of Q₁₀ to continuous warming and the underlying mechanisms are still under debate, especially considering the complex interactions between Q₁₀ , SOM quality, and soil microorganisms. We examined the Q₁₀ of SOM decomposition across a mean annual temperature (MAT) gradient from -1.9 to 5.1°C in temperate mixed forest ecosystems in parallel with SOM quality and bioavailability, microbial taxonomic composition, and functional genes responsible for organic carbon decomposition. Within this temperature gradient of 7.0°C, the Q₁₀ values increased with MAT, but decreased with SOM bioavailability. The Q₁₀ values increased with the prevalence of K-strategy of soil microbial community, which was characterized by: (i) high ratios of oligotrophic to copiotrophic taxa, (ii) ectomycorrhizal to saprotrophic fungi, (iii) functional genes responsible for degradation of recalcitrant to that of labile C, and (iv) low average 16S rRNA operon copy number. Because the recalcitrant organic matter was mainly utilized by the K-strategists, these findings independently support the carbon quality-temperature theory from the perspective of microbial taxonomic composition and functions. A year-long incubation experiment was performed to determine the response of labile and recalcitrant C pools to warming based on the two-pool model. The decomposition of recalcitrant SOM was more sensitive to increased temperature in southern warm regions, which might attribute to the dominance of K-selected microbial communities. It implies that climate warming would mobilize the larger recalcitrant pools in warm regions, exacerbating the positive feedback between increased MAT and CO₂ efflux. This is the first attempt to link temperature sensitivity of SOM decomposition with microbial eco-strategies by incorporating the genetic information and disentangling the complex relationship between Q₁₀ and soil microorganisms.

Divergent above- and below-ground biodiversity pathways mediate disturbance impacts on temperate forest multifunctionality

Biodiversity plays a fundamental role in provisioning and regulating forest ecosystem functions and services. Above-ground (plants) and below-ground (soil microbes) biodiversity could have asynchronous change paces to human-driven land-use impacts. Yet, we know very little how they affect the provision of multiple forest functions related to carbon accumulation, water retention capacity and nutrient cycling simultaneously (i.e. ecosystem multifunctionality; EMF). We used a dataset of 22,000 temperate forest trees from 260 plots within 11 permanent forest sites in Northeastern China, which are recovering from three post-logging disturbances. We assessed the direct and mediating effects of multiple attributes of plant biodiversity (taxonomic, phylogenetic, functional and stand structure) and soil biodiversity (bacteria and fungi) on EMF under the three disturbance levels. We found the highest EMF in highly disturbed rather than undisturbed mature forests. Plant taxonomic, phylogenetic, functional and stand structural diversity had both positive and negative effects on EMF, depending on how the EMF index was quantified, whereas soil microbial diversity exhibited a consistent positive impact. Biodiversity indices explained on average 45% (26%-58%) of the variation in EMF, whereas climate and disturbance together explained on average 7% (0.4%-15%). Our result highlighted that the tremendous effect of biodiversity on EMF, largely overpassing those of both climate and disturbance. While above- (β = 0.02-0.19) and below-ground (β = 0.16-0.26) biodiversity had direct positive effects on EMF, their opposite mediating effects (β = -0.22 vs. β = 0.35 respectively) played as divergent pathways to human disturbance impacts on EMF. Our study sheds light on the need for integrative frameworks simultaneously considering above- and below-ground attributes to grasp the global picture of biodiversity effects on ecosystem functioning and services. Suitable management interventions could maintain both plant and soil microbial biodiversity, and thus guarantee a long-term functioning and provisioning of ecosystem services in an increasing disturbance frequency world.

Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 10¹⁵ electronvolts) indicates the existence of the so-called PeVatrons-cosmic-ray factories that accelerate particles to PeV energies. We need to locate and identify such objects to find the origin of Galactic cosmic rays¹. The principal signature of both electron and proton PeVatrons is ultrahigh-energy (exceeding 100 TeV) γ radiation. Evidence of the presence of a proton PeVatron has been found in the Galactic Centre, according to the detection of a hard-spectrum radiation extending to 0.04 PeV (ref. ²). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane^3-6, unbiased identification and in-depth exploration of PeVatrons requires detection of γ-rays with energies well above 0.1 PeV. Here we report the detection of more than 530 photons at energies above 100 teraelectronvolts and up to 1.4 PeV from 12 ultrahigh-energy γ-ray sources with a statistical significance greater than seven standard deviations. Despite having several potential counterparts in their proximity, including pulsar wind nebulae, supernova remnants and star-forming regions, the PeVatrons responsible for the ultrahigh-energy γ-rays have not yet been firmly localized and identified (except for the Crab Nebula), leaving open the origin of these extreme accelerators.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.