Genomics

Nucleobindin 1 (NUCB1) is a multifunctional conserved protein located in Golgi luminal, nucleus, extracellular and cytosolic pools. NUCB1 is multidomain protein comprised of a signal peptide, a DNA-binding domain, a leucine zipper and Ca2+ -binding domain. The multiple domains and localization of NUCB1 potentiates its interactions with various partners, such as DNA, Gi3 protein, cyclooxygenase 2, LRP10 and RNA suggests its importance in the regulation of many cellular events. We revealed that NUCB1 contains three RNA-binding regions and able to interact with two RNA fragments. It was suggested possible variants of the participation of NUCB1 in the interaction of the two partially complementary RNAs. The RNA-binding properties of the NUCB1 were also confirmed in vivo experiments.

We recently reported the identification of endogenous viral elements (EVEs) originating from the Caulimoviridae family within the alfalfa (Medicago sativa L.) genome. Our subsequent identification of ubiquitous rhabdoviral elements in infected and healthy alfalfa tissues by high throughput sequencing prompted us to suggest that the alfalfa genome might be populated with integrated rhabdoviruses as well. Bioinformatics analysis using 26 publicly available alfalfa genomes proved the suggestion accurate. We found multiple non-retroviral segments of the Rhabdoviridae family belonging to the genera Betanucleorhabdovirus and Betacytorhabdovirus that appeared to be stable constituents of the host genome. In that capacity they could potentially acquire functional roles in alfalfas development and response to environmental stresses. We believe this study reveals the first documented case of rhabdoviruses integrated into the alfalfa genome.

