ABL1-mediated phosphorylation promotes FOXM1-related tumorigenicity by Increasing FOXM1 stability.

The transcription factor FOXM1, which plays critical roles in cell cycle progression and tumorigenesis, is highly expressed in rapidly proliferating cells and various tumor tissues, and high FOXM1 expression is related to a poor prognosis. However, the mechanism responsible for FOXM1 dysregulation is not fully understood. Here, we show that ABL1, a nonreceptor tyrosine kinase, contributes to the high expression of FOXM1 and FOXM1-dependent tumor development. Mechanistically, ABL1 directly binds FOXM1 and mediates FOXM1 phosphorylation at multiple tyrosine (Y) residues. Among these phospho-Y sites, pY575 is indispensable for FOXM1 stability as phosphorylation at this site protects FOXM1 from ubiquitin-proteasomal degradation. The interaction of FOXM1 with CDH1, a coactivator of the E3 ubiquitin ligase anaphase-promoting complex/cyclosome (APC/C), which is responsible for FOXM1 degradation, is significantly inhibited by Y575 phosphorylation. The phospho-deficient FOXM1(Y575F) mutant exhibited increased ubiquitination, a shortened half-life, and consequently a substantially decreased abundance. Compared to wild-type cells, a homozygous Cr-Y575F cell line expressing endogenous FOXM1(Y575F) that was generated by CRISPR/Cas9 showed obviously delayed mitosis progression, impeded colony formation and inhibited xenotransplanted tumor growth. Overall, our study demonstrates that ABL1 kinase is involved in high FOXM1 expression, providing clear evidence that ABL1 may act as a therapeutic target for the treatment of tumors with high FOXM1 expression.

Detecting genetic gain and loss events in terms of protein domain: Method and implementation.

Continuous gain and loss of genes are the primary driving forces of bacterial evolution and environmental adaptation. Studying bacterial evolution in terms of protein domain, which is the fundamental function and evolutionary unit of proteins, can provide a more comprehensive understanding of bacterial differentiation and phenotypic adaptation processes. Therefore, we proposed a phylogenetic tree-based method for detecting genetic gain and loss events in terms of protein domains. Specifically, the method focuses on a single domain to trace its evolution process or on multiple domains to investigate their co-evolution principles. This novel method was validated using 122 Shigella isolates. We found that the loss of a significant number of domains was likely the main driving force behind the evolution of Shigella, which could reduce energy expenditure and preserve only the most essential functions. Additionally, we observed that simultaneously gained and lost domains were often functionally related, which can facilitate and accelerate phenotypic evolutionary adaptation to the environment. All results obtained using our method agree with those of previous studies, which validates our proposed method.

Comparative genomic analysis provides insights into the genetic diversity and pathogenicity of the genus Brucella.

Some Brucella spp. are important pathogens. According to the latest prokaryotic taxonomy, the Brucella genus consists of facultative intracellular parasitic Brucella species and extracellular opportunistic or environmental Brucella species. Intracellular Brucella species include classical and nonclassical types, with different species generally exhibiting host preferences. Some classical intracellular Brucella species can cause zoonotic brucellosis, including B. melitensis, B. abortus, B. suis, and B. canis. Extracellular Brucella species comprise opportunistic or environmental species which belonged formerly to the genus Ochrobactrum and thus nowadays renamed as for example Brucella intermedia or Brucella anthropi, which are the most frequent opportunistic human pathogens within the recently expanded genus Brucella. The cause of the diverse phenotypic characteristics of different Brucella species is still unclear. To further investigate the genetic evolutionary characteristics of the Brucella genus and elucidate the relationship between its genomic composition and prediction of phenotypic traits, we collected the genomic data of Brucella from the NCBI Genome database and conducted a comparative genomics study. We found that classical and nonclassical intracellular Brucella species and extracellular Brucella species exhibited differences in phylogenetic relationships, horizontal gene transfer and distribution patterns of mobile genetic elements, virulence factor genes, and antibiotic resistance genes, showing the close relationship between the genetic variations and prediction of phenotypic traits of different Brucella species. Furthermore, we found significant differences in horizontal gene transfer and the distribution patterns of mobile genetic elements, virulence factor genes, and antibiotic resistance genes between the two chromosomes of Brucella, indicating that the two chromosomes had distinct dynamics and plasticity and played different roles in the survival and evolution of Brucella. These findings provide new directions for exploring the genetic evolutionary characteristics of the Brucella genus and could offer new clues to elucidate the factors influencing the phenotypic diversity of the Brucella genus.

