Prevotella melaninogenica disrupted oral epithelial barrier function via myosin light chain kinase.

OBJECTIVE: Our previous studies have found that the composition ratio of Prevotella melaninogenica (Pm) on buccal mucosa surface of oral lichen planus (OLP) patients increased significantly compared with control. Furthermore, Pm could invade the epithelium of OLP patients. This study aimed to further explore the impact of Pm on oral keratinocytes. MATERIALS AND METHODS: The Pm-human oral keratinocyte (HOK) co-culture model was established to detect monolayer permeability, zona occludens-1 (ZO-1) expression, and intracellular survival of Pm. We performed RNA-seq followed by identification of differentially expressed genes (DEGs) and Gene Ontology (GO) analysis, with a particular focus on myosin light chain kinase (MLCK). An MLCK inhibitor ML-7 was utilized in Pm-HOK co-culture model to assess its effects on monolayer permeability and ZO-1 expression. RESULTS: HOK monolayer permeability was increased, and ZO-1 expression was decreased after co-culture (p < 0.05). Pm could survive in HOK cells. RNA-seq revealed MLCK was an upregulated common DEG. The expression of MLCK in the Pm-HOK co-culture model was upregulated. Inhibition of MLCK rescued the increased epithelial permeability, and ZO-1 expression was upregulated (p < 0.05). CONCLUSION: MLCK may be involved in disrupting epithelial barrier function by Pm.

Predictors for unfavorable prognosis after stroke with perforator artery disease.

Background and purpose: Perforator artery disease (PAD) is an important subtype of ischemic stroke. The risk factors affecting the prognosis of patients with PAD are unclear. This study aimed to investigate the risk factors affecting the unfavorable prognosis of PAD. Methods: Patients with PAD were enrolled from Dushu Lake Hospital Affiliated to Soochow University and diagnosed as stroke with PAD during the period from September 2021 to July 2023 and followed up with a modified Rankin Scale (mRS) after 90 days, defining the mRS of 0-2 as a group with favorable prognosis, and 3-6 as a group with unfavorable functional outcome. Logistic regression was used to identify predictors for PAD. Multiple logistic regression analysis and receiver operating characteristics (ROC) were used to identify predictors of unfavorable prognosis. Results: Of the 181 enrolled patients, 48 (26.5%) were identified with unfavorable prognosis. On multivariate analysis, increased age (OR = 1.076, 95% CI: 1.012 ~ 1.144, p = 0.019), higher National Institutes of Health Stroke Scale (NIHSS) score at admission (OR = 2.930, 95% CI: 1. 905 ~ 4.508, p < 0.001), and increased neutrophil-to-lymphocyte ratio (NLR) (OR = 3.028, 95% CI: 1.615 ~ 5.675, p = 0.001) were independent risk factors for unfavorable prognosis in patients with PAD, and the area under the receiver operating characteristic curve was 0.590, 0.905, and 0.798, and the multi-factor diagnostic model (Model 2) showed reliable diagnostic specificity and sensitivity (area under the curve = 0.956, p < 0.001, specificity 0.805, sensitivity 0.958, accuracy 0.845). Conclusion: Increased baseline NLR and NIHSS score and aging may be independent risk factors for unfavorable prognosis of patients with PAD. NLR can be used as a potential biological indicator to predict the prognosis of stroke with PAD.

Semi-embedded slot waveguide electro-optic modulator.

Electro-optic modulators are essential devices on silicon photonic chips in modern optical communication networks. This paper presents a compact, low-loss electro-optic modulator. The modulation efficiency is greatly improved by embedding the lower half of the slot waveguide into the buried oxide layer and inserting graphene at the junction. The interaction of graphene with an optical field in a waveguide is studied using the finite element method. The functions of phase modulation and absorption modulation are realized by changing the gate voltage to change the chemical potential of graphene. The semi-embedded slot waveguide optical modulator has a length of 50 µm. After simulation verification, it can be used as an electro-absorption modulator and can achieve a modulation depth of 26.38 dB and an insertion loss of 0.60 dB. When used as an electro-refractive modulator, it can be realized with a linear change of phase from zero to π; the total insertion loss is only 0.59 dB. The modulator has a modulation bandwidth of 79.6 GHz, and the energy consumption as electro-absorption and electro-refraction modulation are 0.51 and 1.92 pj/bit, respectively. Compared with common electro-optic modulators, the electro-optic modulator designed in this paper has a higher modulation effect and also takes into account the advantages of low insertion loss and low energy consumption. This research is helpful for the design of higher-performance optical communication network devices.

