Genomic Features Research Articles

Genomic insights of Lactiplantibacillus plantarum CNPC024: a potential probiotic strain producing immune-boosting tryptophan-derived metabolites.

Probiogenomics can provide important insights bout probiotic candidate bacteria. This study aimed to perform an in-depth genomic characterization of the probiotic candidate Lactiplantibacillus plantarum CNPC024 to investigate its probiosis mechanisms, identify metabolic pathways that might benefit the host, and improve the safety assessment for this strain to be effectively used as a probiotic. After whole-genome sequencing in Illumina MiSeq platform, the de novo genome assembly resulted in a 3.2Mb draft genome. According to the Average Nucleotide Identity (ANI) analysis with 46 randomly validated probiotic LAB belonging to the Lactobacillaceae family, the strain showed a 99% nucleotide identity with other L. plantarum probiotic species. We identified a set of determinants conferring tolerance to bile salts and low pH conditions, as well as temperature, oxidative and osmotic stressors via the glutathione-glutaredoxin system (Grxs). As a β‑galactosidase‑producing strain, it has the potential to be used in fermented dairy products for lactose-intolerant individuals. There were no significant hits for transferable antibiotic-resistance genes. We also identified gene clusters associated with production of bacteriocins (plantaricins E, F and K). Lastly, we detected metabolic pathways associated with the production of tryptophan-derived metabolites that could potentially modulate the host's immune responses.

Genomic Insights into Vietnamese Extended-Spectrum β-Lactamase-9-Producing Extensively Drug-Resistant Pseudomonas aeruginosa Isolates Belonging to High-Risk Clone ST357 Obtained from Bulgarian Intensive Care Unit Patients

Extensively drug-resistant P. aeruginosa (XDR-PA) has been highlighted as a serious public health threat. The present study aimed to explore the genomic characteristics of two Vietnamese extended-spectrum β-lactamase-9 (VEB-9)-producing XDR-PA isolates from Bulgaria in comparison to all blaVEB-9-positive strains with available genomes. The isolates designated Pae51 and Pae52 were obtained from tracheobronchial aspirates of intensive care unit (ICU) patients. Antimicrobial susceptibility testing, whole-genome sequencing, RT-qPCR, and phylogenomic analysis were performed. Pae51 and Pae52 were resistant to most antipseudomonal β-lactams including carbapenems, aminoglycosides, and fluoroquinolones but remained susceptible to colistin and cefiderocol. Numerous resistance determinants were detected: blaVEB-9, blaPDC-3, blaOXA-10, blaOXA-50, aac(6′)-II, ant(2″)-Ia, ant(3″)-IIa, aph(3′)-IIb, cprP, catB7, dfrB2, sul1, fosA, and tet(A). Both isolates carried complex integrons with blaVEB-9 and tet(A) embedded next to the conservative 3′ end sequences. A variety of virulence factors were also identified, including the type III secretion system exotoxin U. Pae51 and Pae52 differed by only four SNPs and belonged to the high-risk clone ST357. To our knowledge, this is the first report of blaVEB-9-positive XDR-PA isolates in Bulgaria presenting a detailed genomic analysis. The development of novel antimicrobial strategies for such pathogens should be an essential part of infection control stewardship practices in ICU wards.

Relationship between virulence and carbapenem resistance phenotype of Klebsiella pneumoniae from blood infection: identification of a carbapenem-resistant and hypervirulent strain.

