A high-performance computational workflow to accelerate GATK SNP detection across a 25-genome dataset.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) are the most widely used form of molecular genetic variation studies. As reference genomes and resequencing data sets expand exponentially, tools must be in place to call SNPs at a similar pace. The genome analysis toolkit (GATK) is one of the most widely used SNP calling software tools publicly available, but unfortunately, high-performance computing versions of this tool have yet to become widely available and affordable. RESULTS: Here we report an open-source high-performance computing genome variant calling workflow (HPC-GVCW) for GATK that can run on multiple computing platforms from supercomputers to desktop machines. We benchmarked HPC-GVCW on multiple crop species for performance and accuracy with comparable results with previously published reports (using GATK alone). Finally, we used HPC-GVCW in production mode to call SNPs on a "subpopulation aware" 16-genome rice reference panel with ~ 3000 resequenced rice accessions. The entire process took ~ 16 weeks and resulted in the identification of an average of 27.3 M SNPs/genome and the discovery of ~ 2.3 million novel SNPs that were not present in the flagship reference genome for rice (i.e., IRGSP RefSeq). CONCLUSIONS: This study developed an open-source pipeline (HPC-GVCW) to run GATK on HPC platforms, which significantly improved the speed at which SNPs can be called. The workflow is widely applicable as demonstrated successfully for four major crop species with genomes ranging in size from 400 Mb to 2.4 Gb. Using HPC-GVCW in production mode to call SNPs on a 25 multi-crop-reference genome data set produced over 1.1 billion SNPs that were publicly released for functional and breeding studies. For rice, many novel SNPs were identified and were found to reside within genes and open chromatin regions that are predicted to have functional consequences. Combined, our results demonstrate the usefulness of combining a high-performance SNP calling architecture solution with a subpopulation-aware reference genome panel for rapid SNP discovery and public deployment.

sRNAdeep: a novel tool for bacterial sRNA prediction based on DistilBERT encoding mode and deep learning algorithms.

BACKGROUND: Bacterial small regulatory RNA (sRNA) plays a crucial role in cell metabolism and could be used as a new potential drug target in the treatment of pathogen-induced disease. However, experimental methods for identifying sRNAs still require a large investment of human and material resources. METHODS: In this study, we propose a novel sRNA prediction model called sRNAdeep based on the DistilBERT feature extraction and TextCNN methods. The sRNA and non-sRNA sequences of bacteria were considered as sentences and then fed into a composite model consisting of deep learning models to evaluate classification performance. RESULTS: By filtering sRNAs from BSRD database, we obtained a validation dataset comprised of 2438 positive and 4730 negative samples. The benchmark experiments showed that sRNAdeep displayed better performance in the various indexes compared to previous sRNA prediction tools. By applying our tool to Mycobacterium tuberculosis (MTB) genome, we have identified 21 sRNAs within the intergenic and intron regions. A set of 272 targeted genes regulated by these sRNAs were also captured in MTB. The coding proteins of two genes (lysX and icd1) are implicated in drug response, with significant active sites related to drug resistance mechanisms of MTB. CONCLUSION: In conclusion, our newly developed sRNAdeep can help researchers identify bacterial sRNAs more precisely and can be freely available from https://github.com/pyajagod/sRNAdeep.git .

Convergent dwarfism consequences of minipigs under independent artificial selections.

BACKGROUND: Currently, diverse minipigs have acquired a common dwarfism phenotype through independent artificial selections. Characterizing the population and genetic diversity in minipigs is important to unveil genetic mechanisms regulating their body sizes and effects of independent artificial selections on those genetic mechanisms. However, full understanding for the genetic mechanisms and phenotypic consequences in minipigs still lag behind. RESULTS: Here, using whole genome sequencing data of 41 pig breeds, including eight minipigs, we identified a large genomic diversity in a minipig population compared to other pig populations in terms of population structure, demographic signatures, and selective signatures. Selective signatures reveal diverse biological mechanisms related to body size in minipigs. We also found evidence for neural development mechanism as a minipig-specific body size regulator. Interestingly, selection signatures within those mechanisms containing neural development are also highly different among minipig breeds. Despite those large genetic variances, PLAG1, CHM, and ESR1 are candidate key genes regulating body size which experience different differentiation directions in different pig populations. CONCLUSIONS: These findings present large variances of genetic structures, demographic signatures, and selective signatures in the minipig population. They also highlight how different artificial selections with large genomic diversity have shaped the convergent dwarfism.