Prediction of Potential Suitable Areas and Priority Protection for Cupressus gigantea on the Tibetan Plateau.

Cupressus gigantea (C. gigantea) is an endemic endangered species on the Tibetan Plateau; its potential suitable areas and priority protection in the context of global climate change remain poorly predicted. This study utilized Biomod2 and Marxan to assess the potential suitable areas and priority protection for C. gigantea. Our study revealed that the suitable areas of C. gigantea were concentrated in the southeastern Tibetan Plateau, with the center in Lang County. Temperature was identified as a crucial environmental factor influencing the distribution of C. gigantea. Over the coming decades, the suitable range of C. gigantea expanded modestly, while its overall distribution remained relatively stable. Moreover, the center of the highly suitable areas tended to migrate towards Milin County in the northeast. Presently, significant areas for improvement are needed to establish protected areas for C. gigantea. The most feasible priority protected areas were located between the Lang and Milin counties in Tibet, which have more concentrated and undisturbed habitats. These results provide scientific guidance for the conservation and planning of C. gigantea, contributing to the stability and sustainability of ecosystems.

Spatio-temporal spread and evolution of Lassa virus in West Africa.

BACKGROUND: Lassa fever is a hemorrhagic disease caused by Lassa virus (LASV), which has been classified by the World Health Organization as one of the top infectious diseases requiring prioritized research. Previous studies have provided insights into the classification and geographic characteristics of LASV lineages. However, the factor of the distribution and evolution characteristics and phylodynamics of the virus was still limited. METHODS: To enhance comprehensive understanding of LASV, we employed phylogenetic analysis, reassortment and recombination detection, and variation evaluation utilizing publicly available viral genome sequences. RESULTS: The results showed the estimated the root of time of the most recent common ancestor (TMRCA) for large (L) segment was approximately 634 (95% HPD: [385879]), whereas the TMRCA for small (S) segment was around 1224 (95% HPD: [10301401]). LASV primarily spread from east to west in West Africa through two routes, and in route 2, the virus independently spread to surrounding countries through Liberia, resulting in a wider spread of LASV. From 1969 to 2018, the effective population size experienced two significant increased, indicating the enhanced genetic diversity of LASV. We also found the evolution rate of L segment was faster than S segment, further results showed zinc-binding protein had the fastest evolution rate. Reassortment events were detected in multiple lineages including sub-lineage IIg, while recombination events were observed within lineage V. Significant amino acid changes in the glycoprotein precursor of LASV were identified, demonstrating sequence diversity among lineages in LASV. CONCLUSION: This study comprehensively elucidated the transmission and evolution of LASV in West Africa, providing detailed insights into reassortment events, recombination events, and amino acid variations.

Predicting environmental risks of pharmaceutical residues by wastewater surveillance: An analysis based on pharmaceutical sales and their excretion data.

Pharmaceutical residues are increasingly becoming a significant source of environmental water pollution and ecological risk. This study, leveraging official national pharmaceutical sales statistics, predicts the environmental concentrations of 33 typical pharmaceuticals in the Tianjin area. The results show that 52 % of the drugs have a PEC/MEC (Predicted Environmental Concentration/Measured Environmental Concentration) ratio within the acceptable range of 0.5-2, including atenolol (1.21), carbamazepine (1.22), and sulfamethoxazole (0.91). Among the selected drugs, tetracycline, ciprofloxacin, and acetaminophen had the highest predicted concentrations. The EPI (Estimation Programs Interface) biodegradation model, a tool from the US Environmental Protection Agency, is used to predict the removal efficiency of compounds in wastewater treatment plants. The results indicate that the EPI predictions are acceptable for macrolide antibiotics and ß-blockers, with removal rates of roxithromycin, spiramycin, acetaminophen, and carbamazepine being 14.1 %, 61.2 %, 75.1 %, and 44.5 %, respectively. However, the model proved to be less effective for fluoroquinolone antibiotics. The ECOSAR (Ecological Structure-Activity Relationships) model was used to supplement the assessment of the potential impacts of drugs on aquatic ecosystems, further refining the analysis of pharmaceutical environmental risks. By combining the concentration and detection frequency of pharmaceutical wastewater, this study identified 9 drugs with significant toxicological risks and marked another 24 drugs as substances of potential concern. Additionally, this study provides data support for addressing pharmaceutical residues of priority concern in subsequent research.

