Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing.

Whole-genome reconstruction of bacterial pathogens has become an important tool for tracking transmission and antimicrobial resistance gene spread, but highly accurate and complete assemblies have largely only historically been achievable using hybrid long- and short-read sequencing. We previously found the Oxford Nanopore Technologies (ONT) R10.4/kit12 flowcell/chemistry produced improved assemblies over the R9.4.1/kit10 combination, however long-read only assemblies contained more errors compared to Illumina-ONT hybrid assemblies. ONT have since released an R10.4.1/kit14 flowcell/chemistry upgrade and recommended the use of Bovine Serum Albumin (BSA) during library preparation, both of which reportedly increase accuracy and yield. They have also released updated basecallers trained using native bacterial DNA containing methylation sites intended to fix systematic basecalling errors, including common adenosine (A) to guanine (G) and cytosine (C) to thymine (T) substitutions. To evaluate these improvements, we successfully sequenced four bacterial reference strains, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, and nine genetically diverse E. coli bloodstream infection-associated isolates from different phylogroups and sequence types, both with and without BSA. These sequences were de novo assembled and compared against Illumina-corrected reference genomes. In this small evaluation of 13 isolates we found that nanopore long-read-only R10.4.1/kit 14 assemblies with updated basecallers trained using bacterial methylated DNA produce accurate assemblies with ≥40×depth, sufficient to be cost-effective compared with hybrid ONT/Illumina sequencing in our setting.

Complete genome sequence of a novel mitovirus identified in the phytopathogenic fungus Fusarium oxysporum f. sp. melonis strain T-SD3.

A novel mitovirus was identified in Fusarium oxysporum f. sp. melonis strain T-SD3 and designated as "Fusarium oxysporum mitovirus 3" (FoMV3). The virus was isolated from diseased muskmelon plants with the typical symptom of fusarium wilt. The complete genome of FoMV3 is 2269 nt in length with a predicted AU content of 61.40% and contains a single open reading frame (ORF) using the fungal mitochondrial genetic code. The ORF was predicted to encode a polypeptide of 679 amino acids (aa) containing a conserved RNA-dependent RNA polymerase (RdRp) domain with a molecular mass of 77.39 kDa, which contains six conserved motifs with the highly conserved GDD tripeptide in motif IV. The 5'-untranslated region (UTR) and 3'-UTR of FoMV3 were predicted to fold into stem-loop structures. BLASTp analysis revealed that the RdRp of FoMV3 shared the highest aa sequence identity (83.85%) with that of Fusarium asiaticum mitovirus 5 (FaMV5, a member of the family Mitoviridae) infecting F. asiaticum, the causal agent of wheat fusarium head blight. Phylogenetic analysis further suggested that FoMV3 is a new member of the genus Unuamitovirus within the family Mitoviridae. This is the first report of a new mitovirus associated with F. oxysporum f. sp. melonis.

Multimodal analysis of cfDNA methylomes for early detecting esophageal squamous cell carcinoma and precancerous lesions.

Detecting early-stage esophageal squamous cell carcinoma (ESCC) and precancerous lesions is critical for improving survival. Here, we conduct whole-genome bisulfite sequencing (WGBS) on 460 cfDNA samples from patients with non-metastatic ESCC or precancerous lesions and matched healthy controls. We develop an expanded multimodal analysis (EMMA) framework to simultaneously identify cfDNA methylation, copy number variants (CNVs), and fragmentation markers in cfDNA WGBS data. cfDNA methylation markers are the earliest and most sensitive, detectable in 70% of ESCCs and 50% of precancerous lesions, and associated with molecular subtypes and tumor microenvironments. CNVs and fragmentation features show high specificity but are linked to late-stage disease. EMMA significantly improves detection rates, increasing AUCs from 0.90 to 0.99, and detects 87% of ESCCs and 62% of precancerous lesions with >95% specificity in validation cohorts. Our findings demonstrate the potential of multimodal analysis of cfDNA methylome for early detection and monitoring of molecular characteristics in ESCC.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.

Molnupiravir increases SARS-CoV-2 genome diversity and complexity: A case-control cohort study.