First report and genomic characterization of Escherichia coli O111:H12 serotype from raw mussels in Türkiye.

BACKGROUND: This study aimed to assess the prevalence and genomic characteristics of Shiga-toxigenic (STEC) and Enteroaggregative E. coli (EAEC) strains in raw mussels and ready-to-eat (RTE)-stuffed mussels, focusing on potential public health implications for identifying virulence and antimicrobial resistance genes. RESULTS: The genome sequence analysis identified the E. coli strain named 23EM as serotype O111:H12, with adhesion (fimH-54) and fumarate hydratase (fumC-11) genes. The draft genome (4.9 Mb, 50.6% GC content, 111 contigs, 4,688 genes) is available in NCBI GenBank (accession JAWXVJ000000000). The strain, classified as ST292 and CC ST10, showed high similarity to nonpathogenic E. coli MG1655 but was distinct from pathogenic strains such as EAEC and ExPEC. In silico serotyping revealed the presence of O111-antigen flippase (wzx) and H12-antigen flagellin (fliC) genes. The strain harbors an IncFII (pCoo) plasmid with 96.95% identity. PathogenFinder predicted a 92% probability of being a human pathogen, supported by 720 pathogenic protein families. CRISPR analysis identified one high-evidence sequence with nine spacers and six low-evidence sequences. Phylogenetic analysis using RAxML positioned 23EM close to nonpathogenic E. coli but distant from other pathogenic strains. Antimicrobial resistance genes across multiple classes, including macrolides, fluoroquinolones, and aminoglycosides, were identified. The strain also contains several virulence factors, such as adhesins (e.g., ECP, ELF, TIF, type IV pili), and autotransporter genes (espP, pic), highlighting its significant pathogenic potential and public health risk. CONCLUSIONS: This study highlights the ability of the detection of E. coli strains harboring virulence and antimicrobial resistance genes in mussels, thus emphasizing the importance of ongoing surveillance and careful consideration of the potential risks associated with the consumption of these shellfish.

Pangenome analysis of Shewanella xiamenensis revealed important genetic traits concerning genetic diversity, pathogenicity and antibiotic resistance.

BACKGROUND: Shewanella xiamenensis, widely distributed in natural environments, has long been considered as opportunistic pathogen. Recently, significant changes in the resistance spectrum have been observed in S. xiamenensis, due to acquired antibiotic resistance genes. Therefore, a pan-genome analysis was conducted to illuminate the genomic changes in S. xiamenensis. RESULTS: Phylogenetic analysis revealed three major clusters and three singletons, among which close relationship between several strains was discovered, regardless of their host and niches. The "open" genomes with diversity of accessory and strain-specific genomes took advantage towards diversity environments. The purifying selection pressure was the main force on genome evolution, especially in conservative genes. Only 53 gene families were under positive selection pressure. Phenotypic resistance analysis revealed 21 strains were classified as multi-drug resistance (MDR). Ten types of antibiotic resistance genes and two heavy metal resistance operons were discovered in S. xiamenensis. Mobile genetic elements and horizontal gene transfer increased genome diversity and were closely related to MDR strains. S. xiamenensis carried a variety of virulence genes and macromolecular secretion systems, indicating their important roles in pathogenicity and adaptability. Type IV secretion system was discovered in 15 genomes with various sequence structures, indicating it was originated from different donors through horizontal gene transfer. CONCLUSIONS: This study provided with a detailed insight into the changes in the pan-genome of S. xiamenensis, highlighting its capability to acquire new mobile genetic elements and resistance genes for its adaptation to environment and pathogenicity to human and animals.

