Systematic Transcriptome Profiling of hPSC-derived Osteoblasts Unveils CORIN's Mastery in Governing Osteogenesis Through CEBPD Modulation.

The commitment of stem cells to differentiate into osteoblasts is a highly regulated and complex process that involves the coordination of extrinsic signals and intrinsic transcriptional machinery. While rodent osteoblastic differentiation has been extensively studied, research on human osteogenesis has been limited by cell sources and existing models. Here, we systematically dissect hPSC-derived osteoblasts to identify functional membrane proteins and their downstream transcriptional networks involved in human osteogenesis. Our results reveal an enrichment of type II transmembrane serine protease CORIN in humans but not rodent osteoblasts. Functional analyses demonstrated that CORIN depletion significantly impairs osteogenesis. Genome-wide ChIP enrichment and mechanistic studies show that p38 MAPK-mediated CEBPD upregulation is required for CORIN-modulated osteogenesis. Contrastingly, the type I transmembrane heparan sulfate proteoglycan SDC1 enriched in MSCs exerts a negative regulatory effect on osteogenesis through a similar mechanism. ChIP-seq, bulk and single-cell transcriptomes, and functional validations indicated that CEBPD plays a critical role in controlling osteogenesis. In summary, our findings uncover previously unrecognized CORIN-mediated CEBPD transcriptomic networks in driving human osteoblast lineage commitment.

Individualised neoantigen therapy mRNA-4157 (V940) plus pembrolizumab versus pembrolizumab monotherapy in resected melanoma (KEYNOTE-942): a randomised, phase 2b study.

BACKGROUND: Checkpoint inhibitors are standard adjuvant treatment for stage IIB-IV resected melanoma, but many patients recur. Our study aimed to evaluate whether mRNA-4157 (V940), a novel mRNA-based individualised neoantigen therapy, combined with pembrolizumab, improved recurrence-free survival and distant metastasis-free survival versus pembrolizumab monotherapy in resected high-risk melanoma. METHODS: We did an open-label, randomised, phase 2b, adjuvant study of mRNA-4157 plus pembrolizumab versus pembrolizumab monotherapy in patients, enrolled from sites in the USA and Australia, with completely resected high-risk cutaneous melanoma. Patients with completely resected melanoma (stage IIIB-IV) were assigned 2:1 to receive open-label mRNA-4157 plus pembrolizumab or pembrolizumab monotherapy. mRNA-4157 was administered intramuscularly (maximum nine doses) and pembrolizumab intravenously (maximum 18 doses) in 3-week cycles. The primary endpoint was recurrence-free survival in the intention-to-treat population. This ongoing trial is registered at ClinicalTrials.gov, NCT03897881. FINDINGS: From July 18, 2019, to Sept 30, 2021, 157 patients were assigned to mRNA-4157 plus pembrolizumab combination therapy (n=107) or pembrolizumab monotherapy (n=50); median follow-up was 23 months and 24 months, respectively. Recurrence-free survival was longer with combination versus monotherapy (hazard ratio [HR] for recurrence or death, 0·561 [95% CI 0·309-1·017]; two-sided p=0·053), with lower recurrence or death event rate (24 [22%] of 107 vs 20 [40%] of 50); 18-month recurrence-free survival was 79% (95% CI 69·0-85·6) versus 62% (46·9-74·3). Most treatment-related adverse events were grade 1-2. Grade ≥3 treatment-related adverse events occurred in 25% of patients in the combination group and 18% of patients in the monotherapy group, with no mRNA-4157-related grade 4-5 events. Immune-mediated adverse event frequency was similar for the combination (37 [36%]) and monotherapy (18 [36%]) groups. INTERPRETATION: Adjuvant mRNA-4157 plus pembrolizumab prolonged recurrence-free survival versus pembrolizumab monotherapy in patients with resected high-risk melanoma and showed a manageable safety profile. These results provide evidence that an mRNA-based individualised neoantigen therapy might be beneficial in the adjuvant setting. FUNDING: Moderna in collaboration with Merck Sharp & Dohme, a subsidiary of Merck & Co, Rahway, NJ, USA.

