Novel Gram-positive, catalase-negative, α-haemolytic cocci were isolated from breast milk samples of healthy mothers living in Hanoi, Vietnam. The 16S rRNA gene sequences of these strains varied by 0-2 nucleotide polymorphisms. The 16S rRNA gene sequence of one strain, designated as BME SL 6.1T, showed the highest similarity to those of Streptococcus salivarius NCTC 8618T (99.4 %), Streptococcus vestibularis ATCC 49124T (99.4 %), and Streptococcus thermophilus ATCC 19258T (99.3 %) in the salivarius group. Whole genome sequencing was performed on three selected strains. Phylogeny based on 631 core genes clustered the three strains into the salivarius group, and the strains were clearly distinct from the other species in this group. The average nucleotide identity (ANI) value of strain BME SL 6.1T exhibited the highest identity with S. salivarius NCTC 8618T (88.4 %), followed by S. vestibularis ATCC 49124T (88.3 %) and S. thermophilus ATCC 19258T (87.4 %). The ANI and digital DNA-DNA hybridization values between strain BME SL 6.1T and other species were below the cut-off value (95 and 70 %, respectively), indicating that it represents a novel species of the genus Streptococcus. The strains were able to produce α-galactosidase and acid from raffinose and melibiose. Therefore, we propose to assign the strains to a new species of the genus Streptococcus as Streptococcus raffinosi sp. nov. The type strain is BME SL 6.1T (=VTCC 12812T=NBRC 116368T).

UNASSIGNED: Direct whole genome sequencing (WGS) of Mycobacterium tuberculosis (Mtb) can be used as a tool to study drug resistance, mixed infections, and within-host diversity. However, WGS is challenging to obtain from clinical samples due to low number of bacilli against a high background.
UNASSIGNED: We prospectively collected 34 samples (sputum, n = 17; bronchoalveolar lavage, n = 13; and pus, n = 4) from patients with active tuberculosis (TB). Prior to DNA extraction, we used a ligand-mediated magnetic bead method to enrich Mtb from clinical samples and performed WGS on Illumina platform.
UNASSIGNED: Mtb was definitively identified based on WGS from 88.2% (30/34) of the samples, of which 35.3% (12/34) were smear negative. The overall median genome coverage was 15.2% (interquartile range [IQR], 7.7%-28.2%). There was a positive correlation between load of bacilli on smears and genome coverage (P < .001). We detected 58 genes listed in the World Health Organization mutation catalogue in each positive sample (median coverage, 85% [IQR, 61%-94%]), enabling the identification of mutations missed by routine diagnostics. Mutations causing resistance to rifampicin, isoniazid, streptomycin, and ethambutol were detected in 5 of 34 (14.7%) samples, including the rpoB S441A mutation that confers resistance to rifampicin, which is not covered by Xpert MTB/RIF.
UNASSIGNED: We demonstrate the feasibility of magnetic bead-based enrichment for culture-free WGS of Mtb from clinical specimens, including smear-negative samples. This approach can also be integrated with low-cost sequencing workflows such as targeted sequencing for rapid detection of Mtb and drug resistance.

BACKGROUND: Esophageal carcinoma (EC) and gastric cardiac adenocarcinoma (GCA) have high incidence rates in the Chaoshan region of South China. Multifocal esophageal and cardiac cancer (MECC) is commonly observed in this region in clinical practice. However, the genomic characteristics of MECC remains unclear.
METHODS: In this study, a total of 2123 clinical samples of EC and GCA were analyzed to determine the frequency of multifocal tumors, as well as their occurrence sites and pathological types. Cox proportional hazards regression was used to model the relationship between age, sex, and tumor state concerning survival in our analysis of the cohort of 541 patients with available follow-up data. We performed whole-genome sequencing on 20 tumor foci and 10 normal samples from 10 MECC patients to infer clonal structure on 6 MECC patients to explore genome characteristics.
RESULTS: The MECC rate of EC and GCA was 5.65% (121 of 2123). Age and sex were potential factors that may influence the risk of MECC (p < 0.001). Furthermore, MECC patients showed worse survival compared with single tumor patients. We found that 12 foci from 6 patients were multicentric origin model (MC), which exhibited significant heterogeneity of variations in paired foci and had an increased number of germline mutations in immune genes compared to metastatic model. In MC cases, different lesions in the same patient were driven by distinct mutation and copy number variation (CNV) events. Although TP53 and other driver mutation genes have a high frequency in the samples, their mutation sites show significant heterogeneity in paired tumor specimens. On the other hand, CNV genes exhibited higher concordance in paired samples, especially in the amplification of oncogenes and the deletion of tumor suppressor genes.
CONCLUSIONS: The extent of inter-tumor heterogeneity suggests both monoclonal and polyclonal origins of MECC, which could provide insight into the genome diversity of MECC and guide clinical implementation.

