Genomics

Background: Pancreatic ductal adenocarcinoma (PDAC) is projected to become the second leading cause of cancer-related mortality, with an increasing incidence and dismal 10% survival at 5 years despite best available care. PDAC is characterized by low tumour cellularity, recurrent driver mutation alterations, and substantial inter and intra-tumoural heterogeneity. The influence of individual genetic features such as sex and ancestry on tumour heterogeneity remains poorly understood. Methods: We assembled a cohort of PDAC across five datasets for a total of 7,527 tumours. Measures of genomic instability were assessed across sex and ancestry on a global level (SNV density and percent genome altered by copy number events), and at a gene level for recurrent events. KRAS variant frequencies, ploidy and mRNA abundance were also evaluated. Ancestry was assigned using ADMIXTURE from germline genotypes where available, and patient reported otherwise. Results: We observed ancestry-specific mutation patterns in global point mutation density, with higher SNV burden in African PDAC tumours and lower SNV burden in Asian PDAC tumours as compared to European tumours. We also identified an ancestry bias in the driver gene, GNAS, in which driver SNVs in GNAS was detected more frequently in those of Asian ancestry, and less frequently in those of African ancestry. No sex- or ancestry-associations were found for KRAS, but individual KRAS variants showed a trend towards ancestry-related variation. Conclusions: These findings demonstrate that ancestry contributes to the genomic heterogeneity of PDAC, shaping both global and gene-specific mutation patterns. This work underscores the importance of including diverse ancestries in PDAC molecular profiling efforts which may inform the development of tailored molecular diagnostics and targeted therapies.

ABSTRACT Objective: Skeletal muscle growth rate differs substantially among diverse pig breeds. In the present study, RNA-seq was performed to identify differentially expressed genes (DEGs) across pig breeds with distinct growth rates. CCND2 showed obvious differential expression among pig breeds. Although CCND2 is essential for cell cycle control, its biological function in skeletal muscle differentiation has rarely been reported. This work aims to identify DEGs in pig breeds with different growth rates and explore CCND2 role in porcine skeletal muscle development. Methods: RNA-seq was performed to identify differentially expressed genes (DEGs) across pig breeds with distinct growth rates. Weighted gene co-expression network analysis (WGCNA) was used to identify gene modules potentially associated with skeletal muscle growth and development. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis was performed to validate the role of CCND2 in skeletal muscle differentiation. Results We identified 5,416, 4,846, and 2,777 DEGs in the Duroc pig (DU) vs. Jiangquan black pig (JQ), DU vs. Laiwu pigs (LW), and JQ vs. LW comparison groups, respectively. WGCNA identified two gene modules (lightgreen and darkmagenta) potentially associated with skeletal muscle growth and development. CCND2 exhibited differential expression in comparisons between DU vs JQ and DU vs LW, with higher expression observed in Duroc pigs. qRT-PCR analysis revealed that overexpression of CCND2 at day 6 (D6) of cell differentiation significantly increased the expression levels of CCND2 gene and markers of myogenic differentiation (p < 0.05). Conversely, knockdown of CCND2 significantly decreased the expression levels of CCND2 gene and myogenic differentiation markers (p < 0.05). Conclusions This study identified breed-specific DEGs and key gene modules in porcine skeletal muscle and verified that CCND2 promotes myogenic differentiation. These findings provide a genetic theoretical basis for further investigations into skeletal muscle growth and development. Keywords pig, skeletal muscle, transcriptome, WGCNA, CCND2

