Incorporating prior information in gene expression network-based cancer heterogeneity analysis.

Cancer is molecularly heterogeneous, with seemingly similar patients having different molecular landscapes and accordingly different clinical behaviors. In recent studies, gene expression networks have been shown as more effective/informative for cancer heterogeneity analysis than some simpler measures. Gene interconnections can be classified as "direct" and "indirect," where the latter can be caused by shared genomic regulators (such as transcription factors, microRNAs, and other regulatory molecules) and other mechanisms. It has been suggested that incorporating the regulators of gene expressions in network analysis and focusing on the direct interconnections can lead to a deeper understanding of the more essential gene interconnections. Such analysis can be seriously challenged by the large number of parameters (jointly caused by network analysis, incorporation of regulators, and heterogeneity) and often weak signals. To effectively tackle this problem, we propose incorporating prior information contained in the published literature. A key challenge is that such prior information can be partial or even wrong. We develop a two-step procedure that can flexibly accommodate different levels of prior information quality. Simulation demonstrates the effectiveness of the proposed approach and its superiority over relevant competitors. In the analysis of a breast cancer dataset, findings different from the alternatives are made, and the identified sample subgroups have important clinical differences.

Controlled Noise: Evidence of epigenetic regulation of Single-Cell expression variability.

MOTIVATION: Understanding single-cell expression variability (scEV) or gene expression noise among cells of the same type and state is crucial for delineating population-level cellular function. While epigenetic mechanisms are widely implicated in gene expression regulation, a definitive link between chromatin accessibility and scEV remains elusive. Recent advances in single-cell techniques enable the study of single-cell multiomics data that include the simultaneous measurement of scATAC-seq and scRNA-seq within individual cells, presenting an unprecedented opportunity to address this gap. RESULTS: This paper introduces an innovative testing pipeline to investigate the association between chromatin accessibility and scEV. With single-cell multiomics data of scATAC-seq and scRNA-seq, the pipeline hinges on comparing the prediction performance of scATAC-seq data on gene expression levels between highly variable genes (HVGs) and non-highly variable genes (non-HVGs). Applying this pipeline to paired scATAC-seq and scRNA-seq data from human hematopoietic stem and progenitor cells, we observed a significantly superior prediction performance of scATAC-seq data for HVGs compared to non-HVGs. Notably, there was substantial overlap between well-predicted genes and HVGs. The gene pathways enriched from well-predicted genes are highly pertinent to cell type-specific functions. Our findings support the notion that scEV largely stems from cell-to-cell variability in chromatin accessibility, providing compelling evidence for the epigenetic regulation of scEV and offering promising avenues for investigating gene regulation mechanisms at the single-cell level. AVAILABILITY: The source code and data used in this paper can be found at https://github.com/SiweiCui/EpigeneticControlOfSingle-CellExpressionVariability. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of a copper-responsive small molecule inhibitor of uropathogenic Escherichia coli.

Urinary tract infections (UTIs) are a major global health problem and are caused predominantly by uropathogenic Escherichia coli (UPEC). UTIs are a leading cause of prescription antimicrobial use. Incessant increase in antimicrobial resistance in UPEC and other uropathogens poses a serious threat to the current treatment practices. Copper is an effector of nutritional immunity that impedes the growth of pathogens during infection. We hypothesized that copper would augment the toxicity of select small molecules against bacterial pathogens. We conducted a small molecule screening campaign with a library of 51,098 molecules to detect hits that inhibit a UPEC ΔtolC mutant in a copper-dependent manner. A molecule, denoted as E. coli inhibitor or ECIN, was identified as a copper-responsive inhibitor of wild-type UPEC strains. Our gene expression and metal content analysis results demonstrate that ECIN works in concert with copper to exacerbate Cu toxicity in UPEC. ECIN has a broad spectrum of activity against pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus. Subinhibitory levels of ECIN eliminate UPEC biofilm formation. Transcriptome analysis of UPEC treated with ECIN reveals induction of multiple stress response systems. Furthermore, we demonstrate that L-cysteine rescues the growth of UPEC exposed to ECIN. In summary, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC.IMPORTANCEUrinary tract infection (UTI) is a ubiquitous infectious condition affecting millions of people annually. Uropathogenic Escherichia coli (UPEC) is the predominant etiological agent of UTI. However, UTIs are becoming increasingly difficult to resolve with antimicrobials due to increased antimicrobial resistance in UPEC and other uropathogens. Here, we report the identification and characterization of a novel copper-responsive small molecule inhibitor of UPEC. In addition to E. coli, this small molecule also inhibits pathogens of medical and veterinary significance including Acinetobacter baumannii, Pseudomonas aeruginosa, and methicillin-resistant Staphylococcus aureus.