Currently, the genetic architecture of Middle Eastern populations is underrepresented in global genomic databases. This gap increases the rate of Variants of Uncertain Significance (VUSs) and clinical misinterpretations of genomic data especially in Middle Eastern populations. Whole exome sequencing was conducted on 90 healthy individuals from Jordan and the data were analysed using Principal Component Analysis (PCA) and multi-computational filtering. PCA revealed a double ancestry (EUR-AFR) admixture rather than a triple admixture (EUR-AFR-AMR). More than 3,500 populations-specific variants (PSVs) were identified, of which 72% were singletons. Additionally, 19 variants were significantly enriched compared to the maximum allele frequencies in public global databases (Fisher's exact test with Benjamini-Hochberg false discovery rate correction, p-value < 0.05). Consequently, the results suggest the reclassification of variants of Uncertain Significance (VUS) which reside in the ECE2 gene to likely benign and the variants of Conflicting Classification of Pathogenicity in the genes IL1RN and THPO to benign based on the significant allele frequency (AF=0.0389, p-value < 0.05). Furthermore, a pathogenic ClinVar variant was identified in a healthy individual, warranting careful interpretation. The findings underscore the importance of identifying PSVs in order to minimize or even prevent clinical misdiagnosis and highlight the unique genetic signature in Jordan. The study serves as a foundational resource for precision medicine in the region.

IncC plasmids played an important role in driving antimicrobial drug resistance development in the seventh pandemic Vibrio cholerae O1 El Tor (7PET) lineage during the 1970s and these have begun to re-emerge. In this study, we investigated a comprehensive dataset for 28 complete IncC plasmids -- including 17 newly sequenced plasmid genomes -- from 7PET isolates collected on various continents between 1979 and 2024 and distributed across the global phylogenetic tree for 7PET isolates. IncC type 2 predominated among the V. cholerae plasmids studied, and five new core genome sequence types (cgSTs) were identified. The antimicrobial resistance genes (ARGs) were arranged in islands inserted at specific hotspots within the common IncC backbone and were significantly associated with IncC types or islands. The IncC plasmid backbone has remained stable over the last 50 years, but the ARGs and their associated genomic islands displayed remarkable diversity, underscoring the complex evolution patterns of the 7PET lineage of V. cholerae and its considerable adaptability under selective pressure.

BackgroundMetazoan adenosine-to-inosine (A-to-I) mRNA editing temporospatially diversifies the neuronal transcriptome and proteome. The limited read length from next-generation sequencing (NGS) constrains the quantification of the potentially differential editing levels across different splicing isoforms, restricting our understanding of the extent to which RNA editing contributes to molecular diversity and its interplay with splicing. MethodsWe employed reverse transcription nested PCR (RT-nPCR) and developed a novel interfering-Primer PCR (iPrimer PCR) technique to distinguish different transcripts of any gene. We selected multiple essential genes exhibiting RNA editing in coding sequences (CDSs) or untranslated regions (UTRs) for isoform-specific amplification and Sanger sequencing. ResultsNine different Adar isoforms together with pre-mRNA had distinct editing levels at the S>G auto-recoding site, which was predicted to have isoform-specific effects on catalytic activities. Although pre-mRNA editing might exert isoform-dependent promotion/suppression of splicing, closely located editing sites, such as those in neuronal genes qvr and stj, still exhibited high correlation in editing levels due to co-editing. iPrimer strategy further discovered differential recoding levels between the long/short 3UTR isoforms of gene jef. ConclusionsWe provide the first comprehensive solution for isoform-specific PCR amplification of any gene, enabling quantification of RNA editing level of different isoforms. Our results offer insights into how RNA editing interplays with splicing, and highlight its complicated role in expanding molecular diversity. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=79 SRC="FIGDIR/small/725286v1_ufig1.gif" ALT="Figure 1"> View larger version (17K): org.highwire.dtl.DTLVardef@1ebc82org.highwire.dtl.DTLVardef@1ea365dorg.highwire.dtl.DTLVardef@1971aceorg.highwire.dtl.DTLVardef@160d053_HPS_FORMAT_FIGEXP M_FIG C_FIG We developed isoform-specific PCR followed by Sanger sequencing, and achieved the quantification of differential RNA editing levels in different transcripts of a gene.

BackgroundThe Pacific oyster (Crassostrea/Magallana gigas) is increasingly recognised as a model marine invertebrate. Valued for both ecological and commercial importance, Pacific oysters are farmed widely, supporting global food security by providing a sustainable nutrient-rich source of protein. Despite the significant and recurring economic losses caused by Ostreid herpesvirus (OsHV-1) outbreaks, only a limited number of studies have examined host-pathogen interplay at single-cell resolution. The few available studies largely focus on circulating immune cells (haemocytes), thereby overlooking the complexity of host responses across different tissues and organs. ResultsWe present a detailed single-nucleus transcriptomic atlas of the whole Pacific oysters, including during OsHV-1 infection. A total of 18 distinct transcriptomic clusters were resolved, capturing major cell populations from the gill, mantle, hepatopancreas, adductor muscle, and haemocytes. Notably, three populations- gill ciliary cells, hepatopancreas cells, and an immune-enriched cluster 1- exhibited pronounced transcriptomic responses to OsHV-1 infection. Across the 6, 24, 72, and 96 hours post-infection (hpi) time course, viral transcripts were detected almost exclusively at 72 hpi, with enrichment primarily in adductor muscle cells and two immune cell populations- immature haemocytes, and hyalinocytes. ConclusionsOur findings suggest potential entry portals and tissue-specific replication sites for the OsHV-1 virus in Pacific oysters. This atlas resource provides a high-resolution cellular framework for understanding host-virus interactions and establishes a foundation for future investigations into herpesvirus pathogenesis in marine invertebrates.

Magnaporthe oryzae, the rice blast fungus, plays a role as a model organism for molecular plant-microbe interaction research. Studies on the pathogenic mechanism of this fungus revealed many genes involved in signaling pathways. As multi-omics data are being available, genomic-level researches have been conducted to uncover the underlying biological processes during the pathogenesis of M. oryzae. Identifying the genome-wide protein-protein interaction (PPI) network is one of the omics-level approaches, which helps to understand signaling and regulatory pathways. However, existing biological network resources of M. oryzae are not sufficient to decipher pathogenesis mechanisms due to the abundance of false positives/negatives. In this study, a reliable PPI network database of M. oryzae, MagNet, was constructed with three methods, including homology-based Interolog search, co-expression network construction, and domain-domain interaction (DDI)-based prediction. With three approaches altogether, the pan-network with 5,600,976 interactions was generated, including 217,531 highly confident interactions supported by all three methods. Experimental data on M. oryzae PPIs supported that our PPI network can predict PPIs with higher accuracy compared to the previously constructed databases. MagNet would provide integrated biological network data, which can help to understand the molecular mechanisms of the rice blast fungus. The PPI data can be accessed via https:/magnet.scnu.ac.kr.

Introduction Recent respiratory illness, especially influenza, may trigger acute cardiac events via elevated inflammatory mediators. During the 2018 influenza season in Bangladesh, this study examined whether recent acute clinical respiratory illness (CRI) or laboratory-confirmed influenza was associated with elevated hs-CRP and IL-6, linked to acute cardiac events. Methods A total of 139 participants aged [&ge;]40 were recruited from a Dhaka cardiac hospital: 70 with acute myocardial infarction (AMI), 30 with other acute cardiac events, and 39 healthy individuals. CRI was defined as fever with cough and/or respiratory symptoms within seven days. Respiratory swabs were tested for influenza, and blood was analyzed for hs-CRP and IL-6. Results Median hs-CRP and IL-6 were higher in participants with CRI or influenza but not significantly. Cardiac patients had elevated hs-CRP (9.98 mg/L in other cardiac; 4.86 mg/L in AMI vs. 1.73 mg/L in healthy) and IL-6 (0.1 pg/mL in other cardiac; 0.145 pg/mL in AMI vs. 0.08 pg/mL in healthy) (p<0.001). CRI was not significantly associated with elevated hs-CRP or IL-6, though influenza in healthy participants was linked to higher IL-6. Cardiac patients had a higher risk of hs-CRP [&ge;]3 mg/L and elevated IL-6. Conclusion Cardiac patients showed significantly increased inflammatory markers, but CRI was not clearly linked to inflammation. Further research should assess biomarker utility for early cardiac risk.

Here, we present a comprehensive multiomics analysis of anthocyanin biosynthesis in Rubus armeniacus, known for its dark fruits. A phased genome sequence of the tetraploid blackberry was generated, achieving an N50 of 34 Mb with an assembly size of 1.2 Gbp based on Oxford Nanopore Technology sequencing (ONT). The BUSCO score for the total assembly shows a high completeness of 99.1%. The assembly was separated into 4 pseudohaplophases, with the pseudohaplophase A representing the R. armeniacus genome in 7 chromosome scale contigs, with an N50 of 46 Mbp and 98.8% conserved BUSCO genes. A total of 118,183 protein coding genes were annotated within the genome assembly and all relevant genes encoding enzymes and transcriptional regulators of the anthocyanin biosynthesis pathway were identified within each pseudohaplophase. To further understand the underlying cause of dark pigmentation, the gene expression was analysed during different stages of berry development revealing a strong induction of anthocyanin biosynthesis genes including the anthocyanin activating subgroup 6 MYB transcriptions during the berry ripening process. Further, a quantification of cyanidin-3-O-glucoside in methanolic berry extract, utilizing a UHPLC-HRAM-MS analysis, revealed an approximately 500-fold increase of cyanidin-3-O-glucoside from green to black fruit, indicating that dark pigmentation in R. armeniacus results from high anthocyanin accumulation. Significance statementThis study provides a multiomics analysis of the dark pigmentation of Rubus armeniacus, including a high quality phased assembly and an in-depth analysis of the anthocyanin biosynthesis pathway. A transcriptional and metabolomic analysis revealed that dark berry pigmentation is caused by a high accumulation of cyanidin-3-O-glucoside during fruit ripening.

Vitamin D signaling has recognized roles in female reproductive physiology, but its effects at the chromatin level in endometrial stromal cells are still unclear. Here, we investigated how the active form of vitamin D, 1,25-dihydroxyvitamin D3, or calcitriol, influences the accessible chromatin landscape of human endometrial stromal cells. Assay for transposase-accessible chromatin using sequencing (ATAC-seq) was performed on T-HESCs treated with either a vehicle or 1,25(OH)2D3. Ligand treatment increased overall chromatin accessibility, shown by higher ATAC-seq signal intensity, while causing only minor changes in the total number of called peaks. Peak annotation revealed that accessible regions were spread across both promoter-proximal and distal genomic areas. Integrating this data with CUT&RUN and RNA sequencing showed that most vitamin D-responsive cistromic modifications and transcripts were linked to nearby open chromatin, though fewer were associated with regions that were significantly differentially accessible. These results suggest that 1,25(OH)2D3-dependent transcription mainly occurs within a permissive, pre-accessible chromatin environment. This study offers new evidence that active vitamin D influences the epigenomic landscape of human endometrial stromal cells, establishing the chromatin-based molecular response to a chemically-defined VDR ligand, 1,25(OH)2D3, relevant to stromal differentiation and preparation for decidualization. HighlightsO_LIFirst evidence suggesting the direct impact of active vitamin D, 1,25-dihydroxyvitamin D3, 1,25(OH)2D3, enhanced the signal intensity of chromatin accessibility in human endometrial stromal cells C_LIO_LIMost accessible chromatin regions were shared between vehicle and ligand-treated human endometrial stromal cells C_LIO_LI1,25(OH)2D3-responsive transcription occurs largely within pre-accessible chromatin in human endometrial stromal cells C_LIO_LIAssay for transposase-accessible chromatin sequencing (ATAC-seq) defines a chromatin-level pharmacologic response to a chemically defined VDR ligand in human endometrial stromal cells C_LI

The high efficiency of genome editing presents a challenge when modifying genes associated with viability, welfare, or fertility issues, as implementation of the technology frequently results in mosaic animals with bi-allelic mutations. Combining deactivated Cas9 (dCas9) with Cas9 has been proposed as a strategy to protect one of the two target alleles from editing. We piloted this strategy with 11 genes that are reported as homozygous lethal or associated with welfare issues. We showed that the viability of founders was significantly increased when using 80:20 or 90:10 dCas9:Cas9 ratios, whereas the 70:30 ratio did not yield an equivalent protective effect. The associated overall production rate of mutated founder per manipulated embryo was significantly higher for the 80:20 ratio. Concomitantly, an increased proportion of dCas9 was associated with a significant increase in retention of unedited target alleles but, importantly, did not hinder germline transmission. In addition, editing genes in a paralog cluster with a combination of dCas9 and Cas9 reduced unwanted off-target editing, illustrating a further potential applicability of this approach. This study defines the optimal ratio between dCas9 and Cas9 for strategies aimed at achieving mono-allelic mutations within mosaic founders and proposes a means to reduce the incidence of off-target effects in experiments with limited gRNA options.

Pseudomonas aeruginosa, a Gram-negative opportunistic pathogen is well known for life-threatening acute infections among the human population. The bacterium can withstand most antibiotics by using their high levels of inherent and acquired resistance mechanisms such as Biofilm-EPS, Persistence, and Quorum sensing (QS). Owing to the importance of adaptive antibiotic multi-drug resistance of P. aeruginosa, the current investigation is aimed to explore the phytochemicals derived from mangrove plants as potential agents to control biofilm and drug resistance mechanisms through a multi-mechanistic computational approach. For identifying potential compounds and target, In-silico drug repurposing technique is implemented by docking/virtual screening of 49 phytochemical compounds against 18 proteins involved in the Persister Cell formation, QS, and EPS synthesis in P. aeruginosa which resulted the proteins RelA and SpoT (persistence), PqsA, and PqSR (QS), and PelA and PelB (EPS synthesis) and compounds Taraxerone and Taraxerol to be potential. The results of docking were well corroborated with MD simulations. These targets and compounds explored through in-silico approach, are found to target potential antimicrobial pathways involving EPS synthesis, persistence genes, and QS, aiming to enhance antibiotic efficacy. Further, this study could be reference for in-vivo and in-vitro investigations to evaluate the further effectiveness of the compounds and potentiality of the proteins for MDR therapeutics of P. aeruginosa.

Selective deployment of multiple transcription start sites is a major regulatory feature of human transcriptomes. FANTOM CAGE data exhibit a near-universal TSS deployment parsimony which is disrupted in cancers. We have recently shown that TSS deployment is sensitive to gene function, futile upstream transcription, and cellular biosynthetic states. Patterns in FANTOM CAGE data can reveal mechanisms underlying TSS co-deployments. We propose and test the possibility that some TSSs act like epromoters and act as co-varying hubs of transcriptional activities for multiple other promoters. Using deep analysis of CAGE data implemented through neural networks we show that non-cancers implement transcription co-deployments through cores of epromoter-like TSSs which are generally proximal to their start codons. These TSSs show enhancer-like TFBSs profiles. A comparison with cancer CAGE data shows that the concentrated epromoter core is disrupted in cancers with multiple distal TSSs replacing the proximal TSS cores. We provide evidence that the core TSSs are rich in YY1 and CTCF binding sites and associated with genes coding for transcription factors. Our findings show that covariance of TSS deployment is sensitive to transcriptional resource cost and a parsimonic design of TSS co-deployments depends on proximal TSSs in non-cancers, a mechanism grossly disrupted in cancers. HighlightsO_LIHeterogeneous FANTOM CAGE data contains universal patterns of TSSs co-deployments. C_LIO_LITSS co-deployments exhibit a parsimonious "core-covariant" scheme which is disrupted in cancers. C_LIO_LICore TSSs are enriched in transcription factor binding sites and gene functions which justify biological features of the samples. C_LIO_LIThe DL pipeline we present identifies the core-covariant TSS sets in an unbiased manner. C_LI

Upstream open reading frames (uORFs) in the 5' leader of bacterial mRNAs can modulate gene expression, yet genome-wide identification remains limited. We combined bioinformatic prediction of ribosome-binding sites (RBSs) - a Shine-Dalgarno sequence and a start codon - with experimental validation to uncover new uORFs in Sinorhizobium meliloti 2011. From totally 1106 predicted upstream RBSs (uRBSs), we first examined 15 candidates using eGFP reporters and integrating existing RNA-seq and Ribo-seq data. Translation was detected at 13 sites, with fluorescence intensity broadly correlating with predicted initiation rates. Two uRBSs correspond to gene start sites, thereby refining gene annotations. In nine cases, uRBS mutations affected downstream gene expression in reporter fusions. Among others, the data suggests that a Type I secretion system operon, the RNA chaperone gene hfq, and metabolic genes are regulated by uORFs. Four uORFs acted through translational coupling. We also identified uRBSs that were ribosome-occupied yet (nearly) silent in eGFP assays, and closely spaced to the downstream main RBS (mRBS). These uRBSs probably mediate ribosomal occlusion downregulating lacR and SM2011_RS36230. A re-screen of the prediction set revealed 335 close uRBS/mRBS pairs. Three of them were analyzed, supporting the proposed ribosomal occlusion mechanism for SM2011_RS03630 and SM2011_RS22110, while for glnK translational coupling to an uORF was suggested. These results indicate that uORFs are more widespread in bacteria than previously recognized and suggest that direct ribosomal occlusion of the mRBS is a novel mechanism for down-regulating protein synthesis.

Biological sequences are known to be not random. Thus, the comparison of in silico restriction fragment distributions of random and biological sequences may be an indicator of this non-randomness. Our analyses show that for most of the tested combinations of restriction enzyme and genome sequence the fragments per Megabase of the biological sequence deviate at least more then 10% from the corresponding random sequence. This deviation goes into both directions, i.e. clearly increased values are as common as clearly decreased values. Although there is no species- or restriction-enzyme-specific effect, a clear impact of the GC content both of the restriction site and of the genome sequence can be seen. In contrast to the random sequences, the genome sequences show distinct peaks in their fragment length distributions, hinting to repetitive elements such as transposons.

Cadmium, being a highly toxic metal, perturbs cellular homeostasis by forming stable complexes with numerous thiol-active proteins, ultimately leading to severe liver and lung damage. Despite its well-documented toxicity, the molecular mechanisms governing cadmium export remain poorly understood. Given the chemical similarity between cadmium and copper, we investigated whether the canonical copper-exporting ATPases, ATP7A and ATP7B participate in cadmium handling. Upon Cd treatment in hepatocytes, ATP7B undergoes trafficking to lysosomes via the retromer complex, as also observed in the case of elevated copper, accompanied by the upregulation of acidic lysosomal populations. In contrast, ATP7A expressed in lung adenocarcinoma cells, though exhibit vesicular redistribution upon Cd exposure, does not mediate lysosomal sequestration, suggesting distinct deployment of late secretory pathways by the two copper ATPases in response to cadmium. We have also observed that ATP7B-/- hepatocytes exhibit increased sensitivity to Cd exposure compared to wild-type cells. Whereas, overexpressing the ATP7B amino-terminal copper-binding domain in bacteria alleviates cadmium-induced stress, indicating its capacity to sequester Cd. Caenorhabditis elegans lacking copper-ATPase cua-1, displayed increased Cd sensitivity, while mutants (glo-1-/-), deficient in lysosome-related organelles (LRO), and (lmp-1-/-), deficient in lysosomal membrane glycoprotein, showed reduced resistance to cadmium toxicity. Treatment of the worm with cadmium increases the abundance of lysosomes marked by elevation in lysosomal biogenesis and functional genes, reinforcing the importance of lysosomal pathways in cadmium detoxification. To summarise, we delineated the non-canonical role of copper ATPases and lysosomes in cadmium-induced cellular toxicity.

Prenylated isoflavonoids are widely distributed specialized metabolites within the Fabaceae and contribute to various characteristic biological activities for both plants and humans. Several aromatic prenyltransferases (PTs) have been identified in Glycyrrhiza species, which are the most widely consumed crude drugs in traditional Chinese medicine. However, these enzymes do not sufficiently explain the structural diversity of prenylated flavonoids produced in the Glycyrrhiza genus. To identify additional novel PTs, we used elicited cultured Glycyrrhiza glabra roots as source material, in which elicitor treatment of cultured roots increased the accumulation of multiple prenylated flavonoids. To identify the responsible enzyme, PT candidates were screened using G. uralensis transcriptomes, currently the sole publicly available transcriptomic resource within the genus, and a homolog designated GgBSPT1 (BSPT; a broad-substrate prenyltransferase) was subsequently isolated from elicited cultured G. glabra roots. GgBSPT1 differed from previously identified Glycyrrhiza PTs in both amino acid sequence and enzymatic properties. GgBSPT1 catalyzed 3'-prenylation of isoliquiritigenin and 6-prenylation of five flavonoids, i.e., this PT displayed broad substrate acceptance across 20 distinct flavonoid structures. Overall, elicited cultured G. glabra roots enabled the identification of a previously unrecognized PT that is functionally distinct from earlier reported Glycyrrhiza PTs. This study provides a new insight into the metabolic plasticity of Glycyrrhiza species and expands the enzymatic toolkit for future metabolic engineering of prenylated phytochemicals by the unusually broad substrate specificity of GgBSPT1. Main conclusionUsing cultured Glycyrrhiza glabra roots, we identified a new prenyltransferase involved in the formation of a variety of flavonoids, thereby revealing novel prenylated isoflavonoid pathways in licorice.

Gene model for the ortholog of Lst8 (Lst8) in the May 2011 (WUGSC dyak_caf1/DyakCAF1) Genome Assembly (GenBank Accession: GCA_000005975.1) of Drosophila yakuba. This ortholog was characterized as part of a developing dataset to study the evolution of the Insulin/insulin-like growth factor signaling pathway (IIS) across the genus Drosophila using the Genomics Education Partnership gene annotation protocol for Course-based Undergraduate Research Experiences.

Cryptosporidium parvum is a protozoan parasite responsible for cryptosporidiosis, significantly threatening immunocompromised individuals, particularly HIV/AIDS patients, by causing severe diarrhea and potential mortality. Current treatments are largely ineffective, prompting investigations into new therapeutic options. This study evaluated two antiparasitic drugs: Mebendazole, used for helminth infections, and Artemisinin, used for malaria. The SKSR gene family encodes virulence factors in C. parvum, and Calcium-dependent protein kinase1 (CpCDPK1) regulates the life cycle of C. parvum; targeting these proteins may reduce growth and infection in hosts. In the current study, molecular docking was conducted taking Mebendazole and Artemisinin drugs as ligands, SKSR gene family and CpCDPK1 proteins as drug targets. Results with SKSR showed binding energy of -4.9 kcal/mol, -6.72 kcal/mol for Mebendazole and Artemisinin, respectively. Whereas, with CpCDPK1, the binding energies were -6.44 kcal/mol, -9.18 kcal/mol for Mebendazole and Artemisinin, respectively. Docking of Nitazoxanide (an in-use drug for C. parvum) with SKSR and CpCDPK1 revealed binding energies -4.2 kcal/mol, -4.81 kcal/mol, respectively. The stability of the proteins (targets) upon binding to the ligands was assessed by performing all-atom MD simulations for 100ns using the GROMACS package. No major variations were observed upon binding of Artemisinin and Mebendazole to SKSR and CpCDPK1. The findings of MD simulations imply that both proteins maintain their stability upon binding of Artemisinin and Mebendazole. Molecular Docking and MD simulation studies suggest that Artemisinin and Mebendazole are potential candidates for repurposing in the treatment of C. parvum infections, with recommendations for in vitro studies to validate these findings.

The incidence of vancomycin-resistant Enterococcus (VRE) is rising in hospitals in Canada, and resistance to last-resort antimicrobials including linezolid complicates treatment options for multidrug-resistant isolates. Recent reports from around the globe indicate that both linezolid and vancomycin resistance genes can be co-carried and mobilized by linear plasmids (named pELF) in Enterococcus species, often on the same backbone. We aimed to investigate linezolid resistance and linear plasmid prevalence in VRE bloodstream infection isolates collected by the Canadian Nosocomial Infection Surveillance Program from 2009 to 2024. We found that screening for pELF linear plasmid ends in short reads was a reliable way to predict linear plasmid presence in large-scale surveillance data (100 % accuracy on 85 reference samples). Almost half of the isolates in our collection were predicted to carry pELF plasmids (45.4 %, 941/2071) and we found that this proportion has increased from 2018 (32.2 %, 59/183) to 72 % of isolates between 2021 and 2024 (2021: 68.5 % (115/168); 2022: 71.6 % (146/204); 2023: 72.8 % (166/228); 2024: 71.6 % (235/328)). This trend of increasing linear plasmid carriage is evident from 2018 to 2024 across the dominant emerging sequence types (ST80, ST17, ST117). Linezolid resistance based on phenotypic antimicrobial susceptibility testing was low (1.0 %, 21/2071). Using long read sequencing, we characterized the linezolid resistant isolates and confirmed pELF plasmid presence in 13/21 (61.9 %) isolates. Six isolates harboured pELF plasmids encoding linezolid resistance genes (optrA, cfr(D), poxtA) and five of these also encoded vancomycin resistance genes (vanA). We compared these six plasmids to 39 public plasmid sequences and clustered them using MOB-suite and pling. Overall, this study provides further examples of the co-carriage of vancomycin and linezolid resistance genes on mobile linear plasmids and shows that linear plasmid prevalence is detectable and increasing across VRE in Canada. IMPACT STATEMENTGiven the increasing prevalence of multidrug-resistant hospital-acquired pathogens, resistance to last-resort antibiotics is a global public health threat. Linezolid is a last-resort antibiotic used to treat vancomycin-resistant Enterococcus isolates, and the dissemination of linezolid resistance genes is significantly facilitated by mobile elements that can transfer between unrelated strains and species. Linezolid resistance genes have recently been described on linear plasmids and are often co-localized with other resistance genes on the same plasmid backbone. Consequently, understanding the features and distribution of linear plasmids and those harbouring linezolid resistance genes is crucial for pathogen surveillance and mitigation of resistance. In this work, we used long-read and short-read sequencing to characterize genomic epidemiology of linear plasmids across 16 years of Enterococcus surveillance data in Canada. This study furthers knowledge of linear plasmids by demonstrating that they are relatively common across vancomycin-resistant Enterococcus blood isolates and by providing more examples of co-localized vancomycin and linezolid resistance genes on the same linear plasmid backbone. DATA SUMMARYSequencing data and genome sequences were deposited in National Centre for Biotechnology BioProject PRJNA1279082, and accessions are listed in Table S1. Supplementary materials for this study are available at the Figshare portal through DOI: XXX.