Environmental exposures can induce epigenetic modifications that persist across generations, potentially contributing to the transmission of environmentally-induced phenotypes. However, the extent to which such persistence of molecular changes are maintained independently of genetic mechanisms remains unclear. In this study, we investigated the evolution and genetic determinism of DNA methylation (DNAm) variation across three successive generations of Japanese quails (G0 to G2) belonging to two epilines, defined according to whether their female ancestor had received genistein supplementation (epi+) or not (epi-). Using reduced representation bisulfite sequencing (RRBS), we characterised CpG methylation patterns and identified differentially methylated cytosines (DMCs) and regions (DMRs) between control and genistein-supplemented epilines. Out of 112,745 CpG sites analysed, no DMCs or DMRs were detected between epilines in the first generation following exposure (G0), whereas 621 DMCs (36 DMRs) were identified in G1, and 1,381 DMCs (101 DMRs) in G2. Despite this progressive increase of differential methylation sites, mean differences in methylation rate between epilines were globally small. Similarly, genome-wide genetic differentiation between epilines was very limited. Overall, SNP-based heritability of DNA methylation was 0.19 indicating that genetic variation accounted for only a small part of the global variation of DNAm levels. By contrast, DNAm at differentially methylated CpG sites was extensively under genetic control, with a mean SNP-based heritability > 0.5, and meQTL analysis further identified significant associations between SNPs and methylation levels. Together, these results indicate that local DNAm differentiation in quail is under substantial genetic regulation, while the limited genetic differentiation between groups suggests that genetic variation alone can not fully explain the progressive accumulation of methylation differences following an ancestral exposure to genistein. These results support DNAm as a plausible molecular candidate mechanism in the multigenerational response to environmental exposure and highlight the complex interplay between genetic and epigenetic regulations in the inheritance of phenotypes.

The goat research community benefits from one of the most comprehensive livestock genome assemblies, ARS1. However, the rich diversity of goat breeds suggests that a single individual genome is not enough to represent the species full genomic content. Third-generation sequencing technologies have recently made it possible to generate very high quality assemblies that should enable to elucidate unique and/or complex phenotype determinisms depending on genetic structural variations. In the SeqOccin project, nine goat breeds were selected to represent French goat genetic diversity: Alpine, Boer, Creole, Fosses, Poitevine, Provencale, Pyreneenne, Rove and Saanen. High-molecular-weight DNA was extracted from blood samples and subjected to long read whole-genome sequencing. There reads were used to generate de novo assemblies and were aligned onto the public reference to call structural variants. PacBio HiFi and parental Illumina reads were used to assemble phased haplotypes of the Alpine breed. For the other breeds, PacBio raw (CLR) or ONT reads were generated. Our work demonstrates the feasibility of assembling genomes with a reasonable amount of data and appropriate software packages, delivering one high-quality assembly and 8 standard goat assemblies that pave the way for pangenome construction and structural variant functional studies.

The sheep research community benefits from a genome assembly that has been refined over the last 15 years. However, the rich diversity of sheep breeds suggests that a single individual genome is not enough to represent the species full genomic content. Third-generation sequencing technologies have recently made it possible to generate high quality assemblies that should enable to elucidate unique and/or complex phenotype determinisms depending on genetic structural variations. In the SeqOccin project, twelve breeds were selected to represent (mostly French) sheep genetic diversity: Basco-Bearnaise, Black Belly, Blanche du Massif Central, Corse, Ile de France, Lacaune, Manech Tete Noire, Noire du Velay, Ouessant, Rambouillet, Vendeen, Soay. High-molecular-weight DNA was extracted from blood samples and subjected to long read whole-genome sequencing. The reads were used to generate de novo assemblies and were aligned onto the public reference to call structural variants. PacBio HiFi and parental Illumina reads were used to assemble phased haplotypes of the Lacaune breed. For the other breeds, PacBio raw (CLR) or ONT reads were generated. Our work demonstrates the feasibility of assembling genomes with a reasonable amount of data and appropriate software packages, delivering one high-quality assembly and 11 standard sheep assemblies that pave the way for pangenome construction and structural variant functional studies.

Translation initiation is a critical stage of protein synthesis in prokaryotes, where ribosomes assemble on mRNA to begin decoding genetic information. The interactions between ribosomal RNA (rRNA) and mRNA significantly influence this process, with features like ribosomal binding site (RBS) affecting ribosomal dynamics. However, direct experimental validation of these types of interaction is yet to be provided. In this study, we developed an optimized small sub-unit sequencing (SSU-seq) protocol for Escherichia coli, enabling the detailed examination of ribosomal complexes and subunits dynamics. Unlike Ribo-seq, which primarily maps fully assembled ribosomes, SSU-seq captures both ribosomal subunits and intact ribosomes, offering a comprehensive view of translation initiation and elongation events. By integrating experimental data with predictive models, our study advances the understanding of translation regulation in prokaryotes and underscores the utility of SSU-seq in decoding the complexities of ribosome-mRNA interactions. Our analysis reveals distinct patterns of small subunit distribution across the entire mRNAs with higher density near positions with stronger rRNA-mRNA hybridization. The results may suggest that rRNA-mRNA hybridization positions along the entire transcript facilitate efficient start codon recognition by maintaining proximity between the small subunit and the mRNA.