To investigate the relationship between the virulence and the carbapenem resistance phenotype of Klebsiella pneumoniae from blood infection, and to identify carbapenem-resistant and hypervirulent Klebsiella pneumoniae (CR-HVKP)strains. A total of 192 Klebsiella pneumoniae strains were isolated from blood culture of patients with bloodstream infections from 2016 to 2019, of which 96 isolates were carbapenem-resistant Klebsiella pneumoniae (CRKP) and 96 were carbapenem-sensitive Klebsiella pneumoniae (CSKP). The drug susceptibility was detected by VITEK-2 automatic microbial analyzer; carbapenemase genes, virulence genes and capsule typing were detected by polymerase chain reaction; the high viscosity phenotype of strains was detected by string test, and the genome characteristics of CR-HVKP were detected by whole genome sequencing. Serum killing and biofilm formation test were used to further verify the virulence of CR-HVKP. There were significant differences in drug resistance to common antibiotics, except for minocycline between CSKP and CRKP isolates (all P<0.05). 92 out of 96 CRKP isolates carried carbapenemase genes, mainly blaKPC-2. The string tests were positive in 4 isolates of CRKP and 36 isolates of CSKP (P<0.05). The detection rates of virulence genes Kfu, aerobictin, iutA, ybtS, rmpA, magA, allS, and capsule antigen K1 and K2 in CSKP group were significantly higher than those in CRKP group (all P<0.05). One HVKP strain was detected in the CRKP group (CR-HVKP) and 36 HVKP was detected in the CSKP group (P<0.05). The CR-HVKP strain belonged to the MLST412, serotype K57, expressed iutA, entB, mrkD, fimH, and rmpA virulence genes, and showed strong biofilm formation and significantly increased serum resistance. Whole genome sequencing results showed that this CR-HVKP isolate carried blaSHV-145, blaTEM-1, blaCTX-M-3, fosA6, oqxA5, oqxB26, and aac(3)-IId resistance genes, accompanied by abnormalities in outer membrane protein K (OmpK) 35 and OmpK36. The drug resistance of CRKP is significantly higher than that of CSKP, while CRKP carrying fewer virulence genes in both number and types compared to CSKP. A new MLST type of carbapenem-resistant and hypervirulent Klebsiella pneumoniae strain has been detected, which requires clinical awareness and epidemiological monitoring.

A New Medical Evaluation for Gastric Cancer Patients to Increase the Success Rate of Immunotherapy: A 2024 Update

Researchers have performed numerous studies on immunotherapy because of the high death rate associated with gastric cancer (GC). GC immunotherapy research has made tremendous progress, and we wanted to provide an update on this topic. On the basis of this update, we suggest performing a new medical evaluation before initiating immunotherapy in patients with GC to increase the success rate of immunotherapies. We propose that before patients start GC immunotherapy, they should be evaluated and given a score of one to two points for the following factors: immunopathological features, molecular and genomic features, potential consequences for bacterial pathogens, potential immunotherapeutic resistance and hyperprogressive illness, and the potential to use biomarkers to gauge their prognosis and immunotherapy responses to optimize immunotherapy following surgery. The proposed scoring system could also help in the diagnosis of GC. With all the advances in genetics, immunology, and microbiology, the diagnosis of GC could be improved, not changed. Currently, patients diagnosed with GC undergo surgical resection as the only permanent solution. Patients who meet the maximum score from the presented proposal could be eligible immediately after diagnosis for immunotherapy. Therefore, immunotherapy could be a first-line option for clinicians.

Lynch syndrome-associated and sporadic microsatellite unstable colorectal cancers: different patterns of clonal evolution yield highly similar tumours.

Microsatellite unstable colorectal cancer (MSI-CRC) can arise through germline mutations in mismatch repair (MMR) genes in individuals with Lynch syndrome (LS), or sporadically through promoter methylation of the MMR gene MLH1. Despite the different origins of hereditary and sporadic MSI tumours, their genomic features have not been extensively compared. A prominent feature of MMR-deficient genomes is the occurrence of many indels in short repeat sequences, an understudied mutation type due to the technical challenges of variant calling in these regions. In this study, we performed whole genome sequencing and RNA-sequencing on 29 sporadic and 14 hereditary MSI-CRCs. We compared the tumour groups by analysing genome-wide mutation densities, microsatellite repeat indels, recurrent protein-coding variants, signatures of single base, doublet base, and indel mutations, and changes in gene expression. We show that the mutational landscapes of hereditary and sporadic MSI-CRCs, including mutational signatures and mutation densities genome-wide and in microsatellites, are highly similar. Only a low number of differentially expressed genes were found, enriched to interferon-γ regulated immune response pathways. Analysis of the variance in allelic fractions of somatic variants in each tumour group revealed higher clonal heterogeneity in sporadic MSI-CRCs. Our results suggest that the differing molecular origins of MMR deficiency in hereditary and sporadic MSI-CRCs do not result in substantial differences in the mutational landscapes of these tumours. The divergent patterns of clonal evolution between the tumour groups may have clinical implications, as high clonal heterogeneity has been associated with decreased tumour immunosurveillance and reduced responsiveness to immunotherapy.

