Efficient SARS-CoV-2 variant detection and monitoring with Spike Screen next-generation sequencing.

The emergence and rapid spread of SARS-CoV-2 prompted the global community to identify innovative approaches to diagnose infection and sequence the viral genome because at several points in the pandemic positive case numbers exceeded the laboratory capacity to characterize sufficient samples to adequately respond to the spread of emerging variants. From week 10, 2020, to week 13, 2023, Slovenian routine complete genome sequencing (CGS) surveillance network yielded 41 537 complete genomes and revealed a typical molecular epidemiology with early lineages gradually being replaced by Alpha, Delta, and finally Omicron. We developed a targeted next-generation sequencing based variant surveillance strategy dubbed Spike Screen through sample pooling and selective SARS-CoV-2 spike gene amplification in conjunction with CGS of individual cases to increase throughput and cost-effectiveness. Spike Screen identifies variant of concern (VOC) and variant of interest (VOI) signature mutations, analyses their frequencies in sample pools, and calculates the number of VOCs/VOIs at the population level. The strategy was successfully applied for detection of specific VOC/VOI mutations prior to their confirmation by CGS. Spike Screen complemented CGS efforts with an additional 22 897 samples sequenced in two time periods: between week 42, 2020, and week 24, 2021, and between week 37, 2021, and week 2, 2022. The results showed that Spike Screen can be applied to monitor VOC/VOI mutations among large volumes of samples in settings with limited sequencing capacity through reliable and rapid detection of novel variants at the population level and can serve as a basis for public health policy planning.

Comparative genomic analyses provide insight into the pathogenicity of three Pseudomonas syringae pv. actinidiae strains from Anhui Province, China.

BACKGROUND: Pseudomonas syringae pv. actinidiae (Psa) is an important bacterial plant pathogen that causes severe damage to the kiwifruit industry worldwide. Three Psa strains were recently obtained from different kiwifruit orchards in Anhui Province, China. The present study mainly focused on the variations in virulence and genome characteristics of these strains based on the pathogenicity assays and comparative genomic analyses. RESULTS: Three strains were identified as biovar 3 (Psa3), along with strain QSY6 showing higher virulence than JZY2 and YXH1 in pathogenicity assays. The whole genome assembly revealed that each of the three strains had a circular chromosome and a complete plasmid. The chromosome sizes ranged from 6.5 to 6.6 Mb with a GC content of approximately 58.39 to 58.46%, and a predicted number of protein-coding sequences ranging from 5,884 to 6,019. The three strains clustered tightly with 8 Psa3 reference strains in terms of average nucleotide identity (ANI), whole-genome-based phylogenetic analysis, and pangenome analysis, while they were evolutionarily distinct from other biovars (Psa1 and Psa5). Variations were observed in the repertoire of effectors of the type III secretion system among all 15 strains. Moreover, synteny analysis of the three sequenced strains revealed eight genomic regions containing 308 genes exclusively present in the highly virulent strain QSY6. Further investigation of these genes showed that 16 virulence-related genes highlight several key factors, such as effector delivery systems (type III secretion systems) and adherence (type IV pilus), which might be crucial for the virulence of QSY6. CONCLUSION: Three Psa strains were identified and showed variant virulence in kiwifruit plant. Complete genome sequences and comparative genomic analyses further provided a theoretical basis for the potential pathogenic factors responsible for kiwifruit bacterial canker.

New insight into protective effect against oxidative stress and biosynthesis of exopolysaccharides produced by Lacticaseibacillus paracasei NC4 from fermented eggplant.

