Whole Genome Sequencing Research Articles

Emergence of ECSA-IOL E1-K211E/E2-V264A Lineage of Chikungunya virus during Malaysian 2021 outbreak

BackgroundChikungunya cases was reported to be on the rise in Malaysia from 2019 to 2021. Although potential endemicity was described previously, genotype shift during 2008 outbreak originating from the 2004 Indian Ocean Islands outbreak presents the probability of current CHIKV spread from neighboring countries. This is due to the prevalence of the new IOL sub-lineage which consists of E1-226A wildtype or reverted strains that are circulating in the Indian subcontinent before spreading to neighboring Thailand during 2018–2019 outbreak.MethodIn this study, samples received mostly from the Tangkak, Johor were analyzed. A total 56 CHIKV positive serum samples received in 2021 by Institute of Medical Research Malaysia (IMR), were collected based on sample selection criteria. Selected samples were subjected to total RNA extraction, whole-genome sequencing as well as bioinformatic analysis such as phylogenetic, variant and mutation analysis.ResultsBased on the genomic and phylogenetic analysis, the CHIKV samples from 2021 outbreak were of ECSA-IOL genotype. Genome similarity analysis also revealed that these CHIKVs were highly similar to 2018–2019 outbreak strain from Thailand. In comparison to the 2008 outbreak CHIKV isolate, the current CHIKVs lacked the E1-A226V mutation and harbored the new E1-K211E/E2-V264A sub-linage mutation. Since the E1-K211E/E2-V264A mutation facilitates adaptation to Ae. aegypti as opposed to the E1-A226V mutation which improves adaptation to Ae. albopictus, the emergence 2021 CHIKV outbreak in Malaysia can be postulated due to vector shift. Interestingly, a novel nsP3-T441A/V mutation detected in this study, may also play a role in virus transmission, pathogenicity, fitness and vector adaptation.ConclusionIn summary, the current CHIKV outbreak are strains originated from the Indian subcontinent through Thailand which may have capitalized on vector shifting by adapting to Ae. aegypti. The presence of novel nsP3-T441A/V mutation may also contribute to the spread of this virus across peninsular Malaysia.

Insights into the DNA methylation of Portunus trituberculatus in response to Vibrio parahaemolyticus infection

Vibrio parahaemolyticus is the main pathogen causing acute hepatopancreatic necrotic disease in crustaceans. To elucidate the epigenetic regulatory mechanism of crustacean resistance to V. parahaemolyticus infection, we conducted artificial infection studies on Portunus trituberculatus. The results showed that the mortality rate reached the highest at 12 h of artificial infection, which was 23.69 %. At 72 h after V parahaemolyticus infection, the expression level of DNA demethylase (ten-eleven-translocation protein) Tet was significantly decreased, the expression of DNA methyltransferase Dnmt3B fluctuated significantly. Based on the differential expression levels of Tet and Dnmt3B. We depict for DNA methylation profiles of the whole genome of P. trituberculatus at single-base resolution by using whole-genome bisulfite sequencing (WGBS) on hemolymph tissues. The overall DNA methylation level was low at 2.16% in P. trituberculatus hemolymph. A total of 2,590 differentially methylated regions (DMRs) were identified, of which 1,329 were hypermethylated and 1,261 were hypomethylated, and 1,389 genes were annotated in these DMRs. Differently methylated genes (DMGs) were significantly enriched in ribosomes (KO03010), protein kinases (KO01001), cell cycle (HSA04110), endocrine resistance (HSA01522) and FoxO signaling pathway (KO04068). Finally, we selected six differentially methylated genes for quantitative analysis. The results showed that DNA methylation not only has a negative regulatory effect on gene expression, but also has a positive regulatory effect. These results indicated that DNA methylation in the regulation of genes involved in immune responses contributes to the resistance of P. trituberculatus to V. parahaemolyticus, which is valuable for understanding how crustaceans regulate the innate immune system to defend against bacterial infections.

Genomes of Microtus rodents highlight the importance of olfactory and immune systems in their fast radiation.