Sequence types of enteroaggregative Escherichia coli strains recovered from human, animal, and environmental sources: India

In the current study, we report whole genome sequencing (WGS) data on EAEC strains from India to identify lineages and different sequence types (STs) in our geographical regions across North India. We performed WGS comparative genomics characterization to examine the diversity of 122 EAEC strains collected from a large geographic area from clinical (Human sources) and non-clinical sources (animal and environmental sources). M-PCR for 21 virulence genes was performed. A triplex PCR detected phylogenetic groups A, B1, B2, and D was done. All strains were genome-sequenced, and bioinformatics analysis was performed. EAEC isolates belonged to 29 sequence types, further clustered into 11 clonal complexes, among which CC38 was the largest, containing 38 isolates mainly belonging to two ST types (ST38 and ST315). CC10 was the most diverse group, comprising 8 STs (ST43, ST2706, ST1286, ST 10, ST167, ST34, ST227, and ST4305). The most frequently detected virulence gene among the 96 clinical EAEC isolates was astA (87.5%), followed by ORF3 (62.5%), and aap (54.1 %). ST131, known for multidrug resistance and causing various diseases, was detected in acute diarrhea cases and animal sources, underscoring its clinical significance. These findings indicate the high diversity of EAEC and different sources of unique ST types of EAEC. The study reinforces the One Health approach, highlighting the interconnection between human health, animal reservoirs, and environmental sources in EAEC transmission dynamics. The identification of EAEC strains in animals suggests the zoonotic potential of these pathogens, necessitating comprehensive surveillance and control measures. The study emphasizes the need for continued surveillance and a One Health approach to address the complex transmission dynamics of EAEC.

Genomic characterization of third-generation cephalosporin-resistant Escherichia coli strains isolated from diseased dogs and cats: Report from Japanese Veterinary Antimicrobial Resistance Monitoring

This study investigates the genomic characteristics of canine and feline cefotaxime (CTX, a third-generation cephalosporin)-resistant Escherichia coli using the JVARM, Japanese Veterinary Antimicrobial Resistance Monitoring System, a nationwide monitoring.In this study, whole-genome sequencing (WGS) was performed on 51 canine and 45 feline CTX-resistant E. coli isolates, with certain isolates subjected to pulsed-field gel electrophoresis with S1 nuclease for plasmid–chromosome separation.The most common blaCTX-M genes were blaCTX-M-27 (dogs: 11/51 [21.6 %]; cat: 10/45 [22.2 %]), followed by blaCTX-M-14 (dogs: 10/51 [19.6 %]; cats: 10/45 [22.2 %]), and blaCTX-M-15 (dogs: 9/51 [17.6 %]; cats: 5/45 [11.1 %]). Besides β-lactamase genes, all isolates harbored mdf(A), a multidrug efflux pump, with resistance genes for aminoglycosides, sulfonamides, trimethoprims, macrolides and tetracyclines. None of the isolates had carbapenemase genes, such as blaOXA-48, blaNDM, and blaIMP, whereas most of the isolates showed double mutations in gyrA and parC, which affected quinolone resistance. For the isolates separately analyzed for plasmid and chromosomal DNA via WGS, the majority of CTX-M genes were present on the plasmids. Some plasmids also harbored the same combination of resistance genes and plasmid replicon type, although they differed from isolates derived from different areas of Japan. The predominant plasmids were blaCTX-M-27,aadA5, aph(6)-Id, aph(3”)-Ib, sul1, sul2, tet(A), dfrA17, and mph(A) on IncF. The predominant combination of ST131, O25:H4, and B2 isolates comprised the largest cluster in the minimum spanning tree and the ST131 E. coli harboring blaCTX-M-27 from human in Japan was closely related to these isolates.The results indicated that CTX-resistant canine and feline E. coli harbored multiple plasmids carrying the same combination of resistance genes and emphasizes the need to prevent the spread. Data AvailabilityAll raw short-read sequence data have been deposited in the DNA Data Bank of Japan. (DRR Run No, DRR335726–335821).

