PubMed HealthSearch

SEARCH · PubMed Health

Results for “virtual protein profiling”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

DPAS-Graph: adaptive spatial-feature relation learning for spatial RNA-to-protein prediction and virtual protein profiling.

Paired spatial multi-omics provides a supervised basis for learning RNA-protein correspondence in situ, but predicting protein abundance from spatial transcriptomic data alone remains challenging across tissue contexts and protein panels. Here, we present DPAS-Graph, an adaptive relation-learning framework for spatial RNA-to-protein prediction. Rather than directly merging spatial proximity and transcriptomic similarity as fixed graph priors, DPAS-Graph represents them as two relation channels on a shared edge support and updates their contributions during representation learning for protein prediction. Its Niche-Coupled Field Encoder combines layer-wise edge-relation modeling, intra-branch relation refinement, and cross-branch residual correction to learn spot representations for protein abundance prediction. In a leave-one-dataset-out benchmark across seven paired spatial multi-omics datasets, DPAS-Graph achieved lower aggregate prediction errors and improved spot-level agreement of protein expression profiles, with gains mainly reflected in error-based metrics and PCC-Spot. Spatial autocorrelation and protein-derived domain agreement analyses were further used to characterize the spatial behavior of the predicted protein maps. When applied to external RNA-only spatial sections, DPAS-Graph generated qualitatively interpretable marker-level virtual protein maps, illustrating its use as a complementary tool for protein-level interpretation of transcriptomics-only spatial data.

RNA

Evidence for a structurally homologous Rh-like polypeptide in Rhnull erythrocytes.

Human Rhnull red blood cells fail to react with Rh antibodies, indicating that these cells are either devoid of Rh protein or, like other species, possess antigenically distinct variants. To determine whether Rhnull cells possess an Rh-like polypeptide, 32-kDa proteins from D--, rr, and Rhnull cells were labeled with the cysteine-specific probe, 125I-labeled pyridyldithioethylamine. Size comparisons of labeled proteins in Triton X-100-solubilized membranes from Rh-bearing and Rhnull cells showed similar sedimentation coefficients and Stoke's radii. Immunoprecipitated Rh(D) from D-- cells, Rh(c) from rr cells, and purified 32-kDa proteins from Rhnull cells were digested with alpha-chymotrypsin and examined by high-performance liquid chromatography and by two-dimensional iodopeptide mapping. Analysis of 125I-labeled chymotryptic fragments from immunoprecipitated Rh(D) and Rh(c) showed the labeled peptides from both phenotypes to be virtually identical. High-performance liquid chromatography profiles and iodopeptide maps of 32-kDa Rhnull proteins yielded patterns identical to 32-kDa proteins isolated from D-- cells and rr cells with the exception of one missing 125I-labeled peptide. Further analysis of the Rh-related fragments from Rhnull cells showed significant homology with immunoprecipitated Rh(D) and Rh(c). DNA sequence analysis of cysteine-encoding regions from Rh-bearing and Rhnull cells showed complete identity. These data suggest that Rhnull red blood cells, although serologically distinct, possess an Rh-like protein that is structurally very similar to Rh(D) and Rh(c).

Amino Acid Sequence

Isolation and partial characterization of rat CNS axolemma enriched fractions.

Axolemma-enriched fractions were prepared from rat brain by osmotic shock of a purified preparation of myelinated axons and subsequent separation of myelin, two axolemma-enriched fractions and myelin-free axons by density gradient centrifugation. Compared with the starting whole homogenate, the fractions were enriched in specific activity of Na+K+ ATPase, acetylcholinesterase, 5'nucleotidase as well as 2',3'-cyclic nucleotide 3'phosphohydrolase. Compared with myelin, the axolemmal fractions are greatly enriched in high molecular weight proteins. The 1.0/1.2 fraction has a predominant peak of fucose-labeled glycoprotein with a molecular weight between that of the myelin associated glycoprotein and the Wolfgram protein which is absent from the myelin glycoprotein profile. Polyacrylamide gel electrophoresis showed that the protein profile of myelin isolated by this procedure was similar to that of myelin isolated by other procedures and that the myelin specific basic and proteolipid proteins were virtually absent in the axolemma-enriched fractions. Both axolemma fractions were enriched in higher MW proteins, some of which resembled proteins in the myelin protein profile. Both axolemma-enriched fractions specifically bind between 2 and 3 pmoles of [3H]tetrodotoxin per mg protein. The axolemma-enriched fractions incorporated [3H]leucine and [14C]fucose exclusively into high molecular weight proteins and glycoproteins. In contrast myelin concomitantly isolated with the axolemma-enriched fractions had a significant amount of [3H]leucine labeled protein in myelin proteolipid and basic proteins. In addition to the myelin associated g-ycoprotein the [14C]fucose labeled a glycoprotein of slightly larger apparent molecular weight than proteolipid protein was found in the myelin fraction while the comparable labeled glycoprotein was absent in the axolemma-enriched fractions. The possible extent of contamination of these fractions by myelin or myelin subfractions and relationship of these axolemma-enriched fractions to other axolemma preparations are discussed.

