The CD4 Versus CD8 T Cell Fate Decision: A Multiomics-Informed Perspective.

The choice of developing thymocytes to become CD8+ cytotoxic or CD4+ helper T cells has been intensely studied, but many of the underlying mechanisms remain to be elucidated. Recent multiomics approaches have provided much higher resolution analysis of gene expression in developing thymocytes than was previously achievable, thereby offering a fresh perspective on this question. Focusing on our recent studies using CITE-seq (cellular indexing of transcriptomes and epitopes) analyses of mouse thymocytes, we present a detailed timeline of RNA and protein expression changes during CD8 versus CD4 T cell differentiation. We also revisit our current understanding of the links between T cell receptor signaling and expression of the lineage-defining transcription factors ThPOK and RUNX3. Finally, we propose a sequential selection model to explain the tight linkage between MHC-I versus MHC-II recognition and T cell lineage choice. This model incorporates key aspects of previously proposed kinetic signaling, instructive, and stochastic/selection models.

Molecular and Spatial Signatures of Mouse Embryonic Endothelial Cells at Single-Cell Resolution.

BACKGROUND: Mature endothelial cells (ECs) are heterogeneous, with subtypes defined by tissue origin and position within the vascular bed (ie, artery, capillary, vein, and lymphatic). How this heterogeneity is established during the development of the vascular system, especially arteriovenous specification of ECs, remains incompletely characterized. METHODS: We used droplet-based single-cell RNA sequencing and multiplexed error-robust fluorescence in situ hybridization to define EC and EC progenitor subtypes from E9.5, E12.5, and E15.5 mouse embryos. We used trajectory inference to analyze the specification of arterial ECs (aECs) and venous ECs (vECs) from EC progenitors. Network analysis identified candidate transcriptional regulators of arteriovenous differentiation, which we tested by CRISPR (clustered regularly interspaced short palindromic repeats) loss of function in human-induced pluripotent stem cells undergoing directed differentiation to aECs or vECs (human-induced pluripotent stem cell-aECs or human-induced pluripotent stem cell-vECs). RESULTS: From the single-cell transcriptomes of 7682 E9.5 to E15.5 ECs, we identified 19 EC subtypes, including Etv2+Bnip3+ EC progenitors. Spatial transcriptomic analysis of 15 448 ECs provided orthogonal validation of these EC subtypes and established their spatial distribution. Most embryonic ECs were grouped by their vascular-bed types, while ECs from the brain, heart, liver, and lung were grouped by their tissue origins. Arterial (Eln, Dkk2, Vegfc, and Egfl8), venous (Fam174b and Clec14a), and capillary (Kcne3) marker genes were identified. Compared with aECs, embryonic vECs and capillary ECs shared fewer markers than their adult counterparts. Early capillary ECs with venous characteristics functioned as a branch point for differentiation of aEC and vEC lineages. CONCLUSIONS: Our results provide a spatiotemporal map of embryonic EC heterogeneity at single-cell resolution and demonstrate that the diversity of ECs in the embryo arises from both tissue origin and vascular-bed position. Developing aECs and vECs share common venous-featured capillary precursors and are regulated by distinct transcriptional regulatory networks.

Directed Differentiation of Human Induced Pluripotent Stem Cells to Heart Valve Cells.