Molnupiravir, an oral direct-acting antiviral effective in vitro against SARS-CoV-2, has been largely employed during the COVID-19 pandemic, since December 2021. After marketing and widespread usage, a progressive increase in SARS-CoV-2 lineages characterized by a higher transition/transversion ratio, a characteristic signature of molnupiravir action, appeared in the Global Initiative on Sharing All Influenza Data (GISAID) and International Nucleotide Sequence Database Collaboration (INSDC) databases. Here, we assessed the drug effects by SARS-CoV-2 whole-genome sequencing on 38 molnupiravir-treated persistently positive COVID-19 outpatients tested before and after treatment. Seventeen tixagevimab/cilgavimab-treated outpatients served as controls. Mutational analyses confirmed that SARS-CoV-2 exhibits an increased transition/transversion ratio seven days after initiation of molnupiravir. Moreover we observed an increased G->A ratio compared to controls, which was not related to apolipoprotein B mRNAediting enzyme, catalytic polypeptide-like (APOBEC) activity. In addition, we demonstrated for the first time an increased diversity and complexity of the viral quasispecies.

Whole-genome sequencing of Ganoderma boninense, the causal agent of basal stem rot disease in oil palm, via combined short- and long-read sequencing.

The hemibiotrophic Basidiomycete pathogen Ganoderma boninense (Gb) is the dominant causal agent of oil palm basal stem rot disease. Here, we report a complete chromosomal genome map of Gb using a combination of short-read Illumina and long-read Pacific Biosciences (PacBio) sequencing platforms combined with chromatin conformation capture data from the Chicago and Hi-C platforms. The genome was 55.87 Mb in length and assembled to a high contiguity (N50: 304.34 kb) of 12 chromosomes built from 112 scaffolds, with a total of only 4.34 Mb (~ 7.77%) remaining unplaced. The final assemblies were evaluated for completeness of the genome by using Benchmarking Universal Single Copy Orthologs (BUSCO) v4.1.4, and based on 4464 total BUSCO polyporales group searches, the assemblies yielded 4264 (95.52%) of the conserved orthologs as complete and only a few fragmented BUSCO of 42 (0.94%) as well as a missing BUSCO of 158 (3.53%). Genome annotation predicted a total of 21,074 coding genes, with a GC content ratio of 59.2%. The genome features were analyzed with different databases, which revealed 2471 Gene Ontology/GO (11.72%), 5418 KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthologous/KO (25.71%), 13,913 Cluster of Orthologous Groups of proteins/COG (66.02%), 60 ABC transporter (0.28%), 1049 Carbohydrate-Active Enzymes/CAZy (4.98%), 4005 pathogen-host interactions/PHI (19%), and 515 fungal transcription factor/FTFD (2.44%) genes. The results obtained in this study provide deep insight for further studies in the future.

A multiomic characterization of the leukemia cell line REH using short- and long-read sequencing.

The B-cell acute lymphoblastic leukemia (ALL) cell line REH, with the t(12;21) ETV6::RUNX1 translocation, is known to have a complex karyotype defined by a series of large-scale chromosomal rearrangements. Taken from a 15-yr-old at relapse, the cell line offers a practical model for the study of pediatric B-ALL. In recent years, short- and long-read DNA and RNA sequencing have emerged as a complement to karyotyping techniques in the resolution of structural variants in an oncological context. Here, we explore the integration of long-read PacBio and Oxford Nanopore whole-genome sequencing, IsoSeq RNA sequencing, and short-read Illumina sequencing to create a detailed genomic and transcriptomic characterization of the REH cell line. Whole-genome sequencing clarified the molecular traits of disrupted ALL-associated genes including CDKN2A, PAX5, BTG1, VPREB1, and TBL1XR1, as well as the glucocorticoid receptor NR3C1 Meanwhile, transcriptome sequencing identified seven fusion genes within the genomic breakpoints. Together, our extensive whole-genome investigation makes high-quality open-source data available to the leukemia genomics community.

Copy number variation and elevated genetic diversity at immune trait loci in Atlantic and Pacific herring.