2',3'-Cyclic-Nucleotide Phosphodiesterases

Identifying JAK2 and ANXA5 as Key Genes Linking Obstructive Sleep Apnea and Oxidative Stress via Machine Learning and Multilayer Transcriptomic Integration With Functional Validation.

Obstructive sleep apnea (OSA) is a common and severe sleep disorder closely associated with oxidative stress (OS). This study aims to identify and validate potential OS-related genes associated with OSA through bioinformatics methods. We successfully identified OS-related differentially expressed genes (OS-DEGs) by combining the limma test, weighted correlation network analysis (WGCNA), and OS-related genes from the GeneCards database. Key genes and potential biological roles were further identified using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), enrichment analysis, protein-protein interaction (PPI) network analysis, Lasso regression analysis, random forest algorithm, and support vector machine recursive feature elimination (SVM-RFE) method. Evaluate and validate the accuracy of key genes through receiver operating characteristic (ROC) curve analysis. The human single-cell RNA sequencing (scRNA-seq) dataset is used for cell classification annotation, analysis of key gene single-cell expression profiles, and virtual gene knockout experiments based on the scTenifoldKnk algorithm. Integrating scRNA-seq sequencing, pseudotime trajectory inference, cell-cell communication analysis, and bulk immune infiltration deconvolution reveals monocyte subtype remodeling in OSA. Finally, the expression levels of key genes in clinical samples were validated using real-time quantitative PCR (RT-qPCR) and Western blotting. A total of 57 common DEGs, indicating significant enrichment in OS, inflammation, and tumor pathways, particularly prominent in the immunometabolism pathway. By integrating DEGs, WGCNA, PPI results, and machine learning methods, key genes Janus kinase 2 (JAK2) and ANXA5 were screened out. JAK2 was significantly upregulated under disease conditions, while ANXA5 was significantly downregulated. ROC curve exhibited high accuracy (area under the curve [AUC] > 0.85). Human scRNA-seq analysis revealed that key genes were predominantly highly expressed in monocytes. Virtual knockout experiments demonstrated that these key genes play a crucial role in regulating immune responses and inflammatory reactions. PPI networks and enrichment analysis verified that downstream genes S100P, ALOX5AP, PROK2, and PADI4 may collaboratively participate in immune response and inflammation regulation. Finally, clinical sample experiment further validated the results of bioinformatics analysis. This study provides new research insights for the diagnosis, mechanism research, and treatment development of OSA in the future by integrating multilayer transcriptomic and machine learning techniques.

Humans

Population differentiation under the charge state model.

The extent of divergence between partially isolated sub-populations for electrophoretically detectable alleles was formulated assuming the island model of migration and the charge state model of mutation. At equilibrium the ratio of the variance of charge between the means of k different islands to the average within-island variance of charge was shown to be approximately 4Nemk2/(k-1)2 where Ne is the effective size of each island population and m is the migration rate. This ratio was calculated from published data for eight polymorphic loci in six island populations of Drosophila willistoni. Under the assumption that all variants are selectively neutral, migration rates of greater than 10 adults per generation per island are required to explain the observed similarity of the allelic profiles in D. willistoni. Since the islands studied appear to be virtually completely isolated it was concluded either that the observed protein variants are adaptive and maintained in populations by some form of balancing selection or that the observed variants themselves are neutral but natural selection acts to restrict the appearance of more extreme variants in the charge carried.

Drosophila