Single-nucleus transcriptomics of epicardial adipose tissue from female pigs reveals effects of exercise training on resident innate and adaptive immune cells.

BACKGROUND: Coronary artery disease (CAD) is a leading cause of death in women. Epicardial adipose tissue (EAT) secretes cytokines to modulate coronary artery function, and the release of fatty acids from EAT serves as a readily available energy source for cardiomyocytes. However, despite having beneficial functions, excessive amounts of EAT can cause the secretion of proinflammatory molecules that increase the instability of atherosclerotic plaques and contribute to CAD progression. Although exercise mitigates CAD, the mechanisms by which exercise impacts EAT are unknown. The Yucatan pig is an excellent translational model for the effects of exercise on cardiac function. Therefore, we sought to determine if chronic aerobic exercise promotes an anti-inflammatory microenvironment in EAT from female Yucatan pigs. METHODS: Sexually mature, female Yucatan pigs (n = 7 total) were assigned to sedentary (Sed, n = 3) or exercise (Ex, n = 4) treatments, and coronary arteries were occluded (O) with an ameroid to mimic CAD or remained non-occluded (N). EAT was collected for bulk (n = 7 total) and single nucleus transcriptomic sequencing (n = 2 total, 1 per exercise treatment). RESULTS: Based on the bulk transcriptomic analysis, exercise upregulated S100 family, G-protein coupled receptor, and CREB signaling in neurons canonical pathways in EAT. The top networks in EAT affected by exercise as measured by bulk RNA sequencing were SRC kinase family, fibroblast growth factor receptor, Jak-Stat, and vascular endothelial growth factor. Single nucleus transcriptomic analysis revealed that exercise increased the interaction between immune, endothelial, and mesenchymal cells in the insulin-like growth factor pathway and between endothelial and other cell types in the platelet endothelial cell adhesion molecule 1 pathway. Sub-clustering revealed nine cell types in EAT, with fibroblast and macrophage populations predominant in O-Ex EAT and T cell populations predominant in N-Ex EAT. Unlike the findings for exercise alone as a treatment, there were not increased interactions between endothelial and mesenchymal cells in O-Ex EAT. Coronary artery occlusion impacted the most genes in T cells and endothelial cells. Genes related to fatty acid metabolism were the most highly upregulated in non-immune cells from O-Ex EAT. Sub-clustering of endothelial cells revealed that N-Ex EAT separated from other treatments. CONCLUSIONS: According to bulk transcriptomics, exercise upregulated pathways and networks related to growth factors and immune cell communication. Based on single nucleus transcriptomics, aerobic exercise increased cell-to-cell interaction amongst immune, mesenchymal, and endothelial cells in female EAT. Yet, exercise was minimally effective at reversing alterations in gene expression in endothelial and mesenchymal cells in EAT surrounding occluded arteries. These findings lay the foundation for future work focused on the impact of exercise on cell types in EAT.

Novel Nuclear Roles for Testis-Specific ACTL7A and ACTL7B Supported by In Vivo Characterizations and AI Facilitated In Silico Mechanistic Modeling with Implications for Epigenetic Regulation in Spermiogenesis.