Clade 2.3.4.4b highly pathogenic avian influenza (HPAI) H5N1 virus was detected in the South American sea lions found dead in Santa Catarina, Brazil, in October 2023. Whole genome sequencing and comparative phylogenetic analysis were conducted to investigate the origin, genetic diversity, and zoonotic potentials of the H5N1 viruses. The H5N1 viruses belonged to the genotype B3.2 of clade 2.3.4.4b H5N1 virus, which was identified in North America and disseminated to South America. They have acquired new amino acid substitutions related to mammalian host affinity. Our study provides insights into the genetic landscape of HPAI H5N1 viruses in Brazil, highlighting the continuous evolutionary processes contributing to their possible adaptation to mammalian hosts.

BACKGROUND: Breeding polled goats is a welfare-friendly approach for horn removal in comparison to invasive methods. To gain a comprehensive understanding of the genetic basis underlying polledness in goats, we conducted whole-genome sequencing of 106 Xinong Saanen dairy goats, including 33 horned individuals, 70 polled individuals, and 3 polled intersexuality syndrome (PIS) individuals.
METHODS: The present study employed a genome-wide association study (GWAS) and linkage disequilibrium (LD) analysis to precisely map the genetic locus underlying the polled phenotype in goats.
RESULTS: The analysis conducted in our study revealed a total of 320 genome-wide significant single nucleotide polymorphisms (SNPs) associated with the horned/polled phenotype in goats. These SNPs exhibited two distinct peaks on chromosome 1, spanning from 128,817,052 to 133,005,441 bp and from 150,336,143 to 150,808,639 bp. The present study identified three genome-wide significant SNPs, namely Chr1:129789816, Chr1:129791507, and Chr1:129791577, as potential markers of PIS-affected goats. The results of our LD analysis suggested a potential association between MRPS22 and infertile intersex individuals, as well as a potential association between ERG and the polled trait in goats.
CONCLUSIONS: We have successfully identified three marker SNPs closely linked to PIS, as well as several candidate genes associated with the polled trait in goats. These results may contribute to the development of SNP chips for early prediction of PIS in goats, thereby facilitating breeding programs aimed at producing fertile herds with polled traits.

OBJECTIVE: The fresh-market tomato (Solanum lycopersicum) is bred for direct human consumption. It is selected for specific traits to meet market demands and production systems, and unique genetic variations underlying fresh-market tomato yields have been recently identified. However, DNA sequence variant-trait associations are not yet fully examined even for major traits. To provide a rich genome sequence resource for various genetics and breeding goals for fresh-market tomato traits, we report whole genome sequence data of a pool of contemporary U.S. fresh-market tomatoes.
METHODS: Eighty-one tomatoes were nominated by academic tomato breeding programs in the U.S. Of the 81 tomatoes, 68 were contemporary fresh-market tomatoes, whereas the remaining 13 were relevant fresh-market tomato breeding and germplasm accessions. Whole genome sequencing (WGS) of the 81 tomatoes was conducted using the Illumina next-generation sequencing technology. The polymerase chain reaction (PCR)-free, paired-end sequencing libraries were sequenced on an average depth per sequenced base of 24 × for each tomato. This data note enhances visibility and potential for use of the more diverse, freely accessible whole genome sequence data of contemporary fresh-market tomatoes.

In order to resolve the key genes for weed control by Trichoderma polysporum at the genomic level, we extracted the genomic DNA and sequenced the whole genome of T. polysporum strain HZ-31 on the Illumina Hiseq platform. The raw data was cleaned up using Trimmomatic and checked for quality using FastQC. The sequencing data was assembled using SPAdes, and GeneMark was used to perform gene prediction on the assembly results. The results showed that the genome size of T. polysporum HZ-31 was 39,325,746 bp, with 48% GC content, and the number of genes encoded was 11,998. A total of 148 tRNAs and 45 rRNAs were predicted. A total of 782 genes were annotated in the Carbohydrase Database, 757 genes were annotated to the Pathogen-Host Interaction Database, and 67 gene clusters were identified. In addition, 1023 genes were predicted to be signal peptide proteins. The annotation and functional analysis of the whole genome sequence of T. polymorpha HZ-31 provide a basis for the in-depth study of the molecular mechanism of its herbicidal action and more effective utilization for weed control.