Fermented eggplant is a traditional fermented food, however lactic acid bacteria capable of producing exopolysaccharide (EPS) have not yet been exploited. The present study focused on the production and protective effects against oxidative stress of an EPS produced by Lacticaseibacillus paracasei NC4 (NC4-EPS), in addition to deciphering its genomic features and EPS biosynthesis pathway. Among 54 isolates tested, strain NC4 showed the highest EPS yield and antioxidant activity. The maximum EPS production (2.04 ± 0.11 g/L) was achieved by culturing in MRS medium containing 60 g/L sucrose at 37 °C for 48 h. Under 2 mM H2O2 stress, the survival of a yeast model Saccharomyces cerevisiae treated with 0.4 mg/mL NC4-EPS was 2.4-fold better than non-treated cells, which was in agreement with the catalase and superoxide dismutase activities measured from cell lysates. The complete genome of NC4 composed of a circular chromosome of 2,888,896 bp and 3 circular plasmids. The NC4 genome comprises more genes with annotated function in nitrogen metabolism, phosphorus metabolism, cell division and cell cycle, and iron acquisition and metabolism as compared to other reported L. paracasei. Of note, the eps gene cluster is not conserved across L. paracasei. Pathways of sugar metabolism for EPS biosynthesis were proposed for the first time, in which gdp pathway only present in few plant-derived bacteria was identified. These findings shed new light on the cell-protective activity and biosynthesis of EPS produced by L. paracasei, paving the way for future efforts to enhance yield and tailor-made EPS production for food and pharmaceutical industries.

In-depth comparative pathogenome, virulome, and resistome analysis of an extensive drug resistant Ralstonia mannitolilytica strain isolated from blood.

INTRODUCTION: Ralstonia mannitolilytica is an global opportunistic pathogen responsible for various diseases. In this study, we reported the genome of a R. mannitolilytica isolate responsible for bacteremia in an acute exacerbation of chronic obstructive pulmonary disease (AECOPD). METHODS: Bacterial identification was performed with a Vitek2™ Automated System and 16S rRNA sequencing with BLASTn against the Non-Redundant Protein Sequence (Nr) database. Genome sequencing and analysis were performed using PacBio RS II sequencer, Hierarchical Genome Assembly Process assembly, as well as multiple annotation databases to better understand the innate features. Antibiotic resistance genes and virulence factors were specifically identified through Antibiotic Resistance Genes database and Virulence Factors of Pathogenic Bacteria databases. RESULTS: The complete genome sequence was assembled into two chromosomes with 3,495,817 bp and 1,342,871 bp in length and GC% of 65.37 % and 66.43 %, respectively. The two chromosomes were fully annotated. In chromosome 1 and 2, 19 and 14 antibiotic resistant genes and 48 and 55 virulence factors were predicted, respectively. Specifically, beta-lactam resistance genes blaOXA-443, blaOXA-444 were acquired. CONCLUSIONS: This study aids in the understanding of the innate features of R. mannitolilytica in AECOPD.

Isolation, identification, and pathogenicity of two novel goose astrovirus from goslings with severe gout in China.

Goose astroviruses (GAstVs) are important pathogens which can cause gout in goslings leading to huge economic losses for the goose farming industry in China. In 2023, an infectious disease characterized by visceral gout broke out in commercial goose farms in Guangxi and Guangdong provinces of China. In this study, two GAstV strains of GXNN and GDCS were successfully isolated from these two disease-ridden goose farms. The complete genomic lengths of these two strains were 7166 bp, and phylogenetic analysis showed that they were both GAstV-2 subtypes. The 3-dimensional structures of the capsid protein were predicted and six characteristic mutation sites at amino acid positions 60, 61, 228, 229, 456 and 523 were found within the strong antigenic regions. A recombination event occurred at 6833-7070 nt between the GAstV TZ03 and Turkey astrovirus CA/00 and this was detected in both the GXNN and GDCS strains. Another recombinant event occurred at 63-2747 nt between the GAstV XT1 and GAstV SDPY and this was detected in the GDCS strain. When 1-day-old goslings were infected with the novel GXNN and GDCS strains, they showed severe visceral gout. This was accompanied by enlarged spleens, liver hemorrhages and urate deposits in the kidneys and ureters and their blood urea nitrogen levels were significantly elevated. The mortality rates of the GXNN- and GDCS-infected groups were pathogenically high at 80 % and 60 %, respectively. These results will promote our understanding of the evolution and epidemic potential of GAstVs in China.