A mechanistic role for nuclear function of testis-specific actin related proteins (ARPs) is proposed here through contributions of ARP subunit swapping in canonical chromatin regulatory complexes. This is significant to our understanding of both mechanisms controlling regulation of spermiogenesis, and the expanding functional roles of the ARPs in cell biology. Among these roles, actins and ARPs are pivotal not only in cytoskeletal regulation, but also in intranuclear chromatin organization, influencing gene regulation and nucleosome remodeling. This study focuses on two testis-specific ARPs, ACTL7A and ACTL7B, exploring their intranuclear activities and broader implications utilizing combined in vivo, in vitro, and in silico approaches. ACTL7A and ACTL7B, previously associated with structural roles, are hypothesized here to serve in chromatin regulation during germline development. This study confirms the intranuclear presence of ACTL7B in spermatocytes and round spermatids, revealing a potential role in intranuclear processes, and identifies a putative nuclear localization sequence conserved across mammalian ACTL7B, indicating a potentially unique mode of nuclear transport which differs from conventional actin. Ablation of ACTL7B leads to varied transcriptional changes reported here. Additionally, in the absence of ACTL7A or ACTL7B there is a loss of intranuclear localization of HDAC1 and HDAC3, which are known regulators of epigenetic associated acetylation changes that in turn regulate gene expression. Thus, these HDACs are implicated as contributors to the aberrant gene expression observed in the KO mouse testis transcriptomic analysis. Furthermore, this study employed and confirmed the accuracy of in silico models to predict ARP interactions with Helicase-SANT-associated (HSA) domains, uncovering putative roles for testis-specific ARPs in nucleosome remodeling complexes. In these models, ACTL7A and ACTL7B were found capable of binding to INO80 and SWI/SNF nucleosome remodeler family members in a manner akin to nuclear actin and ACTL6A. These models thus implicate germline-specific ARP subunit swapping within chromatin regulatory complexes as a potential regulatory mechanism for chromatin and associated molecular machinery adaptations in nuclear reorganizations required during spermiogenesis. These results hold implications for male fertility and epigenetic programing in the male-germline that warrant significant future investigation. In summary, this study reveals that ACTL7A and ACTL7B play intranuclear gene regulation roles in male gametogenesis, adding to the multifaceted roles identified also spanning structural, acrosomal, and flagellar stability. ACTL7A and ACTL7B unique nuclear transport, impact on HDAC nuclear associations, impact on transcriptional processes, and proposed mechanism for involvement in nucleosome remodeling complexes supported by AI facilitated in silico modeling contribute to a more comprehensive understanding of the indispensable functions of ARPs broadly in cell biology, and specifically in male fertility.

Experimental and analytical pipeline for sub-genomic RNA landscape of coronavirus by Nanopore sequencer.

Coronaviruses (CoVs), including severe acute respiratory syndrome coronavirus 2, can infect a variety of mammalian and avian hosts with significant medical and economic consequences. During the life cycle of CoV, a coordinated series of subgenomic RNAs, including canonical subgenomic messenger RNA and non-canonical defective viral genomes (DVGs), are generated with different biological implications. Studies that adopted the Nanopore sequencer (ONT) to investigate the landscape and dynamics of viral RNA subgenomic transcriptomes applied arbitrary bioinformatics parameters without justification or experimental validation. The current study used bovine coronavirus (BCoV), which can be performed under biosafety level 2 for library construction and experimental validation using traditional colony polymerase chain reaction and Sanger sequencing. Four different ONT protocols, including RNA direct and cDNA direct sequencing with or without exonuclease treatment, were used to generate RNA transcriptomic libraries from BCoV-infected cell lysates. Through rigorously examining the k-mer, gap size, segment size, and bin size, the optimal cutoffs for the bioinformatic pipeline were determined to remove the sequence noise while keeping the informative DVG reads. The sensitivity and specificity of identifying DVG reads using the proposed pipeline can reach 82.6% and 99.6% under the k-mer size cutoff of 15. Exonuclease treatment reduced the abundance of RNA transcripts; however, it was not necessary for future library preparation. Additional recovery of clipped BCoV nucleotide sequences with experimental validation expands the landscape of the CoV discontinuous RNA transcriptome, whose biological function requires future investigation. The results of this study provide the benchmarks for library construction and bioinformatic parameters for studying the discontinuous CoV RNA transcriptome.IMPORTANCEFunctional defective viral genomic RNA, containing all the cis-acting elements required for translation or replication, may play different roles in triggering cell innate immune signaling, interfering with the canonical subgenomic messenger RNA transcription/translation or assisting in establishing persistence infection. This study does not only provide benchmarks for library construction and bioinformatic parameters for studying the discontinuous coronavirus RNA transcriptome but also reveals the complexity of the bovine coronavirus transcriptome, whose functional assays will be critical in future studies.

Ephrin Forward Signaling Controls Interspecies Cell Competition in Pluripotent Stem Cells.