The teleost whole-genome duplication (TGD) contributed to functional diversification of opsins. Some TGD paralogs, including those of parapinopsin (PP), Vertebrate Ancient (VA) opsin, and long wavelength-sensitive (LWS) opsin, show different absorption spectra and/or expression patterns. However, our knowledge of detailed evolutionary processes and mechanisms by which TGD contributed to opsin diversification is still limited. Here, we report that TGD paralogs of parietopsin (PT1 and PT2) are retained in extant species. While most species, including the zebrafish, have only PT1, several characins, including the red piranha (Pygocentrus nattereri), retain both PT1 and PT2. The Mexican tetra (Astyanax mexicanus) and catfishes have only PT2. To assess the degree of functional diversification between PT1 and PT2, we characterized spectral properties and expression patterns. Maximum absorption spectra differ slightly among PTs. Those of red piranha PT1, PT2, Mexican tetra PT2, and Japanese catfish (Silurus asotus) PT2 were located at 517 nm, 528 nm, 517 nm, and ~535 nm, respectively. Fluorescence in situ hybridization showed that (1) piranha PT1 and PT2 are expressed in the same pineal cells, and (2) they are also co-expressed with PP1. Based on these findings, functional diversification of opsins after gene duplication is discussed.

Detecting somatic mutations is a cornerstone of cancer genomics and clinical genotyping, traditionally relying on whole exome (WES) or whole genome sequencing (WGS) data. While RNA sequencing (RNA-seq) is primarily used to assess gene expression, it can also be utilized for somatic variant detection. RNA-seq inherently captures expressed variants, thereby reducing the identification of passenger mutations and mitigating the annotation biases associated with WES. We performed somatic variant calling and cancer driver detection from tumor RNA-seq data from 1,359 RNA-seq samples with matched WGS from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). Our approach using RNA-seq data, alone, identified at least one putative driver gene in 15 of 16 cancer types, consistent with driver predictions from WGS. We further observed that the variant allele frequency (VAF) of RNA-seq-derived variants reflects their predicted functional impact, providing information specific to transcriptomic data. Because RNA-seq data can be generated and analyzed at lower computational cost and time than WGS or WES, it represents a resource-efficient alternative for cancer driver detection. This study provides an extensive pan-cancer resource of RNA-seq-derived variants from the PCAWG consortium and supports the utility of RNA-seq as a complementary approach for characterizing expressed somatic variation in cancer.

Chromatin profiling methods, such as ChIP-seq (chromatin immunoprecipitation followed by sequencing), are used to characterize the genomic localization of DNA-associated proteins. While ChIP and ChIP-seq have been used for decades, an orthogonal approach, CUT&Tag (Cleavage Under Targets and Tagmentation), is gaining popularity as an efficient and cost-effective alternative. Although previous comparative studies note discrepancies in signal-to-noise ratios and detection bias at certain genomic regions, many differences between ChIP-seq and CUT&Tag results remain largely unreconciled. Here, we systematically investigate the disagreeing signals captured by these two methods across well annotated genomic regions. We assess multiple histone mark profiles generated by different groups in two cell lines (K562 and MCF-7). Overall, our analysis indicates that compared to ChIP-seq, CUT&Tag may have limited sensitivity in low-GC environments and, as previously observed, increased signal in hyper-accessible chromatin. Notably, within GC-poor regions of active gene bodies and Polycomb-repressed domains, CUT&Tag exhibits a loss of H3K36me3 and H3K27me3 signal, respectively, where occupancy of these histone marks is otherwise expected. Further, active promoters with discrepant H3K4me3 and H3K27ac signal between the two assays differ systematically in GC content and chromatin accessibility. Promoters differentially enriched for CUT&Tag signal relative to ChIP-seq generally show higher, broader GC-content profiles and higher DNase-seq and ATAC-seq signal, while promoters enriched for ChIP-seq signal harbor the opposite characteristics. A local bias for high GC content and/or chromatin accessibility in CUT&Tag may also explain its differing signal patterns at nucleosome-depleted regions (NDRs) compared to ChIP-seq and MNase-seq promoter profiles. Altogether, our results highlight that studies focused on profiling the intensity and structure of histone modification occupancy can be sensitive to potential biases of CUT&Tag to GC content and chromatin accessibility. These characteristics should be accounted for when selecting a chromatin profiling approach, analyzing and interpreting data as well as drawing biological conclusions.