Improving nitrofurantoin resistance prediction in Escherichia coli from whole-genome sequence by integrating NfsA/B enzyme assays.

Nitrofurantoin resistance in Escherichia coli is primarily caused by mutations damaging two enzymes, NfsA and NfsB. Studies based on small isolate collections with defined nitrofurantoin MICs have found significant random genetic drift in nfsA and nfsB, making it extremely difficult to predict nitrofurantoin resistance from whole-genome sequence (WGS) where both genes are not obviously disrupted by nonsense or frameshift mutations or insertional inactivation. Here, we report a WGS survey of 200 oqxAB-negative E. coli from community urine samples, of which 34 were nitrofurantoin resistant. We characterized individual non-synonymous mutations seen in nfsA and nfsB among this collection using complementation cloning and NfsA/B enzyme assays in cell extracts. We definitively identified R203C, H11Y, W212R, A112E, and A112T in NfsA and R121C, Q142H, F84S, P163H, W46R, K57E, and V191G in NfsB as amino acid substitutions that reduce enzyme activity sufficiently to cause resistance. In contrast, E58D, I117T, K141E, L157F, A172S, G187D, and A188V in NfsA and G66D, M75I, V93A, and A174E in NfsB are functionally silent in this context. We identified that 9/166 (5.4%) nitrofurantoin-susceptible isolates were "pre-resistant," defined as having loss of function mutations in nfsA or nfsB. Finally, using NfsA/B enzyme assays and proteomics, we demonstrated that 9/34 (26.5%) ribE wild-type nitrofurantoin-resistant isolates also carried functionally wild-type nfsB or nfsB/nfsA. In these cases, NfsA/B activity was reduced through downregulated gene expression. Our biological understanding of nitrofurantoin resistance is greatly improved by this analysis but is still insufficient to allow its reliable prediction from WGS data.

Hybrid-genome sequence analysis of Enterobacter cloacae FACU and morphological characterization: insights into a highly arsenic-resistant strain.

Many organisms have adapted to survive in environments with high levels of arsenic (As), a naturally occurring metalloid with various oxidation states and a common element in human activities. These organisms employ diverse mechanisms to resist the harmful effects of arsenic compounds. Ten arsenic-resistant bacteria were isolated from contaminated wastewater in this study. The most efficient bacterial isolate able to resist 15,000 ppm Na2HAsO4·7H2O was identified using the 16S rRNA gene and whole genome analysis as Enterobacter cloacae FACU. The arsenic E. cloacae FACU biosorption capability was analyzed. To further unravel the genetic determinants of As stress resistance, the whole genome sequence of E. cloacae FACU was performed. The FACU complete genome sequence consists of one chromosome (5.7 Mb) and two plasmids, pENCL 1 and pENCL 2 (755,058 and 1155666 bp, respectively). 7152 CDSs were identified in the E. cloacae FACU genome. The genome consists of 130 genes for tRNA and 21 for rRNAs. The average G + C content was found to be 54%. Sequencing analysis annotated 58 genes related to resistance to many heavy metals, including 16 genes involved in arsenic efflux transporter and arsenic reduction (five arsRDABC genes) and 42 genes related to lead, zinc, mercury, nickel, silver, copper, cadmium and chromium in FACU. Scanning electron microscopy (SEM) confirmed the difference between the morphological responses of the As-treated FACU compared to the control strain. The study highlights the genes involved in the mechanism of As stress resistance, metabolic pathways, and potential activity of E. cloacae FACU at the genetic level.

Phylogenetic Analysis and Comparative Genomics of Brucella abortus and Brucella melitensis Strains in Egypt.