Hereditary retinoblastoma iPSC model reveals aberrant spliceosome function driving bone malignancies.

The RB1 gene is frequently mutated in human cancers but its role in tumorigenesis remains incompletely defined. Using an induced pluripotent stem cell (iPSC) model of hereditary retinoblastoma (RB), we report that the spliceosome is an up-regulated target responding to oncogenic stress in RB1-mutant cells. By investigating transcriptomes and genome occupancies in RB iPSC­derived osteoblasts (OBs), we discover that both E2F3a, which mediates spliceosomal gene expression, and pRB, which antagonizes E2F3a, coregulate more than one-third of spliceosomal genes by cobinding to their promoters or enhancers. Pharmacological inhibition of the spliceosome in RB1-mutant cells leads to global intron retention, decreased cell proliferation, and impaired tumorigenesis. Tumor specimen studies and genome-wide TCGA (The Cancer Genome Atlas) expression profile analyses support the clinical relevance of pRB and E2F3a in modulating spliceosomal gene expression in multiple cancer types including osteosarcoma (OS). High levels of pRB/E2F3a­regulated spliceosomal genes are associated with poor OS patient survival. Collectively, these findings reveal an undiscovered connection between pRB, E2F3a, the spliceosome, and tumorigenesis, pointing to the spliceosomal machinery as a potentially widespread therapeutic vulnerability of pRB-deficient cancers.

Patient-derived iPSCs link elevated mitochondrial respiratory complex I function to osteosarcoma in Rothmund-Thomson syndrome.

Rothmund-Thomson syndrome (RTS) is an autosomal recessive genetic disorder characterized by poikiloderma, small stature, skeletal anomalies, sparse brows/lashes, cataracts, and predisposition to cancer. Type 2 RTS patients with biallelic RECQL4 pathogenic variants have multiple skeletal anomalies and a significantly increased incidence of osteosarcoma. Here, we generated RTS patient-derived induced pluripotent stem cells (iPSCs) to dissect the pathological signaling leading to RTS patient-associated osteosarcoma. RTS iPSC-derived osteoblasts showed defective osteogenic differentiation and gain of in vitro tumorigenic ability. Transcriptome analysis of RTS osteoblasts validated decreased bone morphogenesis while revealing aberrantly upregulated mitochondrial respiratory complex I gene expression. RTS osteoblast metabolic assays demonstrated elevated mitochondrial respiratory complex I function, increased oxidative phosphorylation (OXPHOS), and increased ATP production. Inhibition of mitochondrial respiratory complex I activity by IACS-010759 selectively suppressed cellular respiration and cell proliferation of RTS osteoblasts. Furthermore, systems analysis of IACS-010759-induced changes in RTS osteoblasts revealed that chemical inhibition of mitochondrial respiratory complex I impaired cell proliferation, induced senescence, and decreased MAPK signaling and cell cycle associated genes, but increased H19 and ribosomal protein genes. In summary, our study suggests that mitochondrial respiratory complex I is a potential therapeutic target for RTS-associated osteosarcoma and provides future insights for clinical treatment strategies.

Phase 1b study of pan-AKT inhibitor vevorisertib alone or with paclitaxel or fulvestrant in PIK3CA/AKT/PTEN-mutated advanced solid tumors.