Exploring the multifaceted role of key lncRNA in glioma: From genetic expression to clinical implications and immunotherapy potential

BackgroundLong non-coding RNAs (lncRNAs) are implicated in a variety of regulatory functions within tumors, yet their specific roles in glioma remain underexplored. MethodsWe extracted glioma patient data from The Cancer Genome Atlas and UCSC Xena database for analysis using R, focusing on genomic characterization, biological enrichment, immune evaluation, and the development of a predictive model employing machine learning techniques. Additionally, we conducted cell culture and proliferation assays. ResultsOur analysis revealed that the lncRNA SLC16A1-AS1 plays a pivotal role in glioma pathogenesis and prognosis. We observed that abnormal expression of SLC16A1-AS1 varied with tumor grade, IDH mutation status, and histological type, correlating with worse survival outcomes. Genomically, SLC16A1-AS1 was associated with Tumor Mutational Burden and other prognostic biomarkers. The expression of this lncRNA was also linked to the activation of critical biological pathways and appeared to modulate the immune microenvironment, enhancing the presence of immune cells and checkpoints, which may be predictive of immunotherapy outcomes. Our predictive model, constructed from genes associated with SLC16A1-AS1, accurately forecasted glioma prognosis, strongly correlating with survival and treatment responses. In vitro experiments further demonstrated that SLC16A1-AS1 significantly influences glioma cell proliferation, invasion, and migration, underscoring its role in tumor aggression and its potential as a therapeutic target. ConclusionsThis study underscores the significant influence of SLC16A1-AS1 on glioma progression and prognosis, with its expression correlating with tumor traits and immune responses. The findings highlight the potential of targeting SLC16A1-AS1 in therapeutic strategies aimed at mitigating glioma aggressiveness.

Genomic characterisation and diversity assessment of eight endangered Belgian sheep breeds

Assessing the genetic diversity of local breeds is essential for conserving these unique breeds, which may possess unique traits. This study provides the genomic characterisation of eight indigenous sheep breeds in Belgium based on pedigree and single nucleotide polymorphism (SNP) analysis. A total of 687 sheep were genotyped and were subjected to a rigorous quality control, resulting in a set of 45 978 autosomal SNPs. Pedigree analysis showed breed-average inbreeding estimates between 3.3% and 11.3%. The genomic analysis included an assessment of runs of homozygosity (ROH) to examine the genomic inbreeding coefficient, with breed-average inbreeding coefficients estimated between 4.1% and 8.5%. Runs of homozygosity islands were identified in six of the eight breeds studied, with some exhibiting an incidence of up to 58%. Interestingly, several ROH islands overlapped with other breeds included in this study, as well as with international sheep breeds. Pedigree-based effective population sizes were estimated below 100 for all breeds, whereas genomic-based effective population sizes were below 24, indicating that these eight local sheep breeds are endangered. Principal component analysis, admixture analyses, and Fst computations were used to study the population structure and genetic differences. A neighbour−joining tree using 95 international sheep breeds positioned the eight local breeds in the group of milksheep, Texel sheep and the Scandinavian breeds. Additionally, the investigation of paternal oY1 genotypes revealed diverse lineage origins within the Belgian sheep population. This study refines and deepens our knowledge about the local sheep breeds in Belgium, thereby improving their management and conservation. Moreover, as these breeds are linked to other international breeds, these insights are significant for the global scientific community.

First identification and genomic features of multidrug-resistant Citrobacter freundii ST669 strain isolated from a domesticated duck in Bangladesh