BACKGROUND: Genome-wide comparisons of populations are widely used to explore the patterns of nucleotide diversity and sequence divergence to provide knowledge on how natural selection and genetic drift affect the genome. In this study we have compared whole-genome sequencing data from Atlantic and Pacific herring, two sister species that diverged about 2 million years ago, to explore the pattern of genetic differentiation between the two species. RESULTS: The genome comparison of the two species revealed high genome-wide differentiation but with islands of remarkably low genetic differentiation, as measured by an FST analysis. However, the low FST observed in these islands is not caused by low interspecies sequence divergence (dxy) but rather by exceptionally high estimated intraspecies nucleotide diversity (π). These regions of low differentiation and elevated nucleotide diversity, termed high-diversity regions in this study, are not enriched for repeats but are highly enriched for immune-related genes. This enrichment includes genes from both the adaptive immune system, such as immunoglobulin, T-cell receptor and major histocompatibility complex genes, as well as a substantial number of genes with a role in the innate immune system, e.g. novel immune-type receptor, tripartite motif and tumor necrosis factor receptor genes. Analysis of long-read based assemblies from two Atlantic herring individuals revealed extensive copy number variation in these genomic regions, indicating that the elevated intraspecies nucleotide diversities were partially due to the cross-mapping of short reads. CONCLUSIONS: This study demonstrates that copy number variation is a characteristic feature of immune trait loci in herring. Another important implication is that these loci are blind spots in classical genome-wide screens for genetic differentiation using short-read data, not only in herring, likely also in other species harboring qualitatively similar variation at immune trait loci. These loci stood out in this study because of the relatively high genome-wide baseline for FST values between Atlantic and Pacific herring.

A multigenerational study can detect the evolutionary response to BaP exposure in the non-biting freshwater midge Chironomus riparius.

The release of polycyclic aromatic hydrocarbons (PAHs) into the environment is posing a threat to ecosystems and human health. Benzo(a)pyrene (BaP) is considered a biomarker of PAH exposure and is classified as a Group 1 carcinogen. However, it was not known whether BaP is mutagenic, i.e. induces inherited germline mutations. In this study, we used a recently established method, which combines short-term mutation accumulation lines (MAL) with whole genome sequencing (WGS) to assess mutagenicity in the non-biting midge Chironomus riparius. The mutagenicity analysis was supplemented by an evaluation of the development of population fitness in three successive generations in the case of chronic exposure to BaP at a high concentration (100 µg/L). In addition, the level of ROS-induced oxidative stress was examined in vivo. Exposure to the higher BaP concentration led to an increase in germline mutations relative to the control, while the lower concentration showed no mentionable effect. Against expectations, BaP exposure decreased ROS-level compared to the control and is thus probably not responsible for the increased mutation rate. Likewise, the higher BaP concentration decreased fitness measured as population growth rate per day (PGR) significantly over all generations, without signs of rapid evolutionary adaptations. Our results thus highlighted that high BaP exposure may influence the evolutionary trajectory of organisms.

Whole genome and transcriptome integrated analyses guide clinical care of pediatric poor prognosis cancers.

The role for routine whole genome and transcriptome analysis (WGTA) for poor prognosis pediatric cancers remains undetermined. Here, we characterize somatic mutations, structural rearrangements, copy number variants, gene expression, immuno-profiles and germline cancer predisposition variants in children and adolescents with relapsed, refractory or poor prognosis malignancies who underwent somatic WGTA and matched germline sequencing. Seventy-nine participants with a median age at enrollment of 8.8 y (range 6 months to 21.2 y) are included. Germline pathogenic/likely pathogenic variants are identified in 12% of participants, of which 60% were not known prior. Therapeutically actionable variants are identified by targeted gene report and whole genome in 32% and 62% of participants, respectively, and increase to 96% after integrating transcriptome analyses. Thirty-two molecularly informed therapies are pursued in 28 participants with 54% achieving a clinical benefit rate; objective response or stable disease ≥6 months. Integrated WGTA identifies therapeutically actionable variants in almost all tumors and are directly translatable to clinical care of children with poor prognosis cancers.

SARS-CoV-2 shedding, infectivity, and evolution in an immunocompromised adult patient.

This study aimed to provide further insight into the evolutionary dynamics of SARS-CoV-2 by analyzing the case of a 40-year-old man who had previously undergone autologous hematopoietic stem cell transplantation due to a diffuse large B-cell lymphoma. He developed a persistent SARS-CoV-2 infection lasting at least 218 days and did not manifest a humoral immune response to the virus during this follow-up period. Whole-genome sequencing and viral cultures confirmed a persistent infection with a replication-positive virus that had undergone genetic variation for at least 196 days after symptom onset.

Single-cell mtDNA dynamics in tumors is driven by coregulation of nuclear and mitochondrial genomes.