BACKGROUND: In this first-in-human phase 1b study (ClinicalTrials.gov identifier NCT02761694) of advanced solid tumors with PIK3CA/AKT/PTEN mutations, the authors investigated the safety and efficacy of the pan-AKT inhibitor vevorisertib (MK-4440; ARQ 751) as monotherapy or with paclitaxel or fulvestrant. METHODS: Patients with histologically confirmed, advanced or recurrent, PIK3CA/AKT/PTEN-mutated solid tumors, measurable disease according to Response Evaluation Criteria in Solid Tumors, version 1.1, and an Eastern Cooperative Oncology Group performance status ≤1 received vevorisertib (dose range, 5-100 mg) alone or with paclitaxel 80 mg/m2 or fulvestrant 500 mg. The primary end point was safety and tolerability. Secondary end points included pharmacokinetics and the objective response rate according to Response Evaluation Criteria in Solid Tumors, version 1.1. RESULTS: Of 78 patients enrolled, 58 received vevorisertib monotherapy, 10 received vevorisertib plus paclitaxel, and nine received vevorisertib plus fulvestrant. Dose-limiting toxicity occurred in three patients (vevorisertib monotherapy, n = 2 [grade 3 pruritic and maculopapular rashes]; vevorisertib plus paclitaxel, n = 1 [grade 1 asthenia]). Across doses, treatment-related AEs occurred in 46 patients (79%) with vevorisertib monotherapy, in 10 patients (100%) with vevorisertib plus paclitaxel, and in nine patients (100%) with vevorisertib plus fulvestrant; and grade 3 treatment-related AEs occurred in 13 (22%), 7 (70%), and 3 (33%) patients, respectively. No grade 4/5 treatment-related AEs occurred. Maximum vevorisertib concentrations were reached 1-4 hours after dosing; the elimination half-life ranged from 8.8 to 19.3 hours. The objective response rate was 5% with vevorisertib monotherapy (three partial responses), 20% with vevorisertib plus paclitaxel (two partial responses), and 0% with vevorisertib plus fulvestrant. CONCLUSIONS: Vevorisertib alone or with paclitaxel or fulvestrant had a manageable safety profile, and vevorisertib alone or with paclitaxel had minimal to modest antitumor activity in this patient population with PIK3CA/AKT/PTEN-mutated advanced solid tumors. CLINICAL TRIAL REGISTRATION: ClinicalTrials.gov, NCT02761694.

Data denoising with transfer learning in single-cell transcriptomics.

Single-cell RNA sequencing (scRNA-seq) data are noisy and sparse. Here, we show that transfer learning across datasets remarkably improves data quality. By coupling a deep autoencoder with a Bayesian model, SAVER-X extracts transferable gene-gene relationships across data from different labs, varying conditions and divergent species, to denoise new target datasets.

MicrobiomeExplorer: an R package for the analysis and visualization of microbial communities.

SUMMARY: We developed the MicrobiomeExplorer R package to facilitate the analysis and visualization of microbial communities. The MicrobiomeExplorer R package allows a user to perform typical microbiome analytic workflows and visualize their results, either through the command line or an interactive Shiny application included with the package. In addition to applying common analytical workflows, the application enables automated analysis report generation. AVAILABILITY AND IMPLEMENTATION: Available at https://github.com/zoecastillo/microbiomeExplorer. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SAVER: gene expression recovery for single-cell RNA sequencing.

In single-cell RNA sequencing (scRNA-seq) studies, only a small fraction of the transcripts present in each cell are sequenced. This leads to unreliable quantification of genes with low or moderate expression, which hinders downstream analysis. To address this challenge, we developed SAVER (single-cell analysis via expression recovery), an expression recovery method for unique molecule index (UMI)-based scRNA-seq data that borrows information across genes and cells to provide accurate expression estimates for all genes.

Gene expression distribution deconvolution in single-cell RNA sequencing.

Single-cell RNA sequencing (scRNA-seq) enables the quantification of each gene's expression distribution across cells, thus allowing the assessment of the dispersion, nonzero fraction, and other aspects of its distribution beyond the mean. These statistical characterizations of the gene expression distribution are critical for understanding expression variation and for selecting marker genes for population heterogeneity. However, scRNA-seq data are noisy, with each cell typically sequenced at low coverage, thus making it difficult to infer properties of the gene expression distribution from raw counts. Based on a reexamination of nine public datasets, we propose a simple technical noise model for scRNA-seq data with unique molecular identifiers (UMI). We develop deconvolution of single-cell expression distribution (DESCEND), a method that deconvolves the true cross-cell gene expression distribution from observed scRNA-seq counts, leading to improved estimates of properties of the distribution such as dispersion and nonzero fraction. DESCEND can adjust for cell-level covariates such as cell size, cell cycle, and batch effects. DESCEND's noise model and estimation accuracy are further evaluated through comparisons to RNA FISH data, through data splitting and simulations and through its effectiveness in removing known batch effects. We demonstrate how DESCEND can clarify and improve downstream analyses such as finding differentially expressed genes, identifying cell types, and selecting differentiation markers.