Proteus appendicitidis sp. nov., isolated from the appendiceal pus of an appendicitis patient in Yongzhou, China.

A Gram-negative, facultatively anaerobic, short rod-shaped bacterium, designated as strain HZ0627T, was isolated from the appendiceal pus of a patient with appendicitis in Yongzhou, Hunan, China. This strain was subjected to comprehensive phenotypic, phylogenetic, and genomic analyses using polyphasic taxonomic methods. Phylogenetic analysis of the 16S rRNA gene sequence revealed that this strain belonged to the genus Proteus and the family Morganellaceae, whereas that based on the rpoB gene sequence and phylogenomic analysis demonstrated that this strain was distinctly separated from other type strains of Proteus species. Moreover, whole-genome-based analyses, including in silico DNA-DNA hybridization (isDDH) and average nucleotide identity (ANI), revealed that strain HZ0627T had much lower isDDH rates (24.5-55.6%) and ANI (82.04-93.90%) than those of the thresholds (i.e., 70% and 95%, respectively) for species delineation, when compared to the type strains of other Proteus species. The cellular fatty acid profile of strain HZ0627T was dominated by C16:0 (34.5%), cyclo C17:0 (25.8%), C14:0 (12.6%), C16:1 iso I/14:0 3-OH (7.7%), C18:1ω7c/18:1ω6c (6.5%), and C16:1ω7c/16:1ω6c (4.9%), which clearly differentiated it from the documented type strains of Proteus species. In addition, several specific physiological traits, including optimal growth temperature, tolerance to sodium chloride, and carbon source utilization, differed from those of other Proteus species. Therefore, we propose the name Proteus appendicitidis sp. nov. for strain HZ0627T (= CCTCC AB 2022380T = KCTC 92986T), which represents the type strain of this novel Proteus species.

Analysis of antimicrobial biological activity of a marine Bacillus velezensis NDB.

A new strain of Bacillus velezensis NDB was isolated from Xiangshan Harbor and antibacterial test revealed antibacterial activity of this strain against 12 major pathogenic bacteria. The whole genome of the bacterium was sequenced and found to consist of a 4,214,838 bp circular chromosome and a 7410 bp circular plasmid. Furthermore, it was predicted by AntiSMASH and BAGEL4 to have 12 clusters of secondary metabolism genes for the synthesis of the inhibitors, fengycin, bacillomycin, macrolactin H, bacillaene, and difficidin, and there were also five clusters encoding potentially novel antimicrobial substances, as well as three bacteriocin biosynthesis gene clusters of amylocyclicin, ComX1, and LCI. qRT-PCR revealed significant up-regulation of antimicrobial secondary metabolite synthesis genes after 24 h of antagonism with pathogenic bacteria. Furthermore, MALDI-TOF mass spectrometry revealed that it can secrete surfactin non-ribosomal peptide synthase and polyketide synthase to exert antibacterial effects. GC-MS was used to analyze methanol extract of B. velezensis NDB, a total of 68 compounds were identified and these metabolites include 16 amino acids, 17 acids, 3 amines, 11 sugars, 11 alcohols, 1 ester, and 9 other compounds which can inhibit pathogenic bacteria by initiating the antibiotic secretion pathway. A comparative genomic analysis of gene families showed that the specificity of B. velezensis NDB was mainly reflected in environmental adaptability. Overall, this research on B. velezensis NDB provides the basis for elucidating its biocontrol effect and promotes its future application as a probiotic.

Microbulbifer bruguierae sp. nov., isolated from sediment of mangrove plant Bruguiera sexangula, and comparative genomic analyses of the genus Microbulbifer.