The extent of cell-to-cell variation in tumor mitochondrial DNA (mtDNA) copy number and genotype, and the phenotypic and evolutionary consequences of such variation, are poorly characterized. Here we use amplification-free single-cell whole-genome sequencing (Direct Library Prep (DLP+)) to simultaneously assay mtDNA copy number and nuclear DNA (nuDNA) in 72,275 single cells derived from immortalized cell lines, patient-derived xenografts and primary human tumors. Cells typically contained thousands of mtDNA copies, but variation in mtDNA copy number was extensive and strongly associated with cell size. Pervasive whole-genome doubling events in nuDNA associated with stoichiometrically balanced adaptations in mtDNA copy number, implying that mtDNA-to-nuDNA ratio, rather than mtDNA copy number itself, mediated downstream phenotypes. Finally, multimodal analysis of DLP+ and single-cell RNA sequencing identified both somatic loss-of-function and germline noncoding variants in mtDNA linked to heteroplasmy-dependent changes in mtDNA copy number and mitochondrial transcription, revealing phenotypic adaptations to disrupted nuclear/mitochondrial balance.

Accelerated somatic mutation calling for whole-genome and whole-exome sequencing data from heterogenous tumor samples.

Accurate detection of somatic mutations in DNA sequencing data is a fundamental prerequisite for cancer research. Previous analytical challenges were overcome by consensus mutation calling from four to five popular callers. This, however, increases the already nontrivial computing time from individual callers. Here, we launch MuSE 2, powered by multistep parallelization and efficient memory allocation, to resolve the computing time bottleneck. MuSE 2 speeds up 50 times more than MuSE 1 and eight to 80 times more than other popular callers. Our benchmark study suggests combining MuSE 2 and the recently accelerated Strelka2 achieves high efficiency and accuracy in analyzing large cancer genomic data sets.

Exploring the tumor genomic landscape of aggressive prostate cancer by whole-genome sequencing of tissue or liquid biopsies.

Treatment resistance remains a major issue in aggressive prostate cancer (PC), and novel genomic biomarkers may guide better treatment selection. Circulating tumor DNA (ctDNA) can provide minimally invasive information about tumor genomes, but the genomic landscape of aggressive PC based on whole-genome sequencing (WGS) of ctDNA remains incompletely characterized. Thus, we here performed WGS of tumor tissue (n = 31) or plasma ctDNA (n = 10) from a total of 41 aggressive PC patients, including 11 hormone-naïve, 15 hormone-sensitive, and 15 castration-resistant patients. Across all variant types, we found progressively more altered tumor genomic profiles in later stages of aggressive PC. The potential driver genes most frequently affected by single-nucleotide variants or insertions/deletions included the known PC-related genes TP53, CDK12, and PTEN and the novel genes COL13A1, KCNH3, and SENP3. Etiologically, aggressive PC was associated with age-related and DNA repair-related mutational signatures. Copy number variants most frequently affected 14q11.2 and 8p21.2, where no well-recognized PC-related genes are located, and also frequently affected regions near the known PC-related genes MYC, AR, TP53, PTEN, and BRCA1. Structural variants most frequently involved not only the known PC-related genes TMPRSS2 and ERG but also the less extensively studied gene in this context, PTPRD. Finally, clinically actionable variants were detected throughout all stages of aggressive PC and in both plasma and tissue samples, emphasizing the potential clinical applicability of WGS of minimally invasive plasma samples. Overall, our study highlights the feasibility of using liquid biopsies for comprehensive genomic characterization as an alternative to tissue biopsies in advanced/aggressive PC.

Outbreak of Streptococcus pyogenes emm89 ST646 in a head and neck surgical oncology ward.

Streptococcus pyogenes causes a variety of human infections, and hospital outbreaks with this pathogen have also been reported. The purpose of this study is to describe the clinical characteristics of an outbreak of S. pyogenes involving 15 patients and four healthcare workers (HCWs), as well as the molecular characteristics of the causative isolates. The course and response to the outbreak were reviewed, and information on the characteristics of the patients was extracted retrospectively from the medical records. Whole-genome sequencing of the 16 causative isolates (14 from patients and two from HCWs) was also performed. All 15 patients were postoperative of head and neck cancer with tracheotomy, and 12 had invasive infections, primarily surgical site infections, all of which resolved without causing serious illness. All but the first case was detected more than 7 days after admission. S. pyogenes was detected in two patients after empiric antimicrobial administration was performed on all inpatients and HCWs, and the outbreak was finally contained in approximately 2 months. All isolates detected in patients and HCWs belonged to emm89/clade 3, a hypervirulent clone that has emerged worldwide and was classified as sequence type 646. These isolates had single nucleotide polymorphism (SNP) differences of zero to one, indicating clonal transmission. This study demonstrated an outbreak of S. pyogenes emm89/clade 3 in a ward of patients with head and neck cancer. The global emergence of hypervirulent isolates may increase the risk of outbreaks among high-risk patients. IMPORTANCE: This study describes an outbreak of Streptococcus pyogenes that occurred in a ward caring for patients with head and neck cancer and tracheostomies. Many cases of invasive infections occurred in a short period, and extensive empiric antimicrobial administration on patients and healthcare workers was performed to control the outbreak. Whole-genome sequencing analysis of the causative strains confirmed that it was a monoclonal transmission of strains belonging to emm89/clade 3. The epidemiology and clinical characteristics of S. pyogenes infections have changed with the replacement of the prevalent clones worldwide. In the 1980s, there was a reemergence of S. pyogenes infections in high-income countries due to the spread of hypervirulent emm1 strains. emm89/clade 3 has recently been spreading worldwide and shares common features with emm1, including increased production of two toxins, NADase, and streptolysin O. The outbreak reported here may reflect the high spreading potential and virulence of emm89/clade 3.