Brucellosis is a notifiable disease induced by a facultative intracellular Brucella pathogen. In this study, eight Brucella abortus and eighteen Brucella melitensis strains from Egypt were annotated and compared with RB51 and REV1 vaccines respectively. RAST toolkit in the BV-BRC server was used for annotation, revealing genome length of 3,250,377 bp and 3,285,803 bp, 3289 and 3323 CDS, 48 and 49 tRNA genes, the same number of rRNA (3) genes, 583 and 586 hypothetical proteins, 2697 and 2726 functional proteins for B. abortus and B. melitensis respectively. B. abortus strains exhibit a similar number of candidate genes, while B. melitensis strains showed some differences, especially in the SRR19520422 Faiyum strain. Also, B. melitensis clarified differences in antimicrobial resistance genes (KatG, FabL, MtrA, MtrB, OxyR, and VanO-type) in SRR19520319 Faiyum and (Erm C and Tet K) in SRR19520422 Faiyum strain. Additionally, the whole genome phylogeny analysis proved that all B. abortus strains were related to vaccinated animals and all B. melitensis strains of Menoufia clustered together and closely related to Gharbia, Dameitta, and Kafr Elshiek. The Bowtie2 tool identified 338 (eight B. abortus) and 4271 (eighteen B. melitensis) single nucleotide polymorphisms (SNPs) along the genomes. These variants had been annotated according to type and impact. Moreover, thirty candidate genes were predicted and submitted at GenBank (24 in B. abortus) and (6 in B. melitensis). This study contributes significant insights into genetic variation, virulence factors, and vaccine-related associations of Brucella pathogens, enhancing our knowledge of brucellosis epidemiology and evolution in Egypt.

An analysis pipeline for understanding 6-thioguanine effects on a mouse tumour genome.

Mouse tumour models are extensively used as a pre-clinical research tool in the field of oncology, playing an important role in anticancer drugs discovery. Accordingly, in cancer genomics research, the demand for next-generation sequencing (NGS) is increasing, and consequently, the need for data analysis pipelines is likewise growing. Most NGS data analysis solutions to date do not support mouse data or require highly specific configuration for their use. Here, we present a genome analysis pipeline for mouse tumour NGS data including the whole-genome sequence (WGS) data analysis flow for somatic variant discovery, and the RNA-seq data flow for differential expression, functional analysis and neoantigen prediction. The pipeline is based on standards and best practices and integrates mouse genome references and annotations. In a recent study, the pipeline was applied to demonstrate the efficacy of low dose 6-thioguanine (6TG) treatment on low-mutation melanoma in a pre-clinical mouse model. Here, we further this study and describe in detail the pipeline and the results obtained in terms of tumour mutational burden (TMB) and number of predicted neoantigens, and correlate these with 6TG effects on tumour volume. Our pipeline was expanded to include a neoantigen analysis, resulting in neopeptide prediction and MHC class I antigen presentation evaluation. We observed that the number of predicted neoepitopes were more accurate indicators of tumour immune control than TMB. In conclusion, this study demonstrates the usability of the proposed pipeline, and suggests it could be an essential robust genome analysis platform for future mouse genomic analysis.

Deciphering the possible role of RNA-helicase genes mechanism in response to abiotic stresses in rapeseed (Brassica napus L.).

BACKGROUND: Plants mediate several defense mechanisms to withstand abiotic stresses. Several gene families respond to stress as well as multiple transcription factors to minimize abiotic stresses without minimizing their effects on performance potential. RNA helicase (RH) is one of the foremost critical gene families that can play an influential role in tolerating abiotic stresses in plants. However, little knowledge is present about this protein family in rapeseed (canola). Here, we performed a comprehensive survey analysis of the RH protein family in rapeseed (Brassica napus L.). RESULTS: A total of 133 BnRHs genes have been discovered in this study. By phylogenetic analysis, RHs genes were divided into one main group and a subgroup. Examination of the chromosomal position of the identified genes showed that most of the genes (27%) were located on chromosome 3. All 133 identified sequences contained the main DEXDC domain, the HELICC domain, and a number of sub-domains. The results of biological process studies showed that about 17% of the proteins acted as RHs, 22% as ATP binding, and 14% as mRNA binding. Each part of the conserved motifs, communication network, and three-dimensional structure of the proteins were examined separately. The results showed that the RWC in leaf tissue decreased with higher levels of drought stress and in both root and leaf tissues sodium concentration was increased upon increased levels of salt stress treatments. The proline content were found to be increased in leaf and root with the increased level of stress treatment. Finally, the expression patterns of eight selected RHs genes that have been exposed to drought, salinity, cold, heat and cadmium stresses were investigated by qPCR. The results showed the effect of genes under stress. Examination of gene expression in the Hayola #4815 cultivar showed that all primers except primer #79 had less expression in both leaves and roots than the control level. CONCLUSIONS: New finding from the study have been presented new insights for better understanding the function and possible mechanism of RH in response to abiotic stress in rapeseed.