The characterization of genes and biological functions underlying functional diversification and the formation of species is a major goal of evolutionary biology. In this study, we investigated the fast radiation of Microtus voles, one of the most speciose group of mammals, which shows strong genetic divergence despite few readily observable morphological differences. We produced an annotated reference genome for the common vole, Microtus arvalis, and resequenced the genomes of 10 different species and evolutionary lineages spanning the Microtus speciation continuum. Our full genome sequences illustrate the recent and fast diversification of this group, and we identified genes in highly divergent genomic windows that have likely particular roles in their radiation. We found three biological functions enriched for highly divergent genes in most Microtus species and lineages: olfaction, immunity and metabolism. In particular, olfaction-related genes (mostly olfactory receptors and vomeronasal receptors) are fast evolving in all Microtus species indicating the exceptional importance of the olfactory system in the evolution of these rodents. Of note is e.g. the shared signature among vole species on Olfr1019 which has been associated with fear responses against predator odours in rodents. Our analyses provide a genome-wide basis for the further characterization of the ecological factors and processes of natural and sexual selection that have contributed to the fast radiation of Microtus voles.

A first report of rotavirus B from Zambian pigs leading to the discovery of a novel VP4 genotype P[9]

BackgroundRotavirus B (RVB) causes diarrhea in humans and pigs. Although various RVB strains were identified in humans and various animals globally, little is known about the epidemiology RVB infection in Africa. In this study, we attempted to examine the prevalence of RVB infection in pig populations in Zambia.MethodsMetagenomic analyses were conducted on pig feces collected in Zambia to detect double stranded RNA viruses, including RVB. To clarify the prevalence of RVB infection in pig populations in Zambia, 147 fecal samples were screened for the RVB detection by RT-qPCR. Full genome sequence of a detected RVB was determined by Sanger sequencing and genetically analyzed.ResultsThe metagenomic analyses revealed that RVB sequence reads and contigs of RVB were detected from one fecal sample collected from pigs in Zambia. RT-qPCR screening detected RVB genomes in 36.7% (54/147) of fecal samples. Among 54 positive samples, 13 were positive in non-diarrheal samples (n = 48, 27.1%) and 41 in diarrheal samples (n = 99, 41.4%). Genetic analyses demonstrated that all the segments of ZP18-18, except for VP4, had high nucleotide sequence identities (80.6–92.6%) with all other known RVB strains detected in pigs. In contrast, the VP4 sequence of ZP18-18 was highly divergent from other RVB strains (< 64.6% identities) and formed a distinct lineage in the phylogenetic tree. Notably, the VP8 subunit of the VP4 showed remarkably low amino acid identities (33.3%) to those of known RVB strains, indicating that the VP8 subunit of ZP18-18 was unique among RVB strains. According to the whole genome classification for RVB, ZP18-18 was assigned to a genotype constellation, G18-P[9]-I12-R4-C4-M4-A8-N10-T5-E4-H7 with the newly established VP4 genotype P[9].ConclusionsThis current study updates the geographical distribution and the genetic diversity of RVB. Given the lack of information regarding RVB in Africa, further RVB surveillance is required to assess the potential risk to humans and animals.

Dengue virus surveillance in Nepal yields the first on-site whole genome sequences of isolates from the 2022 outbreak

BackgroundThe 4 serotypes of dengue virus (DENV1-4) can each cause potentially deadly dengue disease, and are spreading globally from tropical and subtropical areas to more temperate ones. Nepal provides a microcosm of this global phenomenon, having met each of these grim benchmarks. To better understand DENV transmission dynamics and spread into new areas, we chose to study dengue in Nepal and, in so doing, to build the onsite infrastructure needed to manage future, larger studies.Methods and resultsDuring the 2022 dengue season, we enrolled 384 patients presenting at a hospital in Kathmandu with dengue-like symptoms; 79% of the study participants had active or recent DENV infection (NS1 antigen and IgM). To identify circulating serotypes, we screened serum from 50 of the NS1+ participants by RT-PCR and identified DENV1, 2, and 3 – with DENV1 and 3 codominant. We also performed whole-genome sequencing of DENV, for the first time in Nepal, using our new on-site capacity. Sequencing analysis demonstrated the DENV1 and 3 genomes clustered with sequences reported from India in 2019, and the DENV2 genome clustered with a sequence reported from China in 2018.ConclusionThese findings highlight DENV’s geographic expansion from neighboring countries, identify China and India as the likely origin of the 2022 DENV cases in Nepal, and demonstrate the feasibility of building onsite capacity for more rapid genomic surveillance of circulating DENV. These ongoing efforts promise to protect populations in Nepal and beyond by informing the development and deployment of DENV drugs and vaccines in real time.