Genomic characterization of extended-spectrum ß-lactamase-producing Enterobacterales isolated from abdominal surgical patients.

Rectal swabs of 104 patients who underwent abdominal surgery were screened for ESBL producers. Sequence types (STs) and resistance genes were identified by whole-genome sequencing of 46 isolates from 17 patients. All but seven isolates were assigned to recognized STs. While 18 ESBL-producing E. coli (EPEC) strains were of unique STs, ESBL-producing K. pneumoniae (EPKP) strains were mainly ST14 or ST15. Eight patients harboured strains of the same ST before and after abdominal surgery. The most prevalent resistant genes in E. coli were blaEC (69.57%), blaCTX-M (65.22%), and blaTEM (36.95%), while blaSHV was present in only K. pneumoniae (41.30%). Overall, genes encoding ß-lactamases of classes A (blaCTX-M, blaTEM, blaZ), C (blaSHV, blaMIR, and blaDHA), and D (blaOXA) were identified, the most prevalent variants being blaCTX-M-15, blaTEM-1B, blaSHV-28, and blaOXA-1. Interestingly, blaCMY-2, the most common pAmpC ß-lactamase genes reported worldwide, and mobile colistin resistance genes, mcr-10-1, were also identified. The presence of blaCMY-2 and mcr-10-1 is concerning as they may constitute a potentially high risk of pan-resistant post-surgical infections. It is imperative that healthcare professionals monitor intra-abdominal surgical site infections rigorously to prevent transmission of faecal ESBL carriage in high-risk patients.

Genomic characterization of human respiratory syncytial virus circulating in Islamabad, Pakistan, during an outbreak in 2022-2023.

In this study, conducted at the National Institute of Health, Islamabad, during an outbreak of human respiratory syncytial virus (hRSV) from December 2022 to January 2023, the first whole-genome sequences of hRSV isolates from Islamabad, Pakistan, were determined. Out of 10 positive samples, five were sequenced, revealing the presence of two genotypes: RSV-A (GA2.3.5, ON1 strain) and RSV-B (GB5.0.5.a, BA-10 strain). A rare non-synonymous substitution (E232G) in G the protein and N276S in the F protein were found in RSV-A. In RSV-B, the unique mutations K191R, Q209R, and I206M were found in the F protein. These mutations could potentially influence vaccine efficacy and viral pathogenicity. This research underscores the importance of genomic surveillance for understanding RSV diversity and guiding public health responses in Pakistan.

A novel fast hybrid capture sequencing method for high-efficiency common human coronavirus whole-genome acquisition.

Rapid and accurate sequencing of the entire viral genome, coupled with continuous monitoring of genetic changes, is crucial for understanding the epidemiology of coronaviruses. We designed a novel method called micro target hybrid capture system (MT-Capture) to enable whole-genome sequencing in a timely manner. The novel design of probes used in target binding exhibits a unique and synergistic "hand-in-hand" conjugation effect. The entire hybrid capture process is within 2.5 hours, overcoming the time-consuming and complex operation characteristics of the traditional liquid-phase hybrid capture (T-Capture) system. By designing specific probes for these coronaviruses, MT-Capture effectively enriched isolated strains and 112 clinical samples of coronaviruses with cycle threshold values below 37. Compared to multiplex PCR sequencing, it does not require frequent primer updates and has higher compatibility. MT-Capture is highly sensitive and capable of tracking variants.IMPORTANCEMT-Capture is meticulously designed to enable the efficient acquisition of the target genome of the common human coronavirus. Coronavirus is a kind of virus that people are generally susceptible to and is epidemic and infectious, and it is the virus with the longest genome among known RNA viruses. Therefore, common human coronavirus samples are selected to evaluate the accuracy and sensitivity of MT-Capture. This method utilizes innovative probe designs optimized through probe conjugation techniques, greatly shortening the time and simplifying the handwork compared with traditional hybridization capture processes. Our results demonstrate that MT-Capture surpasses multiplex PCR in terms of sensitivity, exhibiting a thousandfold increase. Moreover, MT-Capture excels in the identification of mutation sites. This method not only is used to target the coronaviruses but also may be used to diagnose other diseases, including various infectious diseases, genetic diseases, or tumors.

