Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

Zero-shot-capable identification of phage-host relationships with whole-genome sequence representation by contrastive learning.

Accurately identifying phage-host relationships from their genome sequences is still challenging, especially for those phages and hosts with less homologous sequences. In this work, focusing on identifying the phage-host relationships at the species and genus level, we propose a contrastive learning based approach to learn whole-genome sequence embeddings that can take account of phage-host interactions (PHIs). Contrastive learning is used to make phages infecting the same hosts close to each other in the new representation space. Specifically, we rephrase whole-genome sequences with frequency chaos game representation (FCGR) and learn latent embeddings that 'encapsulate' phages and host relationships through contrastive learning. The contrastive learning method works well on the imbalanced dataset. Based on the learned embeddings, a proposed pipeline named CL4PHI can predict known hosts and unseen hosts in training. We compare our method with two recently proposed state-of-the-art learning-based methods on their benchmark datasets. The experiment results demonstrate that the proposed method using contrastive learning improves the prediction accuracy on known hosts and demonstrates a zero-shot prediction capability on unseen hosts. In terms of potential applications, the rapid pace of genome sequencing across different species has resulted in a vast amount of whole-genome sequencing data that require efficient computational methods for identifying phage-host interactions. The proposed approach is expected to address this need by efficiently processing whole-genome sequences of phages and prokaryotic hosts and capturing features related to phage-host relationships for genome sequence representation. This approach can be used to accelerate the discovery of phage-host interactions and aid in the development of phage-based therapies for infectious diseases.

An increased prevalence of carbapenem-resistant hypervirulent Klebsiella pneumoniae associated with the COVID-19 pandemic.

BACKGROUND: Klebsiella pneumoniae (Kp) is a common community-acquired and nosocomial pathogen. Carbapenem-resistant and hypervirulent (CR-hvKp) variants can emerge rapidly within healthcare facilities and impacted by other infectious agents such as COVID-19 virus. METHODS: To understand the impact of COVID-19 virus on the prevalence of CR-hvKp, we accessed Kp genomes with corresponding metadata from GenBank. Sequence types (STs), antimicrobial resistance genes, and virulence genes, and those scores and CR-hvKp were identified. We analyzed population diversity and phylogenetic characteristics of five most common STs, measured the prevalence of CR-hvKp, identified CR-hvKp subtypes, and determined associations between carbapenem resistance gene subtypes with STs and plasmid types. These variables were compared pre- and during the COVID-19 pandemic. FINDINGS: The proportion of CR-hvKp isolates increased within multiple STs in different continents during the COVID-19 pandemic and persistent CR-hvKp subtypes were found in common STs. blaKPC was dominant in CG258, blaKPC-2 was detected in 97 % of the ST11 CR-hvKp, blaNDM subtypes were prominent in ST147 (87.4 %) and ST307 (70.8 %); blaOXA-48 and its subtypes were prevalent in ST15 (80.5 %). The possession of carbapenemase genes was different among subclades from different origins in different periods of time within each ST. IncFIB/IncHI1B hybrid plasmids contained virulence genes and carbapenemase genes and were predominant in ST147 (67.37 %) and ST307 (56.25 %). INTERPRETATION: The prevalence of CR-hvKp increased during the COVID-19 pandemic, which was evident by an increase in local endemic clones. This process was facilitated by the convergence of plasmids containing carbapenemase genes and virulence genes. These findings have implications for the appropriate use of antimicrobials and infection prevention and control during outbreaks of respiratory viruses and pandemic management.

Spatial scale of tuberculosis transmission in Lima, Peru.