In the animal kingdom, evolutionarily conserved mechanisms known as cell competition eliminate unfit cells during development. Interestingly, cell competition also leads to apoptosis of donor cells upon direct contact with host cells from a different species during interspecies chimera formation. The mechanisms underlying how host animal cells recognize and transmit cell death signals to adjacent xenogeneic human cells remain incompletely understood. In this study, we developed an interspecies cell contact reporter system to dissect the mechanisms underlying competitive interactions between mouse and human pluripotent stem cells (PSCs). Through single-cell RNA-seq analyses, we discovered that Ephrin A ligands in mouse cells play a crucial role in signaling cell death to adjacent human cells that express EPHA receptors during interspecies PSC co-culture. We also demonstrated that blocking the Ephrin A-EPHA receptor interaction pharmacologically, and inhibiting Ephrin forward signaling genetically in the mouse cells, enhances the survival of human PSCs and promotes chimera formation both in vitro and in vivo . Our findings elucidate key mechanisms of interspecies PSC competition during early embryogenesis and open new avenues for generating humanized tissues or organs in animals, potentially revolutionizing regenerative medicine.

Intranasally Administered EVs from hiPSC-derived NSCs Alter the Transcriptomic Profile of Activated Microglia and Conserve Brain Function in an Alzheimer's Model.

Antiinflammatory extracellular vesicles (EVs) derived from human induced pluripotent stem cell (hiPSC)-derived neural stem cells (NSCs) hold promise as a disease-modifying biologic for Alzheimer's disease (AD). This study directly addressed this issue by examining the effects of intranasal administrations of hiPSC-NSC-EVs to 3-month-old 5xFAD mice. The EVs were internalized by all microglia, which led to reduced expression of multiple genes associated with disease-associated microglia, inflammasome, and interferon-1 signaling. Furthermore, the effects of hiPSC-NSC-EVs persisted for two months post-treatment in the hippocampus, evident from reduced microglial clusters, inflammasome complexes, and expression of proteins and/or genes linked to the activation of inflammasomes, p38/mitogen-activated protein kinase, and interferon-1 signaling. The amyloid-beta (Aß) plaques, Aß-42, and phosphorylated-tau concentrations were also diminished, leading to better cognitive and mood function in 5xFAD mice. Thus, early intervention with hiPSC-NSC-EVs in AD may help maintain better brain function by restraining the progression of adverse neuroinflammatory signaling cascades.

Parallel degradome-seq and DMS-MaPseq substantially revise the miRNA biogenesis atlas in Arabidopsis.

MicroRNAs (miRNAs) are produced from highly structured primary transcripts (pri-miRNAs) and regulate numerous biological processes in eukaryotes. Due to the extreme heterogeneity of these structures, the initial processing sites of plant pri-miRNAs and the structural rules that determine their processing have been predicted for many miRNAs but remain elusive for others. Here we used semi-active DCL1 mutants and advanced degradome-sequencing strategies to accurately identify the initial processing sites for 147 of 326 previously annotated Arabidopsis miRNAs and to illustrate their associated pri-miRNA cleavage patterns. Elucidating the in vivo RNA secondary structures of 73 pri-miRNAs revealed that about 95% of them differ from in silico predictions, and that the revised structures offer clearer interpretation of the processing sites and patterns. Finally, DCL1 partners Serrate and HYL1 could synergistically and independently impact processing patterns and in vivo RNA secondary structures of pri-miRNAs. Together, our work sheds light on the precise processing mechanisms of plant pri-miRNAs.

Transcriptomic signature, bioactivity and safety of a non-hepatotoxic analgesic generating AM404 in the midbrain PAG region.