Whole-genome sequence of Bacillus paralicheniformis strain R243 isolated from rice (Oryza sativa L.) roots.

We isolated an endophytic strain R243 and sequenced its genome using the NextSeq 1000 (Illumina) sequencing platform. The sequenced genome contains 4,478 protein-coding genes with a mean GC content of 45.56%. The strain R243 was identified as Bacillus paralicheniformis. Gene clusters involved in the synthesis of biocontrol substances were identified.

Production and analysis of synthesized bacterial cellulose by Enterococcus faecalis strain AEF using Phoenix dactylifera and Musa acuminata fruit extracts.

Bacterial cellulose (BC) is a highly versatile biopolymer renowned for its exceptional mechanical strength, water retention, and biocompatibility. These properties make it a valuable material for various industrial and biomedical applications. In this study, Enterococcus faecalis synthesized extracellular BC, utilizing Phoenix dactylifera and Musa acuminata fruit extracts as sustainable carbon sources. LC-MS analysis identified glucose as the primary carbohydrate in these extracts, providing a suitable substrate for BC production. Scanning Electron Microscopy (SEM) revealed a network of BC nanofibers on Congo red agar plates. ATR-FTIR spectroscopy confirmed the presence of characteristic cellulose functional groups, further supporting BC synthesis. X-ray diffraction (XRD) analysis indicated a high crystallinity index of 71%, consistent with the cellulose I structure, as evidenced by peaks at 16.22°, 21.46°, 22.52°, and 34.70°. Whole-genome sequencing of E. faecalis identified vital genes involved in BC biosynthesis, including bcsA, bcsB, diguanylate cyclase (DGC), and 6-phosphofructokinase (pfkA). Antibiotic susceptibility tests revealed resistance to cefotaxime, ceftazidime, and ceftriaxone, while susceptibility to imipenem was observed. Quantitative assessment demonstrated that higher concentrations of fruit extracts (5.0-20mg/mL) significantly enhanced BC production. Cytotoxicity testing via the MTT assay confirmed excellent biocompatibility with NIH/3T3 fibroblast cells, showing high cell viability (97-105%). Unlike commonly studied Gram-negative bacteria like Acetobacter xylinum for BC production, this research focuses on Gram-positive Enterococcus faecalis and utilizes Phoenix dactylifera and Musa acuminata fruit extracts as carbon sources. This approach offers a sustainable and promising avenue for BC production.

Genomic and phenotypic characterization of Pseudomonas aeruginosa isolates from two Mexican cystic fibrosis attention centers.