Comparative and Subtractive Genomics Analysis of Multidrug-Resistant Klebsiella pneumoniae Strains for Novel Target Identification and Drug Repurposing Strategies.

The rapid rise of multidrug-resistant (MDR) Klebsiella pneumoniae has created a major global health challenge due to the limited availability of conserved therapeutic targets effective across diverse resistant strains. In this study, an integrative computational target-discovery and drug-repurposing framework was applied to six clinically relevant K. pneumoniae strains. Comparative genomic analysis identified 3012 conserved genes, which were subsequently filtered to nine essential, non-host homologous proteins. Among these, three conserved cytoplasmic proteins (accD, cpxR, and mraZ) were prioritized for functional analysis, with acetyl-CoA carboxylase subunit beta (accD) emerging as the most promising therapeutic target based on sequence conservation, predicted essentiality, subcellular localization, and pathway association. Structural assessment supported the reliability of the predicted accD model, whereas consensus binding-site analysis identified key residues suitable for ligand interaction. Virtual screening of FDA-approved drugs followed by molecular docking identified several compounds with favorable binding profiles toward accD. Subsequent molecular dynamics simulations, including root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), hydrogen-bond occupancy, principal component analysis (PCA), and PCA-based free energy landscape (FEL) analyses, consistently identified tenapanor, micafungin, deferoxamine, and cobicistat as the most stable protein-ligand complexes, with tenapanor exhibiting the most favorable overall structural and thermodynamic stability profile. These findings identify accD as a promising therapeutic target in MDR K. pneumoniae and suggest several FDA-approved compounds as potential candidates for drug repurposing. Although experimental validation is needed to confirm their biological activity and therapeutic potential, this study demonstrates the potential of integrating comparative genomics with molecular dynamics analyses to support antimicrobial target identification and drug repurposing against MDR bacterial pathogens.

Klebsiella pneumoniae

The rate of access to the organic ligand-binding region of serum albumin is entropy controlled.

The technique of real time dielectric relaxation measurement coupled with a conventional stopped-flow device has made it possible to measure the rates of association and dissociation of the complex of human serum albumin with its most prevalent ligands, the long-chain fatty acids. This association was previously shown to proceed in two steps: a fast, probably diffusion-controlled, nonspecific association, followed by a slower (approximately 3 sec--1) rearrangement of the intermediate protein--ligand configuration, whose kinetics is first order. By use of the Arrhenius relation and standard theory of rate processes it is determined that there is virtually no activation enthalpy in the forward binding reaction and that the rate of access to the interior hydrophobic binding region of serum albumin is controlled by a negative entropy of activation, reflecting a high degree of ordering in the transition state. A complete thermodynamic and kinetic profile of the association reaction is given.

Binding Sites

A procedure for the quantitative recovery of homogeneous populations of undegraded free and bound polysomes from rat liver.

A procedure is described for the preparation of free and bound polysomes from whole homogenates of rat liver tissue. Liver is homogenized in a conventional medium containing glutathione; then after a 12-min centrifugation at 131000g, the free polysomes in the supernatant are saved, while the membrane-bound polysomes in the pellet are suspended in a mixture of ribonuclease inhibitors (cell sap, 250 mM KCl, and glutathione), homogenized in the presence of detergent (Triton X-100), centirfuged for 5 min at 1470g, decanted, and treated with deoxycholate; the polysomes in the two supernatants are harvested by centrifugation through sucrose gradients containing 250 mM KCl and cell sap. Free and bound polysomes prepared in this manner are undegraded, equally active in cell-free protein synthesis, and virtually free of ribonuclease, membranous material, glycogen, deoxycholate, completed protein, and cross-contamination. The recovery of polysomes is approximately 95% and the distribution between the free and membrane-bound state is 25 and 75%, respectively. The molecular weight profiles after sodium dodecyl sulfate-acrylamide gel electrophoresis of the polypeptides completed and released by free and bound polysomes in vitro are different, indicating that there are quantitative differences in the synthesis of various size polypeptides between the two polysome classes. The differential centrifugation procedure is rapid and reproducible, requires much less ultracentrifugation than the isopycnic technique, and provides a nearly quantitative means of separating free and bound polysomes.

Animals

Activation of lipoprotein lipase by native and acylated peptides of apolipoprotein C-II.