Targetable BET proteins- and E2F1-dependent transcriptional program maintains the malignancy of glioblastoma.

Competitive BET bromodomain inhibitors (BBIs) targeting BET proteins (BRD2, BRD3, BRD4, and BRDT) show promising preclinical activities against brain cancers. However, the BET protein-dependent glioblastoma (GBM)-promoting transcriptional network remains elusive. Here, with mechanistic exploration of a next-generation chemical degrader of BET proteins (dBET6), we reveal a profound and consistent impact of BET proteins on E2F1- dependent transcriptional program in both differentiated GBM cells and brain tumor-initiating cells. dBET6 treatment drastically reduces BET protein genomic occupancy, RNA-Pol2 activity, and permissive chromatin marks. Subsequently, dBET6 represses the proliferation, self-renewal, and tumorigenic ability of GBM cells. Moreover, dBET6-induced degradation of BET proteins exerts superior antiproliferation effects compared to conventional BBIs and overcomes both intrinsic and acquired resistance to BBIs in GBM cells. Our study reveals crucial functions of BET proteins and provides the rationale and therapeutic merits of targeted degradation of BET proteins in GBM.

Two-Stage Fast DOA Estimation Based on Directional Antennas in Conformal Uniform Circular Array.

In conformal array radar, due to the directivity of antennas, the responses of the echo signals between different antennas are distinct, and some antennas cannot even receive the target echo signal. These phenomena significantly affect the accuracy of direction-of-arrival (DOA) estimation. To implement accurate DOA estimation in a conformal uniform circular array (UCA) composed of directional antennas, the two-stage fast DOA estimation algorithm is proposed. In the pre-processing stage, multi-target decoupling and target detection are mainly used to obtain the targets' range bin indexes set; in the rough-precise DOA estimation stage, the amplitude and phase information of each antenna are used for rough DOA estimation and precise DOA estimation, respectively. Based on simulation and actual anechoic chamber radar experiments, and through quantitative analyses of the accuracy, validity and elapsed time of the two-stage fast DOA estimation algorithm compared with the directional antenna MUSIC (DA-MUSIC), sub-array MUSIC (S-MUSIC) and Capon-like algorithms, results indicate that the two-stage fast DOA estimation algorithm can rapidly and accurately estimate DOAs in a multi-target scenario without the range-angle pair-matching procedure. Lower computational complexity and superior estimation accuracy provide the two-stage fast DOA estimation algorithm a broader application prospect in the practical engineering field.

Classification Method of Uniform Circular Array Radar Ground Clutter Data Based on Chaotic Genetic Algorithm.

The classification and recognition of radar clutter is helpful to improve the efficiency of radar signal processing and target detection. In order to realize the effective classification of uniform circular array (UCA) radar clutter data, a classification method of ground clutter data based on the chaotic genetic algorithm is proposed. In this paper, the characteristics of UCA radar ground clutter data are studied, and then the statistical characteristic factors of correlation, non-stationery and range-Doppler maps are extracted, which can be used to classify ground clutter data. Based on the clustering analysis, results of characteristic factors of radar clutter data under different wave-controlled modes in multiple scenarios, we can see: in radar clutter clustering of different scenes, the chaotic genetic algorithm can save 34.61% of clustering time and improve the classification accuracy by 42.82% compared with the standard genetic algorithm. In radar clutter clustering of different wave-controlled modes, the timeliness and accuracy of the chaotic genetic algorithm are improved by 42.69% and 20.79%, respectively, compared to standard genetic algorithm clustering. The clustering experiment results show that the chaotic genetic algorithm can effectively classify UCA radar's ground clutter data.

Influence of gene position on the expression divergence of oxidative response genes in intraspecific yeast.