ObjectivesCitrobacter freundii is a prevalent source of nosocomial infections and a well-known cause of diarrheal diseases. In recent years, it has also become increasingly resistant to various antimicrobials. In this study, we screened and characterized a multidrug-resistant (MDR) C. freundii isolate obtained from a domesticated diseased duck to better understand the genetic features, molecular epidemiology, and underlying factors linked to the antimicrobial resistance genes (ARGs) and virulence factor genes (VFGs) of the isolate. MethodsThe C. freundii BAU_TM8 strain was isolated using culturing, staining, biochemical, polymerase chain reaction, and Matrix-assisted laser desorption/ionization-time of flight methods. The MDR properties of the strain were determined by a disk diffusion test. The genomic sequence of C. freundii BAU_TM8 was performed using the Illumina NextSeq2000 platform. The ARGs, VFGs, and genomic functional characteristics of the C. freundii BAU_TM8 strain were identified using several open-source databases. ResultsThe sequence type of this strain was ST669, and the pathogenicity index of the strain was 0.919. Moreover, the strain had an estimated genome length of 5,797,806 bp, harboring 62 contigs, a G + C content of 54.32 %, and five contig L50s with an N50 value of 443,947 bp. Using phylogenetic analysis, this strain was closely related to two strains isolated from human and environmental samples in the USA and China despite huge geographical distances. The C. freundii BAU_TM8 strain consisted of 40 AGRs encoding resistance to 19 antimicrobial categories, e.g., fluoroquinolones, macrolides, folate pathway antagonists, aminoglycosides, tetracyclines, cephalosporins, and others. According to the phenotypic assay and genome sequence, the sensitivity and specificity of resistance profiles of the strain were 100 % and 20 %, respectively. Moreover, the virulence factor database detected 66 VFGs in this strain. This strain contained 1581 subsystems, having 33 % subsystem coverage and 2275 genes encoding amino acid derivatives, carbohydrate metabolism, protein metabolism, cofactors, vitamins, prosthetic groups, pigments, respiration, motility and chemotaxis, stress response, DNA metabolism, nucleosides and nucleotides, and others. ConclusionsTo the best of our knowledge, this is the first WGS report of C. freundii from a domesticated duck in Bangladesh. The ubiquitous occurrence of ARGs and VFGs in the C. freundii BAU_TM8 strain detected in this study highlights the growing concern about antimicrobial resistance in humans, animals, and environments.

A Whole-Genome Scan Revealed Genomic Features and Selection Footprints of Mengshan Cattle

A Whole-Genome Scan Revealed Genomic Features and Selection Footprints of Mengshan Cattle

Tracking clonal evolution of drug resistance in ovarian cancer patients by exploiting structural variants in cfDNA.

Drug resistance is the major cause of therapeutic failure in high-grade serous ovarian cancer (HGSOC). Yet, the mechanisms by which tumors evolve to drug resistant states remains largely unknown. To address this, we aimed to exploit clone-specific genomic structural variations by combining scaled single-cell whole genome sequencing with longitudinally collected cell-free DNA (cfDNA), enabling clonal tracking before, during and after treatment. We developed a cfDNA hybrid capture, deep sequencing approach based on leveraging clone-specific structural variants as endogenous barcodes, with orders of magnitude lower error rates than single nucleotide variants in ctDNA (circulating tumor DNA) detection, demonstrated on 19 patients at baseline. We then applied this to monitor and model clonal evolution over several years in ten HGSOC patients treated with systemic therapy from diagnosis through recurrence. We found drug resistance to be polyclonal in most cases, but frequently dominated by a single high-fitness and expanding clone, reducing clonal diversity in the relapsed disease state in most patients. Drug-resistant clones frequently displayed notable genomic features, including high-level amplifications of oncogenes such as CCNE1 , RAB25 , NOTCH3 , and ERBB2 . Using a population genetics Wright-Fisher model, we found evolutionary trajectories of these features were consistent with drug-induced positive selection. In select cases, these alterations impacted selection of secondary lines of therapy with positive patient outcomes. For cases with matched single-cell RNA sequencing data, pre-existing and genomically encoded phenotypic states such as upregulation of EMT and VEGF were linked to drug resistance. Together, our findings indicate that drug resistant states in HGSOC pre-exist at diagnosis and lead to dramatic clonal expansions that alter clonal composition at the time of relapse. We suggest that combining tumor single cell sequencing with cfDNA enables clonal tracking in patients and harbors potential for evolution-informed adaptive treatment decisions.

Genome characterization of Rift Valley fever virus isolated from cattle, goats and sheep during interepidemic periods in Kenya