Isolation, characterization, and application of a novel Pseudomonas fluorescens phage vB_PF_Y1-MI in contaminated milk.

The food industry has incurred substantial losses from contamination by Pseudomonas fluorescens, emphasizing the critical importance of implementing effective control strategies. Phages are potential sterilizers due to their specific killing abilities and the difficulty bacteria face in developing resistance. However, a significant barrier to their development is the lack of diversity among phage types. In this study, we characterized a novel lytic P. fluorescens phage, named vB_PF_Y1-MI. Phage vB_PF_Y1-MI displayed a latent period of nearly 10 min and a high burst size of 1493 PFU/cell. This phage showed good activity over a wide range of temperature (up to 70 °C) and pH (3-12). The genome of phage vB_PF_Y1-MI spans 93,233 bp with a GC content of 45%. It encompasses 174 open-reading frames and 19 tRNA genes, while no lysogeny or virulence-associated genes were detected. Phylogenetic analysis positions it as a novel unassigned evolutionary lineage within the Caudoviricetes class among related dsDNA phages. Our study provides foundational insights into vB_PF_Y1-MI and emphasizes its potential as an effective biological control agent against P. fluorescens. This research offers crucial theoretical groundwork and technical support for subsequent efforts in preventing and controlling P. fluorescens contamination.

Host genetic variation at a locus near CHD1L impacts HIV sequence diversity in a South African population.

IMPORTANCE: It has been previously shown that genetic variants near CHD1L on chromosome 1 are associated with reduced HIV VL in African populations. However, the impact of these variants on viral diversity and how they restrict viral replication are unknown. We report on a regional association analysis in a South African population and show evidence of selective pressure by variants near CHD1L on HIV RT and gag. Our findings provide further insight into how genetic variability at this locus contributes to host control of HIV in a South African population.

High-quality pan-genome of Escherichia coli generated by excluding confounding and highly similar strains reveals an association between unique gene clusters and genomic islands.

The pan-genome analysis of bacteria provides detailed insight into the diversity and evolution of a bacterial population. However, the genomes involved in the pan-genome analysis should be checked carefully, as the inclusion of confounding strains would have unfavorable effects on the identification of core genes, and the highly similar strains could bias the results of the pan-genome state (open versus closed). In this study, we found that the inclusion of highly similar strains also affects the results of unique genes in pan-genome analysis, which leads to a significant underestimation of the number of unique genes in the pan-genome. Therefore, these strains should be excluded from pan-genome analysis at the early stage of data processing. Currently, tens of thousands of genomes have been sequenced for Escherichia coli, which provides an unprecedented opportunity as well as a challenge for pan-genome analysis of this classical model organism. Using the proposed strategies, a high-quality E. coli pan-genome was obtained, and the unique genes was extracted and analyzed, revealing an association between the unique gene clusters and genomic islands from a pan-genome perspective, which may facilitate the identification of genomic islands.

Benchmarking of long-read sequencing, assemblers and polishers for yeast genome.