Apolipoprotein C-II, a protein found associated with all major classes of plasma lipoproteins, is a potent activator of the enzyme lipoprotein lipase. We have prepared the maleyl, citraconyl and succinyl derivatives of apolipoprotein C-II, and compared the capacities of the intact and tryptically cleaved proteins to activate lipoprotein lipase. The NH2-terminal 50 residue peptide proved virtually inactive, even after removal of the masking groups from the citraconyl derivative. The COOH-terminal 29 residue peptides of maleyl and citraconyl apolipoprotein C-II were more active than the corresponding succinylated peptide. After deacylation of the citraconyl derivative, the COOH-terminal peptide had maximal activity as great as apolipoprotein C-II, although the profile of activation remained dissimilar at low activator concentrations.

Acylation

Peripheral nerve glycoproteins and myelin fine structure during development of rat sciatic nerve.

Developmental changes in relative amounts of peripheral nerve proteins and glycoproteins have been correlated with the degree of morphological myelination at various ages during the first 25 postnatal days in rat sciatic nerve. At birth there is virtually no major myelin glycoprotein (P0), but there is a protein which migrates to the same point on sodium dodecyl sulphate (SDS) polyacrylamide gels as the small myelin basic protein (P2). During the time myelin is being formed in the nerve, the P0 protein increases and the P2 protein appears to decrease in relative amount in the nerve. The accumulation of P0 protein in the nerve correlates extremely well with the degree of myelination in sciatic nerve. At 4-6 days postnatal, smooth membrane profiles are observed which are located within axons and in the inner Schwann cell cytoplasm. Such profiles are also observed to fuse with the axolemma-Schwann cell interface. The profiles may represent membrane material being added to or deleted from the axolemma or myelin during myelination.

Animals

RP-REP Ribosomal Profiling Reports: an open-source cloud-enabled framework for reproducible ribosomal profiling data processing, analysis, and result reporting.

Ribosomal profiling is an emerging experimental technology to measure protein synthesis by sequencing short mRNA fragments undergoing translation in ribosomes. Applied on the genome wide scale, this is a powerful tool to profile global protein synthesis within cell populations of interest. Such information can be utilized for biomarker discovery and detection of treatment-responsive genes. However, analysis of ribosomal profiling data requires careful preprocessing to reduce the impact of artifacts and dedicated statistical methods for visualizing and modeling the high-dimensional discrete read count data. Here we present Ribosomal Profiling Reports (RP-REP), a new open-source cloud-enabled software that allows users to execute start-to-end gene-level ribosomal profiling and RNA-Seq analysis on a pre-configured Amazon Virtual Machine Image (AMI) hosted on AWS or on the user's own Ubuntu Linux server. The software works with FASTQ files stored locally, on AWS S3, or at the Sequence Read Archive (SRA). RP-REP automatically executes a series of customizable steps including filtering of contaminant RNA, enrichment of true ribosomal footprints, reference alignment and gene translation quantification, gene body coverage, CRAM compression, reference alignment QC, data normalization, multivariate data visualization, identification of differentially translated genes, and generation of heatmaps, co-translated gene clusters, enriched pathways, and other custom visualizations. RP-REP provides functionality to contrast RNA-SEQ and ribosomal profiling results, and calculates translational efficiency per gene. The software outputs a PDF report and publication-ready table and figure files. As a use case, we provide RP-REP results for a dengue virus study that tested cytosol and endoplasmic reticulum cellular fractions of human Huh7 cells pre-infection and at 6 h, 12 h, 24 h, and 40 h post-infection. Case study results, Ubuntu installation scripts, and the most recent RP-REP source code are accessible at GitHub. The cloud-ready AMI is available at AWS (AMI ID: RPREP RSEQREP (Ribosome Profiling and RNA-Seq Reports) v2.1 (ami-00b92f52d763145d3)).

AMI

Decreased proteolysis and increased amino acid efflux in aging human fibroblasts.