Introduction. Histoplasma capsulatum is an important fungal pathogen of humans, especially in immunocompromised patients, and the World Health Organisation has listed Histoplasma as a High Priority fungal pathogen. The taxonomy of the genus Histoplasma is still under development, but its only species H. capsulatum has at least of seven lineages, and a range of sublineages. Gap Statement. Allelic multilocus sequence typing (MLST) has been very successfully applied for typing and epidemiology of pathogenic bacteria, but its adoption for fungal pathogens is more limited. Genome sequence-based MLST approaches promise much a higher resolution for comparative phylogenetics, but there are only a small number of genome assemblies available for H. capsulatum. Aim. To develop and validate genome-based MLST typing for H. capsulatum. Methodology. A total of 406 genome assemblies were generated from publicly available Illumina short-read sequencing samples. Selected genome assemblies were subsequently used to generate MLST schemes with chewBBACA, based on whole genome sequences (wgMLST), NCBI-obtained annotations (annMLST) and annotations obtained with Funannotate (fanMLST), and used for lineage identification, and comparative analyses of population structure. Results. Genome assemblies were successfully generated for 406 H. capsulatum samples, with BUSCO completeness scores of over 98%. A subset of 60 genome assemblies was used to generate and validate MLST schemes based on the 60 genome assemblies, 5 genome annotations obtained from Genbank, and 13 genome annotations generated using Funannotate, resulting in schemes with > 8,000 markers. Each of the MLST schemes allowed reproducible separation of all 7 lineages and 6 H. capsulatum Suramericanum sublineages, without the need for high-performance clusters or powerful computing equipment. Comparison of the different schemes using tanglegrams and scatter plots showed good reproducibility (r2 values of 0.8-1.0). Comparison of total sizes of genome assemblies showed significant lineage-dependent size differences, with the Suramericanum, LAmB, India, Capsulatum and Mississippiense lineages averaging a total genome size of 30-33 million bp, and the Africa and Ohiense lineages averaging at 39 million bp. Conclusion. Genome sequence- and especially genome annotation-based MLST analyses are a promising approach that can be adapted as a rapid, high-throughput method for genomic clustering and epidemiological tracing of pathogenic fungi.

Chromosome replication initiation in bacterial genomes typically occurs at a single origin, called oriC, through the binding of the DnaA initiator protein to specific chromosomal repeats known as DnaA-boxes. Here, relying on comparative genomics, we show that the family Fusobacteriaceae, which belongs to the order Fusobacteriales alongside the family Leptotrichiaceae, deviates from this paradigm. Our analysis revealed a stark contrast between the two lineages: while dnaA flanks oriC in Leptotrichiaceae, it is notably missing from Fusobacteriaceae genomes. Adjacent to oriC in Fusobacteriaceae, we identified a gene encoding replication initiator protein A (RPA). We characterized RPA as a Rep-type initiator involved in plasmid replication that likely originated from an iteron-containing plasmid, iterons being the plasmid analogues of DnaA-boxes. Phylogenetic reconstruction indicates that the gene encoding RPA was acquired together with its associated iterons and subsequently domesticated by the Fusobacteriaceae ancestor, successfully displacing the ancestral chromosomal dnaA gene at the origin of replication. This shift to a plasmid-derived replication system may represent an adaptive trait of the Fusobacteriaceae, potentially offering new insights into the biology of Fusobacterium nucleatum, a member of this family known to be a major driver of colorectal cancer progression within the tumor microenvironment.