Whole genome sequencing of recombinant viruses obtained from co-infection and superinfection of Vero cells with modified vaccinia virus ankara vectored influenza vaccine and a naturally occurring cowpox virus.

Modified vaccinia virus Ankara (MVA) has been widely tested in clinical trials as recombinant vector vaccine against infectious diseases and cancers in humans and animals. However, one biosafety concern about the use of MVA vectored vaccine is the potential for MVA to recombine with naturally occurring orthopoxviruses in cells and hosts in which it multiplies poorly and, therefore, producing viruses with mosaic genomes with altered genetic and phenotypic properties. We previously conducted co-infection and superinfection experiments with MVA vectored influenza vaccine (MVA-HANP) and a feline Cowpox virus (CPXV-No-F1) in Vero cells (that were semi-permissive to MVA infection) and showed that recombination occurred in both co-infected and superinfected cells. In this study, we selected the putative recombinant viruses and performed genomic characterization of these viruses. Some putative recombinant viruses displayed plaque morphology distinct of that of the parental viruses. Our analysis demonstrated that they had mosaic genomes of different lengths. The recombinant viruses, with a genome more similar to MVA-HANP (>50%), rescued deleted and/or fragmented genes in MVA and gained new host ranges genes. Our analysis also revealed that some MVA-HANP contained a partially deleted transgene expression cassette and one recombinant virus contained part of the transgene expression cassette similar to that incomplete MVA-HANP. The recombination in co-infected and superinfected Vero cells resulted in recombinant viruses with unpredictable biological and genetic properties as well as recovery of delete/fragmented genes in MVA and transfer of the transgene into replication competent CPXV. These results are relevant to hazard characterization and risk assessment of MVA vectored biologicals.

Anaerolentibacter hominis gen. nov. sp. nov., Diplocloster hominis sp. nov. and Pilosibacter fragilis gen. nov. sp. nov., isolated from human faeces.

Three Gram-positive, obligately anaerobic bacterial strains, namely CSJ-1T, CSJ-3T, and CSJ-4T, were isolated from faeces of healthy persons. They were characterized through a combination of whole-genome sequencing, phenotypic traits, and metabolomic analysis. The genome sizes of CSJ-1T, CSJ-4T, and CSJ-3T were 3.3, 3.8, and 6.1 Mbp, with DNA G+C contents of 47.2, 48.3, and 48.8 mol%, respectively. Strain CSJ-3T was identified as representing a novel species, Diplocloster hominis (type strain CSJ-3T=CGMCC 1.18033T=JCM 36512T) of the genus Diplocloster. The 16S rRNA gene sequence similarity and whole genome average nucleotide identity (gANI) of CSJ-4T to its closest related species, Diplocloster modestus ASD 4241T, were 98.3 and 91.4 %, respectively. Comparative analysis of 16S rRNA gene sequences showed 91.6 % similarity between CSJ-1T and its closest phylogenetic neighbour, Catenibacillus scindens DSM 106146T, and 93.3 % similarity between CSJ-4T and its closest relative strain, Clostridium fessum SNUG30386T. Based on the polyphasic taxonomic results, we proposed two novel genera and three novel species. Strain CSJ-1T was identified as representing a novel species of novel genus, Anaerolentibacter hominis gen. nov. sp. nov. (type strain CSJ-1T=CGMCC 1.18046T=JCM 36511T) of the family Lachnospiraceae, and strain CSJ-4T was identified as representing a novel species of novel genus Pilosibacter fragilis gen. nov. sp. nov. (type strain CSJ-4T=CGMCC 1.18026T= JCM 36513T) of the family Clostridiaceae.