Rift Valley fever virus (RVFV) is a mosquito-borne RNA virus of the Phlebovirus genus in the phenuviridae family. Its genome is trisegmented with small (S), medium (M) and large (L) fragments. In nature, the virus exists as a single serotype that is responsible for outbreaks of Rift Valley fever (RVF), a zoonotic disease that often occurs in Africa and the Middle East. RVFV genomes are thought to undergo both recombination and reassortment and investigations of these events is important for monitoring the emergence of virulent strains and understanding the evolutionary characteristics of this virus. The aim of this study was to characterize the genomes of RVFV isolates from cattle, sheep, and goats collected during an interepidemic period in Kenya between June 2016 and November 2021. A total of 691 serum samples from cattle (n = 144), goats (n = 185) and sheep (n = 362) were analysed at the Central Veterinary Laboratories. The competitive IgM-capture ELISA, was used to screen the samples; 205 samples (29.67%) tested positive for RVFV. Of the 205 positive samples, 42 (20.5%) were from cattle, 57 (27.8%) from goats, and 106 (51.7%) from sheep. All the IgM-positive samples were further analyzed by qPCR, and 24 (11.71%) tested positive with Ct values ranging from 14.788 to 38.286. Two samples, 201808HABDVS from sheep and 201810CML3DVS from cattle, had Ct values of less than 20.0 and yielded whole genome sequences with 96.8 and 96.4 coverage, respectively. There was no statistically significant evidence of recombination in any of the three segments and also phylogenetic analysis showed no evidence of reassortment in the two isolated RVFV segments when compared with other isolates of different lineages from previous outbreaks whose genomes are deposited in the GenBank. No evidence of reassortment leaves room for other factors to be the most probable contributors of change in virulence, pathogenicity and emergence of highly virulent strains of the RVFV.

Comprehensive physiological and genomic characterization of a potential probiotic strain, Lactiplantibacillus plantarum ILSF15, isolated from the gut of tribes of Odisha, India

Characterizing probiotic features of organisms isolated from diverse environments can lead to the discovery of novel strains with promising functional features and health attributes. The present study attempts to characterize a novel probiotic strain isolated from the gut of the tribal population of Odisha, India. Based on 16S rRNA-based phylogeny, the strain was identified as a species of the Lactiplantibacillus genus and was named Lactiplantibacillus plantarum strain ILSF15. The current investigation focuses on elucidating this strain’s genetic and physiological properties associated with probiotic attributes such as biosafety risk, host adaptation/survival traits, and beneficial functional features. The novel strain was observed, in vitro, exhibiting features such as acid/bile tolerance, adhesion to the host enteric epithelial cells, cholesterol assimilation, and pathogen exclusion, indicating its ability to survive the harsh environment of the human GIT and resist the growth of harmful microorganisms. Additionally, the L. plantarum ILSF15 strain was found to harbor genes associated with the metabolism and synthesis of various bioactive molecules, including amino acids, carbohydrates, lipids, and vitamins, highlighting the organism’s ability to efficiently utilize diverse resources and contribute to the host’s nutrition and health. Several genes involved in host adaptation/survival strategies and host-microbe interactions were also identified from the ILSF15 genome. Moreover, L. plantarum strains, in general, were found to have an open pangenome characterized by high genetic diversity and the absence of specific lineages associated with particular habitats, signifying its versatile nature and potential applications in probiotic and functional food industries.

Histopathologic, genomic, transcriptomic, and functional characteristics of eight melanocytic tumors with BRAF fusions showing stronger MAPK pathway activation compared to BRAF V600E tumors.

Activating BRAF gene alterations are central to melanocytic tumor pathogenesis. A small, emerging subset of melanocytic tumors driven by BRAF fusions has distinct therapeutic implications and has been described to have Spitzoid morphology patterns. However, such morphological patterns do not encompass all cases, and little is known about the functional molecular events. We conducted a retrospective search through our molecular archives to identify melanocytic tumors with BRAF fusions. We reviewed clinical, histopathological, and genomic features. We further explored transcriptomic and protein-level findings. Histopathologic patterns varied, with many cases without a distinctive pattern. We identified novel and diverse BRAF gene fusion partners. Differential transcriptomic analysis between low-risk BRAF fusion tumors and reference BRAF V600E tumors showed no differentially expressed genes. However, quantitatively stronger MAPK pathway activation of BRAF fusion tumors over BRAF V600E tumors was demonstrated by statistically significant stronger staining of p-ERK immunohistochemistry. Gene-specific RNA analysis shows comparable BRAF transcript levels between the two groups. The quantitatively stronger activation of the MAPK pathway of BRAF fusion tumors, instead of qualitatively different transcriptomes, may account for the morphology difference from conventional BRAF V600E tumors. BRAF fusions likely act through dysregulated protein function rather than RNA upregulation related to the characteristics of the fusion partners.