Nested micro-ring refractive index sensor based on a subwavelength grating waveguide and the Vernier effect.

In this paper, a nested micro-ring refractive index sensor based on a subwavelength grating waveguide and the Vernier effect is proposed. In this scheme, the nested micro-ring structure is combined with a subwavelength grating structure to enhance the contact area between the optical field and the analyte, and the wavelength offset is doubled through the Vernier effect. The proposed sensor can effectively increase sensing sensitivity, taking into account the improvement of the free spectral range. This structure enables the device to reach a sensitivity of 8030 nm/RIU near 1550 nm wavelength in a deionized water environment, with a detection limit of 5.659×10-5 RIU and free spectral range of 41.956 nm. The device suggested in this study has a greater reduced footprint than the conventional micro-ring resonant sensor, measuring just 35µm×25µm. Due to its high integration, high sensitivity, and large free spectral range compared to conventional micro-ring resonant sensors, such structures are of great value in biosensing and environmental monitoring.

Identifying critical kinematic features of animate motion and contribution to animacy perception.

Humans can distinguish flying birds from drones based solely on motion features when no image information is available. However, it remains unclear which motion features of animate motion induce our animacy perception. To address this, we first analyzed the differences in centroid motion between birds and drones, and discovered that birds exhibit greater acceleration, angular speed, and trajectory fluctuations. We further determined the order of their importance in evoking animacy perception was trajectory fluctuations, acceleration, and speed. More interestingly, people judge whether a moving object is alive using a feature-matching strategy, implying that animacy perception is induced in a key feature-triggered way rather than relying on the accumulation of evidence. Our findings not only shed light on the critical motion features that induce animacy perception and their relative contributions but also have important implications for developing target classification algorithms based on motion features.

Multidimensional tunable graphene chiral metasurface based on coherent control.

The deep application of chiral metasurfaces requires higher flexibility. Herein, we propose a multidimensional tunable chiral graphene metasurface, which uses coherent control to obtain more than 0.8 circular conversion dichroism (CCD) at 2.4 THz as a transmission structure. Its operating frequency can be changed in the 1.3-2.4 THz range, while the amplitude has almost perfect modulation depth in the range of 0-0.8. The mechanism of differential absorption was analyzed through numerical simulation. The device designed is easy to obtain reverse CCD, which is used for unit layout and proves its advantages in near-field imaging. Our work has broadened the path for the development of chiral metasurfaces towards higher degrees of freedom.

Terahertz broadband tunable chiral metamirror based on VO2-metal hybrid structure.

Aiming at the problems of narrow working bandwidth, low efficiency, and complex structure of existing terahertz chiral absorption, we propose a chiral metamirror composed of C-shaped metal split ring and L-shaped vanadium dioxide (VO2). This chiral metamirror is composed of three layers of structure, a gold substrate at the bottom, the first polyethylene cyclic olefin copolymer (Topas) dielectric layer and VO2-metal hybrid structure as the top. Our theoretical results led us to show that this chiral metamirror has a circular dichroism (CD) value greater than 0.9 at 5.70 to 8.55 THz and has a maximum value of 0.942 at f = 7.18 THz. In addition, by adjusting the conductivity of VO2, the CD value can be continuously adjustable from 0 to 0.942, which means that the proposed chiral metamirror supports the free switching of the CD response between the on and off states, and the CD modulation depth exceeds 0.99 in the range of 3 to 10 THz. Moreover, we discuss the influence of structural parameters and the change of incident angle on the performance of the metamirror. Finally, we believe that the proposed chiral metamirror has important reference value in the terahertz range for constructing chiral light detectors, CD metamirrors, switchable chiral absorbers and spin-related systems. This work will provide a new idea for improving the terahertz chiral metamirror operating bandwidth and promote the development of terahertz broadband tunable chiral optical devices.

Boosting charge separation in graphdiyne quantum dots/hollow tubular carbon nitride heterojunction for water pollutant degradation.