BACKGROUND: A main obstacle in current valvular heart disease research is the lack of high-quality homogeneous functional heart valve cells. Human induced pluripotent stem cells (hiPSCs)-derived heart valve cells may help with this dilemma. However, there are no well-established protocols to induce hiPSCs to differentiate into functional heart valve cells, and the networks that mediate the differentiation have not been fully elucidated. METHODS: To generate heart valve cells from hiPSCs, we sequentially activated the Wnt, BMP4, VEGF (vascular endothelial growth factor), and NFATc1 signaling pathways using CHIR-99021, BMP4, VEGF-165, and forskolin, respectively. The transcriptional and functional similarity of hiPSC-derived heart valve cells compared with primary heart valve cells were characterized. Longitudinal single-cell RNA sequencing was used to uncover the trajectory, switch genes, pathways, and transcription factors of the differentiation. RESULTS: An efficient protocol was developed to induce hiPSCs to differentiate into functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells. After 6-day differentiation and CD144 magnetic bead sorting, ≈70% CD144+ cells and 30% CD144- cells were obtained. On the basis of single-cell RNA sequencing data, the CD144+ cells and CD144- cells were found to be highly similar to primary heart valve endothelial cells and primary heart valve interstitial cells in gene expression profile. Furthermore, CD144+ cells had the typical function of primary heart valve endothelial cells, including tube formation, uptake of low-density lipoprotein, generation of endothelial nitric oxide synthase, and response to shear stress. Meanwhile, CD144- cells could secret collagen and matrix metalloproteinases, and differentiate into osteogenic or adipogenic lineages like primary heart valve interstitial cells. Therefore, we identified CD144+ cells and CD144- cells as hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells, respectively. Using single-cell RNA sequencing analysis, we demonstrated that the trajectory of heart valve cell differentiation was consistent with embryonic valve development. We identified the main switch genes (NOTCH1, HEY1, and MEF2C), signaling pathways (TGF-ß, Wnt, and NOTCH), and transcription factors (MSX1, SP5, and MECOM) that mediated the differentiation. Finally, we found that hiPSC-derived valve interstitial-like cells might derive from hiPSC-derived valve endothelial-like cells undergoing endocardial-mesenchymal transition. CONCLUSIONS: In summary, this is the first study to report an efficient strategy to generate functional hiPSC-derived valve endothelial-like cells and hiPSC-derived valve interstitial-like cells from hiPSCs, as well as to elucidate the differentiation trajectory and transcriptional dynamics of hiPSCs differentiated into heart valve cells.

Interferon-γ inducible factor 16 (IFI16) restricts adeno-associated virus type 2 (AAV2) transduction in an immune-modulatory independent way.

We determined the transcription profile of adeno-associated virus type 2 (AAV2)-infected primary human fibroblasts. Subsequent analysis revealed that cells respond to AAV infection through changes in several significantly affected pathways, including cell cycle regulation, chromatin modulation, and innate immune responses. Various assays were performed to validate selected differentially expressed genes and to confirm not only the quality but also the robustness of the raw data. One of the genes upregulated in AAV2-infected cells was interferon-γ inducible factor 16 (IFI16). IFI16 is known as a multifunctional cytosolic and nuclear innate immune sensor for double-stranded as well as single-stranded DNA, exerting its effects through various mechanisms, such as interferon response, epigenetic modifications, or transcriptional regulation. IFI16 thereby constitutes a restriction factor for many different viruses among them, as shown here, AAV2 and thereof derived vectors. Indeed, the post-transcriptional silencing of IFI16 significantly increased AAV2 transduction efficiency, independent of the structure of the virus/vector genome. We also show that IFI16 exerts its inhibitory effect on AAV2 transduction in an immune-modulatory independent way by interfering with Sp1-dependent transactivation of wild-type AAV2 and AAV2 vector promoters. IMPORTANCE: Adeno-associated virus (AAV) vectors are among the most frequently used viral vectors for gene therapy. The lack of pathogenicity of the parental virus, the long-term persistence as episomes in non-proliferating cells, and the availability of a variety of AAV serotypes differing in their cellular tropism are advantageous features of this biological nanoparticle. To deepen our understanding of virus-host interactions, especially in terms of antiviral responses, we present here the first transcriptome analysis of AAV serotype 2 (AAV2)-infected human primary fibroblasts. Our findings indicate that interferon-γ inducible factor 16 acts as an antiviral factor in AAV2 infection and AAV2 vector-mediated cell transduction in an immune-modulatory independent way by interrupting the Sp1-dependent gene expression from viral or vector genomes.

UNMF: a unified nonnegative matrix factorization for multi-dimensional omics data.