Safe and effective pain management is a critical healthcare and societal need. The potential for acute liver injury from paracetamol (ApAP) overdose; nephrotoxicity and gastrointestinal damage from chronic non-steroidal anti-inflammatory drug (NSAID) use; and opioids' addiction are unresolved challenges. We developed SRP-001, a non-opioid and non-hepatotoxic small molecule that, unlike ApAP, does not produce the hepatotoxic metabolite N-acetyl-p-benzoquinone-imine (NAPQI) and preserves hepatic tight junction integrity at high doses. CD-1 mice exposed to SRP-001 showed no mortality, unlike a 70% mortality observed with increasing equimolar doses of ApAP within 72 h. SRP-001 and ApAP have comparable antinociceptive effects, including the complete Freund's adjuvant-induced inflammatory von Frey model. Both induce analgesia via N-arachidonoylphenolamine (AM404) formation in the midbrain periaqueductal grey (PAG) nociception region, with SRP-001 generating higher amounts of AM404 than ApAP. Single-cell transcriptomics of PAG uncovered that SRP-001 and ApAP also share modulation of pain-related gene expression and cell signaling pathways/networks, including endocannabinoid signaling, genes pertaining to mechanical nociception, and fatty acid amide hydrolase (FAAH). Both regulate the expression of key genes encoding FAAH, 2-arachidonoylglycerol (2-AG), cannabinoid receptor 1 (CNR1), CNR2, transient receptor potential vanilloid type 4 (TRPV4), and voltage-gated Ca2+ channel. Phase 1 trial (NCT05484414) (02/08/2022) demonstrates SRP-001's safety, tolerability, and favorable pharmacokinetics, including a half-life from 4.9 to 9.8 h. Given its non-hepatotoxicity and clinically validated analgesic mechanisms, SRP-001 offers a promising alternative to ApAP, NSAIDs, and opioids for safer pain treatment.

A Unified Bayesian Framework for Bi-overlapping-Clustering Multi-omics Data via Sparse Matrix Factorization.

The advances of modern sequencing techniques have generated an unprecedented amount of multi-omics data which provide great opportunities to quantitatively explore functional genomes from different but complementary perspectives. However, distinct modalities/sequencing technologies generate diverse types of data which greatly complicate statistical modeling because uniquely optimized methods are required for handling each type of data. In this paper, we propose a unified framework for Bayesian nonparametric matrix factorization that infers overlapping bi-clusters for multi-omics data. The proposed method adaptively discretizes different types of observations into common latent states on which cluster structures are built hierarchically. The proposed Bayesian nonparametric method is able to automatically determine the number of clusters. We demonstrate the utility of the proposed method using simulation studies and applications to a single-cell RNA-sequencing dataset, a combination of single-cell RNA-sequencing and single-cell ATAC-sequencing dataset, a bulk RNA-sequencing dataset, and a DNA methylation dataset which reveal several interesting findings that are consistent with biological literature.

Aryl Hydrocarbon Receptor (AhR) Signaling in Colonic Cells and Tumors.

The aryl hydrocarbon receptor (AhR) is overexpressed in many tumor types and exhibits tumor-specific tumor promoter and tumor suppressor-like activity. In colon cancer, most but not all studies suggest that the AhR exhibits tumor suppressor activity which is enhanced by AhR ligands acting as agonists. Our studies investigated the role of the AhR in colon tumorigenesis using wild-type and AhR-knockout mice, the inflammation model of colon tumorigenesis using mice treated with azoxymethane (AOM)/dextran sodium sulfate (DSS) and APCS580/+; KrasG12D/+ mice all of which form intestinal tumors. The effects of tissue-specific AhR loss in the intestine of the tumor-forming mice on colonic stem cells, organoid-initiating capacity, colon tumor formation and mechanisms of AhR-mediated effects were investigated. Loss of AhR enhanced stem cell and tumor growth and in the AOM/DSS model AhR-dependent suppression of FOXM1 and downstream genes was important for AhR-dependent anticancer activity. Furthermore, the effectiveness of interleukin-22 (IL22) in colonic epithelial cells was also dependent on AhR expression. IL22 induced phosphorylation of STAT3, inhibited colonic organoid growth, promoted colonic cell proliferation in vivo and enhanced DNA repair in AOM/DSS-induced tumors. In this mouse model, the AhR suppressed SOCS3 expression and enhanced IL22-mediated activation of STAT3, whereas the loss of the AhR increased levels of SOCS3 which in turn inhibited IL22-induced STAT3 activation. In the APCS580/+; KrasG12D/+ mouse model, the loss of the AhR enhanced Wnt signaling and colon carcinogenesis. Results in both mouse models of colon carcinogenesis were complemented by single cell transcriptomics on colonic intestinal crypts which also showed that AhR deletion promoted expression of FOXM1-regulated genes in multiple colonic cell subtypes. These results support the role of the AhR as a tumor suppressor-like gene in the colon.