Non-desirable solar energy absorption and poor charge transfer efficiency are two problems that limit the peroxymonosulfate (PMS) photocatalytic techniques. Herein, a metal-free boron-doped graphdiyne quantum dot (BGDs) modified hollow tubular g-C3N4 photocatalyst (BGD/TCN) was synthesized to activate PMS and achieved effective space separation of carriers for degradation of bisphenol A. With 0.5 mM PMS, the degradation rate of bisphenol A (20 ppm) was 0.0634 min-1, 3.7-fold higher than that of TCN itself. The roles of BGDs in the distribution of electrons and photocatalytic property were well identified by experiments and density functional theory (DFT) calculations. The possible degradation intermediate products of bisphenol A were monitored by mass spectrometer and demonstrated to be nontoxic using ecological structure activity relationship modeling (ECOSAR). Finally, this newly-designed material was successfully applied in actual water bodies, which further renders its promising prospect for actual water remediation.

Thermodynamic restrictions determine ammonia tolerance of methanogenic pathways in Methanosarcina barkeri.

Ammonia is a ubiquitous potential inhibitor of anaerobic digestion processes, mainly exhibiting inhibition towards methanogenic activity. However, knowledge as to how ammonia affects the methanogens is still limited. In this study, we cultured a multitrophic methanogen, Methanosarcina barkeri DSM 800, with acetate, H2/CO2, and methanol to evaluate the influence of ammonia on different methanogenic pathways. Aceticlastic methanogenesis was more sensitive to increased ammonia concentrations than hydrogenotrophic and methylotrophic methanogenesis. Theoretical maximum NH3 tolerances of M. barkeri fed with acetate, H2/CO2, and methanol were calculated to be 39.1 ± 9.0, 104.3 ± 7.4, and 85.7 ± 1.0 mg/L, respectively. The order of the ΔG range of M. barkeri under three methanogenic pathways reflected the order of ammonia tolerance of M. barkeri. Our results provide insights into the role of the thermodynamic potential of methanogenesis on the tolerance of ammonia stress; and shed light on the mechanism of ammonia inhibition on anaerobic digestion.

Genome-wide analysis of MYB family in Nicotiana benthamiana and the functional role of the key members in resistance to Bemisia tabaci.

MYB transcription factors (TFs) play a key role in plant resistance to abiotic and biotical stresses. However, little is currently known about their involvement in the plant defense to piercing-sucking insects. Here, we studied the MYB TFs that responded to and resisted Bemisia tabaci whitefly in the model plant Nicotiana benthamiana. Firstly, a total of 453 NbMYB TFs in N. benthamiana genome were identified and 182 R2R3-MYB TFs were analyzed for molecular characteristics, phylogenetic analysis, genetic structure, motif composition, and cis-elements. Then, six stress-related NbMYB genes were selected for further study. The expression pattern shows they were highly expressed in mature leaves and intensively induced upon whitefly attack. Combined with bioinformatic analysis, overexpression, ß-Glucuronidase (GUS) assay, and virus-induced silencing tests, we determined the transcriptional regulation of these NbMYBs on the genes in lignin biosynthesis and SA-signaling pathways. Meanwhile, we tested the performance of whitefly on plants with increased or silenced NbMYB genes expression and found that NbMYB42, NbMYB107, NbMYB163, and NbMYB423 were resistant to whitefly. Our results contribute to a comprehensive understanding of the MYB TFs in N. benthamiana. Furthermore, our findings will facilitate further studies on the role of MYB TFs in the interaction between plants and piercing-sucking insects.

Targeting Myocardial Mitochondria-STING-Polyamine Axis Prevents Cardiac Hypertrophy in Chronic Kidney Disease.

Chronic kidney disease (CKD) is well recognized as a distinct contributor to cardiac hypertrophy, while the underlying mechanism remains incompletely understood. Here, the authors show that myocardial mitochondrial oxidative damage is early and prominent in CKD and distinctively stimulates the STING-NFκB pathway by releasing mitochondrial DNA to drive cardiac hypertrophy. Furthermore, the authors reveal that ornithine decarboxylase (ODC1)-putrescine metabolic flux is transactivated by NFκB and is required for the STING-NFκB pathway to drive cardiac hypertrophy. Finally, genetic or pharmacologic inhibition of the myocardial mitochondria-STING-NFκB-ODC1 axis significantly prevents CKD-associated cardiac hypertrophy. Therefore, targeting the myocardial mitochoandria-STING-NFκB-ODC1 axis is a promising therapeutic strategy for cardiac hypertrophy in patients with CKD.

Synthesis and Characterization of Fucoidan-Chitosan Nanoparticles Targeting P-Selectin for Effective Atherosclerosis Therapy.