Thirty-nine clinical isolates of Pseudomonas aeruginosa collected from 11 cystic fibrosis (CF) patients at two CF attention centers over 10 years were subjected to whole genome sequencing (WGS). Phenotypic tests (i.e., elastase, motility, biofilm, growth rate, and antibiotic susceptibility) were performed to correlate results. A single strain of P. aeruginosa was found to persist over time in longitudinal isolates. No transmission between patients or centers was observed. A tendency to lack genes related to pyoverdine, flagellum, pili, and O-antigen was observed, whereas those related to biofilm, phenazine, and pyochelin were conserved among isolates. In a patient with a 10-year follow-up, a single strain of P. aeruginosa persisted and showed a gradual decrease in elastase activity and growth rate, demonstrating an adaptive phenotype.IMPORTANCEThis study investigates the genomic and phenotypic characteristics of Pseudomonas aeruginosa isolates from Mexican cystic fibrosis (CF) patients, an underrepresented group in CF research. To our knowledge, it is the first to use whole genome sequencing (WGS) to study longitudinally collected P. aeruginosa isolates from this population, evaluating both genomic features and clonal relationships. Remarkably, the study includes samples from one patient over 10 years, offering an extended observation time compared to existing literature. Unlike similar studies, which often lack phenotypic testing, this research incorporates various virulence-related phenotypic assays, enhancing our understanding of gene-to-phenotype correlations. Two potential mechanisms for the loss of elastolytic activity were identified. Furthermore, we conduct an in-depth mobilome analysis, an area that remains largely unexplored in CF contexts. Whole genome sequencing data are publicly available through the NCBI SRA database, facilitating further re-analysis for studies on P. aeruginosa in CF, as well as epidemiological and population structure research.

Whole-genome sequencing analyses and antibiotic resistance situation of 48 Helicobacter pylori strains isolated in Zhejiang, China

PurposeIn the Zhejiang region, research on Helicobacter pylori is lacking. The purpose of this study was to assess the extent of antibiotic resistance in H. pylori in this region, explore alternative methods for predicting the resistance patterns of H. pylori, and investigate the colonization of native gastric mucosa by other clades of H. pylori in the structure population of this bacterium.MethodsStrains were cultured under microaerobic conditions, and antimicrobial susceptibility testing (AST) was performed via agar dilution. Whole-genome sequencing (WGS) was performed via next-generation sequencing (NGS) technology. Epidemiological data including data from this study and reported articles from Zhejiang, China, were included. Further analyses based on AST, WGS, and epidemiological date include virulence genes, antibiotic resistance-related mutations, and phylogenetic trees based on 7 housekeeping genes and core-genome single nucleotide polymorphisms (SNPs).ResultsThe bacterial isolates in this study presented higher antibiotic resistance rates than previously reported, especially against levofloxacin and clarithromycin. The point mutation A2147G in 23 S rRNA is specific to clarithromycin resistance. Mutations at position/s 87 and/or 91 of the gyrA gene amino acid sequence are highly consistent with levofloxacin resistance highly. The point mutations C1707T in 23 S rRNA and E463K in the gyrB gene have not been previously documented in China. All the bacterial isolates belong to Asian branches in the structure population. The resistance rate to clarithromycin of isolates from hosts born after January 1, 1977 is statistically higher than that of hosts born before 1977.ConclusionEradication therapy based on AST results is urgently needed in Zhejiang. The point mutation A2147G in 23 S rRNA and point mutations in the gyrA gene at amino acid/s 87 and/or 91 are sufficient for predicting resistance to clarithromycin and levofloxacin, respectively. The isolate with the mutation E463K in the gyrB gene represents a significant contribution to the field. Mutations in 23 S rRNA may offer valuable insights into the dynamics of H. pylori transmission among hosts.

Frameshift Mutations in Genes Encoding PBP3 and PBP4 Trigger an Unusual, Extreme β-Lactam Resistance Phenotype in Burkholderia multivorans.

In our curated panel of Burkholderia cepacia complex isolates, Burkholderia multivorans strain AU28442 was unusually highly β-lactam resistant. To explore the molecular mechanisms leading to this phenotype, we performed whole genome sequencing (WGS) and microbiological and biochemical assays. WGS analysis revealed that strain AU28442 produced two β-lactamases, AmpC22 and a novel PenA-like β-lactamase denominated PenA39. Additionally, the strain presented frame-shift mutations in the genes encoding penicillin binding proteins 3 (PBP3) and 4 (PBP4). The antibiotic susceptibilities of the parent AU28442 strain carrying blaPenA39 vs the isogenic E. colistrain producing blaPenA39 were discrepant with ceftazidime MICs of >512 and 1 μg/mL, respectively. Accordingly, PenA39 was found to poorly hydrolyze β-lactams with kcat values of ≤8.8 s-1. An overlay of the crystal structure of PenA39 with PenA1 revealed a shift in the SDN loop in the variant, which may affect the catalytic efficiency of PenA39 toward substrates and inhibitors. Moreover, microscopic examination of AU28442 revealed shortened rod-shaped cells compared to B. multivoransATCC 17616, which carries a full complement of intact PBPs. Further complementation assays confirmed that the loss of PBP3 and PBP4 was the main factor contributing to the high-level β-lactam resistance observed in B. multivoransAU28442. This information allowed us to revert susceptibility by pairing a potent β-lactamase inhibitor with a β-lactam with promiscuous PBP binding. This detailed characterization of B. multivoransprovides an illustration of the myriad ways in which bacteria under antibiotic selection can develop resistance and demonstrates a mechanism to overcome it.