Genetic loci influencing multiple complex traits often remain difficult to interpret because the underlying functional variants, target genes and tissue contexts are unresolved. Here, we investigated the regulatory basis of a shared signal for average daily gain (ADG) and backfat thickness (BF) in Yorkshire pigs. By integrating GWAS in 7,972 pigs, targeted capture sequencing, LD-based fine-mapping, porcine epigenomic annotations, three-dimensional chromatin interaction maps and functional assays, we narrowed the shared SSC1 signal to a 1.59-Mb interval and functionally prioritized rs334769214 (T>G) as the leading candidate variant. The T>G substitution produced marked allele-specific differences in enhancer activity and nuclear-protein binding, while transcription-factor co-transfection assays identified FOXL1-responsive, allele-dependent regulatory activity. Chromatin interaction and CRISPRa analyses supported a regulatory connection between the rs334769214-containing enhancer and PMAIP1, and PMAIP1 overexpression promoted apoptosis in porcine intestinal epithelial cells. In the capture-sequenced subset, the G allele was associated with higher ADG and BF, consistent with a relative loss-of-activation model. Together, these findings delineate a candidate FOXL1-enhancer-PMAIP1 regulatory pathway and demonstrate how multi-omics fine-mapping can advance a shared quantitative-trait locus toward a testable molecular mechanism underlying variation in growth and fatness.

BackgroundTuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), remains a leading global public health challenge. Malnutrition is a major, modifiable driver of this vulnerability to the disease, creating a lethal, bidirectional link with TB. This s review aims to synthesize available evidence to estimate the pooled prevalence of TB and identify associated factors among malnourished children in Ethiopia, providing critical data to guide targeted, integrated public health interventions. MethodsA systematic review and meta-analysis were conducted on studies reporting the prevalence of TB and associated factors among malnourished children in Ethiopia. Electronic databases including PubMed/Medline, EMBASE, CINAHL, ScienceDirect, Academic Search Complete, and Google Scholar were searched. In addition, bibliographies of relevant articles were reviewed, and a search for relevant conference abstracts was performed. Eligible studies were those published in peer-reviewed journals that specifically addressed TB prevalence among malnourished children. A random-effects model was employed to estimate the pooled prevalence of TB and its associated factors. ResultsA total of 746 studies were initially identified through database and catalog searches. Following screening and eligibility assessment, 12 studies were included in the final systematic review and meta-analysis. These studies, published between 2017 and 2024, collectively involved 5,749 participants. The smallest sample size (133) was reported in a study from the Oromia region, while the largest (1,201) was from the Amhara region. The pooled prevalence of TB among malnourished children was estimated at 16%. The highest prevalence (95%) was reported in a study conducted in Addis Ababa, while the lowest (6%) was from Harar. Considerable variation in prevalence was noted across the included studies. ConclusionThis review resulted in the substantial prevalence of TB among malnourished children in Ethiopia, with a pooled prevalence of 16%. This finding highlights the urgent need for integrated healthcare interventions targeting both TB and malnutrition

Background: Sex differences in asthma are well recognized, but sex-associated airway DNA methylation (DNAm) patterns and their relationships with gene expression remain incompletely characterized. We therefore conducted cis- and trans-expression quantitative trait methylation (eQTM) analyses using sex-stratified and CpG-by-sex interaction models. Methods: We analyzed nasal airway epithelial DNAm and transcriptomic data from 298 participants in the CAAPA discovery cohort. Asthma-associated differentially methylated CpGs (DMCs) were identified using sex-combined and sex-stratified models. Cis- and trans-eQTM analyses were then performed in participants with asthma using sex-combined, sex-stratified, and CpG-by-sex interaction models. To prioritize asthma-associated regulatory relationships, eQTMs that were also significant in non-asthmatic controls were excluded from downstream analyses. Replication in an independent bronchial epithelial cohort was defined by concordant effect direction and P < 0.05. Replicated targets underwent pathway and network analyses. Findings: Asthma-associated DMCs were predominantly hypomethylated in females and hypermethylated in males. Among trans-eQTMs, 379 associations were shared across all four models, whereas the majority were model-specific, indicating extensive sex-dependent regulatory architecture. We replicated 3, 1, 9, and 2 cis-eQTMs and 96, 84, 216, and 41 trans-eQTMs in the sex-combined, female-stratified, male-stratified, and interaction models, respectively. Sex-combined targets implicated innate immune and inflammatory signaling (TNF, IRF8); female-stratified targets implicated cytoskeletal remodeling, focal adhesion, and PI3K-Akt signaling (ITGA10); male-stratified targets implicated transcriptional regulation and post-translational modification (MOV10, HDAC7, CSF1, KAT2A); and interaction targets implicated lipid-mediator metabolism and inflammatory signaling and resolution (CYP2C9, LTB4R, CCL19). Our sex-dependent analyses showed that the majority of eQTMs resided within the trans-regulatory architecture, suggesting that sex-dependent gene expression differences are associated predominantly with trans rather than cis regulatory mechanisms. Pathway analyses revealed enrichment of innate immunity in the sex-combined model, structural remodeling in females, transcriptional control in males, and inflammatory resolution in the sex-interaction model. Interpretation: Sex-specific and sex-dependent airway methylation-expression associations were largely undetected by sex-combined analysis. Sex-aware analytical models may help resolve epigenetic mechanisms underlying sex differences in asthma biology.