Atherosclerosis is the key pathogenesis of cardiovascular diseases; oxidative stress, which is induced by the generated excess reactive oxygen species (ROS), has been a crucial mechanism underlying this pathology. Nanoparticles (NPs) represent a novel strategy for the development of potential therapies against atherosclerosis, and multifunctional NPs possessing antioxidative capacities hold promise for amelioration of vascular injury caused by ROS and for evading off-target effects; materials that are currently used for NP synthesis often serve as vehicles that do not possess intrinsic biological activities; however, they may affect the surrounding healthy environment due to decomposition of products. Herein, we used nontoxic fucoidan, a sulfated polysaccharide derived from a marine organism, to develop chitosan-fucoidan nanoparticles (CFNs). Then, by binding to P-selectin, an inflammatory adhesion exhibited molecule expression on the endothelial cells and activated platelets, blocking leukocyte recruitment and rolling on platelets and endothelium. CFNs exhibit antioxidant and anti-inflammatory properties. Nevertheless, by now, the application of CFNs for the target delivery regarding therapeutics specific to atherosclerotic plaques is not well investigated. The produced CFNs were physicochemically characterized using transmission electron microscopy (TEM), together with Fourier transform infrared spectroscopy (FTIR). Evaluations of the in vitro antioxidant as well as anti-inflammatory activities exhibited by CFNs were based on the measurement of their ROS scavenging abilities and investigating inflammatory mediator levels. The in vivo pharmacokinetics and binding efficiency of the CFNs to atherosclerotic plaques were also evaluated. The therapeutic effects indicated that CFNs effectively suppressed local oxidative stress and inflammation by targeting P-selectin in atheromatous plaques and thereby preventing the progression of atherosclerosis.

Function of normal oral mucosa revealed by single-cell RNA sequencing.

The functions of oral mucosa include barrier, sensation, and secretion. The barrier protection function is particularly important, which includes physical barrier and immunological barrier. Few studies have revealed the function of oral mucosa by displaying the map of normal oral mucosal cells from the perspective of single cells. Here, single-cell transcriptome sequencing was used to bring a relatively comprehensive map of the normal oral mucosal cells. In total, 26,398 cells from three cases of normal oral mucosa were analyzed by single-cell RNA-sequencing and 14 distinct cell groups were defined, 7 of which were immune cells. We performed subgroup classification and heterogeneity analysis of epithelial cells, T cells, and macrophagocytes, which found a subpopulation of epithelial cells with high expression of major histocompatibility complex class II molecules, a subpopulation CD8+ GZMK+ T cells, and two kinds of active macrophagocytes. Meanwhile, we identified ligand-receptor pairs among the major cell types to explore the interactions and how they maintain the homeostasis of normal oral mucosa. Based on these results, the epithelial barrier function, immunological barrier function, and potential maintenance function of stromal cells in the oral mucosa were described at the single-cell level, which provides basic data resources for further studies of oral mucosal diseases.

DUSP2-mediated inhibition of tubular epithelial cell pyroptosis confers nephroprotection in acute kidney injury.

Rationale: Acute kidney injury (AKI) is pathologically characterized by renal tubular epithelial cell (RTEC) death and interstitial inflammation, while their pathogenesis remains incompletely understood. Dual-specificity phosphatase 2 (DUSP2) recently emerges as a crucial regulator of cell death and inflammation in a wide range of diseases, but its roles in renal pathophysiology are largely unknown. Methods: The expression of DUSP2 in the kidney was characterized by histological analysis in renal tissues from patients and mice with AKI. The role and mechanism of DUSP2-mediated inhibition of tubular epithelial cell pyroptosis in AKI were evaluated both in vivo and in vitro, and confirmed in RTEC-specific deletion of DUSP2 mice. Results: Here, we show that DUSP2 is enriched in RTECs in the renal tissue of both human and mouse and mainly positions in the nucleus. Further, we reveal that loss-of-DUSP2 in RTECs not only is a common feature of human and murine AKI but also positively contributes to AKI pathogenesis. Especially, RTEC-specific deletion of DUSP2 sensitizes mice to AKI by promoting RTEC pyroptosis and the resultant interstitial inflammation. Mechanistic studies show that gasdermin D (GSDMD), which mediates RTEC pyroptosis, is identified as a transcriptional target of activated STAT1 during AKI, whereas DUSP2 as a nuclear phosphatase deactivates STAT1 to restrict GSDMD-mediated RTEC pyroptosis. Importantly, DUSP2 overexpression in RTECs via adeno-associated virus-mediated gene transfer significantly ameliorates AKI. Conclusion: Our findings demonstrate a hitherto unrecognized role of DUSP2-STAT1 axis in regulating RTEC pyroptosis in AKI, highlighting that DUSP2-STAT1 axis is an attractive therapeutic target for AKI.