The release of radioactivity was measured following variable times of dual labeling of normal human fibroblasts at early and late passage and late-passage cell strains derived from the progeria and Werner's syndrome of premature aging. In the rapid first phase to 20 min, all 3 late-passage (aged) cells released more acid-soluble radioactivity into the medium than early-passage normal cells in the order Werner's greater than progeria greater than late-passage normal greater than early-passage normal, virtually all of this radioactivity emanating from intracellular material that was acid-soluble after prelabeling (zero time). In the slower second phase from 20 min to 3 h, all 4 cell types showed approximately parallel profiles of release. Following preloading with the non-utilizable amino acid alpha-amino[14C] isobutyric acid, progeria and Werner's cells released radioactivity more rapidly than did early- and late-passage normal cells in that order. In contrast, direct measurements on the net loss of counts from acid-insoluble material, i.e. true proteolysis, revealed that all 3 aged cells degraded proteins of short half-life more slowly but proteins of long half-life were degraded at the same rate as those of young normal cells. The results indicate that aged cells have a reduced proteolytic capacity and increased amino acid efflux. The latter process is probably due to the higher proportion of labeled amino acids in the cell water of aged cells at zero time and possibly increased membrane leakiness. The reason for the decreased proteolysis is not clear but it may relate to the rising proportion of defective proteins, the increased protein content and the loss of replicative capacity in aging cells.

Aging

Partially hepatectomized rats: a model for the study of the effect of toxins on the plasma protein profiles of nascent hepatocytes.

A useful framework is proposed for unifying the synthesis of plasma proteins and their degradation by, or release from, liver cells of intact and partially hepatectomized rats, in which synthesis and release of acute-phase plasma proteins occur in synchrony with the internalization and catabolism of plasma and extracellular proteins. The catabolism of proteins and other hepato-intracellular glycoproteins during sepsis or trauma is essential to provide constituent amino acids and carbohydrates for the synthesis of acute-phase plasma proteins. Increases in the plasma levels of acute-phase response proteins in sham-operated rats reached a maximum between 1 and 2 d after mock surgery, and had returned virtually to control levels within 6 d. By contrast, acute-phase proteins in the plasma of partially hepatectomized rats were decreased by 10-20% of their initial values after 24 h. A maximum acute-phase response on d 7 after the operation was characterized by an increase of 181, 445, and 19% for alpha-1-acid glycoprotein, hepatoglobin, and hemopexin, whereas other acute-phase proteins remained below control levels, for example, by 11, 25, and 38% for albumin, transferrin, and prealbumin, respectively. This delayed response suggests that the nascent liver cells had inherited the capacity of the parent cells to respond to inflammatory signal and had synthesized acute-phase plasma proteins. Accordingly, a time frame for the application of toxin to nascent hepatocytes is suggested. An increased activity (300 +/- 50%) for both bound and free neuraminidase in remnant liver tissue 19 h post partial hepatectomy suggested that hepatic regenerating factor(s) were produced in liver tissue via the hepatic bound and/or free neuraminidase-mediated desialylation of humoral substrates. By contrast, circulating levels of lysosomal enzymes alpha-fucosidase and beta-N-acetyl-D-glucosaminidase were increased marginally after 24 h but had returned nearly to control levels after 7 d, suggesting that lysosomal acid hydrolases do not play a major role in regenerative DNA synthesis, mitosis, or in the synthesis of acute-phase plasma proteins.

Acute-Phase Proteins

Food-processing approaches to altering allergenic potential of milk-based formula.

All the major cow milk proteins in their native states are potential allergens in infants with milk allergy. Heat treatment can reduce the antigenicity of whey proteins considerably, but it has virtually no effect on the antigenicity of casein. Infants allergic to milk still react to heat-denatured whey proteins. Therefore heat denaturation alone cannot produce a formula with low allergenicity. Enzyme hydrolysis reduces the antigenicity and allergenicity of protein. Partial hydrolysis produces hydrolysate consisting mainly of large peptides, whereas extensive hydrolysis produces hydrolysate containing a mixture of large and small peptides and free amino acids. Enzyme hydrolysis often produces bitter peptides and destroys the physical functionality of protein. When hydrolysate formula is made, casein or whey protein hydrolysate is ultrafiltered to remove large residual peptides. Certain amino acids are fortified to provide a balanced amino acid profile. The final formulation must comply with the recommendations of the Codex Alimentarius Commission (FAO/WHO) and the U.S. Infant Formula Act to provide adequate infant nutrition. Protein hydrolysate-based formula can be improved by further reducing the residual allergenicity, increasing the small peptide content, and improving the taste.

Food Handling