Genomic Characterization of a Dog-Mediated Rabies Outbreak in El Pedregal, Arequipa, Peru.

Rabies, a re-emerging zoonosis with the highest known human case fatality rate, has been largely absent from Peru, except for endemic circulation in the Puno region on the Bolivian border and re-emergence in Arequipa City in 2015, where it has persisted. In 2021, an outbreak occurred in the rapidly expanding city of El Pedregal near Arequipa, followed by more cases in 2022 after nearly a year of epidemiological silence. While currently under control, questions persist regarding the origin of the El Pedregal outbreak and implications for maintaining rabies control in Peru. We sequenced 25 dog rabies virus (RABV) genomes from the El Pedregal outbreak (n=11) and Arequipa City (n=14) from 2021-2023 using Nanopore sequencing in Peru. Historical genomes from Puno (n=4, 2010-2012) and Arequipa (n=5, 2015-2019), were sequenced using an Illumina approach in the UK. In total, 34 RABV genomes were analyzed, including archived and newly obtained samples. The genomes were analyzed phylogenetically to understand the outbreak's context and origins. Phylogenomic analysis identified two genetic clusters in El Pedregal: 2021 cases stemmed from a single introduction unrelated to Arequipa cases, while the 2022 sequence suggested a new introduction from Arequipa rather than persistence. In relation to canine RABV diversity in Latin America, all new sequences belonged to a new minor clade, Cosmopolitan Am5, sharing relatives from Bolivia, Argentina, and Brazil. Genomic insights into the El Pedregal outbreak revealed multiple introductions over a 2-year window. Eco-epidemiological conditions, including migratory worker patterns, suggest human-mediated movement drove introductions. Despite outbreak containment, El Pedregal remains at risk of dog-mediated rabies due to ongoing circulation in Arequipa, Puno, and Bolivia. Human-mediated movement of dogs presents a major risk for rabies re-emergence in Peru, jeopardizing regional dog-mediated rabies control. Additional sequence data is needed for comprehensive phylogenetic analyses.

Quantum mechanical electronic and geometric parameters for DNA k-mers as features for machine learning

We are witnessing a steep increase in model development initiatives in genomics that employ high-end machine learning methodologies. Of particular interest are models that predict certain genomic characteristics based solely on DNA sequence. These models, however, treat the DNA as a mere collection of four, A, T, G and C, letters, dismissing the past advancements in science that can enable the use of more intricate information from nucleic acid sequences. Here, we provide a comprehensive database of quantum mechanical (QM) and geometric features for all the permutations of 7-meric DNA in their representative B, A and Z conformations. The database is generated by employing the applicable high-cost and time-consuming QM methodologies. This can thus make it seamless to associate a wealth of novel molecular features to any DNA sequence, by scanning it with a matching k-meric window and pulling the pre-computed values from our database for further use in modelling. We demonstrate the usefulness of our deposited features through their exclusive use in developing a model for A->C mutation rates.

Recurrent selection shapes the genomic landscape of differentiation between a pair of host-specialized haplodiploids that diverged with gene flow.