MetNC: Predicting Metabolites in vivo for Natural Compounds.

Natural compounds (NCs) undergo complicated biotransformation in vivo to produce diverse forms of metabolites dynamically, many of which are of high medicinal value. Predicting the profiles of chemical products may help to narrow down possible candidates, yet current computational methods for predicting biotransformation largely focus on synthetic compounds. Here, we proposed a method of MetNC, a tailor-made method for NC biotransformation prediction, after exploring the overall patterns of NC in vivo metabolism. Based on 850 pairs of the biotransformation dataset validated by comprehensive in vivo experiments with sourcing compounds from medicinal plants, MetNC was designed to produce a list of potential metabolites through simulating in vivo biotransformation and then prioritize true metabolites into the top list according to the functional groups in compound structures and steric hindrance around the reaction sites. Among the well-known peers of GLORYx and BioTransformer, MetNC gave the highest performance in both the metabolite coverage and the ability to short-list true products. More importantly, MetNC seemed to display an extra advantage in recommending the microbiota-transformed metabolites, suggesting its potential usefulness in the overall metabolism estimation. In summary, complemented to those techniques focusing on synthetic compounds, MetNC may help to fill the gap of natural compound metabolism and narrow down those products likely to be identified in vivo.

Small RNA and Degradome Sequencing Reveal Important MicroRNA Function in Nicotiana tabacum Response to Bemisia tabaci.

MicroRNAs (miRNAs), a class of small non-coding regulatory RNAs, are key molecules in many biological and metabolic processes of plant growth, development and stress response via targeting mRNAs. The phloem-feeding insect whitefly Bemisia tabaci (Hemiptera, Aleyrodidae) is a serious pest that causes devastating harm to agricultural production worldwide. However, the function of host miRNAs in the response to whitefly infestation remains unclear. Here, we sequenced the small RNA and degradome of tobacco (Nicotiana tabacum L.), after and before infestation by B. tabaci. We identified 1291 miRNAs belonging to 138 miRNA families including 706 known miRNAs and 585 novel miRNAs. A total of 47 miRNAs were differentially expressed, of which 30 were upregulated and 17 were downregulated by whitefly exposure. Then, computational analysis showed that the target genes of differential miRNAs were involved in R gene regulation, plant innate immunity, plant pathogen defense, the plant hormone signal pathway and abiotic stress tolerance. Furthermore, degradome analysis demonstrated that 253 mRNAs were cleaved by 66 miRNAs. Among them, the targets cleaved by upregulated miR6025, miR160, miR171, miR166 and miR168 are consistent with our prediction, suggesting that pathogen-related miRNAs may function in plant defense against whitefly. Moreover, our results show that plant miRNA response and miRNA-mediated post-transcriptional regulation for phloem-feeding insect infestation are similar to pathogen invasion. Our study provides additional data to further elucidate how host plants respond and defend the phloem-feeding insects.

Renal Klotho and inorganic phosphate are extrinsic factors that antagonistically regulate hematopoietic stem cell maintenance.

The composition and origin of extrinsic cues required for hematopoietic stem cell (HSC) maintenance are incompletely understood. Here we identify renal Klotho and inorganic phosphate (Pi) as extrinsic factors that antagonistically regulate HSC maintenance in the bone marrow (BM). Disruption of the Klotho-Pi axis by renal Klotho deficiency or Pi excess causes Pi overload in the BM niche and Pi retention in HSCs, leading to alteration of HSC maintenance. Mechanistically, Pi retention is mediated by soluble carrier family 20 member 1 (SLC20A1) and sensed by diphosphoinositol pentakisphosphate kinase 2 (PPIP5K2) to enhance Akt activation, which then upregulates SLC20A1 to aggravate Pi retention and augments GATA2 activity to drive the expansion and megakaryocyte/myeloid-biased differentiation of HSCs. However, kidney-secreted soluble Klotho directly maintains HSC pool size and differentiation by restraining SLC20A1-mediated Pi absorption of HSCs. These findings uncover a regulatory role of the Klotho-Pi axis orchestrated by the kidneys in BM HSC maintenance.

Genome-Wide Analysis of the Peroxidase Gene Family and Verification of Lignin Synthesis-Related Genes in Watermelon.