Resonance Raman studies of Rieske-type proteins.

Resonance Raman (RR) spectra are reported for the [2Fe-2S] Rieske protein from Thermus thermophilus (TRP) and phthalate dioxygenase from Pseudomonas cepacia (PDO) as a function of pH and excitation wavelength. Depolarization ratio measurements are presented for the RR spectra of spinach ferredoxin (SFD), TRP, and PDO at 74 K. By comparison with previously published RR spectra of SFD, we suggest reasonable assignments for the spectra of TRP and PDO. The spectra of PDO exhibit virtually no pH dependence, while significant changes are observed in TRP spectra upon raising the pH from 7.3 to 10.1. One band near 270 cm-1, which consists of components at 266 cm-1 and 274 cm-1, is attributed to Fe(III)-N(His) stretching motions. We suggest that these two components arise from conformers having a protonated-hydrogen-bonded imidazole (266 cm-1) and deprotonated-hydrogen-bonded imidazolate (274 cm-1) coordinated to the Fe/S cluster and that the relative populations of the two species are pH-dependent; a simple structural model is proposed to account for this behavior in the respiratory-type Rieske proteins. In addition, we have identified RR peaks associated with the bridging and terminal sulfur atoms of the Fe-S-N cluster. The RR excitation profiles of peaks associated with these atoms are indistinguishable from each other in TRP (pH 7.3) and PDO and differ greatly from those of [2Fe-2S] ferrodoxins. The profiles are bimodal with maxima near 490 nm and > approx. 550 nm. By contrast, bands associated with the Fe-N stretch show a somewhat different enhancement profile. Upon reduction, RR peaks assigned to Fe-N vibrations are no longer observed, with the resulting spectrum being remarkably similar to that reported for reduced adrenodoxin. This indicates that only modes associated with Fe-S bonds are observed and supports the idea that the reducing electron resides on the iron atom coordinated to the two histidine residues. Taken as a whole, the data are consistent with an St2FeSb2Fe[N(His)]t2 structure for the Rieske-type cluster.

Burkholderia cepacia

Aberrant interleukin (IL)-4 and IL-5 production in vitro by CD4+ helper T cells from atopic subjects.

The cytokine secretion profiles of T cell lines (TCL) specific for purified protein derivative (PPD) or streptokinase (SK), contemporarily derived from nine atopic and nine nonatopic individuals, were compared. Upon stimulation with phorbol myristate acetate (PMA) plus anti-CD3 monoclonal antibody (mAb), all TCL from both atopics and nonatopics produced interleukin (IL)-2 and interferon (IFN)-gamma. The mean IL-2 production by PPD- or SK-specific TCL from both atopics and nonatopics was similar, whereas the mean IFN-gamma production by TCL derived from atopics was significantly lower. In addition, both PPD- and SK-specific TCL from atopics produced detectable amounts of IL-4 and IL-5, whereas the corresponding TCL derived from nonatopics did not. A total number of 107 and 99 PPD-specific CD4+ T cell clones (TCC) were then derived from TCL of 4 atopic and 4 nonatopic donors and assessed for their profile of cytokine production in response to stimulation with either PMA plus anti-CD3 mAb or the specific antigen. Under both these experimental conditions, virtually all PPD-specific TCC from both atopic and nonatopic individuals produced IL-2 and IFN-gamma. In contrast, the great majority of PPD-specific TCC derived from nonatopic individuals did not produce IL-4 and IL-5, whereas high proportions of PPD-specific TCC derived from atopic donors displayed the ability to produce noticeable amounts of IL-4 and IL-5 besides IL-2 and IFN-gamma. These data indicate that CD4+ T cells from atopic individuals are able to produce IL-4 and IL-5 in response to bacterial antigens, such as PPD and SK, that usually evoke responses with a restricted type-1 T helper (Th1)-like cytokine profile in nonatopic individuals. Aberrant IL-4 production by Th cells may represent one of the immune alterations responsible for enhanced IgE antibody production in atopic people.

Adult

Histology-Based Virtual RNA Inference Identifies Pathways Associated With Metastasis Risk in Colorectal Cancer.