Objective: Current guidelines recommend different adjuvant treatment strategies for women with endometrioid early stage endometrial cancer exhibiting a p53 abnormal pattern. We aimed to evaluate recurrence patterns and explore the clinical implications for adjuvant therapy in this population. Materials and methods: This multicenter retrospective study included 86 patients with endometrioid, p53 abnormal, FIGO stage I, II endometrial cancer treated between 2013 and 2023 at six hospitals in the Republic of Korea. Clinicopathologic data were collected through medical record review. Progression free survival (PFS) was estimated using the Kaplan Meier method, and prognostic factors were analyzed using Cox proportional hazards regression models. Results: The 12, 24, and 36 month PFS rates were 91%, 78%, and 76%, respectively. Among the 16 patients who experienced recurrence, 12 (75%) had distant metastases. Adjuvant chemotherapy showed a trend toward improved PFS compared with adjuvant radiation therapy with or without chemotherapy (p=0.10). In multivariable analysis, the p53 overexpression pattern was independently associated with worse PFS compared with the p53 null pattern (HR 8.59, 95% CI 1.61 to 45.79, p=0.01). Conclusion: Given the poor PFS, the predominance of distant recurrence and the observed trend toward improved PFS with chemotherapy, adjuvant chemotherapy may be considered in women with endometrioid early stage endometrial cancer exhibiting a p53 abnormal pattern. Further studies are warranted to clarify the prognostic impact of p53 abnormality subtypes.

The sea urchin is a long-studied model used to study fertilization and various developmental processes. While its genome and transcriptomics data are well documented, the global quantitative proteomic dynamics that complement mRNA expression during development are not available. Here we utilize an "in-cell proteomics" strategy to profile the proteomic changes in 5 developmental time points (egg, 32-cell, blastula, gastrula and larva) of the purple sea urchin, Strongylocentrotus purpuratus. This strategy bypasses cell lysis needed in conventional proteomic experiments and detects over 6,000 protein groups using only 50 eggs/embryos. Over 3,000 proteins have significant changes throughout development that cluster into 6 unique expression patterns. Analysis of signaling proteins and some transcription factors involved in cell specification and differentiation align well with published literature, supporting that proteomes obtained reflect biology of the organism. Furthermore, a majority of mRNA and their corresponding protein levels are discordantly correlated across development. Overall, this study provides a comprehensive analysis of proteomic changes occurring throughout sea urchin development and establishes a powerful and sensitive high-throughput approach suitable for detecting proteomic changes with a small number of embryos.

Drug targets and vaccines for bacterial diseases are challenging to develop. Computational subtractive genomics and reverse vaccinology offer a viable substitute for conventional drug and vaccine development approaches. In recent years, Pasteurella multocida, a Gram-negative pathogen responsible for severe infections in humans along with major diseases in livestock and poultry, has developed resistance to multiple drugs. As a result, there is an urgent demand for research into reliable, efficient, highly potent, and cost-effective strategies to combat the rise of Pasteurella multocida infections. To address this, the present study integrated an in silico subtractive genomics and reverse vaccinology workflow to pinpoint proteins with potential value as both drug and vaccine targets. From the complete Pasteurella multocida Pm70 proteome, we identified 446 proteins that had no similarities with human proteins, out of which 129 were found to be essential for the pathogens survival. These 129 proteins are therefore potential candidates for drug development. Further metabolic pathway analysis, virulence potential, subcellular localization, and structural analyses revealed two outer membrane proteins, namely, 47 kDa outer membrane protein and lipopolysaccharide-assembly protein LptD, as promising candidates for therapeutic development. Subsequent antigenicity analysis revealed both outer membrane proteins as potential vaccine targets, exhibiting highly favorable antigenicity scores. Candidate B- and T-cell epitopes were systematically mapped from these two targets based on antigenicity, immunogenicity, toxicity, and allergenicity potential and only high-confidence epitopes were assembled using adjuvants and class-specific linkers to construct the final multi-epitope vaccine (MEV). The designed MEV construct was demonstrated to be highly antigenic, non-toxic and non-allergenic to humans, while possessing favorable structural and physicochemical properties, and exhibited strong binding interactions with human toll-like receptors with minimal molecular deformability. Although further experimental validation is recommended to fully assess their efficacy, the present study establishes a solid foundation for the development of novel therapeutics and vaccine candidates against P. multocida, ultimately contributing to the reduction of its rising infections.