Monkeypox: Can we count on the current smallpox immunization?

Two vaccines ACAM 2000 and JYNNEOS have obtained approval from the Food and Drug Administration as preventive measures against monkeypox, contributing significantly to the management of the monkeypox epidemic. Nonetheless, research has demonstrated that smallpox vaccination offers approximately 88.8% protection against monkeypox, while immunization with these vaccines generates relatively low levels of neutralizing antibodies. In this work, we performed a comprehensive comparison of antigens between the 2022-2023 monkeypox strains and the smallpox vaccine strains. Our analysis has revealed considerable amino acid changes in all 27 antigens, including core and envelope proteins. Amino acid substitutions within B cell epitopes were observed in 26 of these antigens, with at least half of the antigen substitutions occurring within B cell epitopes in 20 out of the 26 antigens analyzed. These findings may raise potential concerns regarding the efficacy of these vaccines.

Evolution of the ability to evade host innate immune defense by Talaromyces marneffei.

Talaromyces (Penicillium) marneffei is an intracellular pathogenic fungus. Some strains of this fungus have been misidentified due to the similarity between Talaromyces and Penicillium. T. marneffei has mainly been found to afflict immunocompromised individuals, causing respiratory, skin, and systemic mycosis. Mp1p is a key virulence factor that can help T. marneffei evade clearance by the normally functioning immune system. Understanding how novel functions arise is an intriguing question in many fields of biology. Mp1p has two homologous domains (Mp1p-LBD1 and Mp1p-LBD2). Sequence similarity searches with Mp1p-LBD sequences revealed Mp1p homologs in many other pathogenic fungi. Integrated information on the taxonomic distribution, phylogenetic relationships, and sequence similarity of Mp1p domains revealed that the ancestor of Mp1p-LBDs was acquired through horizontal gene transfer (HGT). Additional evidence revealed that Mp1p homologs have undergone extensive gene duplications in T. marneffei. Mp1p might be a result of gene fusion following gene duplication. Furthermore, we propose a new method for identifying Talaromyces and identify 4 strains with misclassification errors. Our results characterize the evolutionary mechanism of T. marneffei evasion of host innate immune defense and clearly demonstrate the role of gene duplication and HGT in the evolution of host immune escape by T. marneffei.

Distinct roles of the two SEC scaffold proteins, AFF1 and AFF4, in regulating RNA Pol II transcription elongation.

The super elongation complex (SEC) containing P-TEFb plays a critical role in regulating transcription elongation. AFF1 and AFF4, members of the AF4/FMR2 family, act as central scaffold proteins of SEC and are associated with various human diseases. However, their precise roles in transcriptional control remain unclear. We here reveal differences in the genomic distribution patterns of AFF1 and AFF4 around transcription start sites (TSSs). AFF1 mainly binds upstream of the TSSs, while AFF4 is enriched downstream of the TSSs. Notably, disruption of AFF4 results in slow elongation and early termination in a subset of AFF4 bound active genes, whereas AFF1 deletion leads to fast elongation and transcriptional readthrough in the same gene subset. Additionally, AFF1 knockdown increases AFF4 levels at chromatin, and vice versa. In summary, these findings demonstrate that AFF1 and AFF4 function antagonistically to regulate Pol II transcription.

Genetic characteristics of classical astroviruses in Shenzhen, China, 2016-2019.

Human astrovirus (HAstV) is a single-stranded, positive-sense RNA virus and is the leading cause of viral gastroenteritis. However, despite its prevalence, astroviruses still remain one of the least studied enteroviruses. In this study, we sequenced 11 classical astrovirus strains from clinical samples collected in Shenzhen, China from 2016 to 2019, analyzed their genetic characteristics, and deposited them into GenBank. We conducted phylogenetic analysis using IQ-TREE software, with references to astrovirus sequences worldwide. The phylogeographic analysis was performed using the Bayesian Evolutionary Analysis Sampling Trees program, through Bayesian Markov Chain Monte Carlo sampling. We also conducted recombination analysis with the Recombination Detection Program. The newly sequenced strains were categorized as HAstV genotype 1, which is the predominant genotype in Shenzhen. Phylogeographic reconstruction indicated that HAstV-1 may have migrated from the United States to China, followed by frequent transmission between China and Japan. The recombination analysis revealed recombination events within and across genotypes, and identified a recombination-prone region that produced relatively uniform recombination breakpoints and fragment lengths. The genetic analysis of HAstV strains in Shenzhen addresses the current lack of astrovirus data in the region of Shenzhen and provides key insights to the evolution and transmission of astroviruses worldwide. These findings highlight the importance of improving surveillance of astroviruses.