BACKGROUND: The long reads of the third-generation sequencing significantly benefit the quality of the de novo genome assembly. However, its relatively high single-base error rate has been criticized. Currently, sequencing accuracy and throughput continue to improve, and many advanced tools are constantly emerging. PacBio HiFi sequencing and Oxford Nanopore Technologies (ONT) PromethION are two up-to-date platforms with low error rates and ultralong high-throughput reads. Therefore, it is urgently needed to select the appropriate sequencing platforms, depths and genome assembly tools for high-quality genomes in the era of explosive data production. METHODS: We performed 455 (7 assemblers with 4 polishing pipelines or without polishing on 13 subsets with different depths) and 88 (4 assemblers with or without polishing on 11 subsets with different depths) de novo assemblies of Yeast S288C on high-coverage ONT and HiFi datasets, respectively. The assembly quality was evaluated by Quality Assessment Tool (QUAST), Benchmarking Universal Single-Copy Orthologs (BUSCO) and the newly proposed Comprehensive_score (C_score). In addition, we applied four preferable pipelines to assemble the genome of nonreference yeast strains. RESULTS: The assembler plays an essential role in genome construction, especially for low-depth datasets. For ONT datasets, Flye is superior to other tools through C_score evaluation. Polishing by Pilon and Medaka improve accuracy and continuity of the preassemblies, respectively, and their combination pipeline worked well in most quality metrics. For HiFi datasets, Flye and NextDenovo performed better than other tools, and polishing is also necessary. Enough data depth is required for high-quality genome construction by ONT (>80X) and HiFi (>20X) datasets.

Increased resistance against tellurite is conferred by a mutation in the promoter region of uncommon tellurite resistance gene tehB in the ter-negative Shiga toxin-producing Escherichia coli O157:H7.

Resistance to potassium tellurite (PT) is an important indicator in isolating Shiga toxin-producing Escherichia coli (STEC) O157:H7 and other major STEC serogroups. Common resistance determinant genes are encoded in the ter gene cluster. We found an O157:H7 isolate that does not harbor ter but is resistant to PT. One nonsynonymous mutation was found in another PT resistance gene, tehA, through whole-genome sequence analyses. To elucidate the contribution of this mutation to PT resistance, complementation of tehA and the related gene tehB in isogenic strains and quantitative RT‒PCR were performed. The results indicated that the point mutation not only changed an amino acid of tehA, but also was positioned on a putative internal promoter of tehB and increased PT resistance by elevating tehB mRNA expression. Meanwhile, the amino acid change in tehA had negligible impact on the PT resistance. Comprehensive screening revealed that 2.3% of O157:H7 isolates in Japan did not harbor the ter gene cluster, but the same mutation in tehA was not found. These results suggested that PT resistance in E. coli can be enhanced through one mutational event even in ter-negative strains. IMPORTANCE: Selective agents are important for isolating Shiga toxin-producing Escherichia coli (STEC) because the undesirable growth of microflora should be inhibited. Potassium tellurite (PT) is a common selective agent for major STEC serotypes. In this study, we found a novel variant of PT resistance genes, tehAB, in STEC O157:H7. Molecular experiments clearly showed that one point mutation in a predicted internal promoter region of tehB upregulated the expression of the gene and consequently led to increased resistance to PT. Because tehAB genes are ubiquitous across E. coli, these results provide universal insight into PT resistance in this species.

Genome characterisation and comparative analysis of Schaalia dentiphila sp. nov. and its subspecies, S. dentiphila subsp. denticola subsp. nov., from the human oral cavity.