Here, we present the whole genome sequence of Bt S2160-1, a potential alternative to the mosquitocidal model strain, Bti. One chromosome genome and four mega-plasmids were contained in Bt S2160-1, and 13 predicted genes encoding predicted insecticidal crystal proteins were identified clustered on one plasmid pS2160-1p2 containing two pathogenic islands (PAIs) designed as PAI-1 (Cry54Ba, Cry30Ea4, Cry69Aa-like, Cry50Ba2-like, Cry4Ca1-like, Cry30Ga2, Cry71Aa-like, Cry72Aa-like, Cry70Aa-like, Cyt1Da2-like and Vpb4C1-like) and PAI-2 (Cyt1Aa-like, and Tpp80Aa1-like). The clusters appear to represent mosquitocidal toxin islands similar to pathogenicity islands. Transcription/translation of 10 of the 13 predicted genes was confirmed by whole-proteome analysis using LTQ-Orbitrap LC-MS/MS. In summary, the present study identified the existence of a mosquitocidal toxin island in Bacillus thuringiensis, and provides important genomic information for understanding the insecticidal mechanism of B. thuringiensis.

BACKGROUND: Missing data is a common challenge in mass spectrometry-based metabolomics, which can lead to biased and incomplete analyses. The integration of whole-genome sequencing (WGS) data with metabolomics data has emerged as a promising approach to enhance the accuracy of data imputation in metabolomics studies.
METHODS: In this study, we propose a novel method that leverages the information from WGS data and reference metabolites to impute unknown metabolites. Our approach utilizes a multi-scale variational autoencoder to jointly model the burden score, polygenetic risk score (PGS), and linkage disequilibrium (LD) pruned single nucleotide polymorphisms (SNPs) for feature extraction and missing metabolomics data imputation. By learning the latent representations of both omics data, our method can effectively impute missing metabolomics values based on genomic information.
RESULTS: We evaluate the performance of our method on empirical metabolomics datasets with missing values and demonstrate its superiority compared to conventional imputation techniques. Using 35 template metabolites derived burden scores, PGS and LD-pruned SNPs, the proposed methods achieved R2-scores > 0.01 for 71.55 % of metabolites.
CONCLUSIONS: The integration of WGS data in metabolomics imputation not only improves data completeness but also enhances downstream analyses, paving the way for more comprehensive and accurate investigations of metabolic pathways and disease associations. Our findings offer valuable insights into the potential benefits of utilizing WGS data for metabolomics data imputation and underscore the importance of leveraging multi-modal data integration in precision medicine research.

Four rod-shaped, non-motile, non-spore-forming, facultative anaerobic, Gram-stain-positive lactic acid bacteria, designated as EB0058T, SCR0080, LD0937T and SCR0063T, were isolated from different corn and grass silage samples. The isolated strains were characterized using a polyphasic approach and EB0058T and SCR0080 were identified as Lacticaseibacillus zeae by 16S rRNA gene sequence analysis. Based on whole-genome sequence-based characterization, EB0058T and SCR0080 were separated into a distinct clade from Lacticaseibacillus zeae DSM 20178T, together with CECT9104 and UD2202, whose genomic sequences are available from NCBI GenBank. The average nucleotide identity (ANI) values within the new subgroup are 99.9 % and the digital DNA-DNA hybridization (dDDH) values are 99.3-99.9 %, respectively. In contrast, comparison of the new subgroup with publicly available genomic sequences of L. zeae strains, including the type strain DSM 20178T, revealed dDDH values of 70.2-72.5 % and ANI values of 96.2-96.6 %. Based on their chemotaxonomic, phenotypic and phylogenetic characteristics, EB0058T and SCR0080 represent a new subspecies of L. zeae. The name Lacticaseibacillus zeae subsp. silagei subsp. nov. is proposed with the type strain EB0058T (=DSM 116376T=NCIMB 15474T). According to the results of 16S rRNA gene sequencing, LD0937T and SCR0063T are members of the Lacticaseibacillus group. The dDDH value between the isolates LD0937T and SCR0063T was 67.6 %, which is below the species threshold of 70 %, clearly showing that these two isolates belong to different species. For both strains, whole genome-sequencing revealed that the closest relatives within the Lacticaseibacillus group were Lacticaseibacillus huelsenbergensis DSM 115425 (dDDH 66.5 and 65.9 %) and Lacticaseibacillus casei DSM 20011T (dDDH 64.1 and 64.9 %). Based on the genomic, chemotaxonomic and morphological data obtained in this study, two novel species, Lacticaseibacillus parahuelsenbergensis sp. nov. and Lacticaseibacillus styriensis sp. nov. are proposed and the type strains are LD0937T (=DSM 116105T=NCIMB 15471T) and SCR0063T (=DSM 116297T=NCIMB 15473T), respectively.