Characterisation of Neisseria gonorrhoeae strain differences in patients with multisite infection

BackgroundNeisseria gonorrhoeae, the aetiological agent of gonorrhoea, is an increasing global health priority due to high levels of antimicrobial resistance (AMR). It is estimated that up to 42% of patients...

Phylogenetic Characterization of Novel Reassortant 2.3.4.4b H5N8 Highly Pathogenic Avian Influenza Viruses Isolated from Domestic Ducks in Egypt During the Winter Season 2021–2022

Avian influenza (AI) is an extremely contagious viral disease of domestic and wild birds that can spread rapidly among bird populations, inducing serious economic losses in the poultry industry. During the winter season 2021–2022, we isolated seventeen highly pathogenic avian influenza (HPAI) H5N8 viruses from outbreaks involving ducks in Egypt, occurring in both backyard and farm settings. The aim of this study was to pinpoint genetic key substitutions (KSs) that could heighten the risk of a human pandemic by influencing the virus’s virulence, replication ability, host specificity, susceptibility to drugs, or transmissibility. To understand their evolution, origin, and potential risks for a human pandemic, whole-genome sequencing and phylogenetic analysis were conducted. Our analysis identified numerous distinctive mutations in the Egyptian H5N8 viruses, suggesting potential enhancements in virulence, resistance to antiviral drugs, and facilitation of transmission in mammals. In this study, at least five genotypes within one genome constellation of H5N8 viruses were identified, raising concerns about the potential emergence of novel viruses with altered characteristics through reassortment between different genotypes and distinct groups. These findings underscore the role of ducks in the virus’s evolutionary process and emphasize the urgent need for enhanced biosecurity measures in domestic duck farms to mitigate pandemic risk.

Polygenic Landscape of Cryptogenic New-Onset Refractory Status Epilepticus: A Comprehensive Whole-Genome Sequencing Study.

Cryptogenic new-onset refractory status epilepticus (cNORSE) is a devastating condition with unclear pathogenesis. Here, we analyzed the genetic underprints of 31 cNORSE patients from an autoimmune encephalitis observational cohort through whole-genome sequencing. Compared to their controls, cNORSE patients exhibited elevated polygenic risk scores (PRS) for traits associated with autoimmune diseases. The individual PRS against these diseases were correlated with specific clinical phenotypes of cNORSE. The variants were enriched in genes expressed in the central nervous system and lymphocytes. These results suggest a shared genetic framework between cNORSE and other autoimmune/autoinflammatory diseases, and its involvement in the disease pathogenesis. ANN NEUROL 2024.

Novel genetic variants of banana streak MY virus and banana streak IM virus naturally infecting banana in Northeast India.