Understanding the genetics of adaptation and speciation is critical for a complete picture of how biodiversity is generated and maintained. Heterogeneous genomic differentiation between diverging taxa is commonly documented, with genomic regions of high differentiation interpreted as resulting from differential gene flow, linked selection and reduced recombination rates. Disentangling the roles of each of these non-exclusive processes in shaping genome-wide patterns of divergence is challenging but will enhance our knowledge of the repeatability of genomic landscapes across taxa. Here, we combine whole-genome resequencing and genome feature data to investigate the processes shaping the genomic landscape of differentiation for a sister-species pair of haplodiploid pine sawflies, Neodiprion lecontei and Neodiprion pinetum. We find genome-wide correlations between genome features and summary statistics are consistent with pervasive linked selection, with patterns of diversity and divergence more consistently predicted by exon density and recombination rate than the neutral mutation rate (approximated by dS). We also find that both global and local patterns of FST, dXY and π provide strong support for recurrent selection as the primary selective process shaping variation across pine sawfly genomes, with some contribution from balancing selection and lineage-specific linked selection. Because inheritance patterns for haplodiploid genomes are analogous to those of sex chromosomes, we hypothesize that haplodiploids may be especially prone to recurrent selection, even if gene flow occurred throughout divergence. Overall, our study helps fill an important taxonomic gap in the genomic landscape literature and contributes to our understanding of the processes that shape genome-wide patterns of genetic variation.

Exploring the functional and genomic features of Cellulophaga lytica NFXS1, a zeaxanthin and lytic enzyme-producing marine bacterium that promotes microalgae growth

Beneficial microalgae-bacteria interactions mediate important marine ecological processes and present potential for exploitation in the industrial production of microalgae. Members of the bacterial genus Cellulophaga are known for their ability to interact with several marine organisms, including microalgae, however, not much is understood about their properties. In this work we present the detailed characterization of the marine bacterium, Cellulophaga lytica NFXS1, including its ability to synthesize valuable compounds such as zeaxanthin, its ability to modulate the growth and productivity of two relevant marine microalgae, Nannochloropsis oceanica and Phaeodactylum tricornutum, and the analysis of its complete genome sequence. The obtained results showed that C. lytica NFXS1 used several complex carbohydrates, indicating strong lytic enzymatic activities. Additionally, co-cultivation studies revealed that C. lytica NFXS1 colonized the microalgae phycosphere, catabolized microalgae exudates, promoted the growth of microalgae, and increased fucoxanthin productivity in P. tricornutum. Genomic analysis revealed the basis for the microalgae growth-promoting activities of strain NFXS1, which are mostly mediated by gliding motility, biofilm formation, biosynthesis of phytohormones (auxins, cytokinins), polyamine biosynthesis, and resistance to oxidative stress induced by zeaxanthin production and ROS detoxifying enzymes. A vast array of genes encoding glycosyl hydrolases and pectate lyase family enzymes were also found in the NFXS1 chromosome. Overall, these findings brought new information regarding the properties of C. lytica NFXS1, which can be used as a model organism to study microalgae-bacteria interactions and is a strong candidate for the development of a wide range of marine biotechnological solutions.

Insights into the global genomic features of Salmonella enterica serovar Gallinarum biovars Gallinarum and Pullorum.

Characterize global genomic features of 86 genomes of Salmonella Gallinarum (SG) and Pullorum (SP), which are important pathogens causing systemic infections in poultry. All genomes harbored efflux pump encoding gene mdsA and gold tolerance genes golS and golT. Aminoglycoside (aac(6')-Ib, aadA5, aph(6)-Id, aph(3'')-Ib, ant(2'')-Ia), beta-lactam (blaTEM-1, blaTEM-135), efflux pump (mdsB), fosfomycin (fosA3), sulfonamide (sul1, sul2), tetracycline [tet(A)], trimethoprim (dfrA17), acid (asr), and disinfectant (qacEdelta1) resistance genes, gyrA, gyrB, and parC quinolone resistance point mutations, and mercury tolerance genes (mer) were found in different frequencies. Additionally, 310 virulence genes, pathogenicity islands (including SPI-1, 2, 3, 4, 5, 6, 9, 10, 12, 13, and 14), plasmids [IncFII(S), ColpVC, IncX1, IncN, IncX2, and IncC], and prophages (Fels-2, ST104, 500465-1, pro483, Gifsy-2, 103203_sal5, Fels-1, RE-2010, vB_SenS-Ent2, and L-413C) were detected. MLST showed biovar-specific sequence types, and core genome MLST showed country-specific and global-related clusters. SG and SP global strains carry many virulence factors and important antimicrobial resistance genes. The diverse plasmids and prophages suggest genetic variability. MLST and cgMLST differentiated biovars and showed profiles occurring locally or worldwide.