Characterization and Diversity of Klebsiella pneumoniae Prophages.

Klebsiella pneumoniae is a common human commensal and opportunistic pathogen. In recent years, the clinical isolation and resistance rates of K. pneumoniae have shown a yearly increase, leading to a special interest in mobile genetic elements. Prophages are a representative class of mobile genetic elements that can carry host-friendly genes, transfer horizontally between strains, and coevolve with the host's genome. In this study, we identified 15,946 prophages from the genomes of 1437 fully assembled K. pneumoniae deposited in the NCBI database, with 9755 prophages on chromosomes and 6191 prophages on plasmids. We found prophages to be notably diverse and widely disseminated in the K. pneumoniae genomes. The K. pneumoniae prophages encoded multiple putative virulence factors and antibiotic resistance genes. The comparison of strain types with prophage types suggests that the two may be related. The differences in GC content between the same type of prophages and the genomic region in which they were located indicates the alien properties of the prophages. The overall distribution of GC content suggests that prophages integrated on chromosomes and plasmids may have different evolutionary characteristics. These results suggest a high prevalence of prophages in the K. pneumoniae genome and highlight the effect of prophages on strain characterization.

Evaluating the effect of SARS-CoV-2 spike mutations with a linear doubly robust learner.

Driven by various mutations on the viral Spike protein, diverse variants of SARS-CoV-2 have emerged and prevailed repeatedly, significantly prolonging the pandemic. This phenomenon necessitates the identification of key Spike mutations for fitness enhancement. To address the need, this manuscript formulates a well-defined framework of causal inference methods for evaluating and identifying key Spike mutations to the viral fitness of SARS-CoV-2. In the context of large-scale genomes of SARS-CoV-2, it estimates the statistical contribution of mutations to viral fitness across lineages and therefore identifies important mutations. Further, identified key mutations are validated by computational methods to possess functional effects, including Spike stability, receptor-binding affinity, and potential for immune escape. Based on the effect score of each mutation, individual key fitness-enhancing mutations such as D614G and T478K are identified and studied. From individual mutations to protein domains, this paper recognizes key protein regions on the Spike protein, including the receptor-binding domain and the N-terminal domain. This research even makes further efforts to investigate viral fitness via mutational effect scores, allowing us to compute the fitness score of different SARS-CoV-2 strains and predict their transmission capacity based solely on their viral sequence. This prediction of viral fitness has been validated using BA.2.12.1, which is not used for regression training but well fits the prediction. To the best of our knowledge, this is the first research to apply causal inference models to mutational analysis on large-scale genomes of SARS-CoV-2. Our findings produce innovative and systematic insights into SARS-CoV-2 and promotes functional studies of its key mutations, serving as reliable guidance about mutations of interest.

A humanized anti-human adenovirus 55 monoclonal antibody with good neutralization ability.

Background: Human adenovirus type 55 (HAdV55) has a re-emerged as pathogen causing an acute respiratory disease presenting as a severe lower respiratory illness that can cause death. To date, there is no HAdV55 vaccine or treatment available for general use. Methods: Herein, a monoclonal antibody specific for HAdV55, mAb 9-8, was isolated from an scFv-phage display library derived from mice immunized with the purified inactived-HAdV55 virions. By using ELISA and a virus micro-neutralization assay, we evaluated the binding and neutralizing activity of mAb 9-8 following humanization. Western blotting analysis and antigen-antibody molecular docking analysis were used to identify the antigenic epitopes that the humanized monoclonal antibody 9-8-h2 recognized. After that, their thermal stability was determined. Results: MAb 9-8 showed potent neutralization activity against HAdV55. After humanization, the humanized neutralizing monoclonal antibody (9-8-h2) was identified to neutralize HAdV55 infection with an IC50 of 0.6050 nM. The mAb 9-8-h2 recognized HAdV55 and HAdV7 virus particles, but not HAdV4 particles. Although mAb 9-8-h2 could recognize HAdV7, it could not neutralize HAdV7. Furthermore, mAb 9-8-h2 recognized a conformational neutralization epitope of the fiber protein and the crucial amino acid residues (Arg 288, Asp 157, and Asn 200) were identified. MAb 9-8-h2 also showed favorable general physicochemical properties, including good thermostability and pH stability. Conclusions: Overall, mAb 9-8-h2 might be a promising molecule for the prevention and treatment of HAdV55.