Spatially targeted interventions may be effective alternatives to individual or population-based prevention strategies against tuberculosis (TB). However, their efficacy may depend on the mechanisms that lead to geographically constrained hotspots. Local TB incidence may reflect high levels of local transmission; conversely, they may point to frequent travel of community members to high-risk areas. We used whole-genome sequencing to explore patterns of TB incidence and transmission in Lima, Peru. Between 2009 and 2012, we recruited incident pulmonary TB patients and their household contacts, whom we followed for the occurrence of TB disease. We used whole-genome sequences of 2,712 Mycobacterial tuberculosis isolates from 2,440 patients to estimate pariwise genomic distances and compared these to the spatial distance between patients' residences. Genomic distances increased rapidly as spatial distances increased and remained high beyond 2 km of separation. Next, we divided the study catchment area into 1 × 1 km grid-cell surface units and used household spatial coordinates to locate each TB patient to a specific cell. We estimated cell-specific transmission by calculating the proportion of patients in each cell with a pairwise genomic distance of 10 or fewer single-nucleotide polymorphisms. We found that cell-specific TB incidence and local transmission varied widely but that cell-specific TB incidence did not correlate closely with our estimates of local transmission (Cohen's k = 0.27). These findings indicate that an understanding of the spatial heterogeneity in the relative proportion of TB due to local transmission may help guide the implementation of spatially targeted interventions.

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats.

BACKGROUND: Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate. RESULTS: To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds. CONCLUSIONS: The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Unveiling biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2.

BACKGROUND: A growing number of studies have demonstrated that the polar regions have the potential to be a significant repository of microbial resources and a potential source of active ingredients. Genome mining strategy plays a key role in the discovery of bioactive secondary metabolites (SMs) from microorganisms. This work highlighted deciphering the biosynthetic potential of an Arctic marine-derived strain Aspergillus sydowii MNP-2 by a combination of whole genome analysis and antiSMASH as well as feature-based molecular networking (MN) in the Global Natural Products Social Molecular Networking (GNPS). RESULTS: In this study, a high-quality whole genome sequence of an Arctic marine strain MNP-2, with a size of 34.9 Mb was successfully obtained. Its total number of genes predicted by BRAKER software was 13,218, and that of non-coding RNAs (rRNA, sRNA, snRNA, and tRNA) predicted by using INFERNAL software was 204. AntiSMASH results indicated that strain MNP-2 harbors 56 biosynthetic gene clusters (BGCs), including 18 NRPS/NRPS-like gene clusters, 10 PKS/PKS-like gene clusters, 8 terpene synthse gene clusters, 5 indole synthase gene clusters, 10 hybrid gene clusters, and 5 fungal-RiPP gene clusters. Metabolic analyses of strain MNP-2 grown on various media using GNPS networking revealed its great potential for the biosynthesis of bioactive SMs containing a variety of heterocyclic and bridge-ring structures. For example, compound G-8 exhibited a potent anti-HIV effect with an IC50 value of 7.2 nM and an EC50 value of 0.9 nM. Compound G-6 had excellent in vitro cytotoxicities against the K562, MCF-7, Hela, DU145, U1975, SGC-7901, A549, MOLT-4, and HL60 cell lines, with IC50 values ranging from 0.10 to 3.3 µM, and showed significant anti-viral (H1N1 and H3N2) activities with IC50 values of 15.9 and 30.0 µM, respectively. CONCLUSIONS: These findings definitely improve our knowledge about the molecular biology of genus A. sydowii and would effectively unveil the biosynthetic potential of strain MNP-2 using genomics and metabolomics techniques.

Genomic and transcriptomic landscape to decipher the genetic basis of hyperpigmentation in Lanping black-boned sheep (Ovis aries).

BACKGROUND: Lanping black-boned sheep (LPB) represent a distinctive mammalian species characterized by hyperpigmentation, resulting in black bone and muscle features, in contrast to their conventional counterparts exhibiting red muscle and white bone. The genetic basis underlying LPB hyperpigmentation has remained enigmatic. METHODS: In this study, we conducted whole-genome sequencing of 100 LPB and 50 Lanping normal sheep (LPN), and integrated this data with 421 sequenced datasets from wild and domestic sheep, shedding light on the genetic backdrop and genomic variations associated with LPB. Furthermore, we performed comparative RNA-Seq analysis using liver sample to pinpoint genes implicated in the pigmentation process. We generated a comprehensive dataset comprising 97,944,357 SNPs from 571 sheep, facilitating an in-depth exploration of genetic factors. RESULTS: Population genetic structure analysis revealed that the LPB breed traces its origin back to LPN, having evolved into a distinct breed. The integration of positively selected genes with differentially expressed genes identified two candidates, ERBB4 and ROR1, potentially linked to LPB hyperpigmentation. Comparative analysis of ERBB4 and ROR1 mRNA relative expression levels in liver, spleen, and kidney tissues of LPB, in comparison to Diqing sheep, revealed significant upregulation, except for ERBB4 in the liver. Gene expression heatmaps further underscored marked allelic frequency disparities in different populations. CONCLUSION: Our findings establish the evolutionary lineage of the LPB breed from LPN and underscore the involvement of ERBB4 and ROR1 genes in melanin synthesis. These results enhance our comprehension of the molecular basis of hyperpigmentation and contribute to a more comprehensive depiction of sheep diversity.