A Gram-negative, non-motile, strictly aerobic, rod-shaped bacterium, designated as H12T, was isolated from the sediments of mangrove plant Bruguiera sexangula taken from Dapeng district, Shenzhen, PR China. The pairwise 16S rRNA gene sequence analysis showed that strain H12T shared high identity levels with species of the genus Microbulbifer, with the highest similarity level of 98.5 % to M. pacificus SPO729T, followed by 98.1 % to M. donghaiensis CN85T. Phylogenetic analysis using core-genome sequences showed that strain H12T formed a cluster with type species of M. pacificus SPO729T and M. harenosus HB161719T. The complete genome of strain H12T was 4 481 396 bp in size and its DNA G+C content was 56.7 mol%. The average nucleotide identity and digital DNA-DNA hybridization values among strain H12T and type species of genus Microbulbifer were below the cut-off levels of 95-96 and 70 %, respectively. The predominant cellular fatty acids of strain H12T were iso-C15 : 0 (22.5 %) and C18 : 1 ω7c (13.9 %). Ubiquinone-8 was detected as the major respiratory quinone. The polar lipids of strain H12T comprised one phosphatidylglycerol, one phosphatidylethanolamine, one unidentified aminoglycophospholipid, one unidentified glycophospholipid, three unidentified glycolipids, two unidentified aminolipids, and one unidentified lipid. Based on polyphasic evidence, strain H12T represents a novel species of the genus Microbulbifer, for which the name Microbulbifer bruguierae sp. nov. is proposed. The type strain is H12T (=KCTC 92859T=MCCC 1K08451T). Comparative genomic analyses of strain H12T with strains of the genus Microbulbifer reveal its potential in degradation of pectin.

Complete genome sequence of a novel mitovirus identified in the phytopathogenic fungus Alternaria tenuissima.

In this study, a novel positive-sense single-stranded RNA (+ ssRNA) mycovirus, Alternaria tenuissima mitovirus 1 (AtMV1), was identified in Alternaria tenuissima strain YQ-2-1, a phytopathogenic fungus causing leaf blight on muskmelon. The genome of AtMV1 is a single RNA molecule that is 3013 nt in length with an A + U content of 66.58% and contains a single open reading frame (ORF) using the fungal mitochondrial genetic code. The ORF was predicted to encode a 313-amino-acid RNA-dependent RNA polymerase (RdRp) with a molecular mass of 35.48 kDa, which contains six conserved motifs with the highly conserved GDD tripeptide in motif IV. The 5' and 3' untranslated regions were predicted to fold into stem-loop and panhandle secondary structures. The results of a BLASTp search revealed that the amino acid (aa) sequence of RdRp of AtMV1 shared the highest sequence similarity (51.04% identity) with that of Sichuan mito-like virus 30, a member of the genus Duamitovirus within the family Mitoviridae. Phylogenetic analysis based on the aa sequence of the RdRp suggested that AtMV1 is a novel member of the genus Duamitovirus. To our knowledge, this is the first report of the complete genome sequence of a new mitovirus infecting A. tenuissima.

The full genome characterization of avian encephalomyelitis virus, Iran: a vertical transmission case.

Avian encephalomyelitis (AE) is an important infectious poultry disease worldwide that is caused by avian encephalomyelitis virus (AEV). The causative virus can be transmitted both horizontally and vertically. In the present study, an AEV suspected outbreak with typical neurological signs occurred in broilers. Histopathological examination, RT-PCR assay and full genome sequencing were applied to confirm the presence of AEV. Phylogenetic analysis of the full genome sequence showed that the detected AEV strain at 7055 nucleotide length is classified in cluster I and is closely related to vaccinal USA and China originated isolates. Although, the outbreaks of AE in progeny of vaccinated breeders have been reported previously, the source of infection was unknown. Based on the results obtained in this study, the outbreaks are vaccine-originated. This study provides the first whole genome analysis of AEV from Iran and reveals that the AEV possesses a hepatitis C virus-like internal ribosome entry site.