BACKGROUND: Schaalia species are primarily found among the oral microbiota of humans and other animals. They have been associated with various infections through their involvement in biofilm formation, modulation of host responses, and interaction with other microorganisms. In this study, two strains previously indicated as Actinomyces spp. were found to be novel members of the genus Schaalia based on their whole genome sequences. RESULTS: Whole-genome sequencing revealed both strains with a genome size of 2.3 Mbp and GC contents of 65.5%. Phylogenetics analysis for taxonomic placement revealed strains NCTC 9931 and C24 as distinct species within the genus Schaalia. Overall genome-relatedness indices including digital DNA-DNA hybridization (dDDH), and average nucleotide/amino acid identity (ANI/AAI) confirmed both strains as distinct species, with values below the species boundary thresholds (dDDH < 70%, and ANI and AAI < 95%) when compared to nearest type strain Schaalia odontolytica NCTC 9935 T. Pangenome and orthologous analyses highlighted their differences in gene properties and biological functions compared to existing type strains. Additionally, the identification of genomic islands (GIs) and virulence-associated factors indicated their genetic diversity and potential adaptive capabilities, as well as potential implications for human health. Notably, CRISPR-Cas systems in strain NCTC 9931 underscore its adaptive immune mechanisms compared to strain C24. CONCLUSIONS: Based on these findings, strain NCTC 9931T (= ATCC 17982T = DSM 43331T = CIP 104728T = CCUG 18309T = NCTC 14978T = CGMCC 1.90328T) represents a novel species, for which the name Schaalia dentiphila subsp. dentiphila sp. nov. subsp. nov. is proposed, while strain C24T (= NCTC 14980T = CGMCC 1.90329T) represents a distinct novel subspecies, for which the name Schaalia dentiphila subsp. denticola. subsp. nov. is proposed. This study enriches our understanding of the genomic diversity of Schaalia species and paves the way for further investigations into their roles in oral health. SIGNIFICANCE: This research reveals two Schaalia strains, NCTC 9931 T and C24T, as novel entities with distinct genomic features. Expanding the taxonomic framework of the genus Schaalia, this study offers a critical resource for probing the metabolic intricacies and resistance patterns of these bacteria. This work stands as a cornerstone for microbial taxonomy, paving the way for significant advances in clinical diagnostics.

Characterization of the genome and cell invasive phenotype of Vibrio diabolicus Cg5 isolated from mass mortality of Pacific oyster, Crassostrea gigas.

Vibrio is an important group of aquatic animal pathogens, which has been identified as the main pathogenic factor causing mass summer mortality of Crassostrea gigas in northern China. This study aims to investigate the potential pathogenic mechanisms of Vibrio Cg5 isolate in C. gigas. We sequenced and annotated the genome of Vibrio Cg5 to analyze potential virulence factors. The gentamicin protection assays were performed with C. gigas primary cells to reveal the cell-invasive behavior of Cg5. The genome analysis showed that Cg5 was a strain of human disease-associated pathogen with multiple antibiotic resistance, and four virulence factors associated with intracellular survival were present in the genome. The gentamicin protection assays showed that Cg5 could potentially invade the cells of C. gigas, indicating that Cg5 could be a facultative intracellular pathogen of C. gigas. These results provide insights into the pathogenic mechanism of V. diabolicus, an emerging pathogenic Vibrio on aquatic animals, which would be valuable in preventing and controlling diseases in oysters.

Isolation and characterization of a G9P[23] porcine rotavirus strain AHFY2022 in China.