Genomic dissection of the correlation between milk yield and various health traits using functional and evolutionary information about imputed sequence variants of 34,497 German Holstein cows.

BACKGROUND: Over the last decades, it was subject of many studies to investigate the genomic connection of milk production and health traits in dairy cattle. Thereby, incorporating functional information in genomic analyses has been shown to improve the understanding of biological and molecular mechanisms shaping complex traits and the accuracies of genomic prediction, especially in small populations and across-breed settings. Still, little is known about the contribution of different functional and evolutionary genome partitioning subsets to milk production and dairy health. Thus, we performed a uni- and a bivariate analysis of milk yield (MY) and eight health traits using a set of ~34,497 German Holstein cows with 50K chip genotypes and ~17 million imputed sequence variants divided into 27 subsets depending on their functional and evolutionary annotation. In the bivariate analysis, eight trait-combinations were observed that contrasted MY with each health trait. Two genomic relationship matrices (GRM) were included, one consisting of the 50K chip variants and one consisting of each set of subset variants, to obtain subset heritabilities and genetic correlations. In addition, 50K chip heritabilities and genetic correlations were estimated applying merely the 50K GRM. RESULTS: In general, 50K chip heritabilities were larger than the subset heritabilities. The largest heritabilities were found for MY, which was 0.4358 for the 50K and 0.2757 for the subset heritabilities. Whereas all 50K genetic correlations were negative, subset genetic correlations were both, positive and negative (ranging from -0.9324 between MY and mastitis to 0.6662 between MY and digital dermatitis). The subsets containing variants which were annotated as noncoding related, splice sites, untranslated regions, metabolic quantitative trait loci, and young variants ranked highest in terms of their contribution to the traits` genetic variance. We were able to show that linkage disequilibrium between subset variants and adjacent variants did not cause these subsets` high effect. CONCLUSION: Our results confirm the connection of milk production and health traits in dairy cattle via the animals` metabolic state. In addition, they highlight the potential of including functional information in genomic analyses, which helps to dissect the extent and direction of the observed traits` connection in more detail.

Variability of macrolide-resistant profile in Mycobacterium avium complex pulmonary disease.

This single-center retrospective study aimed to analyze the variability of macrolide resistance (MR) in 68 patients with Mycobacterium avium complex pulmonary disease. Among 25 patients treated without macrolides, 13 (52%) reverted to macrolide-susceptible (MS) profiles. Only one (2%) of 43 patients who continued macrolide treatment showed this change. We compared 30 MR isolates with recent specimens. Among them, seven shifted to MS (five attributed to clonally related strains; two resulting from reinfection or polyclonal infection).

Stimulating T cell responses against patient-derived breast cancer cells with neoantigen peptide-loaded peripheral blood mononuclear cells.