Divergent banana streak viruses (BSV) were characterized from banana plants exhibiting diverse symptoms in the Northeast region (NER) of India. Using rolling circle amplification (RCA), the complete genome sequences of seven episomal banana streak MY virus (BSMYV) isolates, including two novel variants, and two new banana streak IM virus (BSIMV) isolates were characterized. The novel BSMYV genetic variants were associated with conspicuous necrosis on newly emerged leaves, peduncle distortion, pseudostem internal necrosis, in addition to common streak symptoms. For complete genome nucleotide sequences, BSMYV-IN4 and IN5 shared 77-79% identity with other BSMYVs, while BSMYV-IN7 and IN8 exhibited identities of 77-97%. This study reports for the first time, the complete genomes of two banana streak IM virus (BSIMV-IN1 and -IN2) infecting triploid banana hybrids exhibiting leaf distortion, stunted rosette-like growth, and necrosis, sharing 87% sequence identity with reference BSIMV genome (GenBank accession no. HQ593112). Phylogenetic inference based on complete genomes revealed the distinct and congruent placement of BSMYV-IN4 and IN5 within the BSMYV cluster. Pairwise sequence comparisons of the conserved RT/RNase H nucleotide (nt) sequences revealed that the BSMYV-IN7 and IN4 isolates showed 85% and 97% identity to BSMYV (AY805074), respectively, which shared highest nt identity with BSMYV-IN6, IN9, and IN10, at 100%. The RT/RNase H nt sequences of BSIMV-IN1 and IN2 had 98% identity with the BSIMV (HQ593112), but were characterized as novel variants of BSIMV based on complete genomes. An analysis of relative synonymous codon usage (RSCU) pattern in the ORFIII polyprotein of BSMYV and BSIMV isolates revealed AGA and AGG (arginine) as the most frequently overrepresented codons (>1.5), evolutionary conserved in the genome of both species. A total of 14 recombination events were detected among the 36 BSV genomes, with recombination breakpoints mainly located in the ORFI, III, and IGR genomic regions. A novel phylogenetic cluster, comprised of BSMYV-IN4 and IN5 within the clade I was probably derived from heterologous recombination between parents resembling banana streak VN virus (BSVNV; AY750155) and banana streak GF virus (BSGFV; KJ013507) isolates. The present study conclusively reports the infection of genetically and symptomatically distinct variants of BSMYV and BSIMV infecting banana hybrids in NER India.

Identifying factors that influence the use of pathogen genomics in Australia and New Zealand: a protocol

IntroductionPathogen genomics, where whole genome sequencing technologies are used to produce complete genomic sequences of pathogens, is being increasingly used for infectious disease surveillance and outbreak response. Although proof-of-concept studies have highlighted the viability of using pathogen genomics in public health, few studies have investigated how end-users utilize pathogen genomics in public health. We describe a protocol for a study that aims to identify key factors that influence the use of pathogen genomics to inform public health responses against infectious diseases in Australia and New Zealand.MethodsWe will use qualitative comparative analysis (QCA), a case-oriented methodology that systematically compares and analyses multiple cases (or ‘units of analysis’), to identify multiple pathways leading to the use of pathogen genomics results in public health actions. As part of the process, we will develop a rubric to identify and define the use of pathogen genomics and individual factors affecting this process. Simultaneously, we will identify cases where pathogen genomics has been used in public health across Australia and New Zealand. Data for these cases will be collected from document review of publicly available and confidential documents and semi-structured interviews with technicians and end-users and summarized in a case report. These case reports will form the basis for scoring each case on the extent of the use of pathogen genomics data and the presence or absence of specific factors such as the ease of extracting essential information from pathogen genomics reports and perceptions toward pathogen genomics. Using the scores, cases will be analyzed using QCA techniques to identify pathways leading to the use of pathogen genomics data. These pathways will be interpreted alongside the cases to provide rich explanations of the use of pathogen genomics in public health.DiscussionThis study will improve our understanding of the key factors that facilitate or hinder the use of pathogen genomics to inform public health authorities and end-users. These findings may inform ways to enhance the use of pathogen genomics data in public health.

Phylogenetic Insights into H7Nx Influenza Viruses: Uncovering Reassortment Patterns and Geographic Variability