DsRNA sequencing revealed a previously missed terminal sequence of a +ssRNA virus that infects dinoflagellate Heterocapsa circularisquama.

Heterocapsa circularisquama RNA virus (HcRNAV) is the only dinoflagellate-infecting RNA virus cultured. However, only two strains of HcRNAV have been registered with complete genome sequences (strains 34 and 109 for UA and CY types, respectively). To extend the genomic information of HcRNAV, we performed full-genome sequencing of an unsequenced strain of HcRNAV (strain A8) using the fragmented and primer-ligated double-stranded RNA (dsRNA) sequencing (FLDS) method. The complete genome of HcRNAV A8 with 4457 nucleotides (nt) was successfully determined, and sequence alignment of the major capsid protein gene suggested that A8 was a UA-type strain, consistent with its intraspecific host specificity. The complete sequence was found to be 80 nt longer at the 5' terminus than the registered sequences of HcRNAV strains (34 and 109), suggesting that FLDS is more reliable for determining the terminal sequence than conventional methods (5' Rapid Amplification of cDNA End). Our study contributes to a better understanding of dinoflagellate-infecting viruses with limited sequence data.

Detection and genetic characterization of chicken astrovirus and avian nephritis virus from hatchery to farm.

RESEARCH HIGHLIGHTS: Identified the role of the hatchery in astrovirus transmission.Sequenced the avian nephritis virus complete genome.Investigated tissue distribution of astrovirus from egg to chicks.Demonstrated co-infection of ANV/CAstV.

Genomics and taxonomy of the glyphosate-degrading, copper-tolerant rhizospheric bacterium Achromobacter insolitus LCu2.

A rhizosphere strain, Achromobacter insolitus LCu2, was isolated from alfalfa (Medicago sativa L.) roots. It was able to degrade of 50% glyphosate as the sole phosphorus source, and was found resistant to 10 mM copper (II) chloride, and 5 mM glyphosate-copper complexes. Inoculation of alfalfa seedlings and potato microplants with strain LCu2 promoted plant growth by 30-50%. In inoculated plants, the toxicity of the glyphosate-copper complexes to alfalfa seedlings was decreased, as compared with the noninoculated controls. The genome of A. insolitus LCu2 consisted of one circular chromosome (6,428,890 bp) and encoded 5843 protein genes and 76 RNA genes. Polyphasic taxonomic analysis showed that A. insolitus LCu2 was closely related to A. insolitus DSM23807T on the basis of the average nucleotide identity of the genomes of 22 type strains and the multilocus sequence analysis. Genome analysis revealed genes putatively responsible for (1) plant growth promotion (osmolyte, siderophore, and 1-aminocyclopropane-1-carboxylate deaminase biosynthesis and auxin metabolism); (2) degradation of organophosphonates (glyphosate oxidoreductase and multiple phn clusters responsible for the transport, regulation and C-P lyase cleavage of phosphonates); and (3) tolerance to copper and other heavy metals, effected by the CopAB-CueO system, responsible for the oxidation of copper (I) in the periplasm, and by the efflux Cus system. The putative catabolic pathways involved in the breakdown of phosphonates are predicted. A. insolitus LCu2 is promising in the production of crops and the remediation of soils contaminated with organophosphonates and heavy metals.

Genomic and Comparative Transcriptomic Analyses Reveal Key Genes Associated with the Biosynthesis Regulation of Okaramine B in Penicillium daleae NBP-49626.