Factor analysis, ranging from principal component analysis to nonnegative matrix factorization, represents a foremost approach in analyzing multi-dimensional data to extract valuable patterns, and is increasingly being applied in the context of multi-dimensional omics datasets represented in tensor form. However, traditional analytical methods are heavily dependent on the format and structure of the data itself, and if these change even slightly, the analyst must change their data analysis strategy and techniques and spend a considerable amount of time on data preprocessing. Additionally, many traditional methods cannot be applied as-is in the presence of missing values in the data. We present a new statistical framework, unified nonnegative matrix factorization (UNMF), for finding informative patterns in messy biological data sets. UNMF is designed for tidy data format and structure, making data analysis easier and simplifying the development of data analysis tools. UNMF can handle a wide range of data structures and formats, and works seamlessly with tensor data including missing observations and repeated measurements. The usefulness of UNMF is demonstrated through its application to several multi-dimensional omics data, offering user-friendly and unified features for analysis and integration. Its application holds great potential for the life science community. UNMF is implemented with R and is available from GitHub (https://github.com/abikoushi/moltenNMF).

Model-based evaluation of spatiotemporal data reduction methods with unknown ground truth through optimal visualization and interpretability metrics.

Optimizing and benchmarking data reduction methods for dynamic or spatial visualization and interpretation (DSVI) face challenges due to many factors, including data complexity, lack of ground truth, time-dependent metrics, dimensionality bias and different visual mappings of the same data. Current studies often focus on independent static visualization or interpretability metrics that require ground truth. To overcome this limitation, we propose the MIBCOVIS framework, a comprehensive and interpretable benchmarking and computational approach. MIBCOVIS enhances the visualization and interpretability of high-dimensional data without relying on ground truth by integrating five robust metrics, including a novel time-ordered Markov-based structural metric, into a semi-supervised hierarchical Bayesian model. The framework assesses method accuracy and considers interaction effects among metric features. We apply MIBCOVIS using linear and nonlinear dimensionality reduction methods to evaluate optimal DSVI for four distinct dynamic and spatial biological processes captured by three single-cell data modalities: CyTOF, scRNA-seq and CODEX. These data vary in complexity based on feature dimensionality, unknown cell types and dynamic or spatial differences. Unlike traditional single-summary score approaches, MIBCOVIS compares accuracy distributions across methods. Our findings underscore the joint evaluation of visualization and interpretability, rather than relying on separate metrics. We reveal that prioritizing average performance can obscure method feature performance. Additionally, we explore the impact of data complexity on visualization and interpretability. Specifically, we provide optimal parameters and features and recommend methods, like the optimized variational contractive autoencoder, for targeted DSVI for various data complexities. MIBCOVIS shows promise for evaluating dynamic single-cell atlases and spatiotemporal data reduction models.

Division-Independent Differentiation of Muscle Stem Cells During a Growth Stimulus.

Adult muscle stem cells (MuSCs) are known to replicate upon activation before differentiating and fusing to regenerate myofibers. It is unclear whether MuSC differentiation is intrinsically linked to cell division, which has implications for stem cell population maintenance. We use single-cell RNA-sequencing to identify transcriptionally diverse subpopulations of MuSCs after 5 days of a growth stimulus in adult muscle. Trajectory inference in combination with a novel mouse model for tracking MuSC-derived myonuclei and in vivo labeling of DNA replication revealed an MuSC population that exhibited division-independent differentiation and fusion. These findings demonstrate that in response to a growth stimulus in the presence of intact myofibers, MuSC division is not obligatory.

Single-Cell Transcriptome Reveals Potential Mechanisms for Coronary Artery Lesions in Kawasaki Disease.

BACKGROUND: Coronary artery lesions (CALs) are the most common and major complication of Kawasaki disease (KD) in developed countries. However, the underlying immunologic mechanisms of CAL development in KD remain unclear. METHODS: Here, we conducted single-cell transcriptome analyses of 212 210 peripheral blood mononuclear cells collected from a cross-sectional cohort of 16 children, including 4 patients with KD with CALs, 5 patients with KD without CALs, 4 healthy controls, and 3 febrile controls. RESULTS: KD altered the proportion of peripheral blood mononuclear cells, including an increasing trend in inflammatory cells (megakaryocytes and monocytes) and a decreasing trend in lymphocytes (eg, CD4+ T, CD8+ T, mucosal-associated invariant T, natural killer, and γδ T cells), highlighting the potential presence of lymphopenia phenomenon in KD. Our data indicated the presence of inflammatory cytokine storm in patients with KD with CALs, caused by systemic upregulation of TNFSF13B (tumor necrosis factor superfamily member 13b), CXCL16 (C-X-C motif chemokine ligand 16), TNFSF10 (tumor necrosis factor superfamily member 10), and IL1RN (interleukin 1 receptor antagonist), mainly produced by monocytes (especially for the Mono_CD14-CD16 cluster) and megakaryocytes. We also found that myeloid cells of patients with KD, particularly in those with CALs, might play a role in vascular injury (eg, increased MMP [matrix metalloproteinase] 9, MMP17, and MMP25) and immune cell recruitment. The immune landscape of patients with KD with CALs was featured by lower exhaustion levels in natural killer cells, a high cytotoxic state in the CD8_Pro cluster, and activation of the complement system in monocytes. Additionally, the activation of B cells was more pronounced in the early stage of KD. CONCLUSIONS: Collectively, this study provides a comprehensive understanding of the roles of various immune cells and inflammatory cytokine storms in the development of CALs in KD and offers a valuable resource for identifying novel therapeutic targets for patients with KD with CALs.

Unravelling the immune signature of herpes zoster: Insights into pathophysiology and the HLA risk profile.

The varicella-zoster virus (VZV) infects over 95% of the population. VZV reactivation causes herpes zoster (HZ), known as shingles, primarily affecting the elderly and immunocompromised individuals. However, HZ can also occur in otherwise healthy individuals. We analyzed the immune signature and risk profile in HZ patients using a genome-wide association study across different UK Biobank HZ cohorts. Additionally, we conducted one of the largest HZ HLA association studies to date, coupled with transcriptomic analysis of pathways underlying HZ susceptibility. Our findings highlight the significance of the MHC locus for HZ development, identifying five protective and four risk HLA alleles. This demonstrates that HZ susceptibility is largely governed by variations in the MHC. Furthermore, functional analyses revealed the upregulation of type I interferon and adaptive immune responses. These findings provide fresh molecular insights into the pathophysiology and the activation of innate and adaptive immune responses triggered by symptomatic VZV reactivation.

Targeting Epsins to Inhibit Fibroblast Growth Factor Signaling While Potentiating Transforming Growth Factor-ß Signaling Constrains Endothelial-to-Mesenchymal Transition in Atherosclerosis.

BACKGROUND: Epsin endocytic adaptor proteins are implicated in the progression of atherosclerosis; however, the underlying molecular mechanisms have not yet been fully defined. In this study, we determined how epsins enhance endothelial-to-mesenchymal transition (EndoMT) in atherosclerosis and assessed the efficacy of a therapeutic peptide in a preclinical model of this disease. METHODS: Using single-cell RNA sequencing combined with molecular, cellular, and biochemical analyses, we investigated the role of epsins in stimulating EndoMT using knockout in Apoe-/- and lineage tracing/proprotein convertase subtilisin/kexin type 9 serine protease mutant viral-induced atherosclerotic mouse models. The therapeutic efficacy of a synthetic peptide targeting atherosclerotic plaques was then assessed in Apoe-/- mice. RESULTS: Single-cell RNA sequencing and lineage tracing revealed that epsins 1 and 2 promote EndoMT and that the loss of endothelial epsins inhibits EndoMT marker expression and transforming growth factor-ß signaling in vitro and in atherosclerotic mice, which is associated with smaller lesions in the Apoe-/- mouse model. Mechanistically, the loss of endothelial cell epsins results in increased fibroblast growth factor receptor-1 expression, which inhibits transforming growth factor-ß signaling and EndoMT. Epsins directly bind ubiquitinated fibroblast growth factor receptor-1 through their ubiquitin-interacting motif, which results in endocytosis and degradation of this receptor complex. Consequently, administration of a synthetic ubiquitin-interacting motif-containing peptide atheroma ubiquitin-interacting motif peptide inhibitor significantly attenuates EndoMT and progression of atherosclerosis. CONCLUSIONS: We conclude that epsins potentiate EndoMT during atherogenesis by increasing transforming growth factor-ß signaling through fibroblast growth factor receptor-1 internalization and degradation. Inhibition of EndoMT by reducing epsin-fibroblast growth factor receptor-1 interaction with a therapeutic peptide may represent a novel treatment strategy for atherosclerosis.

Genome-wide identification and expression characterization of the GH3 gene family of tea plant (Camellia sinensis).

To comprehensively understand the characteristics of the GH3 gene family in tea plants (Camellia sinensis), we identified 17 CsGH3 genes and analyzed their physicochemical properties, phylogenetic relationships, gene structures, promoters, and expression patterns in different tissues. The study showed that the 17 CsGH3 genes are distributed on 9 chromosomes, and based on evolutionary analysis, the CsGH3 members were divided into three subgroups. Gene duplication analysis revealed that segmental duplications have a significant impact on the amplification of CsGH3 genes. In addition, we identified and classified cis-elements in the CsGH3 gene promoters and detected elements related to plant hormone responses and non-biotic stress responses. Through expression pattern analysis, we observed tissue-specific expression of CsGH3.3 and CsGH3.10 in flower buds and roots. Moreover, based on predictive analysis of upstream regulatory transcription factors of CsGH3, we identified the potential transcriptional regulatory role of gibberellin response factor CsDELLA in CsGH3.14 and CsGH3.15. In this study, we found that CsGH3 genes are involved in a wide range of activities, such as growth and development, stress response, and transcription. This is the first report on CsGH3 genes and their potential roles in tea plants. In conclusion, these results provide a theoretical basis for elucidating the role of GH3 genes in the development of perennial woody plants and offer new insights into the synergistic effects of multiple hormones on plant growth and development in tea plants.

Prognostic significance of circulating tumor cells and tumor related transcripts in small cell lung cancer: A step further to clinical implementation.

Small-cell lung cancer (SCLC) is a fatal disease with limited treatment options. Circulating tumor cells (CTCs) in liquid biopsy samples may serve as predictive and prognostic biomarkers; but the analysis of CTCs is still challenging. By using microfluidic or density gradient CTC enrichment in combination with immunofluorescent (IF) staining or qPCR of CTC-related transcripts, we achieved a 60.8% to 88.0% positivity in SCLC blood samples. Epithelial and neuroendocrine transcripts including the druggable target DLL3 were associated with shorter overall survival (OS), indicating the clinical value of these markers in terms of differential diagnosis and treatment decisions. High CTC counts and the presence of CTC duplets detected by IF staining were prognostic for OS, and thus may serve as indicators of disease progression or therapy failure. In patient samples with high CTC load detected by IF staining, a concordance of the transcripts positivity in circulating free plasma RNA and CTCs was observed. Our data emphasize the role of CTCs and CTC-related transcripts and underline the clinical value of liquid biopsy analysis in SCLC.

Implementable assay for monitoring minimum residual disease after radical treatment for colorectal cancer.

Considering the cost and invasiveness of monitoring postoperative minimal residual disease (MRD) of colorectal cancer (CRC) after adjuvant chemoradiotherapy (ACT), we developed a favorable approach based on methylated circulating tumor DNA to detect MRD after radical resection. Analyzing the public database, we identified the methylated promoter regions of the genes FGD5, GPC6, and MSC. Using digital polymerase chain reaction (dPCR), we termed the "amplicon of methylated sites using a specific enzyme" assay as "AMUSE." We examined 180 and 114 pre- and postoperative serial plasma samples from 28 recurrent and 19 recurrence-free pathological stage III CRC patients, respectively. The results showed 22 AMUSE-positive of 28 recurrent patients (sensitivity, 78.6%) and 17 AMUSE-negative of 19 recurrence-free patients (specificity, 89.5%). AMUSE predicted recurrence 208 days before conventional diagnosis using radiological imaging. Regarding ACT evaluation by the reactive response, 19 AMUSE-positive patients during their second or third blood samples showed a significantly poorer prognosis than the other patients (p = 9E-04). The AMUSE assay stratified four groups by the altered patterns of tumor burden postoperatively. Interestingly, only 34.8% of cases tested AMUSE-negative during ACT treatment, indicating eligibility for ACT. The AMUSE assay addresses the clinical need for accurate MRD monitoring with universal applicability, minimal invasiveness, and cost-effectiveness, thereby enabling the timely detection of recurrences. This assay can effectively evaluate the efficacy of ACT in patients with stage III CRC following curative resection. Our study strongly recommends reevaluating the clinical application of ACT using the AMUSE assay.

De novo transcriptome assembly and global analysis of differential gene expression of aphid tolerant wild mustard Rorippa indica (L.) Hiern infested by mustard aphid Lipaphis Erysimi (L.) Kaltenbach.

Rapeseed-mustard, the oleiferous Brassica species are important oilseed crops cultivated all over the globe. Mustard aphid Lipaphis erysimi (L.) Kaltenbach is a major threat to the cultivation of rapeseed-mustard. Wild mustard Rorippa indica (L.) Hiern shows tolerance to mustard aphids as a nonhost and hence is an important source for the bioprospecting of potential resistance genes and defense measures to manage mustard aphids sustainably. We performed mRNA sequencing of the R. indica plant uninfested and infested by the mustard aphids, harvested at 24 hours post-infestation. Following quality control, the high-quality reads were subjected to de novo assembly of the transcriptome. As there is no genomic information available for this potential wild plant, the raw reads will be useful for further bioinformatics analysis and the sequence information of the assembled transcripts will be helpful to design primers for the characterization of specific gene sequences. In this study, we also used the generated resource to comprehensively analyse the global profile of differential gene expression in R. indica in response to infestation by mustard aphids. The functional enrichment analysis of the differentially expressed genes reveals a significant immune response and suggests the possibility of chitin-induced defense signaling.

Gene expression analysis of drought tolerance and cuticular wax biosynthesis in diploid and tetraploid induced wallflowers.

Whole-genome doubling leads to cell reprogramming, upregulation of stress genes, and establishment of new pathways of drought stress responses in plants. This study investigated the molecular mechanisms of drought tolerance and cuticular wax characteristics in diploid and tetraploid-induced Erysimum cheiri. According to real-time PCR analysis, tetraploid induced wallflowers exhibited increased expression of several genes encoding transcription factors (TFs), including AREB1 and AREB3; the stress response genes RD29A and ERD1 under drought stress conditions. Furthermore, two cuticular wax biosynthetic pathway genes, CER1 and SHN1, were upregulated in tetraploid plants under drought conditions. Leaf morphological studies revealed that tetraploid leaves were covered with unique cuticular wax crystalloids, which produced a white fluffy appearance, while the diploid leaves were green and smooth. The greater content of epicuticular wax in tetraploid leaves than in diploid leaves can explain the decrease in cuticle permeability as well as the decrease in water loss and improvement in drought tolerance in wallflowers. GC‒MS analysis revealed that the wax components included alkanes, alcohols, aldehydes, and fatty acids. The most abundant wax compound in this plant was alkanes (50%), the most predominant of which was C29. The relative abundance of these compounds increased significantly in tetraploid plants under drought stress conditions. These findings revealed that tetraploid-induced wallflowers presented upregulation of multiple drought-related and wax biosynthesis genes; therefore, polyploidization has proved useful for improving plant drought tolerance.

Comprehensive genome-wide analysis of wheat xylanase inhibitor protein (XIP) genes: unveiling their role in Fusarium head blight resistance and plant immune mechanisms.

In this comprehensive genome-wide study, we identified and classified 83 Xylanase Inhibitor Protein (XIP) genes in wheat, grouped into five distinct categories, to enhance understanding of wheat's resistance to Fusarium head blight (FHB), a significant fungal threat to global wheat production. Our analysis reveals the unique distribution of XIP genes across wheat chromosomes, particularly at terminal regions, suggesting their role in the evolutionary expansion of the gene family. Several XIP genes lack signal peptides, indicating potential alternative secretion pathways that could be pivotal in plant defense against FHB. The study also uncovers the sequence homology between XIPs and chitinases, hinting at a functional diversification within the XIP gene family. Additionally, the research explores the association of XIP genes with plant immune mechanisms, particularly their linkage with plant hormone signaling pathways like abscisic acid and jasmonic acid. XIP-7A3, in particular, demonstrates a significant increase in expression upon FHB infection, highlighting its potential as a key candidate gene for enhancing wheat's resistance to this disease. This research not only enriches our understanding of the XIP gene family in wheat but also provides a foundation for future investigations into their role in developing FHB-resistant wheat cultivars. The findings offer significant implications for wheat genomics and breeding, contributing to the development of more resilient crops against fungal diseases.

Single-Cell RNA Sequencing of Coronary Perivascular Adipose Tissue From End-Stage Heart Failure Patients Identifies SPP1+ Macrophage Subpopulation as a Target for Alleviating Fibrosis.

BACKGROUND: Perivascular adipose tissue (PVAT) is vital for vascular homeostasis, and PVAT dysfunction is associated with increased atherosclerotic plaque burden. But the mechanisms underlining coronary PVAT dysfunction in coronary atherosclerosis remain elusive. METHODS: We performed single-cell RNA sequencing of the stromal vascular fraction of coronary PVAT from 3 groups of heart transplant recipients with end-stage heart failure, including 3 patients with nonobstructive coronary atherosclerosis, 3 patients with obstructive coronary artery atherosclerosis, and 4 nonatherosclerosis control subjects. Bioinformatics was used to annotate the cellular populations, depict the cellular developmental trajectories and interactions, and explore the differences among 3 groups of coronary PVAT at the cellular and molecular levels. Pathological staining, quantitative real-time polymerase chain reaction, and in vitro studies were performed to validate the key findings. RESULTS: Ten cell types were identified among 67 936 cells from human coronary PVAT. Several cellular subpopulations, including SPP1+ (secreted phosphoprotein 1) macrophages and profibrotic fibroadipogenic progenitor cells, were accumulated in PVAT surrounding atherosclerotic coronary arteries compared with nonatherosclerosis coronary arteries. The fibrosis percentage was increased in PVAT surrounding atherosclerotic coronary arteries, and it was positively associated with the grade of coronary artery stenosis. Cellular interaction analysis suggested OPN (osteopontin) secreted by SPP1+ macrophages interacted with CD44 (cluster of differentiation 44)/integrin on fibroadipogenic progenitor cells. Strikingly, correlation analyses uncovered that higher level of SPP1 in PVAT correlates with a more severe fibrosis degree and a higher coronary stenosis grade. In vitro studies showed that conditioned medium from atherosclerotic coronary PVAT promoted the migration and proliferation of fibroadipogenic progenitor cells, while such effect was prevented by blocking CD44 or integrin. CONCLUSIONS: SPP1+ macrophages accumulated in the PVAT surrounding atherosclerotic coronary arteries, and they promoted the migration and proliferation of fibroadipogenic progenitor cells via OPN-CD44/integrin interaction and thus aggravated the fibrosis of coronary PVAT, which was positively correlated to the coronary stenosis burden. Therefore, SPP1+ macrophages in coronary PVAT may participate in the progression of coronary atherosclerosis.

Single-Molecule Spatial Transcriptomics of Human Thoracic Aortic Aneurysms Uncovers Calcification-Related CARTPT-Expressing Smooth Muscle Cells.

BACKGROUND: Although single-cell RNA-sequencing is commonly applied to dissect the heterogeneity in human tissues, it involves the preparation of single-cell suspensions via cell dissociation, causing loss of spatial information. In this study, we employed high-resolution single-cell transcriptome imaging to reveal rare smooth muscle cell (SMC) types in human thoracic aortic aneurysm (TAA) tissue samples. METHODS: Single-molecule spatial distribution of transcripts from 140 genes was analyzed in fresh-frozen human TAA samples with region and sex-matched controls. In vitro studies and tissue staining were performed to examine human CART prepropeptide (CARTPT) regulation and function. RESULTS: We captured thousands of cells per sample including a spatially distinct CARTPT-expressing SMC subtype enriched in male TAA samples. Immunoassays confirmed human CART (cocaine- and amphetamine-regulated transcript) protein enrichment in male TAA tissue and truncated CARTPT secretion into cell culture medium. Oxidized low-density lipoprotein, a cardiovascular risk factor, induced CARTPT expression, whereas CARTPT overexpression in human aortic SMCs increased the expression of key osteochondrogenic transcription factors and reduced contractile gene expression. Recombinant human CART treatment of human SMCs further confirmed this phenotype. Alizarin red staining revealed calcium deposition in male TAA samples showing similar localization with human CART staining. CONCLUSIONS: Here, we demonstrate the feasibility of single-molecule imaging in uncovering rare SMC subtypes in the diseased human aorta, a difficult tissue to dissociate. We identified a spatially distinct CARTPT-expressing SMC subtype enriched in male human TAA samples. Our functional studies suggest that human CART promotes osteochondrogenic switch of aortic SMCs, potentially leading to medial calcification of the thoracic aorta.

From correlation to causation using directed topological overlap matrix: Applications in genomics.

Most causal discovery tools assume the local causal Markov condition. However, the theoretical assumptions that underlie the local causal Markov condition are often not met in practice. This is especially marked in genomics, where the unwanted presence of measurement errors, averaging effects, and feedback loops significantly undermine the legitimacy of the local causal Markov condition. Furthermore, these causal discovery algorithms require very large samples, orders above what is often available. In this paper, relaxing the local causal Markov condition and using Reichenbach's common cause principle instead, we present a more flexible approach to causal discovery, the directed topological overlap matrix (DTOM). DTOM is robust w.r.t. the presence of measurement errors, averaging effects, feedback loops, and is significantly more sample efficient. We study the utility of DTOM for discovering causal relations in biological data using three real gene expression data-sets. We first examine if DTOM can help distinguish the Myostatin mutation in the Piedmontese cattle by contrasting the muscle transcriptomes of the Piedmontese and Wagyu crosses: the Myostatin mutation is the cause of the double-muscling the Piedmontese cattle are famous for. We then consider a large-scale gene deletion study in yeast. We show that DTOM allows us to distinguish the deleted gene in a sample knowing only the set of differentially expressed genes in that sample. We then examine the progression of Alzheimer's disease (AD) under the lens of DTOM. The genes implicated as having a causal role in the progression of AD by our DTOM analysis were significantly enriched in cellular components that had been repeatedly implicated in the progression of AD.

Genome-Wide Identification, Evolution, and Expression Analysis of the DIR Gene Family in Schima superba.

Schima superba, commonly known as the Chinese guger tree, is highly adaptable and tolerant of poor soil conditions. It is one of the primary species forming the evergreen broad-leaved forests in southern China. Dirigent proteins (DIRs) play crucial roles in the synthesis of plant lignin and lignans, secondary metabolism, and response to adversity stress. However, research on the DIR gene family in S. superba is currently limited. This study identified 24 SsDIR genes, categorizing them into three subfamilies. These genes are unevenly distributed across 13 chromosomes, with 83% being intronless. Collinearity analysis indicated that tandem duplication played a more significant role in the expansion of the gene family compared to segmental duplication. Additionally, we analyzed the expression patterns of SsDIRs in different tissues of S. superba. The SsDIR genes exhibited distinct expression patterns across various tissues, with most being specifically expressed in the roots. Further screening identified SsDIR genes that may regulate drought stress, with many showing differential expression under drought stress conditions. In the promoter regions of SsDIRs, various cis-regulatory elements involved in developmental regulation, hormone response, and stress response were identified, which may be closely related to their diverse regulatory functions. This study will contribute to the further functional identification of SsDIR genes, providing insights into the biosynthetic pathways of lignin and lignans and the mechanisms of plant stress resistance.