Downregulation of miRNA-155-5p contributes to the adipogenic activity of 2-ethylhexyl diphenyl phosphate in 3T3-L1 preadipocytes.

2-Ethylhexyl diphenyl phosphate (EHDPP) is a commonly used organophosphorus flame retardant and food packaging material. Because of its high lipophilic and bioaccumulative properties, adipocytes are the primary target of EHDPP. However, the toxicity of EHDPP on preadipocytes and the potential mechanism have not been fully elucidated. MicroRNAs (miRNAs) are thought to be an important mediator that contribute to the toxicity of environmental contaminants. To identify the miRNAs specifically responsible for EHDPP exposure and their role in EGDPP's toxicity in preadipocytes, the adipogenic effects and miRNA expression profiling were performed on 3T3-L1 preadipocytes exposed to EHDPP. EHDPP at concentrations of 1-10 µM promoted adipocyte differentiation, as evidenced by lipid staining, triglyceride content, and expression of adipogenesis markers. MiRNA-seq analysis revealed that 7 differentially expressed miRNAs were recognized under EHDPP exposure, with miR-155-5p being the top down-regulated miRNA. Quantitative reverse transcription PCR (RT-qPCR) analysis showed that miR-155-5p level fell sharply during the first 2 days and continued to fall dose-dependently throughout the EHDPP exposure period. MiR-155-5p inhibition promotes adipocyte differentiation, whereas its overexpression counteracted EHDPP-induced adipogenesis. Luciferase reporter assay identified CCAAT/enhancer-binding protein beta (C/EBPß) as a target of miR-155-5p in 3T3-L1 preadipocytes in response to EHDPP. Taken together, EHDPP exposure down-regulated miR-155-5p, which then increased C/EBPß and peroxisome proliferator-activated receptor γ (PPARγ) expression and promoted adipogenesis in preadipocytes.

Triphenyl phosphate proved more potent than its metabolite diphenyl phosphate in inducing hepatic insulin resistance through endoplasmic reticulum stress.

Triphenyl phosphate (TPHP) is a widely used flame retardant and plasticizer and has been detected extensively in environmental media, wildlife and human bodies. Several epidemiological and animal studies have revealed that TPHP exposure is positively associated with glucose homeostasis disruption and diabetes. However, the effects of TPHP on hepatic glucose homeostasis and the underlying mechanisms remain unclear. The present work aimed to investigate the cytotoxicity and glucose metabolism disruption of TPHP and its metabolite diphenyl phosphate (DPHP) within hepatocytes. The cell viability assay undertaken on human normal liver (L02) cells showed that TPHP exhibited more potent hepatotoxicity than DPHP. RNA sequencing (RNA-seq) data showed that TPHP and DPHP presented different modes of toxic action. Insulin resistance is one of the predominant toxicities for TPHP, but not for DPHP. The insulin-stimulated glucose uptake and glycogen synthesis were impaired by TPHP, while DPHP exhibited no significant impairment on these factors. TPHP exposure induced endoplasmic reticulum (ER) stress, and the ER stress antagonist 4-PBA restored the impairment of insulin-stimulated glucose uptake and glycogen synthesis induced by TPHP. TPHP could also induce liver ER stress and insulin resistance in mice. Taken together, the results suggested that TPHP induces more potent insulin resistance through ER stress than its metabolite DPHP.

Comparative genomics reveals new insights into the evolution of the IncA and IncC family of plasmids.