Restricted production of fungal secondary metabolites hinders the ability to conduct comprehensive research and development of novel biopesticides. Okaramine B from Penicillium demonstrates remarkable insecticidal efficacy; however, its biosynthetic yield is low, and its regulatory mechanism remains unknown. The present study found that the yield difference was influenced by fermentation modes in okaramine-producing strains and performed genomic and comparative transcriptome analysis of P. daleae strain NBP-49626, which exhibits significant features. The NBP-49626 genome is 37.4 Mb, and it encodes 10,131 protein-encoding genes. Up to 5097 differentially expressed genes (DEGs) were identified during the submerged and semi-solid fermentation processes. The oka gene cluster, lacking regulatory and transport genes, displayed distinct transcriptional patterns in response to the fermentation modes and yield of Okaramine B. Although transcription trends of most known global regulatory genes are inconsistent with those of oka, this study identified five potential regulatory genes, including two novel Zn(II)2Cys6 transcription factors, Reg2 and Reg19. A significant correlation was also observed between tryptophan metabolism and Okaramine B yields. In addition, several transporter genes were identified as DEGs. These results were confirmed using real-time quantitative PCR. This study provides comprehensive information regarding the regulatory mechanism of Okaramine B biosynthesis in Penicillium and is critical to the further yield improvement for the development of insecticides.

Does Salmonella diarizonae 58:r:z53 Isolated from a Mallard Duck Pose a Threat to Human Health?

Salmonella diarizonae (IIIb) is frequently isolated from reptiles and less frequently from birds and mammals. However, its isolation from invasive human infections has not been widely reported. Migratory mallard ducks are excellent bioindicators of pathogen presence and pathogen antibiotic resistance (AMR). We present the first isolation from a mallard duck in central Europe of the antibiotic-resistant Salmonella enterica subsp. diarizonae with the unique antigenic pattern 58:r:z53 and report its whole-genome sequencing, serosequencing, and genotyping, which enabled the prediction of its pathogenicity and comparison with phenotypic AMR. The isolated strain was highly similar to S. diarizonae isolated from humans and food. Twenty-four AMR genes were detected, including those encoding aminoglycoside, fluoroquinolone, macrolide, carbapenem, tetracycline, cephalosporin, nitroimidazole, peptide antibiotic, and disinfecting agent/antiseptic resistance. Six Salmonella pathogenicity islands were found (SPI-1, SPI-2, SPI-3, SPI-5, SPI-9, and SPI-13). An iron transport system was detected in SPI-1 centisome C63PI. Plasmid profile analyses showed three to be present. Sequence mutations in the invA and invF genes were noted, which truncated and elongated the proteins, respectively. The strain also harbored genes encoding type-III secretion-system effector proteins and many virulence factors found in S. diarizonae associated with human infections. This study aims to elucidate the AMR and virulence genes in S. enterica subsp. diarizonae that may most seriously threaten human health.

Analysis of codon usage patterns in 48 Aconitum species.

BACKGROUND: The Aconitum genus is a crucial member of the Ranunculaceae family. There are 350 Aconitum species worldwide, with about 170 species found in China. These species are known for their various pharmacological effects and are commonly used to treat joint pain, cold abdominal pain, and other ailments. Codon usage bias (CUB) analysis contributes to evolutionary relationships and phylogeny. Based on protein-coding sequences (PCGs), we selected 48 species of Aconitum for CUB analysis. RESULTS: The results revealed that Aconitum species had less than 50% GC content. Furthermore, the distribution of GC content was irregular and followed a trend of GC1 > GC2 > GC3, indicating a bias towards A/T bases. The relative synonymous codon usage (RSCU) heat map revealed the presence of conservative codons with slight variations within the genus. The effective number of codons (ENC)-Plot and the parity rule 2 (PR2)-bias plot analysis indicate that natural selection is the primary factor influencing the variation in codon usage. As a result, we screened various optimal codons and found that A/T bases were preferred as the last codon. Furthermore, our Maximum Likelihood (ML) analysis based on PCGs among 48 Aconitum species yielded results consistent with those obtained from complete chloroplast (cp.) genome data. This suggests that analyzing mutation in PCGs is an efficient method for demonstrating the phylogeny of species at the genus level. CONCLUSIONS: The CUB analysis of 48 species of Aconitum was mainly influenced by natural selection. This study reveals the CUB pattern of Aconitum and lays the foundation for future genetic modification and phylogenetic analyses.