Influenza A viruses (IAVs), which belong to the Orthomyxoviridae family, are RNA viruses characterized by a segmented genome that allows them to evolve and adapt rapidly. These viruses are mainly transmitted by wild waterfowl. In this study, we investigated the evolutionary processes of H7Nx (H7N1, H7N2, H7N3, H7N4, H7N5, H7N6, H7N7, H7N8, H7N9) viruses, which pose a significant pandemic risk due to the known cases of human infection and their potential for rapid genetic evolution and reassortment. The complete genome sequences of H7Nx influenza viruses (n = 3239) were compared between each other to investigate their phylogenetic relationships and reassortment patterns. For the selected viruses, phylogenetic trees were constructed for eight genome segments (PB2, PB1, PA, HA, NP, NA, M, NS) to assess the genetic diversity and geographic distribution of these viruses. Distinct phylogenetic clades with remarkable geographic patterns were found for the different segments. While the viruses were consistently grouped by subtype based on the NA segment sequences, the phylogeny of the other segment sequences, with the exception of the NS segment, showed distinct grouping patterns based on geographic origin rather than formal subtype assignment. Reassortment events leading to complex phylogenetic relationships were frequently observed. In addition, multiple cases of previously undescribed reassortments between subtypes were detected, emphasizing the fluidity of H7Nx virus populations. These results indicate a high degree of genetic diversity and reassortment within H7Nx influenza viruses. In other words, H7Nx viruses exist as constantly changing combinations of gene pools rather than stable genetic lineages.

First insights into the prevalence, genetic characteristics, and pathogenicity of Bacillus cereus from generations worldwide.

Bacillus cereus, a global threat, is one of the major causes of toxin-induced foodborne diseases. However, a comprehensive assessment of the prevalence and characteristics of B. cereus worldwide is still lacking. Here, we applied whole-genome sequence analysis to 191 B. cereus collected in Africa, America, Asia, Europe, and Oceania from the 1900s to 2022, finding that CC142 dominated the global B. cereus clonal complex. The results provided direct evidence that B. cereus could spread through the food chain and intercontinentally. B. cereus from different generations worldwide showed coherence in the antibiotic-resistant gene and virulence and biofilm gene profiles, although with high genomic heterogeneity. The BCI-BCII-vanZF-fosB profiles and virulence and biofilm genes were detected at high rates, and we emphasized that B. cereus would pose a serious challenge to global public health and clinical treatment.IMPORTANCEThis study first emphasized the prevalence, genetic characteristics, and pathogenicity of Bacillus cereus worldwide from the 1900s to 2022 using whole-genome sequence analysis. The CC142 dominated the global Bacillus cereus clonal complex. Moreover, we revealed a close evolutionary relationship between the isolates from different sources. B. cereus isolates from different generations worldwide showed coherence in potential pathogenicity, although with high genomic heterogeneity. The BCI-BCII-vanZF-fosB profiles and virulence and biofilm genes were detected at high rates, and we emphasized that B. cereus would pose a serious challenge to global public health and clinical treatment.

Flavobacterium capsici sp. nov., isolated from the rhizospheric soils of bell pepper (Capsicum annuum L.).

Two separate bacterial strains, PMTSA4T and PMR2A8, were isolated from the rhizospheric soils of bell pepper plants grown in a plant nursery. These strains are Gram-negative, non-motile and rod-shaped and grow in aerobic conditions. They exhibit a positive reaction for catalase activity but negative results for oxidase activity. Phylogenetic analysis of the 16S rRNA gene sequences revealed that the strains PMTSA4T and PMR2A8 are closely related to Flavobacterium piscinae ICH-30T (95.6%, respectively), Flavobacterium ahnfeltiae 10Alg 130T (95.5%) and Flavobacterium maris KMM 9535T (95.3%), aligning them within the genus Flavobacterium. Digital DNA-DNA hybridization (dDDH) values and average nucleotide identities (ANIs) of the whole-genome sequences for the two strains and related Flavobacterium species were significantly below the established thresholds for prokaryotic species delineation (<70% for dDDH and <95% for ANI). The observed values were as follows: Flavobacterium aquatile LMG 4008T (dDDH: 19.8% and ANI: 75.5%), F. piscinae ICH-30T (dDDH: 18.6% and ANI: 73.3%) and F. stagni WWJ 16T (dDDH: 18.5% and ANI: 72.0%). The strains have genome sizes of 3 068 185 bp and 3 068 330 bp, with a G+C content of 32.5 mol%. In phenotypic characterization, the new strains grew at 10-35 °C and tolerated up to 4% NaCl at pH 5-9 (optimum pH 8). The predominant cellular fatty acids were observed to be iso-C15:0, iso-C17:0 3-OH and iso-C15:0 3-OH. Menaquinone-6 was the predominant quinone. Considering the results from phenotypic, chemotaxonomic, phylogenetic and genomic analyses, it is proposed that the strains PMTSA4T and PMR2A8 represent a novel species within the genus Flavobacterium.