Incompatibility groups IncA and IncC plasmids are of great concern due to their ability to disseminate antibiotic resistance in bacteria via conjugative transfer. A deep understanding of their genomic structures and evolutionary characteristics is of great significance for improving our knowledge about its multidrug-resistance evolution and dissemination. However, current knowledge of their backbone structure, features of core functional modules and the characteristics of variable regions is based on a few plasmids, which highlights the need for a comprehensive systematic study. The present study thoroughly compared and analysed 678 IncA and IncC plasmid genomes. We found that their core functional genes were occasionally deficient and sometimes existed as multiple functional copies/multiple families, which resulted in much diversity. The phylogeny of 13 core functional genes corresponded well to the plasmid subtypes. The conjugative transfer system gained diverse complexity and exhibited many previously unnoticed types with multiple combinations. The insertion of mobile genetic elements (MGEs) in plasmids varied between types and was present in 4 insertion spots in different types of plasmids with certain types of transposons, integrons and insertion sequences. The impact of gene duplication, deletion, the insertion of MGEs, genome rearrangement and recombination resulted in the complex dynamic variable backbone of IncA and IncC plasmids. And IncA and IncC plasmids were more complex than their closest relative SXT/R391 integrative conjugative elements (ICEs), which included nearly all of the diversity of SXT/R391 in key systems. Our work demonstrated a global and systematic view of the IncA and IncC plasmids and provides many new insights into their genome evolution.

A Vaccine with Multiple Receptor-Binding Domain Subunit Mutations Induces Broad-Spectrum Immune Response against SARS-CoV-2 Variants of Concern.

With the emergence of more variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the immune evasion of these variants from existing vaccines, the development of broad-spectrum vaccines is urgently needed. In this study, we designed a novel SARS-CoV-2 receptor-binding domain (RBD) subunit (RBD5m) by integrating five important mutations from SARS-CoV-2 variants of concern (VOCs). The neutralization activities of antibodies induced by the RBD5m candidate vaccine are more balanced and effective for neutralizing different SARS-CoV-2 VOCs in comparison with those induced by the SARS-CoV-2 prototype strain RBD. Our results suggest that the RBD5m vaccine is a good broad-spectrum vaccine candidate able to prevent disease from several different SARS-CoV-2 VOCs.

A Standardized Framework for Better Understanding of Phenotypic Differences within Bacterial Phyla Based on Protein Domain.

We propose a standardized framework to classify target species based on their protein domains, which can be utilized in different contexts, like eukaryotes and prokaryotes. In this study, by applying the framework to the bacterial kingdom as an implementation example and comparing the results with the current taxonomy standards at the phylum level, we came to the conclusion that the sequence of domains rather than the content of domains in a protein and the presence of one domain rather than the number of occurrences of one domain play more important roles in deciding bacterial phenotypes as well as matching the current taxonomy. In addition, the comparison also helps us to better focus on the species that conflict with the current phylum category, as well as to further investigate their phenotypic or genotypic differences. IMPORTANCE A 3-step framework was designed which can be applied to clustering species based on their protein domains, and different candidate models are proposed in each step for better adaptation of various scenarios. We show its implementation for the bacterial kingdom as an example, which helps us to find the most appropriate model combination that will best reflect the relationship between domains and phenotypes in this context. In addition, identifying species that are distant in the results but should be closely related phylogenetically can help us to focus on the mismatch for better understanding of their key phenotypic or genotypic differences.

Comparative Genomic Analysis of Vibrio cincinnatiensis Provides Insights into Genetic Diversity, Evolutionary Dynamics, and Pathogenic Traits of the Species.

Vibrio cincinnatiensis is a poorly understood pathogenic Vibrio species, and the underlying mechanisms of its genetic diversity, genomic plasticity, evolutionary dynamics, and pathogenicity have not yet been comprehensively investigated. Here, a comparative genomic analysis of V. cincinnatiensis was constructed. The open pan-genome with a flexible gene repertoire exhibited genetic diversity. The genomic plasticity and stability were characterized by the determinations of diverse mobile genetic elements (MGEs) and barriers to horizontal gene transfer (HGT), respectively. Evolutionary divergences were exhibited by the difference in functional enrichment and selective pressure between the different components of the pan-genome. The evolution on the Chr I and Chr II core genomes was mainly driven by purifying selection. Predicted essential genes in V. cincinnatiensis were mainly found in the core gene families on Chr I and were subject to stronger evolutionary constraints. We identified diverse virulence-related elements, including the gene clusters involved in encoding flagella, secretion systems, several pili, and scattered virulence genes. Our results indicated the pathogenic potential of V. cincinnatiensis and highlighted that HGT events from other Vibrio species promoted pathogenicity. This pan-genome study provides comprehensive insights into this poorly understood species from the genomic perspective.