Assessment of beneficial effects and identification of host adaptation-associated genes of Ligilactobacillus salivarius isolated from badgers.

BACKGROUND: Ligilactobacillus salivarius has been frequently isolated from the gut microbiota of humans and domesticated animals and has been studied as a candidate probiotic. Badger (Meles meles) is known as a "generalist" species that consumes complex foods and exhibits tolerance and resistance to certain pathogens, which can be partly attributed to the beneficial microbes such as L. salivarius in the gut microbiota. However, our understanding of the beneficial traits and genomic features of badger-originated L. salivarius remains elusive. RESULTS: In this study, nine L. salivarius strains were isolated from wild badgers' feces, one of which exhibited good probiotic properties. Complete genomes of the nine L. salivarius strains were generated, and comparative genomic analysis was performed with the publicly available complete genomes of L. salivarius obtained from humans and domesticated animals. The strains originating from badgers harbored a larger genome, a higher number of protein-coding sequences, and functionally annotated genes than those originating from humans and chickens. The pan-genome phylogenetic tree demonstrated that the strains originating from badgers formed a separate clade, and totally 412 gene families (12.6% of the total gene families in the pan-genome) were identified as genes gained by the last common ancestor of the badger group. The badger group harbored significantly more gene families responsible for the degradation of complex carbohydrate substrates and production of polysaccharides than strains from other hosts; many of these were acquired by gene gain events. CONCLUSIONS: A candidate probiotic and nine L. salivarius complete genomes were obtained from the badgers' gut microbiome, and several beneficial genes were identified to be specifically present in the badger-originated strains that were gained in the evolution. Our study provides novel insights into the adaptation of L. salivarius to the intestinal habitat of wild badgers and provides valuable strain and genome resources for the development of L. salivarius as a probiotic.

Sequencing and characterization of Helcococcus ovis: a comprehensive comparative genomic analysis of virulence.

BACKGROUND: Helcococcus ovis (H. ovis) is an emerging bacterial pathogen that commonly causes opportunistic respiratory, mammary, and uterine infections across mammalian hosts. This study applied long- and short-read whole genome sequencing technologies to identify virulence factors in five H. ovis isolates with low, medium, and high virulence phenotypes. RESULTS: The resulting assemblies contained one circular chromosome ranging from 1,744,566 to 1,850,083 bp in length and had a mean GC content of 27.6%. Phylogenetic and nucleotide identity analyses found low virulence strain KG38 to be part of a clade that forms an outgroup apart from the rest of the H. ovis taxon. Assembling the first complete genomes of the species revealed major genomic rearrangements in KG38. One to six prophage regions were identified in each genome. A novel pathogenicity island was found exclusively in the two high virulence strains (KG37 and KG104), along with two hypothetical transmembrane proteins designated as putative VFs. Finally, three zinc ABC transporters and three Type-II/IV secretion systems were identified as possible virulence determinants in this species. The low virulence strain KG38 has fewer intact paralogs of these operons in its genome compared to the higher virulence isolates, which strongly suggests a role in virulence. This strain is also missing four putative virulence factors (VFs) found in other isolates associated with adherence (collagen adhesin precursor), immune evasion (choline-binding protein A and a PspA-like hypothetical protein) and cell wall synthesis (glycerol-3-phosphate cytidylyltransferase). CONCLUSIONS: In this study, we assembled reference-quality complete genomes for five H. ovis strains to identify putative virulence factors. Phylogenetic analyses of H. ovis isolates revealed the presence of a clade representing a potentially novel species within the genus Helcococcus. A novel pathogenicity island and two hypothetical transmembrane proteins were found exclusively in high-virulence strains. The identification of Zinc ABC transporters and Type-II/IV secretion systems as possible virulence determinants, along with the differences in operon content between the low and high virulence isolates, strongly suggests they also play a role in the bacterium's pathogenicity. Taken together, these findings are a valuable first step toward deciphering the pathogenesis of H. ovis infections.