Breast cancer stands as a formidable global health challenge for women. While neoantigens exhibit efficacy in activating T cells specific to cancer and instigating anti-tumor immune responses, the accuracy of neoantigen prediction remains suboptimal. In this study, we identified neoantigens from the patient-derived breast cancer cells, PC-B-142CA and PC-B-148CA cells, utilizing whole-genome and RNA sequencing. The pVAC-Seq pipeline was employed, with minor modification incorporating criteria (1) binding affinity of mutant (MT) peptide with HLA (IC50 MT) ≤ 500 nm in 3 of 5 algorithms and (2) IC50 wild type (WT)/MT > 1. Sequencing results unveiled 2513 and 3490 somatic mutations, and 646 and 652 non-synonymous mutations in PC-B-142CA and PC-B-148CA, respectively. We selected the top 3 neoantigens to perform molecular dynamic simulation and synthesized 9-12 amino acid neoantigen peptides, which were then pulsed onto healthy donor peripheral blood mononuclear cells (PBMCs). Results demonstrated that T cells activated by ADGRL1E274K, PARP1E619K, and SEC14L2R43Q peptides identified from PC-B-142CA exhibited significantly increased production of interferon-gamma (IFN-γ), while PARP1E619K and SEC14L2R43Q peptides induced the expression of CD107a on T cells. The % tumor cell lysis was notably enhanced by T cells activated with MT peptides across all three healthy donors. Moreover, ALKBH6V83M and GAAI823T peptides from PC-B-148CA remarkably stimulated IFN-γ- and CD107a-positive T cells, displaying high cell-killing activity against target cancer cells. In summary, our findings underscore the successful identification of neoantigens with anti-tumor T cell functions and highlight the potential of personalized neoantigens as a promising avenue for breast cancer treatment.

Whole-genome sequencing and analysis of Chryseobacterium arthrosphaerae from Rana nigromaculata.

Chryseobacterium arthrosphaerae strain FS91703 was isolated from Rana nigromaculata in our previous study. To investigate the genomic characteristics, pathogenicity-related genes, antimicrobial resistance, and phylogenetic relationship of this strain, PacBio RS II and Illumina HiSeq 2000 platforms were used for the whole genome sequencing. The genome size of strain FS91703 was 5,435,691 bp and GC content was 37.78%. A total of 4,951 coding genes were predicted; 99 potential virulence factors homologs were identified. Analysis of antibiotic resistance genes revealed that strain FS91703 harbored 10 antibiotic resistance genes in 6 categories and 2 multidrug-resistant efflux pump genes, including adeG and farA. Strain FS91703 was sensitive to ß-lactam combination drugs, cephem, monobactam and carbapenems, intermediately resistant to phenicol, and resistant to penicillin, aminoglycosides, tetracycline, fluoroquinolones, and folate pathway inhibitors. Phylogenetic analysis revealed that strain FS91703 and C. arthrosphaerae CC-VM-7T were on the same branch of the phylogenetic tree based on 16 S rRNA; the ANI value between them was 96.99%; and the DDH values were 80.2, 72.2 and 81.6% by three default calculation formulae. These results suggested that strain FS91703 was a species of C. arthrosphaerae. Pan-genome analysis showed FS91703 had 566 unique genes compared with 13 other C. arthrosphaerae strains, and had a distant phylogenetic relationship with the other C. arthrosphaerae strains of the same branch in phylogenetic tree based on orthologous genes. The results of this study suggest that strain FS91703 is a multidrug-resistant and highly virulent bacterium, that differs from other C. arthrosphaerae strains at the genomic level. The knowledge about the genomic characteristics and antimicrobial resistance of strain FS91703 provides valuable insights into this rare species, as well as guidance for the treatment of the disease caused by FS91703 in Rana nigromaculata.

PE/PPE mutations in the transmission of Mycobacterium tuberculosis in China revealed by whole genome sequencing.