Uropathogenic Escherichia coli (UPEC) is the predominant causative agent of urinary tract infections. The UPEC strain, CFT073, has a mosaic genome containing 60 genomic islands enriched in virulence genes. In this study, we explore the small proteome of CFT073 using ribosome profiling (Ribo-seq) coupled with RNA sequencing (RNA-seq), and reveal 338 translated, previously unannotated, putative small open reading frames (sORFs) with less than 100 codons. Among them, 78 sORFs are located in intergenic regions. Fifteen intergenic sORFs encode a predicted transmembrane helix (smemORFs), of which seven are hypothetical or have not been characterized in other E. coli strains. The majority of the remaining putative sORFs encodes uncharacterized or novel soluble proteins (ssolORFs). Considering the evolution of many sORFs is likely to be recent, functional sORFs are more likely to have conserved gene neighborhood in other genomes. Searching into the genomes of 10 other E. coli strains and four enterobacteria species, synteny analysis reveals that 10 of 15 putative smemORFs, and 27 of 63 putative ssolORFs, have a homologous sequence in a similar genomic context in at least one other genome, supporting their likelihood of encoding functional proteins. Further validation of putative sORFs using proteomics, cell-free protein synthesis, and mNeonGreen tagging provide strong evidence for the translation of nine predicted intergenic sORFs. Notably, four previously uncharacterized small membrane proteins have increased abundance under oxidative stress, suggesting a potential role in the stress response. Together, our findings reveal the small proteome of an extraintestinal pathogenic E. coli and provide a catalog of putative sORFs for further investigation.

Ovarian cancer stands out as one of the most lethal gynecological cancers among women in the United States. While patients with ovarian cancer often show positive responses to treatment, it is common for the cancer to reappear and develop resistance to treatment, so it is important to explore and create novel combinations to increase patient survival rates. ATRi (Ataxia telangiectasia and Rad3-related inhibitors) and ONC201 monotherapies have shown efficacy in ovarian cancer in preclinical models. We hypothesized that combining ATRi with ONC201 would enhance this efficacy. ONC201 induces apoptosis through the TRAIL-pathway following activation of the integrated stress response and inactivation of Akt which is upregulated in ovarian cancer. Six human ovarian cancer cell lines were treated alone and with the novel drug combination of ceralasertib (ATRi) and ONC201. We found an IC50 range of 0.88-20.16 microM for ONC201 and IC50 range of 0.55-54.75 microM for ceralasertib. The combination of ONC201 plus ceralasertib resulted in PARP cleavage consistent with apoptosis along with reduction in Bcl-XL, Xiap and p-Akt expression after 72 hours of treatment. Cytokine profiling showed a decrease in markers of inflammation and tumor progression. Our data provide evidence of synergy from the combination of ATRi ceralasertib and ONC201 against Ovarian cancer cells that merits further investigation including clinical translation.

Prospero-related homeobox (Prox) genes encode a family of transcription factors that play essential roles in the development of several organs and systems, including the central nervous system, lymphatic endothelium, musculature, and liver. Despite their developmental importance, the evolutionary history of the vertebrate Prox gene family remains poorly understood. In this study we combined phylogenetic and synteny analysis to characterise the evolution of the Prox family in vertebrates. Our results reveal that two to three Prox subfamilies were already present in the last common ancestor of jawed vertebrates. We clarify the identity and evolutionary relationships of well-studied members of this family and identify multiple independent losses of Prox2 and Prox3 genes in specific vertebrate lineages. Furthermore, we uncover evidence for the existence of a fourth Prox gene in the ancestral vertebrate genome, which was subsequently lost. Overall, this study provides the first comprehensive analysis of the evolutionary history of the vertebrate Prox gene family and establishes a foundations for future studies on the functional roles of these genes.