A near-complete genome assembly of Brassica rapa provides new insights into the evolution of centromeres.

Brassica rapa comprises many important cultivated vegetables and oil crops. However, Chiifu v3.0, the current B. rapa reference genome, still contains hundreds of gaps. Here, we presented a near-complete genome assembly of B. rapa Chiifu v4.0, which was 424.59 Mb with only two gaps, using Oxford Nanopore Technology (ONT) ultralong-read sequencing and Hi-C technologies. The new assembly contains 12 contigs, with a contig N50 of 38.26 Mb. Eight of the ten chromosomes were entirely reconstructed in a single contig from telomere to telomere. We found that the centromeres were mainly invaded by ALE and CRM long terminal repeats (LTRs). Moreover, there is a high divergence of centromere length and sequence among B. rapa genomes. We further found that centromeres are enriched for Copia invaded at 0.14 MYA on average, while pericentromeres are enriched for Gypsy LTRs invaded at 0.51 MYA on average. These results indicated the different invasion mechanisms of LTRs between the two structures. In addition, a novel repetitive sequence PCR630 was identified in the pericentromeres of B. rapa. Overall, the near-complete genome assembly, B. rapa Chiifu v4.0, offers valuable tools for genomic and genetic studies of Brassica species and provides new insights into the evolution of centromeres.

Complete Genome Sequencing Reveals Unusual Equine Rotavirus A of Bat Origin from India.

Rotaviruses are the most common viral agents associated with foal diarrhea. Between 2014 and 2017, the annual prevalence of rotavirus in diarrheic foals ranged between 18 and 28% in Haryana (India). Whole-genome sequencing of two equine rotavirus A (ERVA) isolates (RVA/Horse-wt/IND/ERV4/2017 and RVA/Horse-wt/IND/ERV6/2017) was carried out to determine the genotypic constellations (GCs) of ERVAs. The GCs of both the isolates were G3-P[3]-I8-R3-C3-M3-A9-N3-T3-E3-H6, a unique combination reported for ERVAs so far. Both the isolates carried VP6 of genotype I8, previously unreported from equines. Upon comparison with RVAs of other species, the GC of both isolates was identical to that of a bat rotavirus strain, MSLH14, isolated from China in 2012. The nucleotide sequences of the genes encoding VP3, NSP2, and NSP3 shared >95.81% identity with bat RVA strains isolated from Africa (Gabon). The genes encoding VP1, VP2, VP7, NSP1, and NSP4 shared 94.82% to 97.12% nucleotide identities with the human strains which have zoonotic links to bats (RCH272 and MS2015-1-0001). The VP6 genes of both strains were distinct and had the highest similarity of only 87.08% with that of CMH222, a human strain of bat origin. The phylogenetic analysis and lineage studies revealed that VP7 of both isolates clustered in a new lineage (lineage X) of the G3 genotype with bat, human, and alpaca strains. Similarly, VP4 clustered in a distinct P[3] lineage. These unusual findings highlight the terra incognita of the genomic diversity of equine rotaviruses and support the need for the surveillance of RVAs in animals and humans with a "one health" approach. IMPORTANCE Rotaviruses are globally prevalent diarrheal pathogens in young animals including foals, piglets, calves, goats, sheep, cats, and dogs along with humans. The genome of rotaviruses consists of 11 segments, which enables them to undergo reshuffling by reassortment of segments from multiple species during mixed infections. In this study, the prevalence of equine rotaviruses was 32.11% in organized equine farms of North India. The complete genome analysis of two ERVA isolates revealed an unusual genomic constellation, which was previously reported only in a bat RVA strain. A segment-wise phylogenetic analysis revealed that most segments of both isolates were highly similar either to bat or to bat-like human rotaviruses. The occurrence of unusual bat-like rotaviruses in equines emphasizes the need of extensive surveillance of complete genomes of both animal and human rotaviruses with a "one health" approach.