OBJECTIVE: This study aims to examine the impact of PE/PPE gene mutations on the transmission of Mycobacterium tuberculosis (M. tuberculosis) in China. METHODS: We collected the whole genome sequencing (WGS) data of 3202 M. tuberculosis isolates in China from 2007 to 2018 and investigated the clustering of strains from different lineages. To evaluate the potential role of PE/PPE gene mutations in the dissemination of the pathogen, we employed homoplastic analysis to detect homoplastic single nucleotide polymorphisms (SNPs) within these gene regions. Subsequently, logistic regression analysis was conducted to analyze the statistical association. RESULTS: Based on nationwide M. tuberculosis WGS data, it has been observed that the majority of the M. tuberculosis burden in China is caused by lineage 2 strains, followed by lineage 4. Lineage 2 exhibited a higher number of transmission clusters, totaling 446 clusters, of which 77 were cross-regional clusters. Conversely, there were only 52 transmission clusters in lineage 4, of which 9 were cross-regional clusters. In the analysis of lineage 2 isolates, regression results showed that 4 specific gene mutations, PE4 (position 190,394; c.46G > A), PE_PGRS10 (839,194; c.744 A > G), PE16 (1,607,005; c.620T > G) and PE_PGRS44 (2,921,883; c.333 C > A), were significantly associated with the transmission of M. tuberculosis. Mutations of PE_PGRS10 (839,334; c.884 A > G), PE_PGRS11 (847,613; c.1455G > C), PE_PGRS47 (3,054,724; c.811 A > G) and PPE66 (4,189,930; c.303G > C) exhibited significant associations with the cross-regional clusters. A total of 13 mutation positions showed a positive correlation with clustering size, indicating a positive association. For lineage 4 strains, no mutations were found to enhance transmission, but 2 mutation sites were identified as risk factors for cross-regional clusters. These included PE_PGRS4 (338,100; c.974 A > G) and PPE13 (976,897; c.1307 A > C). CONCLUSION: Our results indicate that some PE/PPE gene mutations can increase the risk of M. tuberculosis transmission, which might provide a basis for controlling the spread of tuberculosis.

Whole-genome sequencing and assessment of a novel protein- and gossypol-degrading Bacillus subtilis strain isolated from intestinal digesta of Tibetan Pigs.

BACKGROUND: With the rapid development of animal husbandry, the demand for protein feed resources is increasing. Cottonseed meal (CSM) and soybean meal (SBM) are rich sources of protein. However, their application is limited due to the existence of anti-nutrients, which can be harmful to the digestion and absorption. A strain of Bacillus subtilis (Mafic-Y7) was isolated from digesta of intestines of Tibetan pigs. The strain showed high protease activity, which helps in degrading proteinic anti-nutritional factors in grain meal and in vitro degradation of free gossypol. In order to better understand this isolated strain, whole genome of Mafic-Y7 strain was sequenced and analyzed. Different effects on various grain meals were identified. RESULT: The GC-depth Poisson distributions showed no bias suggesting high-quality genome assembly of Mafic-Y7. The whole genome sequencing showed that one chromosome with 4,248,845 base pairs(bp)and the genes total length with 3,736,524 bp was predicted in Mafic-Y7. Additionally, Mafic-Y7 possessed 4,254 protein-coding genes, and several protease genes were annotated by aligning them with databases. There are 55 protease genes, one phytase gene and one laccase gene were annotated in the gene sequence of Mafic-Y7. The average nucleotide identity between Mafic-Y7 and the GCA-000009045.1 homologous genome was 0.9938, suggesting a close genetic relationship between them at the species level. Compared with the closest four whole genomes, Mafic-Y7 was annotated the most abundant of protease genes (55 genes). The fermentation supernatant of Mafic-Y7 could increase the content of small peptides, water-soluble proteins, and acid-soluble proteins in vitro by 411%, 281% and 317% in SBM and 420%, 257% and 338% in CSM. After fermentation in grain meal by Mafic-Y7, the degradation rate of anti-nutritional factors in SBM, such as trypsin inhibitor, glycinin, and ß-conglycinin was greater than 70%, and lectin was greater than 30%. The degradation rates of anti-nutritional factors in CSM, such as gossypol and phytic acid, were 82% and 26%, respectively.

Genetic identification of familial hypercholesterolemia within whole genome sequences in 6820 newborns.