Phenotypic variation can arise from differences in the protein coding sequence and in the regulatory elements. However, little is known about the contribution of regulatory difference to the expression divergence, especially the cis and trans regulatory variation to the expression divergence in intraspecific populations. In this study, we used two different yeast strains, BY4743 and RM11-1a/α, to study the regulatory variation to the expression divergence between BY and RM under oxidative stress condition. Our results indicated that the expression divergence of BY and RM is mainly due to trans regulatory variations under both normal and oxidative stress conditions. However, cis regulatory variation seems to play a very important role in oxidative stress response in yeast because 36% of genes showed an increase in cis regulatory variation effect compared with 13% of genes that showed an increase in trans regulatory variation effect after oxidative stress. Our data also indicated that genes located on the longer arm of the chromosomes are more susceptible to cis variation effect under oxidative stress than genes on the shorter arm of the chromosomes.

Synthesis and biological evaluation of novel SIPI-7623 derivatives as farnesoid X receptor (FXR) antagonists.

Most of reported steroidal FXR antagonists are restricted due to low potency. We described the design and synthesis of novel nonsteroidal scaffold SIPI-7623 derivatives as FXR antagonists. The most potent compound A-11 (IC50 = 7.8 ± 1.1 µM) showed better activity compared to SIPI-7623 (IC50 = 40.8 ± 1.7 µM) and guggulsterone (IC50 = 45.9 ± 1.1 µM). Docking of A-11 in FXR's ligand-binding domain was also studied.

BreakPoint Surveyor: a pipeline for structural variant visualization.

SUMMARY: BreakPoint Surveyor (BPS) is a computational pipeline for the discovery, characterization, and visualization of complex genomic rearrangements, such as viral genome integration, in paired-end sequence data. BPS facilitates interpretation of structural variants by merging structural variant breakpoint predictions, gene exon structure, read depth, and RNA-sequencing expression into a single comprehensive figure. AVAILABILITY AND IMPLEMENTATION: Source code and sample data freely available for download at https://github.com/ding-lab/BreakPointSurveyor, distributed under the GNU GPLv3 license, implemented in R, Python and BASH scripts, and supported on Unix/Linux/OS X operating systems. CONTACT: lding@wustl.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

[Metabolites of tanshinone â and tanshinone â ¡A in vivo in rats].

An efficient method of ultra-high performance liquid chromatography coupled with linear ion trap-Orbitrap (UHPLC-LTQ-Orbitrap) mass spectrometer was established to elucidate the in vivo metabolites of tanshinone Ⅰ and tanshinone ⅡA in rats. Urine and plasma samples were collected after oral gavage. After processing biological sample by solid phase extraction, Waters ACQUITY HPLC BEH C18 column (2.1 mm×100 mm, 1.7 µm) was used with 0.1% formic acid (A) - acetonitrile (B) solution as the mobile phase for gradient elution. The plasma, urine and the blank samples were then analyzed by ESI-LTQ-Orbitrap equipped with an ESI ion source under positive ion mode. On the basis of the accurate mass measurements, multiple mass spectra and comparison of data with published literature, a total of 26 metabolites were tentatively identified and characterized in the rat samples. Among them, 7 metabolites were derived from tanshinone Ⅰ through metabolic pathways of glucuronide conjugation, hydroxylation, reduction reaction, demethylation reaction, methylation, sulfate conjugation and their composite reactions. Nineteen metabolites were derived from tanshinone ⅡA through metabolic pathways of hydroxylation, reduction reaction, methylation, sulfate conjugation, glucuronidation, glucosylation and their complicated reactions. The results showed that the metabolism of tanshinone Ⅰ and tanshinone ⅡA in rats could be comprehensively clarified by using UHPLC-LTQ-Orbitrap mass spectrometer, providing material basis for the further research in terms of pharmacodynamics, toxicology, and secondary development of Chinese medicine.

Paper-Based MicroRNA Expression Profiling from Plasma and Circulating Tumor Cells.

