Search

Resistance training is used to combat skeletal muscle function decline in older adults. Few studies have been designed specific for females, resulting in very limited treatment options for skeletal muscle atrophy in aging women. Here, we analyzed the gene expression profiles of skeletal muscle samples from sedentary young women, sedentary older women, and resistance-trained older women, using microarray data from public database. A total of 45 genes that were differentially expressed during female muscle aging and reversed by resistance training were identified. Functional and pathway enrichment analysis, protein–protein interaction network analysis, and receiver operating characteristic analysis were performed to reveal the key genes and pathways involved in the effects of resistance training on female muscle aging. The collagen genes COL1A1, COL3A1, and COL4A1 were identified important regulators of female muscle aging and resistance training, by modulating multiple signaling pathways, such as PI3 kinase-Akt signaling, focal adhesions, extracellular matrix-receptor interactions, and relaxin signaling. Interestingly, the expression of CDKN1A and TP63 were increased during aging, and further upregulated by resistance training in older women, suggesting they may negatively affect resistance training outcomes. Our findings provide novel insights into the molecular mechanisms of resistance training on female muscle aging and identify potential biomarkers and targets for clinical intervention. [ABSTRACT FROM AUTHOR]

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

ObjectivePeople with HIV (PWH) are at greater risk of developing lung diseases even when they are antiretroviral therapy (ART)-adherent and virally suppressed. The most common pulmonary function abnormality in PWH is that of impaired diffusing capacity of the lungs for carbon monoxide (DL CO ), which is an independent risk factor for increased mortality in PWH. Earlier work has identified several plasma biomarkers of inflammation and immune activation to be associated with decreased DL CO . However, the underpinning molecular mechanisms of HIV-associated impaired DL CO are largely unknown.DesignCross-sectional pilot study with PWH with normal DL CO (values greater than or equal to the lower limit of normal, DL CO  ≥ LLN, N = 9) or abnormal DL CO (DL CO

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p 0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

Bone marrow aspirate concentrate (BMAC) and adipose-derived stromal vascular fraction (ADSVF) are the most marketed stem cell therapies to treat a variety of conditions in the general population and elite athletes. Both tissues have been used interchangeably clinically even though their detailed composition, heterogeneity, and mechanisms of action have neither been rigorously inventoried nor compared. This lack of information has prevented investigations into ideal dosages and has facilitated anecdata and misinformation. Here, we analyzed single-cell transcriptomes, proteomes, and flow cytometry profiles from paired clinical-grade BMAC and ADSVF. This comparative transcriptional atlas challenges the prevalent notion that there is one therapeutic cell type present in both tissues. We also provide data of surface markers that may enable isolation and investigation of cell (sub)populations. Furthermore, the proteome atlas highlights intertissue and interpatient heterogeneity of injected proteins with potentially regenerative or immunomodulatory capacities. An interactive webtool is available online.

Cardiac allograft vasculopathy (CAV) causes late graft failure and mortality after heart transplantation. Donor-specific antibodies (DSAs) lead to chronic endothelial cell injury, inflammation, and arterial intimal thickening. In this study, GeoMx digital spatial profiling was used to analyze arterial areas of interest (AOIs) from CAV+DSA+ rejected cardiac allografts (N = 3; 22 AOIs total). AOIs were categorized based on CAV neointimal thickening and underwent whole transcriptome and protein profiling. By comparing our transcriptomic data with that of healthy control vessels of rapid autopsy myocardial tissue, we pinpointed specific pathways and transcripts indicative of heightened inflammatory profiles in CAV lesions. Moreover, we identified protein and transcriptomic signatures distinguishing CAV lesions exhibiting low and high neointimal lesions. AOIs with low neointima showed increased markers for activated inflammatory infiltrates, endothelial cell activation transcripts, and gene modules involved in metalloproteinase activation and TP53 regulation of caspases. Inflammatory and apoptotic proteins correlated with inflammatory modules in low neointima AOIs. High neointima AOIs exhibited elevated TGFβ-regulated transcripts and modules enriched for platelet activation/aggregation. Proteins associated with growth factors/survival correlated with modules enriched for proliferation/repair in high neointima AOIs. Our findings reveal novel insight into immunological mechanisms mediating CAV pathogenesis.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in FSHD-hi and FSHD-lo states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

During pulmonary mucormycosis, inhaled sporangiospores adhere to, germinate, and invade airway epithelial cells to establish infection. We provide evidence that HIF1α plays dual roles in airway epithelial cells during Mucorales infection. We observed an increase in HIF1α protein accumulation and increased expression of many known HIF1α-responsive genes during in vitro infection, indicating that HIF1α signaling is activated by Mucorales infection. Inhibition of HIF1α signaling led to a substantial decrease in the ability of R. delemar to invade cultured airway epithelial cells. Transcriptome analysis revealed that R. delemar infection induces the expression of many pro-inflammatory genes whose expression was significantly reduced by HIF1α inhibition. Importantly, pharmacological inhibition of HIF1α increased survival in a mouse model of pulmonary mucormycosis without reducing fungal burden. These results suggest that HIF1α plays two opposing roles during mucormycosis: one that facilitates the ability of Mucorales to invade the host cells and one that facilitates the ability of the host to mount an innate immune response.

It has proved challenging to quantitatively relate the proteome to the transcriptome on a per-gene basis. Recent advances in data analytics have enabled a biologically meaningful modularization of the bacterial transcriptome. We thus investigate whether matched datasets of transcriptomes and proteomes from bacteria under diverse conditions can be modularized in the same way to reveal novel relationships between their compositions. We find that; (1) the modules of the proteome and the transcriptome are comprised of a similar list of gene products, (2) the modules in the proteome often represent combinations of modules from the transcriptome, (3) known transcriptional and post-translational regulation is reflected in differences between two sets of modules, allowing for knowledge-mapping when interpreting module functions, and (4) through statistical modeling, absolute proteome allocation can be inferred from the transcriptome alone. Quantitative and knowledge-based relationships can thus be found at the genome-scale between the proteome and transcriptome in bacteria.

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

Meningiomas, although mostly benign, can be recurrent and fatal. World Health Organization (WHO) grading of the tumor does not always identify high-risk meningioma, and better characterizations of their aggressive biology are needed. To approach this problem, we combined 13 bulk RNA sequencing (RNA-seq) datasets to create a dimension-reduced reference landscape of 1,298 meningiomas. The clinical and genomic metadata effectively correlated with landscape regions, which led to the identification of meningioma subtypes with specific biological signatures. The time to recurrence also correlated with the map location. Further, we developed an algorithm that maps new patients onto this landscape, where the nearest neighbors predict outcome. This study highlights the utility of combining bulk transcriptomic datasets to visualize the complexity of tumor populations. Further, we provide an interactive tool for understanding the disease and predicting patient outcomes. This resource is accessible via the online tool Oncoscape, where the scientific community can explore the meningioma landscape.

MOTIVATION: Recent advances in spatial transcriptomics allow spatially resolved gene expression measurements with cellular or even sub-cellular resolution, directly characterizing the complex spatiotemporal gene expression landscape and cell-to-cell interactions in their native microenvironments. Due to technology limitations, most spatial transcriptomic technologies still yield incomplete expression measurements with excessive missing values. Therefore, gene imputation is critical to filling in missing data, enhancing resolution, and improving overall interpretability. However, existing methods either require additional matched single-cell RNA-seq data, which is rarely available, or ignore spatial proximity or expression similarity information. RESULTS: To address these issues, we introduce Impeller, a path-based heterogeneous graph learning method for spatial transcriptomic data imputation. Impeller has two unique characteristics distinct from existing approaches. First, it builds a heterogeneous graph with two types of edges representing spatial proximity and expression similarity. Therefore, Impeller can simultaneously model smooth gene expression changes across spatial dimensions and capture similar gene expression signatures of faraway cells from the same type. Moreover, Impeller incorporates both short- and long-range cell-to-cell interactions (e.g. via paracrine and endocrine) by stacking multiple GNN layers. We use a learnable path operator in Impeller to avoid the over-smoothing issue of the traditional Laplacian matrices. Extensive experiments on diverse datasets from three popular platforms and two species demonstrate the superiority of Impeller over various state-of-the-art imputation methods. AVAILABILITY AND IMPLEMENTATION: The code and preprocessed data used in this study are available at https://github.com/aicb-ZhangLabs/Impeller and https://zenodo.org/records/11212604.

With aging skeletal muscle fibers undergo repeating cycles of denervation and reinnervation. In approximately the 8th decade of life reinnervation no longer keeps pace, resulting in the accumulation of persistently denervated muscle fibers that in turn cause an acceleration of muscle dysfunction. The significance of denervation in important clinical outcomes with aging is poorly studied. The Study of Muscle, Mobility, and Aging (SOMMA) is a large cohort study with the primary objective to assess how aging muscle biology impacts clinically important traits. Using transcriptomics data from vastus lateralis muscle biopsies in 575 participants we have selected 49 denervation-responsive genes to provide insights to the burden of denervation in SOMMA, to test the hypothesis that greater expression of denervation-responsive genes negatively associates with SOMMA participant traits that included time to walk 400 meters, fitness (VO2peak), maximal mitochondrial respiration, muscle mass and volume, and leg muscle strength and power. Consistent with our hypothesis, increased transcript levels of: a calciumdependent intercellular adhesion glycoprotein (CDH15), acetylcholine receptor subunits (CHRNA1, CHRND, CHRNE), a glycoprotein promoting reinnervation (NCAM1), a transcription factor regulating aspects of muscle organization (RUNX1), and a sodium channel (SCN5A) were each negatively associated with at least 3 of these traits. VO2peak and maximal respiration had the strongest negative associations with 15 and 19 denervation-responsive genes, respectively. In conclusion, the abundance of denervationresponsive gene transcripts is a significant determinant of muscle and mobility outcomes in aging humans, supporting the imperative to identify new treatment strategies to restore innervation in advanced age.

A quarter of humanity is estimated to have been exposed to Mycobacterium tuberculosis (Mtb) with a 5-10% risk of developing tuberculosis (TB) disease. Variability in responses to Mtb infection could be due to host or pathogen heterogeneity. Here, we focused on host genetic variation in a Peruvian population and its associations with gene regulation in monocyte-derived macrophages and dendritic cells (DCs). We recruited former household contacts of TB patients who previously progressed to TB (cases, n = 63) or did not progress to TB (controls, n = 63). Transcriptomic profiling of monocyte-derived DCs and macrophages measured the impact of genetic variants on gene expression by identifying expression quantitative trait loci (eQTL). We identified 330 and 257 eQTL genes in DCs and macrophages (False Discovery Rate (FDR) < 0.05), respectively. Four genes in DCs showed interaction between eQTL variants and TB progression status. The top eQTL interaction for a protein-coding gene was with FAH, the gene encoding fumarylacetoacetate hydrolase, which mediates the last step in mammalian tyrosine catabolism. FAH expression was associated with genetic regulatory variation in cases but not controls. Using public transcriptomic and epigenomic data of Mtb-infected monocyte-derived dendritic cells, we found that Mtb infection results in FAH downregulation and DNA methylation changes in the locus. Overall, this study demonstrates effects of genetic variation on gene expression levels that are dependent on history of infectious disease and highlights a candidate pathogenic mechanism through pathogen-response genes. Furthermore, our results point to tyrosine metabolism and related candidate TB progression pathways for further investigation.

BackgroundSkeletal muscle dysfunction is a common extrapulmonary manifestation of chronic obstructive pulmonary disease (COPD). Alterations in skeletal muscle myosin heavy chain expression, with reduced type I and increased type II myosin heavy chain expression, are associated with COPD severity when studied in largely male cohorts. The objectives of this study were (1) to define an abnormal myofibre proportion phenotype in both males and females with COPD and (2) to identify transcripts and transcriptional networks associated with abnormal myofibre proportion in COPD.MethodsForty-six participants with COPD were assessed for body composition, strength, endurance and pulmonary function. Skeletal muscle biopsies from the vastus lateralis were assayed for fibre-type distribution and cross-sectional area via immunofluorescence microscopy and RNA-sequenced to generate transcriptome-wide gene expression data. Sex-stratified k-means clustering of type I and IIx/IIax fibre proportions was used to define abnormal myofibre proportion in participants with COPD and contrasted with previously defined criteria. Single transcripts and weighted co-expression network analysis modules were tested for correlation with the abnormal myofibre proportion phenotype.ResultsAbnormal myofibre proportion was defined in males with COPD (n = 29) as 22% type IIx/IIax fibres and in females with COPD (n = 17) as 12% type IIx/IIax fibres. Half of the participants with COPD were classified as having an abnormal myofibre proportion. Participants with COPD and an abnormal myofibre proportion had lower median handgrip strength (26.1 vs. 34.0 kg, P = 0.022), 6-min walk distance (300 vs. 353 m, P = 0.039) and forced expiratory volume in 1 s-to-forced vital capacity ratio (0.42 vs. 0.48, P = 0.041) compared with participants with COPD and normal myofibre proportions. Twenty-nine transcripts were associated with abnormal myofibre proportions in participants with COPD, with the upregulated NEB, TPM1 and TPM2 genes having the largest fold differences. Co-expression network analysis revealed that two transcript modules were significantly positively associated with the presence of abnormal myofibre proportions. One of these co-expression modules contained genes classically associated with muscle atrophy, as well as transcripts associated with both type I and type II myofibres, and was enriched for genetic loci associated with bone mineral density.ConclusionsOur findings indicate that there are significant transcriptional alterations associated with abnormal myofibre proportions in participants with COPD. Transcripts canonically associated with both type I and type IIa fibres were enriched in a co-expression network associated with abnormal myofibre proportion, suggesting altered transcriptional regulation across multiple fibre types.

Correlating transcriptional profiles with imaging-derived phenotypes has the potential to reveal possible molecular architectures associated with cognitive functions, brain development and disorders. Competitive null models built by resampling genes and self-contained null models built by spinning brain regions, along with varying test statistics, have been used to determine the significance of transcriptional associations. However, there has been no systematic evaluation of their performance in imaging transcriptomics analyses. Here, we evaluated the performance of eight different test statistics (mean, mean absolute value, mean squared value, max mean, median, Kolmogorov-Smirnov (KS), Weighted KS and the number of significant correlations) in both competitive null models and self-contained null models. Simulated brain maps (n = 1,000) and gene sets (n = 500) were used to calculate the probability of significance (Psig) for each statistical test. Our results suggested that competitive null models may result in false positive results driven by co-expression within gene sets. Furthermore, we demonstrated that the self-contained null models may fail to account for distribution characteristics (e.g., bimodality) of correlations between all available genes and brain phenotypes, leading to false positives. These two confounding factors interacted differently with test statistics, resulting in varying outcomes. Specifically, the sign-sensitive test statistics (i.e., mean, median, KS, Weighted KS) were influenced by co-expression bias in the competitive null models, while median and sign-insensitive test statistics were sensitive to the bimodality bias in the self-contained null models. Additionally, KS-based statistics produced conservative results in the self-contained null models, which increased the risk of false negatives. Comprehensive supplementary analyses with various configurations, including realistic scenarios, supported the results. These findings suggest utilizing sign-insensitive test statistics such as mean absolute value, max mean in the competitive null models and the mean as the test statistic for the self-contained null models. Additionally, adopting the confounder-matched (e.g., coexpression-matched) null models as an alternative to standard null models can be a viable strategy. Overall, the present study offers insights into the selection of statistical tests for imaging transcriptomics studies, highlighting areas for further investigation and refinement in the evaluation of novel and commonly used tests.

BACKGROUND: Single-cell transcriptome sequencing (scRNA-Seq) has allowed new types of investigations at unprecedented levels of resolution. Among the primary goals of scRNA-Seq is the classification of cells into distinct types. Many approaches build on existing clustering literature to develop tools specific to single-cell. However, almost all of these methods rely on heuristics or user-supplied parameters to control the number of clusters. This affects both the resolution of the clusters within the original dataset as well as their replicability across datasets. While many recommendations exist, in general, there is little assurance that any given set of parameters will represent an optimal choice in the trade-off between cluster resolution and replicability. For instance, another set of parameters may result in more clusters that are also more replicable. RESULTS: Here, we propose Dune, a new method for optimizing the trade-off between the resolution of the clusters and their replicability. Our method takes as input a set of clustering results-or partitions-on a single dataset and iteratively merges clusters within each partitions in order to maximize their concordance between partitions. As demonstrated on multiple datasets from different platforms, Dune outperforms existing techniques, that rely on hierarchical merging for reducing the number of clusters, in terms of replicability of the resultant merged clusters as well as concordance with ground truth. Dune is available as an R package on Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/Dune.html . CONCLUSIONS: Cluster refinement by Dune helps improve the robustness of any clustering analysis and reduces the reliance on tuning parameters. This method provides an objective approach for borrowing information across multiple clusterings to generate replicable clusters most likely to represent common biological features across multiple datasets.

We present the Manatee variational autoencoder model to predict transcription factor (TF) perturbation-induced transcriptomes. We demonstrate that the Manatee in silico perturbation analysis recapitulates target transcriptomic phenotypes in diverse cellular lineage transitions. We further propose the Manatee in silico screening analysis for prioritizing TF combinations targeting desired transcriptomic phenotypes.

The SUVmax is a measure of FDG uptake and is related with tumor aggressiveness in thyroid cancer, however, its association with molecular pathways is unclear. Here, we investigated the relationship between SUVmax and gene expression profiles in 80 papillary thyroid cancer (PTC) patients. We conducted an analysis of DEGs and enriched pathways in relation to SUVmax and tumor size. SUVmax showed a positive correlation with tumor size and correlated with glucose metabolic process. The genes that indicate thyroid differentiation, such as SLC5A5 and TPO, were negatively correlated with SUVmax. Unsupervised analysis revealed that SUVmax positively correlated with DNA replication(r = 0.29, p = 0.009), pyrimidine metabolism(r = 0.50, p

Heat stress exposure in intermittent heat waves and subsequent exposure during war theaters pose a clinical challenge that can lead to multi-organ dysfunction and long-term complications in the elderly. Using an aged mouse model and high-throughput sequencing, this study investigated the molecular dynamics of the liver-brain connection during heat stress exposure. Distinctive gene expression patterns induced by periodic heat stress emerged in both brain and liver tissues. An altered transcriptome profile showed heat stress-induced altered acute phase response pathways, causing neural, hepatic, and systemic inflammation and impaired synaptic plasticity. Results also demonstrated that proinflammatory molecules such as S100B, IL-17, IL-33, and neurological disease signaling pathways were upregulated, while protective pathways like aryl hydrocarbon receptor signaling were downregulated. In parallel, Rantes, IRF7, NOD1/2, TREM1, and hepatic injury signaling pathways were upregulated. Furthermore, current research identified Orosomucoid 2 (ORM2) in the liver as one of the mediators of the liver-brain axis due to heat exposure. In conclusion, the transcriptome profiling in elderly heat-stressed mice revealed a coordinated network of liver-brain axis pathways with increased hepatic ORM2 secretion, possibly due to gut inflammation and dysbiosis. The above secretion of ORM2 may impact the brain through a leaky blood-brain barrier, thus emphasizing intricate multi-organ crosstalk.

BACKGROUND: Both defects in mismatch repair (dMMR) and high microsatellite instability (MSI-H) have been recognised as crucial biomarkers that guide treatment strategies and disease management in colorectal cancer (CRC). As MMR and MSI tests are being widely conducted, an increasing number of MSI-H tumours have been identified in CRCs with mismatch repair proficiency (pMMR). The objective of this study was to assess the clinical features of patients with pMMR/MSI-H CRC and elucidate the underlying molecular mechanism in these cases. METHODS: From January 2015 to December 2018, 1684 cases of pMMR and 401 dMMR CRCs were enrolled. Of those patients, 93 pMMR/MSI-H were identified. The clinical phenotypes and prognosis were analysed. Frozen and paraffin-embedded tissue were available in 35 patients with pMMR/MSI-H, for which comprehensive genomic and transcriptomic analyses were performed. FINDINGS: In comparison to pMMR/MSS CRCs, pMMR/MSI-H CRCs exhibited significantly less tumour progression and better long-term prognosis. The pMMR/MSI-H cohorts displayed a higher presence of CD8+ T cells and NK cells when compared to the pMMR/MSS group. Mutational signature analysis revealed that nearly all samples exhibited deficiencies in MMR genes, and we also identified deleterious mutations in MSH3-K383fs. INTERPRETATION: This study revealed pMMR/MSI-H CRC as a distinct subgroup within CRC, which manifests diverse clinicopathological features and long-term prognostic outcomes. Distinct features in the tumour immune-microenvironment were observed in pMMR/MSI-H CRCs. Pathogenic deleterious mutations in MSH3-K383fs were frequently detected, suggesting another potential biomarker for identifying MSI-H. FUNDING: This work was supported by the Science and Technology Commission of Shanghai Municipality (20DZ1100101).

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

BACKGROUND: Sulfur (S) is a mineral nutrient essential for plant growth and development, which is incorporated into diverse molecules fundamental for primary and secondary metabolism, plant defense, signaling, and maintaining cellular homeostasis. Although, S starvation response is well documented in the dicot model Arabidopsis thaliana, it is not clear if the same transcriptional networks control the response also in the monocots. RESULTS: We performed series of physiological, expression, and metabolite analyses in two model monocot species, one representing the C3 plants, Oryza sativa cv. kitaake, and second representing the C4 plants, Setaria viridis. Our comprehensive transcriptomic analysis revealed twice as many differentially expressed genes (DEGs) in S. viridis than in O. sativa under S-deficiency, consistent with a greater loss of sulfur and S-containing metabolites under these conditions. Surprisingly, most of the DEGs and enriched gene ontology terms were species-specific, with an intersect of only 58 common DEGs. The transcriptional networks were different in roots and shoots of both species, in particular no genes were down-regulated by S-deficiency in the roots of both species. CONCLUSIONS: Our analysis shows that S-deficiency seems to have different physiological consequences in the two monocot species and their nutrient homeostasis might be under distinct control mechanisms.

Bridging molecular information to ecosystem-level processes would provide the capacity to understand system vulnerability and, potentially, a means for assessing ecosystem health. Here, we present an integrated dataset containing environmental and metagenomic information from plant-associated microbial communities, plant transcriptomics, plant and soil metabolomics, and soil chemistry and activity characterization measurements derived from the model tree species Populus trichocarpa. Soil, rhizosphere, root endosphere, and leaf samples were collected from 27 different P. trichocarpa genotypes grown in two different environments leading to an integrated dataset of 318 metagenomes, 98 plant transcriptomes, and 314 metabolomic profiles that are supported by diverse soil measurements. This expansive dataset will provide insights into causal linkages that relate genomic features and molecular level events to system-level properties and their environmental influences.

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an older transcriptomic signature and Alzheimers disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Single-cell and single-nucleus RNA-sequencing technologies capture the expression of plant genes at an unprecedented resolution. Therefore, these technologies are gaining traction in plant molecular and developmental biology for elucidating the transcriptional changes across cell types in a specific tissue or organ, upon treatments, in response to biotic and abiotic stresses, or between genotypes. Despite the rapidly accelerating use of these technologies, collective and standardized experimental and analytical procedures to support the acquisition of high-quality data sets are still missing. In this commentary, we discuss common challenges associated with the use of single-cell transcriptomics in plants and propose general guidelines to improve reproducibility, quality, comparability, and interpretation and to make the data readily available to the community in this fast-developing field of research.

In situ transcriptomic techniques promise a holistic view of tissue organization and cell-cell interactions. There has been a surge of multiplexed RNA in situ mapping techniques but their application to human tissues has been limited due to their large size, general lower tissue quality and high autofluorescence. Here we report DART-FISH, a padlock probe-based technology capable of profiling hundreds to thousands of genes in centimeter-sized human tissue sections. We introduce an omni-cell type cytoplasmic stain that substantially improves the segmentation of cell bodies. Our enzyme-free isothermal decoding procedure allows us to image 121 genes in large sections from the human neocortex in 20 healthy and pathological cell states, and identified diseased niches enriched in transcriptionally altered epithelial cells and myofibroblasts.

UNLABELLED: Limosilactobacillus reuteri, a probiotic microbe instrumental to human health and sustainable food production, adapts to diverse environmental shifts via dynamic gene expression. We applied the independent component analysis (ICA) to 117 RNA-seq data sets to decode its transcriptional regulatory network (TRN), identifying 35 distinct signals that modulate specific gene sets. Our findings indicate that the ICA provides a qualitative advancement and captures nuanced relationships within gene clusters that other methods may miss. This study uncovers the fundamental properties of L. reuteris TRN and deepens our understanding of its arginine metabolism and the co-regulation of riboflavin metabolism and fatty acid conversion. It also sheds light on conditions that regulate genes within a specific biosynthetic gene cluster and allows for the speculation of the potential role of isoprenoid biosynthesis in L. reuteris adaptive response to environmental changes. By integrating transcriptomics and machine learning, we provide a system-level understanding of L. reuteris response mechanism to environmental fluctuations, thus setting the stage for modeling the probiotic transcriptome for applications in microbial food production. IMPORTANCE: We have studied Limosilactobacillus reuteri, a beneficial probiotic microbe that plays a significant role in our health and production of sustainable foods, a type of foods that are nutritionally dense and healthier and have low-carbon emissions compared to traditional foods. Similar to how humans adapt their lifestyles to different environments, this microbe adjusts its behavior by modulating the expression of genes. We applied machine learning to analyze large-scale data sets on how these genes behave across diverse conditions. From this, we identified 35 unique patterns demonstrating how L. reuteri adjusts its genes based on 50 unique environmental conditions (such as various sugars, salts, microbial cocultures, human milk, and fruit juice). This research helps us understand better how L. reuteri functions, especially in processes like breaking down certain nutrients and adapting to stressful changes. More importantly, with our findings, we become closer to using this knowledge to improve how we produce more sustainable and healthier foods with the help of microbes.

Macrophages are abundant immune cells in the microenvironment of diffuse large B-cell lymphoma (DLBCL). Macrophage estimation by immunohistochemistry shows varying prognostic significance across studies in DLBCL, and does not provide a comprehensive analysis of macrophage subtypes. Here, using digital spatial profiling with whole transcriptome analysis of CD68+ cells, we characterize macrophages in distinct spatial niches of reactive lymphoid tissues (RLTs) and DLBCL. We reveal transcriptomic differences between macrophages within RLTs (light zone /dark zone, germinal center/ interfollicular), and between disease states (RLTs/ DLBCL), which we then use to generate six spatially-derived macrophage signatures (MacroSigs). We proceed to interrogate these MacroSigs in macrophage and DLBCL single-cell RNA-sequencing datasets, and in gene-expression data from multiple DLBCL cohorts. We show that specific MacroSigs are associated with cell-of-origin subtypes and overall survival in DLBCL. This study provides a spatially-resolved whole-transcriptome atlas of macrophages in reactive and malignant lymphoid tissues, showing biological and clinical significance.

Without new transcription, gene expression across the oocyte-to-embryo transition (OET) relies instead on regulation of mRNA poly(A) tails to control translation. However, how tail dynamics shape translation across the OET in mammals remains unclear. We perform long-read RNA sequencing to uncover poly(A) tail lengths across the mouse OET and, incorporating published ribosome profiling data, provide an integrated, transcriptome-wide analysis of poly(A) tails and translation across the entire transition. We uncover an extended wave of global deadenylation during fertilization in which short-tailed, oocyte-deposited mRNAs are translationally activated without polyadenylation through resistance to deadenylation. Subsequently, in the embryo, mRNAs are readenylated and translated in a surge of global polyadenylation. We further identify regulation of poly(A) tail length at the isoform level and stage-specific enrichment of mRNA sequence motifs among regulated transcripts. These data provide insight into the stage-specific mechanisms of poly(A) tail regulation that orchestrate gene expression from oocyte to embryo in mammals.

Alternative polyadenylation plays an important role in cancer initiation and progression; however, current transcriptome-wide association studies mostly ignore alternative polyadenylation when identifying putative cancer susceptibility genes. Here, we perform a pan-cancer 3 untranslated region alternative polyadenylation transcriptome-wide association analysis by integrating 55 well-powered (n > 50,000) genome-wide association studies datasets across 22 major cancer types with alternative polyadenylation quantification from 23,955 RNA sequencing samples across 7,574 individuals. We find that genetic variants associated with alternative polyadenylation are co-localized with 28.57% of cancer loci and contribute a significant portion of cancer heritability. We further identify 642 significant cancer susceptibility genes predicted to modulate cancer risk via alternative polyadenylation, 62.46% of which have been overlooked by traditional expression- and splicing- studies. As proof of principle validation, we show that alternative alleles facilitate 3 untranslated region lengthening of CRLS1 gene leading to increased protein abundance and promoted proliferation of breast cancer cells. Together, our study highlights the significant role of alternative polyadenylation in discovering new cancer susceptibility genes and provides a strong foundational framework for enhancing our understanding of the etiology underlying human cancers.

Childhood maltreatment has been repeatedly linked to a higher incidence of health conditions with an underlying proinflammatory component, such as asthma, chronic obstructive pulmonary disease, stroke, and cardiovascular disease. Childhood maltreatment has also been linked to elevated systemic inflammation prior to the onset of disease. However, childhood maltreatment is highly comorbid with other risk factors which have also been linked to inflammation, namely major depression. The present analysis addresses this issue by assessing the association of maltreatment with genome-wide transcriptional profiling of immune cells collected from four orthogonal groups of adolescents (aged 13-17): maltreated and not maltreated in childhood, with and without major depressive disorder. Maltreatment and psychiatric history were determined using semi-structured clinical interviews and cross-validated using self-report questionnaires. Dried whole blood spots were collected from each participant (n = 133) and assayed to determine the extent to which maltreatment in childhood was associated with a higher prevalence of transcriptional activity among differentially expressed genes, specific immune cell subtypes, and up- or down-regulation of genes involved in immune function after accounting for current major depression. Maltreatment was associated with increased interferon regulatory factor (IRF) transcriptional activity (p = 0.03), as well as nuclear factor erythroid-2 related factor 1 (NRF1; p = 0.002) and MAF (p = 0.01) among up-regulated genes, and increased activity of nuclear factor kappa beta (NF-κB) among down-regulated genes (p = 0.01). Non-classical CD16+ monocytes were implicated in both the up- and down-regulated genes among maltreated adolescents. These data provide convergent evidence supporting the role of maltreatment in altering intracellular and molecular markers of immune function, as well as implicate monocyte/macrophage functions as mechanisms through which childhood maltreatment may shape lifelong immune development and function.

BACKGROUND: Diabetes is a spectrum of metabolic diseases affecting millions of people worldwide. The loss of pancreatic β-cell mass by either autoimmune destruction or apoptosis, in type 1-diabetes (T1D) and type 2-diabetes (T2D), respectively, represents a pathophysiological process leading to insulin deficiency. Therefore, therapeutic strategies focusing on restoring β-cell mass and β-cell insulin secretory capacity may impact disease management. This study took advantage of powerful integrative bioinformatic tools to scrutinize publicly available diabetes-associated gene expression data to unveil novel potential molecular targets associated with β-cell dysfunction. METHODS: A comprehensive literature search for human studies on gene expression alterations in the pancreas associated with T1D and T2D was performed. A total of 6 studies were selected for data extraction and for bioinformatic analysis. Pathway enrichment analyses of differentially expressed genes (DEGs) were conducted, together with protein-protein interaction networks and the identification of potential transcription factors (TFs). For noncoding differentially expressed RNAs, microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), which exert regulatory activities associated with diabetes, identifying target genes and pathways regulated by these RNAs is fundamental for establishing a robust regulatory network. RESULTS: Comparisons of DEGs among the 6 studies showed 59 genes in common among 4 or more studies. Besides alterations in mRNA, it was possible to identify differentially expressed miRNA and lncRNA. Among the top transcription factors (TFs), HIPK2, KLF5, STAT1 and STAT3 emerged as potential regulators of the altered gene expression. Integrated analysis of protein-coding genes, miRNAs, and lncRNAs pointed out several pathways involved in metabolism, cell signaling, the immune system, cell adhesion, and interactions. Interestingly, the GABAergic synapse pathway emerged as the only common pathway to all datasets. CONCLUSIONS: This study demonstrated the power of bioinformatics tools in scrutinizing publicly available gene expression data, thereby revealing potential therapeutic targets like the GABAergic synapse pathway, which holds promise in modulating α-cells transdifferentiation into β-cells.

BACKGROUND: Luminal A tumours generally have a favourable prognosis but possess the highest 10-year recurrence risk among breast cancers. Additionally, a quarter of the recurrence cases occur within 5 years post-diagnosis. Identifying such patients is crucial as long-term relapsers could benefit from extended hormone therapy, while early relapsers might require more aggressive treatment. METHODS: We conducted a study to explore non-structural chromosome maintenance condensin I complex subunit Hs (NCAPH) role in luminal A breast cancer pathogenesis, both in vitro and in vivo, aiming to identify an intratumoural gene expression signature, with a focus on elevated NCAPH levels, as a potential marker for unfavourable progression. Our analysis included transgenic mouse models overexpressing NCAPH and a genetically diverse mouse cohort generated by backcrossing. A least absolute shrinkage and selection operator (LASSO) multivariate regression analysis was performed on transcripts associated with elevated intratumoural NCAPH levels. RESULTS: We found that NCAPH contributes to adverse luminal A breast cancer progression. The intratumoural gene expression signature associated with elevated NCAPH levels emerged as a potential risk identifier. Transgenic mice overexpressing NCAPH developed breast tumours with extended latency, and in Mouse Mammary Tumor Virus (MMTV)-NCAPHErbB2 double-transgenic mice, luminal tumours showed increased aggressiveness. High intratumoural Ncaph levels correlated with worse breast cancer outcome and subpar chemotherapy response. A 10-gene risk score, termed Gene Signature for Luminal A 10 (GSLA10), was derived from the LASSO analysis, correlating with adverse luminal A breast cancer progression. CONCLUSIONS: The GSLA10 signature outperformed the Oncotype DX signature in discerning tumours with unfavourable outcomes, previously categorised as luminal A by Prediction Analysis of Microarray 50 (PAM50) across three independent human cohorts. This new signature holds promise for identifying luminal A tumour patients with adverse prognosis, aiding in the development of personalised treatment strategies to significantly improve patient outcomes.

In single-cell RNA sequencing (scRNA-Seq), gene expression is assessed individually for each cell, allowing the investigation of developmental processes, such as embryogenesis and cellular differentiation and regeneration, at unprecedented resolution. In such dynamic biological systems, cellular states form a continuum, e.g., for the differentiation of stem cells into mature cell types. This process is often represented via a trajectory in a reduced-dimensional representation of the scRNA-Seq dataset. While many methods have been suggested for trajectory inference, it is often unclear how to handle multiple biological groups or conditions, e.g., inferring and comparing the differentiation trajectories of wild-type and knock-out stem cell populations. In this manuscript, we present condiments, a method for the inference and downstream interpretation of cell trajectories across multiple conditions. Our framework allows the interpretation of differences between conditions at the trajectory, cell population, and gene expression levels. We start by integrating datasets from multiple conditions into a single trajectory. By comparing the cells conditions along the trajectorys path, we can detect large-scale changes, indicative of differential progression or fate selection. We also demonstrate how to detect subtler changes by finding genes that exhibit different behaviors between these conditions along a differentiation path.

Crown gall disease (Agrobacterium tumefaciens), crown/root rot disease (Phytophthora spp.), root lesion disease (Pratylenchus vulnus) and tree vigor are key traits affecting the productivity and quality of walnuts in California. Unchallenged hybrid rootstocks were analyzed by RNA-seq to examine pre-formed factors affecting these traits. Enrichment analysis of the differentially expressed genes revealed that the increased expression of cell wall biogenesis-related genes plays a key role in susceptibility to A. tumefaciens, susceptibility to Phytophthora spp. and increased vigor. Analysis of the predicted subcellular loci of the encoded proteins revealed that many gene products associated with vigor and susceptibility were targeted to the plasma membrane and extracellular space, connecting these traits to sustaining barrier function. We observed that RNA processing and splicing, along with predicted nuclear targeting, were associated with resistance to A. tumefaciens, resistance to Phytophthora spp. and low vigor. Four genes within the J. microcarpa QTL region for resistance to A. tumefaciens and Phytophthora spp. were represented among our transcripts, with two of the genes being differentially expressed in association with resistance to A. tumefaciens and decreased vigor. No differential expression related to Phytophthora spp. or P. vulnus resistance was observed in this region. Additionally, the J. microcarpa haplotype expressed more transcripts associated with resistance to A. tumefaciens, Phytophthora spp. and low vigor, but not P. vulnus, than the J. regia haplotype. We also report unique and shared hormone and defense responses associated with each trait. This research suggests a link between cell wall biogenesis, vigor and critical root diseases of walnut.

Issatchenkia orientalis is a promising industrial chassis to produce biofuels and bioproducts due to its high tolerance to multiple environmental stresses such as low pH, heat, and other chemicals otherwise toxic for the most widely used microbes. Yet, little is known about specific mechanisms of such tolerance in this organism, hindering our ability to engineer this species to produce valuable biochemicals. Here, we report a comprehensive study of the mechanisms of acidic tolerance in this species via transcriptome profiling across variable pH for 12 different strains with different phenotypes. We found multiple regulatory mechanisms involved in tolerance to low pH in different strains of I. orientalis, marking potential targets for future gene editing and perturbation experiments.

Laryngeal squamous cell carcinoma (LSCC) is a malignant tumor with a high probability of metastasis. The tumor microenvironment (TME) plays a critical role in cancer metastasis. To gain insights into the TME of LSCC, we conducted single-cell RNA-seq (scRNA-seq) on samples collected from LSCC patients with or without lymphatic metastasis. The stem and immune cell signatures in LSCC suggest their roles in tumor invasion and metastasis. Infiltration of a large number of regulatory T cells, dysplastic plasma cells, and macrophages that are at the early development stage in the cancerous tissue indicates an immunosuppressive state. Abundant neutrophils detected at the cancer margins reflect the inflammatory microenvironment. In addition to dynamic ligand-receptor interactions between the stromal and myeloid cells, the enhanced autophagy in endothelial cells and fibroblasts implies a role in nutrient supply. Taken together, the comprehensive atlas of LSCC obtained allowed us to identify a complex yet unique TME of LSCC, which may help identify potential diagnostic biomarkers and therapeutic targets for LSCC.

BackgroundThe accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.ResultsA total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.ConclusionsThese validated results show significant improvement over current bovine genome annotations.

BackgroundStatins lower circulating low-density lipoprotein cholesterol (LDLC) levels and reduce cardiovascular disease risk. Though highly efficacious in general, there is considerable inter-individual variation in statin efficacy that remains largely unexplained.MethodsTo identify novel genes that may modulate statin-induced LDLC lowering, we used RNA-sequencing data from 426 control- and 2 µM simvastatin-treated lymphoblastoid cell lines (LCLs) derived from European and African American ancestry participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial (ClinicalTrials.gov Identifier: NCT00451828). We correlated statin-induced changes in LCL gene expression with plasma LDLC statin response in the corresponding CAP participants. For the most correlated gene identified (ZNF335), we followed up in vivo by comparing plasma cholesterol levels, lipoprotein profiles, and lipid statin response between wild-type mice and carriers of a hypomorphic (partial loss of function) missense mutation in Zfp335 (the mouse homolog of ZNF335).ResultsThe statin-induced expression changes of 147 human LCL genes were significantly correlated to the plasma LDLC statin responses of the corresponding CAP participants in vivo (FDR = 5%). The two genes with the strongest correlations were zinc finger protein 335 (ZNF335 aka NIF-1, rho = 0.237, FDR-adj p = 0.0085) and CCR4-NOT transcription complex subunit 3 (CNOT3, rho = 0.233, FDR-adj p = 0.0085). Chow-fed mice carrying a hypomorphic missense (R1092W; aka bloto) mutation in Zfp335 had significantly lower non-HDL cholesterol levels than wild-type C57BL/6J mice in a sex combined model (p = 0.04). Furthermore, male (but not female) mice carrying the Zfp335R1092W allele had significantly lower total and HDL cholesterol levels than wild-type mice. In a separate experiment, wild-type mice fed a control diet for 4 weeks and a matched simvastatin diet for an additional 4 weeks had significant statin-induced reductions in non-HDLC (-43 ± 18% and -23 ± 19% for males and females, respectively). Wild-type male (but not female) mice experienced significant reductions in plasma LDL particle concentrations, while male mice carrying Zfp335R1092W allele(s) exhibited a significantly blunted LDL statin response.ConclusionsOur in vitro and in vivo studies identified ZNF335 as a novel modulator of plasma cholesterol levels and statin response, suggesting that variation in ZNF335 activity could contribute to inter-individual differences in statin clinical efficacy.

OBJECTIVES: To determine the dysregulated signaling pathways of head and neck squamous cell carcinoma associated with circulating tumor cells (CTCs) via single-cell molecular characterization. INTRODUCTION: Head and neck squamous cell carcinoma (HNSCC) has a significant global burden and is a disease with poor survival. Despite trials exploring new treatment modalities to improve disease control rates, the 5 year survival rate remains low at only 60%. Most cancer malignancies are reported to progress to a fatal phase due to the metastatic activity derived from treatment-resistant cancer cells, regarded as one of the most significant obstacles to develope effective cancer treatment options. However, the molecular profiles of cancer cells have not been thoroughly studied. METHODS: Here, we examined in-situ HNSCC tumors and pairwisely followed up with the downstream circulating tumor cells (CTCs)-based on the surrogate biomarkers to detect metastasis that is established in other cancers - not yet being fully adopted in HNSCC treatment algorithms. RESULTS: Specifically, we revealed metastatic HNSCC patients have complex CTCs that could be defined through gene expression and mutational gene profiling derived from completed single-cell RNASeq (scRNASeq) that served to confirm molecular pathways inherent in these CTCs. To enhance the reliability of our findings, we cross-validated those molecular profiles with results from previously published studies. CONCLUSION: Thus, we identified 5 dysregulated signaling pathways in CTCs to derive HNSCC biomarker panels for screening HNSCC in situ tumors.

Cardiometabolic disorders (CMD) are a growing public health problem across the world. Among the known cardiometabolic risk factors are compounds that induce endocrine and metabolic dysfunctions, such as endocrine disrupting chemicals (EDCs). To date, how EDCs influence molecular programs and cardiometabolic risks has yet to be fully elucidated, especially considering the complexity contributed by species-, chemical-, and dose-specific effects. Moreover, different experimental and analytical methodologies employed by different studies pose challenges when comparing findings across studies. To explore the molecular mechanisms of EDCs in a systematic manner, we established a data-driven computational approach to meta-analyze 30 human, mouse, and rat liver transcriptomic datasets for 4 EDCs, namely bisphenol A (BPA), bis(2-ethylhexyl) phthalate (DEHP), tributyltin (TBT), and perfluorooctanoic acid (PFOA). Our computational pipeline uniformly re-analyzed pre-processed quality-controlled microarray data and raw RNAseq data, derived differentially expressed genes (DEGs) and biological pathways, modeled gene regulatory networks and regulators, and determined CMD associations based on gene overlap analysis. Our approach revealed that DEHP and PFOA shared stable transcriptomic signatures that are enriched for genes associated with CMDs, suggesting similar mechanisms of action such as perturbations of peroxisome proliferator-activated receptor gamma (PPARγ) signaling and liver gene network regulators VNN1 and ACOT2. In contrast, TBT exhibited highly divergent gene signatures, pathways, network regulators, and disease associations from the other EDCs. In addition, we found that the rat, mouse, and human BPA studies showed highly variable transcriptomic patterns, providing molecular support for the variability in BPA responses. Our work offers insights into the commonality and differences in the molecular mechanisms of various EDCs and establishes a streamlined data-driven workflow to compare molecular mechanisms of environmental substances to elucidate the underlying connections between chemical exposure and disease risks.

IntroductionMaternal intervillous monocytes (MIMs) and fetal Hofbauer cells (HBCs) are myeloid-derived immune cells at the maternal-fetal interface. Maternal reproductive history is associated with differential risk of pregnancy complications. The molecular phenotypes and roles of these distinct monocyte/macrophage populations and the influence of gravidity on these phenotypes has not been systematically investigated.MethodsHere, we used RNA sequencing to study the transcriptional profiles of MIMs and HBCs in normal term pregnancies.ResultsOur analyses revealed distinct transcriptomes of MIMs and HBCs. Genes involved in differentiation and cell organization pathways were more highly expressed in MIMs vs. HBCs. In contrast, HBCs had higher expression of genes involved in inflammatory responses and cell surface receptor signaling. Maternal gravidity influenced monocyte programming, as expression of pro-inflammatory molecules was significantly higher in MIMs from multigravidae compared to primigravidae. In HBCs, multigravidae displayed enrichment of gene pathways involved in cell-cell signaling and differentiation.DiscussionOur results demonstrated that MIMs and HBCs have highly divergent transcriptional signatures, reflecting their distinct origins, locations, functions, and roles in inflammatory responses. Furthermore, maternal gravidity influences the gene signatures of MIMs and HBCs, potentially modulating the interplay between tolerance and trained immunity. The phenomenon of reproductive immune memory may play a novel role in the differential susceptibility of primigravidae to pregnancy complications.

This study focuses on understanding the transcriptional heterogeneity of activated platelets and its impact on diseases such as sepsis, COVID-19, and systemic lupus erythematosus (SLE). Recognizing the limited knowledge in this area, our research aims to dissect the complex transcriptional profiles of activated platelets to aid in developing targeted therapies for abnormal and pathogenic platelet subtypes. We analyzed single-cell transcriptional profiles from 47,977 platelets derived from 413 samples of patients with these diseases, utilizing Deep Neural Network (DNN) and eXtreme Gradient Boosting (XGB) to distinguish transcriptomic signatures predictive of fatal or survival outcomes. Our approach included source data annotations and platelet markers, along with SingleR and Seurat for comprehensive profiling. Additionally, we employed Uniform Manifold Approximation and Projection (UMAP) for effective dimensionality reduction and visualization, aiding in the identification of various platelet subtypes and their relation to disease severity and patient outcomes. Our results highlighted distinct platelet subpopulations that correlate with disease severity, revealing that changes in platelet transcription patterns can intensify endotheliopathy, increasing the risk of coagulation in fatal cases. Moreover, these changes may impact lymphocyte function, indicating a more extensive role for platelets in inflammatory and immune responses. This study identifies crucial biomarkers of platelet heterogeneity in serious health conditions, paving the way for innovative therapeutic approaches targeting platelet activation, which could improve patient outcomes in diseases characterized by altered platelet function.

INTRODUCTION: Staphylococcus aureus bacteremia (SAB) is a life-threatening infection particularly involving methicillin-resistant S. aureus (MRSA). In contrast to resolving MRSA bacteremia (RB), persistent MRSA bacteremia (PB) blood cultures remain positive despite appropriate antibiotic treatment. Host immune responses distinguishing PB vs. RB outcomes are poorly understood. Here, integrated transcriptomic, IL-10 cytokine levels, and genomic analyses sought to identify signatures differentiating PB vs. RB outcomes. METHODS: Whole-blood transcriptomes of propensity-matched PB (n=28) versus RB (n=30) patients treated with vancomycin were compared in one independent training patient cohort. Gene expression (GE) modules were analyzed and prioritized relative to host IL-10 cytokine levels and DNA methyltransferase-3A (DNMT3A) genotype. RESULTS: Differential expression of T and B lymphocyte gene expression early in MRSA bacteremia discriminated RB from PB outcomes. Significant increases in effector T and B cell signaling pathways correlated with RB, lower IL-10 cytokine levels and DNMT3A heterozygous A/C genotype. Importantly, a second PB and RB patient cohort analyzed in a masked manner demonstrated high predictive accuracy of differential signatures. DISCUSSION: Collectively, the present findings indicate that human PB involves dysregulated immunity characterized by impaired T and B cell responses associated with excessive IL-10 expression in context of the DNMT3A A/A genotype. These findings reveal distinct immunologic programs in PB vs. RB outcomes, enable future studies to define mechanisms by which host and/or pathogen drive differential signatures and may accelerate prediction of PB outcomes. Such prognostic assessment of host risk could significantly enhance early anti-infective interventions to avert PB and improve patient outcomes.

INTRODUCTION: Hypoxia is a common pathological driver contributing to various forms of pulmonary vascular diseases leading to pulmonary hypertension (PH). Pulmonary interstitial macrophages (IMs) play pivotal roles in immune and vascular dysfunction, leading to inflammation, abnormal remodeling, and fibrosis in PH. However, IMs response to hypoxia and their role in PH progression remain largely unknown. We utilized a murine model of hypoxia-induced PH to investigate the repertoire and functional profiles of IMs in response to acute and prolonged hypoxia, aiming to elucidate their contributions to PH development. METHODS: We conducted single-cell transcriptomic analyses to characterize the repertoire and functional profiles of murine pulmonary IMs following exposure to hypobaric hypoxia for varying durations (0, 1, 3, 7, and 21 days). Hallmark pathways from the mouse Molecular Signatures Database were utilized to characterize the molecular function of the IM subpopulation in response to hypoxia. RESULTS: Our analysis revealed an early acute inflammatory phase during acute hypoxia exposure (Days 1-3), which was resolved by Day 7, followed by a pro-remodeling phase during prolonged hypoxia (Days 7-21). These phases were marked by distinct subpopulations of IMs: MHCIIhiCCR2+EAR2+ cells characterized the acute inflammatory phase, while TLF+VCAM1hi cells dominated the pro-remodeling phase. The acute inflammatory phase exhibited enrichment in interferon-gamma, IL-2, and IL-6 pathways, while the pro-remodeling phase showed dysregulated chemokine production, hemoglobin clearance, and tissue repair profiles, along with activation of distinct complement pathways. DISCUSSION: Our findings demonstrate the existence of distinct populations of pulmonary interstitial macrophages corresponding to acute and prolonged hypoxia exposure, pivotal in regulating the inflammatory and remodeling phases of PH pathogenesis. This understanding offers potential avenues for targeted interventions, tailored to specific populations and distinct phases of the disease. Moreover, further identification of triggers for pro-remodeling IMs holds promise in unveiling novel therapeutic strategies for pulmonary hypertension.

The import of plant-derived small sugars by sugar transporters (STs) has received increasing interest due to its important biological role and great industrial potential. STs are important targets of genetic engineering to improve fungal plant biomass conversion. Comparatively analysis of the genome-wide prevalence and transcriptomics of STs was performed in four filamentous fungi: Aspergillus niger, Aspergillus nidulans, Penicillium subrubescens and Trichoderma reesei. Using phylogenetic analysis and literature mining, their predicted STs were divided into ten subfamilies with putative sugar specificities assigned. In addition, transcriptome analysis revealed complex expression profiles among different STs subfamilies and fungal species, indicating a sophisticated transcriptome regulation and functional diversity of fungal STs. Several STs showed strong co-expression with other genes involved in sugar utilization, encoding CAZymes and sugar catabolic enzymes. This study provides new insights into the diversity of STs at the genomic/transcriptomic level, facilitating their biochemical characterization and metabolic engineering.

Spatially resolved transcriptomics technologies have opened new avenues for understanding gene expression heterogeneity in spatial contexts. However, existing methods for identifying spatially variable genes often focus solely on statistical significance, limiting their ability to capture continuous expression patterns and integrate spot-level covariates. To address these challenges, we introduce spVC, a statistical method based on a generalized Poisson model. spVC seamlessly integrates constant and spatially varying effects of covariates, facilitating comprehensive exploration of gene expression variability and enhancing interpretability. Simulation and real data applications confirm spVC's accuracy in these tasks, highlighting its versatility in spatial transcriptomics analysis.

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.