Familial hypercholesterolemia (FH) is defined as a monogenic disease, characterized by elevated low-density lipoprotein cholesterol (LDL-C) levels. FH remains underdiagnosed and undertreated in Chinese. We whole-genome sequenced 6820 newborns from Qingdao of China to investigate the FH-related gene (LDLR, APOB, PCSK9) mutation types, carrier ratio and genotype-phenotype correlation. In this study, the prevalence of FH in Qingdao of China was 0.47% (95% CI: 0.32%-0.66%). The plasma lipid levels of FH-related gene mutation carriers begin to increase as early as infant. T-CHO and LDL-C of FH infants was higher by 48.1% (p < 0.001) and 42.9% (p < 0.001) relative to non-FH infants. A total of 22 FH infants and their parent participate in further studies. The results indicated that FH infant parent noncarriers have the normal plasma lipid level, while T-CHO and LDL-C increased in FH infants and FH infant parent carriers, but no difference between the groups. This highlights the importance of genetic factors. In conclusion, the spectrum of FH-causing mutations in the newborns of Qingdao, China was described for the first time. These data can serve as a considerable dataset for next-generation sequencing analysis of the Chinese population with FH and potentially helping reform regional policies for early detection and prevention of FH.

Taxonogenomic analysis of the Xanthomonas translucens complex leads to the descriptions of Xanthomonas cerealis sp. nov. and Xanthomonas graminis sp. nov.

The pathovar-based taxonomy of the Xanthomonas translucens group is very confusing due to an overlap of plant host ranges and level of host specificity. Here, whole-genome sequence-based parameters (digital DNA-DNA hybridization and blast-based average nucleotide identity), phylogenomic, biochemical and phenotypical data were used to taxonomically analyse the 11 known pathovars of the X. translucens complex. This polyphasic approach taxonomically assigned the 11 pathovars of X. translucens complex into three distinct species, two of which are new: X. translucens, X. cerealis sp. nov. and X. graminis sp. nov. X. translucens consists of three pathovars: pv. translucens (=pv. hordei), pv. pistaciae strain A ICMP 16316PT and pv. undulosa (=pv. secalis). X. cerealis sp. nov. encompasses the pv. cerealis strain LMG 679PT and pv. pistaciae strain B ICMP 16317PT with genome similarity of 92.7% (dDDH) and 99.0% (ANIb) suggesting taxonomically similar genotypes. The other new species, X. graminis sp. nov., consists of the remaining five designated pathovars (pv. graminis, pv. arrhenatheri, pv. poae, pv. phleipratensis and pv. phlei) with highly variable dDDH and ANIb values ranging from 74.5 to 93.0% and from 96.7 to 99.2%, respectively, an indication of a very divergent taxonomic group. Only strains of pvs. phlei and phleipratensis showed the highest genomic similarities of 93.0% (dDDH) and 99.2% (ANIb), suggesting synonymic pathovars as both infect the same plant hosts. The dDDH and ANI data were corroborated by phylogenomics clustering. The fatty acid contents were similar but the type strain of X. graminis sp. nov. exhibited 20% less C15 : 0 iso and 40% more C17 : 0 iso fatty acids than the other species. Based on phenotypic, biochemical and whole-genome sequence data, we propose two new species, Xanthomonas cerealis sp. nov. and Xanthomonas graminis sp. nov. with type strains LMG 679T (=NCPPB 1944T) and LMG 726T (=NCPPB 2700T), respectively.

Genome-based taxonomic identification and safety assessment of an Enterococcus strain isolated from a homemade dairy product.

The aim of this study was to explore the taxonomic identification and evaluate the safety of a bacterium, Enterococcus lactis IDCC 2105, isolated from homemade cheese in Korea, using whole genome sequence (WGS) analysis. It sought to identify the species level of this Enterococcus spp., assess its antibiotic resistance, and evaluate its virulence potential. WGS analysis confirmed the bacterial strain IDCC 2105 as E. lactis and identified genes responsible for resistance to erythromycin and clindamycin, specifically msrC, and eatAv, which are chromosomally located, indicating a minimal risk for horizontal gene transfer. The absence of plasmids in E. lactis IDCC 2105 further diminishes the likelihood of resistance gene dissemination. Additionally, our investigation into seven virulence factors, including hemolysis, platelet aggregation, biofilm formation, hyaluronidase, gelatinase, ammonia production, and ß-glucuronidase activity, revealed no detectable virulence traits. Although bioinformatic analysis suggested the presence of collagen adhesion genes acm and scm, these were not corroborated by phenotypic virulence assays. Based on these findings, E. lactis IDCC 2105 presents as a safe strain for potential applications, contributing valuable information on its taxonomy, antibiotic resistance profile, and lack of virulence factors, supporting its use in food products.