Blood-derived mitochondrial DNA copy number is associated with Alzheimer disease, Alzheimer-related biomarkers and serum metabolites

BackgroundBlood-derived mitochondrial DNA copy number (mtDNA-CN) is a proxy measurement of mitochondrial function in the peripheral and central systems. Abnormal mtDNA-CN not only indicates impaired mtDNA replication and transcription machinery but also dysregulated biological processes such as energy and lipid metabolism. However, the relationship between mtDNA-CN and Alzheimer disease (AD) is unclear.MethodsWe performed two-sample Mendelian randomization (MR) using publicly available summary statistics from GWAS for mtDNA-CN and AD to investigate the causal relationship between mtDNA-CN and AD. We estimated mtDNA-CN using whole-genome sequence data from blood and brain samples of 13,799 individuals from the Alzheimer’s Disease Sequencing Project. Linear and Cox proportional hazards models adjusting for age, sex, and study phase were used to assess the association of mtDNA-CN with AD. The association of AD biomarkers and serum metabolites with mtDNA-CN in blood was evaluated in Alzheimer’s Disease Neuroimaging Initiative using linear regression. We conducted a causal mediation analysis to test the natural indirect effects of mtDNA-CN change on AD risk through the significantly associated biomarkers and metabolites.ResultsMR analysis suggested a causal relationship between decreased blood-derived mtDNA-CN and increased risk of AD (OR = 0.68; P = 0.013). Survival analysis showed that decreased mtDNA-CN was significantly associated with higher risk of conversion from mild cognitive impairment to AD (HR = 0.80; P = 0.002). We also identified significant associations of mtDNA-CN with brain FDG-PET (β = 0.103; P = 0.022), amyloid-PET (β = 0.117; P = 0.034), CSF amyloid-β (Aβ) 42/40 (β=-0.124; P = 0.017), CSF t-Tau (β = 0.128; P = 0.015), p-Tau (β = 0.140; P = 0.008), and plasma NFL (β=-0.124; P = 0.004) in females. Several lipid species, amino acids, biogenic amines in serum were also significantly associated with mtDNA-CN. Causal mediation analyses showed that about a third of the effect of mtDNA-CN on AD risk was mediated by plasma NFL (P = 0.009), and this effect was more significant in females (P < 0.005).ConclusionsOur study indicates that mtDNA-CN measured in blood is predictive of AD and is associated with AD biomarkers including plasma NFL particularly in females. Further, we illustrate that decreased mtDNA-CN possibly increases AD risk through dysregulation of mitochondrial lipid metabolism and inflammation.

Whole-genome assembly and annotation of the firecracker penstemon (Penstemon eatonii).

The penstemons are ornamental annual flowering plants native to the Intermountain West and Rocky Mountains and commonly used for urban landscaping. Elite commercial penstemons are generally susceptible to abiotic stresses, including drought, root rot, cold, and high salinity. Firecracker penstemon (Penstemon eatonii), however, is much more tolerant to these stresses than most elite cultivars. Importantly, firecracker penstemon has been reported to hybridize with many other penstemons and therefore provides the opportunity to develop more tolerant elite cultivars through strategic crossing. To facilitate the study and utilization of firecracker penstemon, we sequenced and annotated the genome of a P. eatonii accession collected from Utah, USA. We also performed low-coverage, whole-genome sequencing of 26 additional accessions from three different varieties of P. eatonii. This chromosome-scale genome assembly is the most contiguous and complete Penstemon genome sequenced to date.