Colorectal cancer (CRC) remains a major health concern, with >150,000 new diagnoses and >50,000 deaths annually in the United States, underscoring an urgent need for improved screening, prognostication, disease management, and therapeutic approaches. The tumor microenvironment (TME)-comprising cancerous and immune cells interacting within the tumor's spatial architecture-plays a critical role in disease progression and treatment outcomes, reinforcing its importance as a prognostic marker for metastasis and recurrence risk. However, traditional methods for TME characterization, such as bulk transcriptomics and multiplex protein assays, lack sufficient spatial resolution. Although spatial transcriptomics (ST) allows for the high-resolution mapping of whole transcriptomes at near-cellular resolution, current ST technologies (eg, Visium and Xenium) are limited by high costs, low throughput, and issues with reproducibility, preventing their widespread application in large-scale molecular epidemiology studies. In this study, we refined and implemented virtual RNA inference (VRI) to derive ST-level molecular information directly from hematoxylin and eosin (H&E)-stained tissue images. Our VRI models were trained on the largest matched CRC ST data set to date, comprising 45 patients and >300,000 Visium spots from primary tumors. Using state-of-the-art deep learning models (UNI, ResNet-50, Vision Transformer, and Vision Mamba), we achieved a median Spearman's correlation coefficient of 0.546 between predicted and measured spot-level expression. As validation, VRI-derived gene signatures linked to specific tissue regions (tumor, interface, submucosa, stroma, serosa, muscularis, and inflammation) showed strong concordance with signatures generated via direct ST, and VRI performed accurately in estimating cell-type proportions spatially from H&E slides. In an expanded CRC cohort controlling for tumor invasiveness and clinical factors, we further identified VRI-derived gene signatures significantly associated with key prognostic outcomes, including metastasis status. Although certain tumor-related pathways are not fully captured by histology alone, our findings highlight the ability of VRI to infer a wide range of "histology-associated" biological pathways at near-cellular resolution without requiring ST profiling. Future efforts will extend this framework to expand TME phenotyping from standard H&E tissue images, with the potential to accelerate translational CRC research at scale.

Humans

Multistage Genetic, Transcriptomic, and Single-Cell Evidence Prioritizes MAP1LC3A among Ferroptosis-Related Genes in Glioblastoma.

Glioblastoma (GBM) remains a highly aggressive malignancy, and the contribution of ferroptosis-related genes to disease susceptibility remains incompletely understood. A genetically anchored, multistage framework was applied to prioritize ferroptosis-related genes associated with GBM. Among 483 genes curated from FerrDb V2, 315 had candidate cis-expression quantitative trait loci (cis-eQTLs) in eQTLGen, 250 retained at least three independent instruments after linkage disequilibrium clumping, and 226 yielded valid inverse-variance weighted (IVW) Mendelian randomization estimates using a GBM genome-wide association study comprising 6,183 cases and 18,169 controls. Thirty-four genes met the exploratory discovery criteria of P < 0.05 and a Benjamini-Hochberg false discovery rate (BH-FDR) < 0.20, with directionally concordant Bayesian weighted Mendelian randomization (BWMR) estimates. Replication-stage Mendelian randomization using GTEx V10 whole-blood cis-eQTLs supported four genes: ATG7, RPTOR, MAP1LC3A, and CHMP6. Evaluation across three independent tumor-control transcriptomic cohorts demonstrated that MAP1LC3A was consistently downregulated in tumor tissue and showed a significant random-effects pooled estimate (log&#x2082; fold change, -1.273; 95% confidence interval, -1.625 to -0.920; false discovery rate = 0.016), whereas the other three genes lacked comparable cross-cohort statistical support. Single-cell virtual knockout analysis was subsequently performed in a patient-balanced subset of 2,400 malignant cells selected from 4,916 eligible cells across 20 adult IDH-wild-type GBM tumors. Across five independently seeded runs, 3, 15, 4, and 7 robust downstream genes were identified for ATG7, RPTOR, MAP1LC3A, and CHMP6, respectively. The resulting consensus sets comprised 17 unique genes, with RND3 shared across all four targets. Gene Ontology analysis indicated enrichment of cell-adhesion and cell-surface processes, whereas no KEGG or Reactome pathways remained significant after multiple-testing correction. Collectively, these findings prioritize MAP1LC3A for future experimental investigation while distinguishing genetic association, tumor-expression concordance, and computational perturbation from definitive evidence of causality or mechanism.

Humans