Rotavirus group A (RVA) is a main pathogen causing diarrheal diseases in humans and animals. Various genotypes are prevalent in the Chinese pig herd. The genetic diversity of RVA lead to distinctly characteristics. In the present study, a porcine RVA strain, named AHFY2022, was successfully isolated from the small intestine tissue of piglets with severe diarrhea. The AHFY2022 strain was identified by cytopathic effects (CPE) observation, indirect immunofluorescence assay (IFA), electron microscopy (EM), high-throughput sequencing, and pathogenesis to piglets. The genomic investigation using NGS data revealed that AHFY2022 exhibited the genotypes G9-P[23]-I5-R1-C1-M1-A8-N1-T1-E1-H1, using the online platform the Bacterial and Viral Bioinformatics Resource Center (BV-BRC) (https://www.bv-brc.org/). Moreover, experimental inoculation in 5-day-old and 27-day-old piglets demonstrated that AHFY2022 caused severe diarrhea, fecal shedding, small intestinal villi damage, and colonization in all challenged piglets. Taken together, our results detailed the virological features of the porcine rotavirus G9P[23] from China, including the whole-genome sequences, genotypes, growth kinetics in MA104 cells and the pathogenicity in suckling piglets.

Microbacterium aquilitoris sp. nov. and Microbacterium gwkjiense sp. nov., isolated from beach.

Three Gram-reaction-positive bacterial strains, designated KSW-18T, KSW2-22, and KSW4-11T, were isolated from seawater, and two dried seaweed samples collected at Gwakji Beach in Jeju, Republic of Korea, respectively, and their taxonomic positions were examined by a polyphasic approach. The 16S rRNA gene phylogeny showed that strain KSW4-11T was tightly associated with Microbacterium oleivorans NBRC 103075T, while strains KSW-18T and KSW2-22 formed a distinctive subline at the base of a clade including the above two strains. The three isolates showed high sequence similarity with one another (99.7-99.9%; 1-4 nt differences) and Microbacterium oleivorans (99.8-99.9%; 1-3 nt differences). The chemotaxonomic features were typical for the genus Microbacterium; Lysine as the diagnostic diamino acid and N-glycolylated muramic acid of the peptidoglycans, the predominant menaquinones of MK-11, MK-10 and MK-12, the major fatty acids of anteiso-C15:0 and anteiso-C17:0, and the major polar lipids including diphosphatidylglycerol, phosphatidylglycerol, and two or three unidentified glycolipids. In core genome-based phylogenetic tree, strains KSW-18T and KSW2-22 were closely associated with Microbacterium oleivorans NBRC 103075T, while strain KSW4-11T formed a distinctive subline at the base of a clade including the above three strains, in contrast to the 16S rRNA gene tree. Strains KSW-18T and KSW2-22 shared an OrthoANIu of 98.6% and a digital DNA-DNA hybridization of 87.6% with each other, representing that they were strains of a species, while the OrthoANIu and digital DNA-DNA hybridization values between strains KSW-18T and KSW4-11T, and between both of these isolates and all members of the genus Microbacterium were ≤86.5% and ≤30.7%, respectively. The analyses of overall genomic relatedness indices and phenotypic distinctness support that the three isolates represent two new species of the genus Microbacterium. Based on the results obtained here, Microbacterium aquilitoris sp. nov. (type strain KSW-18T = KCTC 49623T = NBRC 115222T) and Microbacterium gwkjiense sp. nov. (type strain KSW4-11T = KACC 23321T = DSM 116380T) are proposed.

Flavobacterium acetivorans sp. nov. and Flavobacterium galactosidilyticum sp. nov., psychrotolerant flavobacteria isolated from cultivated fish.

Two yellow-coloured strains, F-29T and F-340T, were isolated from fish farms in Antalya and Mugla in 2015 and 2017 during surveillance studies. The 16S rRNA gene sequence analysis showed that both strains belong to the genus Flavobacterium. A polyphasic approach involving a comprehensive genome analysis was employed to ascertain the taxonomic provenance of the strains. The overall genome-relatedness indices of digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) between the strains and the other members of the genus Flavobacterium were found to be well below the established thresholds of 70 and 95 %, respectively. The whole-genome-based phylogenetic analysis revealed that strain F-29T is closely related to Flavobacterium granuli (dDDH 39.3 % and ANI 89.4 %), while strain F-340T has a close relationship with the type strain of Flavobacterium pygoscelis (dDDH 25.6 % and ANI 81.5 %). Both strains were psychrotolerant with an optimum growth temperature of 25 °C. The chemotaxonomic characteristics of the strains were typical of the genus Flavobacterium. Both strains had phosphatidylethanolamine, aminolipids and unidentified lipids in their polar lipid profile and MK-6 as the isoprenoid quinone. The major fatty acids were iso-C15 : 0 and anteiso-C15 : 0. The genome size of the strains was 3.5 Mb, while G+C contents were 35.3 mol% for strain F-29T and 33.4 mol% for strain F-340T. Overall, the characterizations confirmed that both strains are representatives of two novel species within the genus Flavobacterium, for which the names Flavobacterium acetivorans sp. nov. and Flavobacterium galactosidilyticum sp. nov. are proposed, with F-29T (JCM 34193T=KCTC 82253T) and F-340T (JCM 34203T=KCTC 82263T) as the type strains, respectively.