BACKGROUND: Molecular characterization of circulating tumor cells (CTCs) holds great promise for monitoring metastatic progression and characterizing metastatic disease. However, leukocyte and red blood cell contamination of routinely isolated CTCs makes CTC-specific molecular characterization extremely challenging. METHODS: Here we report the use of a paper-based medium for efficient extraction of microRNAs (miRNAs) from limited amounts of biological samples such as rare CTCs harvested from cancer patient blood. Specifically, we devised a workflow involving the use of Flinders Technology Associates (FTA)® Elute Card with a digital PCR-inspired "partitioning" method to extract and purify miRNAs from plasma and CTCs. RESULTS: We demonstrated the sensitivity of this method to detect miRNA expression from as few as 3 cancer cells spiked into human blood. Using this method, background miRNA expression was excluded from contaminating blood cells, and CTC-specific miRNA expression profiles were derived from breast and colorectal cancer patients. Plasma separated out during purification of CTCs could likewise be processed using the same paper-based method for miRNA detection, thereby maximizing the amount of patient-specific information that can be derived from a single blood draw. CONCLUSIONS: Overall, this paper-based extraction method enables an efficient, cost-effective workflow for maximized recovery of small RNAs from limited biological samples for downstream molecular analyses.

A Low Cost Bluetooth Low Energy Transceiver for Wireless Sensor Network Applications with a Front-end Receiver-Matching Network-Reusing Power Amplifier Load Inductor.

In this work, a low cost Bluetooth Low Energy (BLE) transceiver for wireless sensor network (WSN) applications, with a receiver (RX)-matching network-reusing power amplifier (PA) load inductor, is presented. In order to decrease the die area, only two inductors were used in this work. Besides the one used in the voltage control oscillator (VCO), the PA load inductor was reused as the RX impedance matching component in the front-end. Proper controls have been applied to achieve high transmitter (TX) input impedance when the transceiver is in the receiving mode, and vice versa. This allows the TRX-switch/matching network integration without significant performance degradation. The RX adopted a low-IF structure and integrated a single-ended low noise amplifier (LNA), a current bleeding mixer, a 4th complex filter and a delta-sigma continuous time (CT) analog-to-digital converter (ADC). The TX employed a two-point PLL-based architecture with a non-linear PA. The RX achieved a sensitivity of -93 dBm and consumes 9.7 mW, while the TX achieved a 2.97% error vector magnitude (EVM) with 9.4 mW at 0 dBm output power. This design was fabricated in a 0.11 µm complementary metal oxide semiconductor (CMOS) technology and the front-end circuit only occupies 0.24 mm². The measurement results verify the effectiveness and applicability of the proposed BLE transceiver for WSN applications.

Highly diversified fungi are associated with the achlorophyllous orchid Gastrodia flavilabella.

BACKGROUND: Mycoheterotrophic orchids are achlorophyllous plants that obtain carbon and nutrients from their mycorrhizal fungi. They often show strong preferential association with certain fungi and may obtain nutrients from surrounding photosynthetic plants through ectomycorrhizal fungi. Gastrodia is a large genus of mycoheterotrophic orchids in Asia, but Gastrodia species' association with fungi has not been well studied. We asked two questions: (1) whether certain fungi were preferentially associated with G. flavilabella, which is an orchid in Taiwan and (2) whether fungal associations of G. flavilabella were affected by the composition of fungi in the environment. RESULTS: Using next-generation sequencing, we studied the fungal communities in the tubers of Gastrodia flavilabella and the surrounding soil. We found (1) highly diversified fungi in the G. flavilabella tubers, (2) that Mycena species were the predominant fungi in the tubers but minor in the surrounding soil, and (3) the fungal communities in the G. flavilabella tubers were clearly distinct from those in the surrounding soil. We also found that the fungal composition in soil can change quickly with distance. CONCLUSIONS: G. flavilabella was associated with many more fungi than previously thought. Among the fungi in the tuber of G. flavilabella, Mycena species were predominant, different from the previous finding that adult G. elata depends on Armillaria species for nutritional supply. Moreover, the preferential fungus association of G. flavilabella was not significantly influenced by the composition of fungi in the environment.

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts.

BACKGROUND: Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants. RESULTS: Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting. CONCLUSION: Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.