Watermelon (Citrullus lanatus) is an important horticultural crop worldwide, but peel cracking caused by peel hardness severely decreases its quality. Lignification is one of the important functions of class III peroxidase (PRX), and its accumulation in the plant cell wall leads to cell thickening and wood hardening. For in-depth physiological and genetical understanding, we studied the relationship between peel hardness and lignin accumulation and the role of PRXs affecting peel lignin biosynthesis using genome-wide bioinformatics analysis. The obtained results showed that lignin accumulation gradually increased to form the peel stone cell structure, and tissue lignification led to peel hardness. A total of 79 ClPRXs (class III) were identified using bioinformatics analysis, which were widely distributed on 11 chromosomes. The constructed phylogenetics indicated that ClPRXs were divided into seven groups and eleven subclasses, and gene members of each group had highly conserved intron structures. Repeated pattern analysis showed that deletion and replication events occurred during the process of ClPRX amplification. However, in the whole-protein sequence alignment analysis, high homology was not observed, although all contained four conserved functional sites. Repeated pattern analysis showed that deletion and replication events occurred during ClPRXs' amplification process. The prediction of the promoter cis-acting element and qRT-PCR analysis in four tissues (leaf, petiole, stem, and peel) showed different expression patterns for tissue specificity, abiotic stress, and hormone response by providing a genetic basis of the ClPRX gene family involved in a variety of physiological processes in plants. To our knowledge, we for the first time report the key roles of two ClPRXs in watermelon peel lignin synthesis. In conclusion, the extensive data collected in this study can be used for additional functional analysis of ClPRXs in watermelon growth and development and hormone and abiotic stress response.

HIT 2.0: an enhanced platform for Herbal Ingredients' Targets.

Literature-described targets of herbal ingredients have been explored to facilitate the mechanistic study of herbs, as well as the new drug discovery. Though several databases provided similar information, the majority of them are limited to literatures before 2010 and need to be updated urgently. HIT 2.0 was here constructed as the latest curated dataset focusing on Herbal Ingredients' Targets covering PubMed literatures 2000-2020. Currently, HIT 2.0 hosts 10 031 compound-target activity pairs with quality indicators between 2208 targets and 1237 ingredients from more than 1250 reputable herbs. The molecular targets cover those genes/proteins being directly/indirectly activated/inhibited, protein binders, and enzymes substrates or products. Also included are those genes regulated under the treatment of individual ingredient. Crosslinks were made to databases of TTD, DrugBank, KEGG, PDB, UniProt, Pfam, NCBI, TCM-ID and others. More importantly, HIT enables automatic Target-mining and My-target curation from daily released PubMed literatures. Thus, users can retrieve and download the latest abstracts containing potential targets for interested compounds, even for those not yet covered in HIT. Further, users can log into 'My-target' system, to curate personal target-profiling on line based on retrieved abstracts. HIT can be accessible at http://hit2.badd-cao.net.

Leaf surface microtopography shaping the bacterial community in the phyllosphere: evidence from 11 tree species.

Phyllosphere bacteria are an important component of environmental microbial communities and are closely related to plant health and ecosystem stability. However, the relationships between the inhabitation and assembly of phyllosphere bacteria and leaf microtopography are still obscure. In this study, the phyllosphere bacterial communities and leaf microtopographic features (vein density, stomatal length, and density) of eleven tree species were fully examined. Both the absolute abundance and diversity of phyllosphere bacterial communities were significantly different among the tree species, and leaf vein density dominated the variation. TITAN analysis showed that leaf vein density also played more important roles in regulating the relative abundance of bacteria than stomatal features, and 6 phyla and 62 genera of phyllosphere bacteria showed significant positive responses to leaf vein density. Moreover, LEfSe analysis showed that the leaves with higher vein density had more bacterial biomarkers. Leaf vein density also changed the co-occurrence pattern of phyllosphere bacteria, and the co-occurrence network demonstrated more negative correlations and more nodes on the leaves with larger leaf vein density, indicating that higher densities of leaf veins improved the stability of the phyllosphere bacterial community. Phylogenetic analysis showed that deterministic processes (especially homogeneous selection) dominated the assembly process of phyllosphere bacterial communities. The leaf vein density increased the degree of bacterial clustering at the phylogenetic level. Therefore, the inhabitation and assembly of the phyllosphere bacterial community are related to leaf microtopography, which provides deeper insight into the interaction between plants and bacteria.