Rare serotype c Haemophilus influenzae invasive isolate: characterization of the first case in Portugal.

We report for the first time in Portugal a serotype c Haemophilus influenzae isolated from an adult, with HIV-1 infection. Whole-genome sequencing characterized the isolate as clonal complex ST-7, albeit with a novel MLST (ST2754) due to a unique atpG profile. Integration of this genome with other available H. influenzae serotype c genomes from PubMLST revealed its overall genetic distinctiveness, with the closest related isolate being identified in France in 2020. This surveillance study, involving collaboration among hospitals and reference laboratory, successfully contributed to the identification and characterization of this rare serotype.

An overview of next generation sequencing strategies and genomics tools used for tuberculosis research.

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

Whole-genome sequencing of Pseudoalteromonas piscicida 2515 revealed its antibacterial potency against Vibrio anguillarum: a preliminary invitro study.

Pseudoalteromonas piscicida 2515, isolated from Litopenaeus vannamei culture water, is a potential marine probiotic with broad anti-Vibrio properties. However, genomic information on P. piscicida 2515 is scarce. In this study, the general genomic characteristics and probiotic properties of the P. piscicida 2515 strain were analysed. In addition, we determined the antibacterial mechanism of this bacterial strain by scanning electron microscopy (SEM). The results indicated that the whole-genome sequence of P. piscicida 2515 contained one chromosome and one plasmid, including a total length of 5,541,406 bp with a G + C content of 43.24%, and 4679 protein-coding genes were predicted. Various adhesion-related genes, amino acid and vitamin metabolism and biosynthesis genes, and stress-responsive genes were found with genome mining tools. The presence of genes encoding chitin, bromocyclic peptides, lantibiotics, and sactipeptides showed the strong antibacterial activity of the P. piscicida 2515 strain. Moreover, in coculture with Vibrio anguillarum, P. piscicida 2515 displayed vesicle/pilus-like structures located on its surface that possibly participated in its bactericidal activity, representing an antibacterial mechanism. Additionally, 16 haemolytic genes and 3 antibiotic resistance genes, including tetracycline, fluoroquinolone, and carbapenem were annotated, but virulence genes encoding enterotoxin FM (entFM), cereulide (ces), and cytotoxin K were not detected. Further tests should be conducted to confirm the safety characteristics of P. piscicida 2515, including long-term toxicology tests, ecotoxicological assessment, and antibiotic resistance transfer risk assessment. Our results here revealed a new understanding of the probiotic properties and antibacterial mechanism of P. piscicida 2515, in addition to theoretical information for its application in aquaculture.

MFSD2A frameshift variant in Kerry Hill sheep with microcephaly.

Microcephaly is a rare neurodevelopmental disorder characterized by reduced skull circumference and brain volume that occurs sporadically in farm animals. We investigated an early-onset neurodegenerative disorder observed in seven lambs of purebred Kerry Hill sheep. Clinical signs included inability to stand or severe ataxia, convulsions, and early death. Diagnostic imaging and brain necropsy confirmed microcephaly. The pedigree of the lambs suggested monogenic autosomal recessive inheritance. We sequenced the genome of one affected lamb, and comparison with 115 control genomes revealed a single private protein-changing variant. This frameshift variant, MFSD2A: c.285dupA, p.(Asp96fs*9), represents a 1-bp duplication predicted to truncate 80% of the open reading frame. MFSD2A is a transmembrane protein that is essential for maintaining blood-brain barrier homeostasis and plays a key role in regulating brain lipogenesis. Human MFSD2A pathogenic variants are associated with a neurodevelopmental disorder with progressive microcephaly, spasticity, and brain imaging abnormalities (NEDMISBA, OMIM 616486). Here we present evidence for the occurrence of a recessively inherited form of microcephaly in sheep due to a loss-of-function variant in MFSD2A (OMIA 002371-9940). To the best of our knowledge, this is the first report of a spontaneous MFSD2A variant in domestic animals.