Spliceosome-targeted therapies trigger an antiviral immune response in triple-negative breast cancer.

Many oncogenic insults deregulate RNA splicing, often leading to hypersensitivity of tumors to spliceosome-targeted therapies (STTs). However, the mechanisms by which STTs selectively kill cancers remain largely unknown. Herein, we discover that mis-spliced RNA itself is a molecular trigger for tumor killing through viral mimicry. In MYC-driven triple-negative breast cancer, STTs cause widespread cytoplasmic accumulation of mis-spliced mRNAs, many of which form double-stranded structures. Double-stranded RNA (dsRNA)-binding proteins recognize these endogenous dsRNAs, triggering antiviral signaling and extrinsic apoptosis. In immune-competent models of breast cancer, STTs cause tumor cell-intrinsic antiviral signaling, downstream adaptive immune signaling, and tumor cell death. Furthermore, RNA mis-splicing in human breast cancers correlates with innate and adaptive immune signatures, especially in MYC-amplified tumors that are typically immune cold. These findings indicate that dsRNA-sensing pathways respond to global aberrations of RNA splicing in cancer and provoke the hypothesis that STTs may provide unexplored strategies to activate anti-tumor immune pathways.

Gene expression imputation and cell-type deconvolution in human brain with spatiotemporal precision and its implications for brain-related disorders.

As the most complex organ of the human body, the brain is composed of diverse regions, each consisting of distinct cell types and their respective cellular interactions. Human brain development involves a finely tuned cascade of interactive events. These include spatiotemporal gene expression changes and dynamic alterations in cell-type composition. However, our understanding of this process is still largely incomplete owing to the difficulty of brain spatiotemporal transcriptome collection. In this study, we developed a tensor-based approach to impute gene expression on a transcriptome-wide level. After rigorous computational benchmarking, we applied our approach to infer missing data points in the widely used BrainSpan resource and completed the entire grid of spatiotemporal transcriptomics. Next, we conducted deconvolutional analyses to comprehensively characterize major cell-type dynamics across the entire BrainSpan resource to estimate the cellular temporal changes and distinct neocortical areas across development. Moreover, integration of these results with GWAS summary statistics for 13 brain-associated traits revealed multiple novel trait-cell-type associations and trait-spatiotemporal relationships. In summary, our imputed BrainSpan transcriptomic data provide a valuable resource for the research community and our findings help further studies of the transcriptional and cellular dynamics of the human brain and related diseases.

Electrochemical Skeletal Indole Editing via Nitrogen Atom Insertion by Sustainable Oxygen Reduction Reaction.

Skeletal molecular editing gained considerable recent momentum and emerged as a uniquely powerful tool for late-stage diversifications. Thus far, superstoichiometric amounts of costly hypervalent iodine(III) reagents were largely required for skeletal indole editing. In contrast, we herein show that electricity enables sustainable nitrogen atom insertion reactions to give bio-relevant quinazoline scaffolds without stoichiometric chemical redox-waste product. The transition metal-free electro-editing was enabled by the oxygen reduction reaction (ORR) and proved robust on scale, while tolerating a variety of valuable functional groups.

Single-Cell RNA Sequencing Uncovers Paracrine Functions of the Epicardial-Derived Cells in Arrhythmogenic Cardiomyopathy.

BACKGROUND: Arrhythmogenic cardiomyopathy (ACM) manifests with sudden death, arrhythmias, heart failure, apoptosis, and myocardial fibro-adipogenesis. The phenotype typically starts at the epicardium and advances transmurally. Mutations in genes encoding desmosome proteins, including DSP (desmoplakin), are major causes of ACM. METHODS: To delineate contributions of the epicardium to the pathogenesis of ACM, the Dsp allele was conditionally deleted in the epicardial cells in mice upon expression of tamoxifen-inducible Cre from the Wt1 locus. Wild type (WT) and Wt1-CreERT2:DspW/F were crossed to Rosa26mT/mG (R26mT/mG) dual reporter mice to tag the epicardial-derived cells with the EGFP (enhanced green fluorescent protein) reporter protein. Tagged epicardial-derived cells from adult Wt1-CreERT2:R26mT/mG and Wt1-CreERT2: R26mT/mG:DspW/F mouse hearts were isolated by fluorescence-activated cell staining and sequenced by single-cell RNA sequencing. RESULTS: WT1 (Wilms tumor 1) expression was progressively restricted postnatally and was exclusive to the epicardium by postnatal day 21. Expression of Dsp was reduced in the epicardial cells but not in cardiac myocytes in the Wt1-CreERT2:DspW/F mice. The Wt1-CreERT2:DspW/F mice exhibited premature death, cardiac dysfunction, arrhythmias, myocardial fibro-adipogenesis, and apoptosis. Single-cell RNA sequencing of ≈18 000 EGFP-tagged epicardial-derived cells identified genotype-independent clusters of endothelial cells, fibroblasts, epithelial cells, and a very small cluster of cardiac myocytes, which were confirmed on coimmunofluorescence staining of the myocardial sections. Differentially expressed genes between the paired clusters in the 2 genotypes predicted activation of the inflammatory and mitotic pathways-including the TGFß1 (transforming growth factor ß1) and fibroblast growth factors-in the epicardial-derived fibroblast and epithelial clusters, but predicted their suppression in the endothelial cell cluster. The findings were corroborated by analysis of gene expression in the pooled RNA-sequencing data, which identified predominant dysregulation of genes involved in epithelial-mesenchymal transition, and dysregulation of 146 genes encoding the secreted proteins (secretome), including genes in the TGFß1 pathway. Activation of the TGFß1 and its colocalization with fibrosis in the Wt1-CreERT2:R26mT/mG:DspW/F mouse heart was validated by complementary methods. CONCLUSIONS: Epicardial-derived cardiac fibroblasts and epithelial cells express paracrine factors, including TGFß1 and fibroblast growth factors, which mediate epithelial-mesenchymal transition, and contribute to the pathogenesis of myocardial fibrosis, apoptosis, arrhythmias, and cardiac dysfunction in a mouse model of ACM. The findings uncover contributions of the epicardial-derived cells to the pathogenesis of ACM.

Longitudinal RNA-Seq Analysis of the Repeatability of Gene Expression and Splicing in Human Platelets Identifies a Platelet SELP Splice QTL.

RATIONALE: Longitudinal studies are required to distinguish within versus between-individual variation and repeatability of gene expression. They are uniquely positioned to decipher genetic signal from environmental noise, with potential application to gene variant and expression studies. However, longitudinal analyses of gene expression in healthy individuals-especially with regards to alternative splicing-are lacking for most primary cell types, including platelets. OBJECTIVE: To assess repeatability of gene expression and splicing in platelets and use repeatability to identify novel platelet expression quantitative trait loci (QTLs) and splice QTLs. METHODS AND RESULTS: We sequenced the transcriptome of platelets isolated repeatedly up to 4 years from healthy individuals. We examined within and between individual variation and repeatability of platelet RNA expression and exon skipping, a readily measured alternative splicing event. We find that platelet gene expression is generally stable between and within-individuals over time-with the exception of a subset of genes enriched for the inflammation gene ontology. We show an enrichment among repeatable genes for associations with heritable traits, including known and novel platelet expression QTLs. Several exon skipping events were also highly repeatable, suggesting heritable patterns of splicing in platelets. One of the most repeatable was exon 14 skipping of SELP. Accordingly, we identify rs6128 as a platelet splice QTL and define an rs6128-dependent association between SELP exon 14 skipping and race. In vitro experiments demonstrate that this single nucleotide variant directly affects exon 14 skipping and changes the ratio of transmembrane versus soluble P-selectin protein production. CONCLUSIONS: We conclude that the platelet transcriptome is generally stable over 4 years. We demonstrate the use of repeatability of gene expression and splicing to identify novel platelet expression QTLs and splice QTLs. rs6128 is a platelet splice QTL that alters SELP exon 14 skipping and soluble versus transmembrane P-selectin protein production.

The spliceosome is a therapeutic vulnerability in MYC-driven cancer.

MYC (also known as c-MYC) overexpression or hyperactivation is one of the most common drivers of human cancer. Despite intensive study, the MYC oncogene remains recalcitrant to therapeutic inhibition. MYC is a transcription factor, and many of its pro-tumorigenic functions have been attributed to its ability to regulate gene expression programs. Notably, oncogenic MYC activation has also been shown to increase total RNA and protein production in many tissue and disease contexts. While such increases in RNA and protein production may endow cancer cells with pro-tumour hallmarks, this increase in synthesis may also generate new or heightened burden on MYC-driven cancer cells to process these macromolecules properly. Here we discover that the spliceosome is a new target of oncogenic stress in MYC-driven cancers. We identify BUD31 as a MYC-synthetic lethal gene in human mammary epithelial cells, and demonstrate that BUD31 is a component of the core spliceosome required for its assembly and catalytic activity. Core spliceosomal factors (such as SF3B1 and U2AF1) associated with BUD31 are also required to tolerate oncogenic MYC. Notably, MYC hyperactivation induces an increase in total precursor messenger RNA synthesis, suggesting an increased burden on the core spliceosome to process pre-mRNA. In contrast to normal cells, partial inhibition of the spliceosome in MYC-hyperactivated cells leads to global intron retention, widespread defects in pre-mRNA maturation, and deregulation of many essential cell processes. Notably, genetic or pharmacological inhibition of the spliceosome in vivo impairs survival, tumorigenicity and metastatic proclivity of MYC-dependent breast cancers. Collectively, these data suggest that oncogenic MYC confers a collateral stress on splicing, and that components of the spliceosome may be therapeutic entry points for aggressive MYC-driven cancers.

GATA2/3-TFAP2A/C transcription factor network couples human pluripotent stem cell differentiation to trophectoderm with repression of pluripotency.

To elucidate the molecular basis of BMP4-induced differentiation of human pluripotent stem cells (PSCs) toward progeny with trophectoderm characteristics, we produced transcriptome, epigenome H3K4me3, H3K27me3, and CpG methylation maps of trophoblast progenitors, purified using the surface marker APA. We combined them with the temporally resolved transcriptome of the preprogenitor phase and of single APA+ cells. This revealed a circuit of bivalent TFAP2A, TFAP2C, GATA2, and GATA3 transcription factors, coined collectively the "trophectoderm four" (TEtra), which are also present in human trophectoderm in vivo. At the onset of differentiation, the TEtra factors occupy multiple sites in epigenetically inactive placental genes and in OCT4 Functional manipulation of GATA3 and TFAP2A indicated that they directly couple trophoblast-specific gene induction with suppression of pluripotency. In accordance, knocking down GATA3 in primate embryos resulted in a failure to form trophectoderm. The discovery of the TEtra circuit indicates how trophectoderm commitment is regulated in human embryogenesis.

A paradox of transcriptional and functional innate interferon responses of human intestinal enteroids to enteric virus infection.

The intestinal epithelium can limit enteric pathogens by producing antiviral cytokines, such as IFNs. Type I IFN (IFN-α/ß) and type III IFN (IFN-λ) function at the epithelial level, and their respective efficacies depend on the specific pathogen and site of infection. However, the roles of type I and type III IFN in restricting human enteric viruses are poorly characterized as a result of the difficulties in cultivating these viruses in vitro and directly obtaining control and infected small intestinal human tissue. We infected nontransformed human intestinal enteroid cultures from multiple individuals with human rotavirus (HRV) and assessed the host epithelial response by using RNA-sequencing and functional assays. The dominant transcriptional pathway induced by HRV infection is a type III IFN-regulated response. Early after HRV infection, low levels of type III IFN protein activate IFN-stimulated genes. However, this endogenous response does not restrict HRV replication because replication-competent HRV antagonizes the type III IFN response at pre- and posttranscriptional levels. In contrast, exogenous IFN treatment restricts HRV replication, with type I IFN being more potent than type III IFN, suggesting that extraepithelial sources of type I IFN may be the critical IFN for limiting enteric virus replication in the human intestine.

The Human Lung Cell Atlas: A High-Resolution Reference Map of the Human Lung in Health and Disease.

Lung disease accounts for every sixth death globally. Profiling the molecular state of all lung cell types in health and disease is currently revolutionizing the identification of disease mechanisms and will aid the design of novel diagnostic and personalized therapeutic regimens. Recent progress in high-throughput techniques for single-cell genomic and transcriptomic analyses has opened up new possibilities to study individual cells within a tissue, classify these into cell types, and characterize variations in their molecular profiles as a function of genetics, environment, cell-cell interactions, developmental processes, aging, or disease. Integration of these cell state definitions with spatial information allows the in-depth molecular description of cellular neighborhoods and tissue microenvironments, including the tissue resident structural and immune cells, the tissue matrix, and the microbiome. The Human Cell Atlas consortium aims to characterize all cells in the healthy human body and has prioritized lung tissue as one of the flagship projects. Here, we present the rationale, the approach, and the expected impact of a Human Lung Cell Atlas.

Integrative Multi-omic Analysis of Human Platelet eQTLs Reveals Alternative Start Site in Mitofusin 2.

Platelets play a central role in ischemic cardiovascular events. Cardiovascular disease (CVD) is a major cause of death worldwide. Numerous genome-wide association studies (GWASs) have identified loci associated with CVD risk. However, our understanding of how these variants contribute to disease is limited. Using data from the platelet RNA and expression 1 (PRAX1) study, we analyzed cis expression quantitative trait loci (eQTLs) in platelets from 154 normal human subjects. We confirmed these results in silico by performing allele-specific expression (ASE) analysis, which demonstrated that the allelic directionality of eQTLs and ASE patterns correlate significantly. Comparison of platelet eQTLs with data from the Genotype-Tissue Expression (GTEx) project revealed that a number of platelet eQTLs are platelet specific and that platelet eQTL peaks localize to the gene body at a higher rate than eQTLs from other tissues. Upon integration with data from previously published GWASs, we found that the trait-associated variant rs1474868 coincides with the eQTL peak for mitofusin 2 (MFN2). Additional experimental and computational analyses revealed that this eQTL is linked to an unannotated alternate MFN2 start site preferentially expressed in platelets. Integration of phenotype data from the PRAX1 study showed that MFN2 expression levels were significantly associated with platelet count. This study links the variant rs1474868 to a platelet-specific regulatory role for MFN2 and demonstrates the utility of integrating multi-omic data with eQTL analysis in disease-relevant tissues for interpreting GWAS results.

Effects of genetic variation in protease activated receptor 4 after an acute coronary syndrome: Analysis from the TRACER trial.

Variation in platelet response to thrombin may affect the safety and efficacy of PAR antagonism. The Thr120 variant of the common single nucleotide polymorphism (SNP) rs773902 in the protease-activated receptor (PAR) 4 gene is associated with higher platelet aggregation compared to the Ala120 variant. We investigated the relationship between the rs773902 SNP with major bleeding and ischemic events, safety, and efficacy of PAR1 inhibition in 6177 NSTE ACS patients in the TRACER trial. There was a lower rate of GUSTO moderate/severe bleeding in patients with the Thr120 variant. The difference was driven by a lower rate in the smaller homozygous group (recessive model, HR 0.13 [0.02-0.92] P = 0.042). No significant differences were observed in the ischemic outcomes. The excess in bleeding observed with PAR1 inhibition was attenuated in patients with the Thr120 variant, but the interactions were not statistically significant. In summary, lower major bleeding rates were observed in the overall TRACER cohort with the hyperreactive PAR4 Thr120 variant. The increase in bleeding with vorapaxar was attenuated with the Thr120 variant, but we could not demonstrate an interaction with PAR1 inhibition. These findings warrant further exploration, including those of African ancestry where the A allele (Thr120) frequency is ~65%.

Common variants in the human platelet PAR4 thrombin receptor alter platelet function and differ by race.

Human platelets express 2 thrombin receptors: protease-activated receptor (PAR)-1 and PAR4. Recently, we reported 3.7-fold increased PAR4-mediated aggregation kinetics in platelets from black subjects compared with white subjects. We now show that platelets from blacks (n = 70) express 14% more PAR4 protein than those from whites (n = 84), but this difference is not associated with platelet PAR4 function. Quantitative trait locus analysis identified 3 common single nucleotide polymorphisms in the PAR4 gene (F2RL3) associated with PAR4-induced platelet aggregation. Among these single nucleotide polymorphisms, rs773902 determines whether residue 120 in transmembrane domain 2 is an alanine (Ala) or threonine (Thr). Compared with the Ala120 variant, Thr120 was more common in black subjects than in white subjects (63% vs 19%), was associated with higher PAR4-induced human platelet aggregation and Ca2+ flux, and generated greater inositol 1,4,5-triphosphate in transfected cells. A second, less frequent F2RL3 variant, Phe296Val, was only observed in blacks and abolished the enhanced PAR4-induced platelet aggregation and 1,4,5-triphosphate generation associated with PAR4-Thr120. PAR4 genotype did not affect vorapaxar inhibition of platelet PAR1 function, but a strong pharmacogenetic effect was observed with the PAR4-specific antagonist YD-3 [1-benzyl-3(ethoxycarbonylphenyl)-indazole]. These findings may have an important pharmacogenetic effect on the development of new PAR antagonists.

Human platelet microRNA-mRNA networks associated with age and gender revealed by integrated plateletomics.

There is little data considering relationships among human RNA, demographic variables, and primary human cell physiology. The platelet RNA and expression-1 study measured platelet aggregation to arachidonic acid, ADP, protease-activated receptor (PAR) 1 activation peptide (PAR1-AP), and PAR4-AP, as well as mRNA and microRNA (miRNA) levels in platelets from 84 white and 70 black healthy subjects. A total of 5911 uniquely mapped mRNAs and 181 miRNAs were commonly expressed and validated in a separate cohort. One hundred twenty-nine mRNAs and 15 miRNAs were differentially expressed (DE) by age, and targets of these miRNAs were over-represented among these mRNAs. Fifty-four mRNAs and 9 miRNAs were DE by gender. Networks of miRNAs targeting mRNAs, both DE by age and gender, were identified. The inverse relationship in these RNA pairs suggests miRNAs regulate mRNA levels on aging and between genders. A simple, interactive public web tool (www.plateletomics.com) was developed that permits queries of RNA levels and associations among RNA, platelet aggregation and demographic variables. Access to these data will facilitate discovery of mechanisms of miRNA regulation of gene expression. These results provide new insights into aging and gender, and future platelet RNA association studies must account for age and gender.

Leveraging Integrated RNA Sequencing to Decipher Adrenomedullin's Protective Mechanisms in Experimental Bronchopulmonary Dysplasia.

Bronchopulmonary dysplasia (BPD) is a chronic lung disease commonly affecting premature infants, with limited therapeutic options and increased long-term consequences. Adrenomedullin (Adm), a proangiogenic peptide hormone, has been found to protect rodents against experimental BPD. This study aims to elucidate the molecular and cellular mechanisms through which Adm influences BPD pathogenesis using a lipopolysaccharide (LPS)-induced model of experimental BPD in mice. Bulk RNA sequencing of Adm-sufficient (wild-type or Adm+/+) and Adm-haplodeficient (Adm+/-) mice lungs, integrated with single-cell RNA sequencing data, revealed distinct gene expression patterns and cell type alterations associated with Adm deficiency and LPS exposure. Notably, computational integration with cell atlas data revealed that Adm-haplodeficient mouse lungs exhibited gene expression signatures characteristic of increased inflammation, natural killer (NK) cell frequency, and decreased endothelial cell and type II pneumocyte frequency. Furthermore, in silico human BPD patient data analysis supported our cell type frequency finding, highlighting elevated NK cells in BPD infants. These results underscore the protective role of Adm in experimental BPD and emphasize that it is a potential therapeutic target for BPD infants with an inflammatory phenotype.

Single-cell multiomics decodes regulatory programs for mouse secondary palate development.

Perturbations in gene regulation during palatogenesis can lead to cleft palate, which is among the most common congenital birth defects. Here, we perform single-cell multiome sequencing and profile chromatin accessibility and gene expression simultaneously within the same cells (n = 36,154) isolated from mouse secondary palate across embryonic days (E) 12.5, E13.5, E14.0, and E14.5. We construct five trajectories representing continuous differentiation of cranial neural crest-derived multipotent cells into distinct lineages. By linking open chromatin signals to gene expression changes, we characterize the underlying lineage-determining transcription factors. In silico perturbation analysis identifies transcription factors SHOX2 and MEOX2 as important regulators of the development of the anterior and posterior palate, respectively. In conclusion, our study charts epigenetic and transcriptional dynamics in palatogenesis, serving as a valuable resource for further cleft palate research.

DeepFace: Deep-learning-based framework to contextualize orofacial-cleft-related variants during human embryonic craniofacial development.

Orofacial clefts (OFCs) are among the most common human congenital birth defects. Previous multiethnic studies have identified dozens of associated loci for both cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP). Although several nearby genes have been highlighted, the "casual" variants are largely unknown. Here, we developed DeepFace, a convolutional neural network model, to assess the functional impact of variants by SNP activity difference (SAD) scores. The DeepFace model is trained with 204 epigenomic assays from crucial human embryonic craniofacial developmental stages of post-conception week (pcw) 4 to pcw 10. The Pearson correlation coefficient between the predicted and actual values for 12 epigenetic features achieved a median range of 0.50-0.83. Specifically, our model revealed that SNPs significantly associated with OFCs tended to exhibit higher SAD scores across various variant categories compared to less related groups, indicating a context-specific impact of OFC-related SNPs. Notably, we identified six SNPs with a significant linear relationship to SAD scores throughout developmental progression, suggesting that these SNPs could play a temporal regulatory role. Furthermore, our cell-type specificity analysis pinpointed the trophoblast cell as having the highest enrichment of risk signals associated with OFCs. Overall, DeepFace can harness distal regulatory signals from extensive epigenomic assays, offering new perspectives for prioritizing OFC variants using contextualized functional genomic features. We expect DeepFace to be instrumental in accessing and predicting the regulatory roles of variants associated with OFCs, and the model can be extended to study other complex diseases or traits.

Single-nucleus multiomics reveals the disrupted regulatory programs in three brain regions of sporadic early-onset Alzheimer's disease.

Sporadic early-onset Alzheimer's disease (sEOAD) represents a significant but less-studied subtype of Alzheimer's disease (AD). Here, we generated a single-nucleus multiome atlas derived from the postmortem prefrontal cortex, entorhinal cortex, and hippocampus of nine individuals with or without sEOAD. Comprehensive analyses were conducted to delineate cell type-specific transcriptomic changes and linked candidate cis- regulatory elements (cCREs) across brain regions. We prioritized seven conservative transcription factors in glial cells in multiple brain regions, including RFX4 in astrocytes and IKZF1 in microglia, which are implicated in regulating sEOAD-associated genes. Moreover, we identified the top 25 altered intercellular signaling between glial cells and neurons, highlighting their regulatory potential on gene expression in receiver cells. We reported 38 cCREs linked to sEOAD-associated genes overlapped with late-onset AD risk loci, and sEOAD cCREs enriched in neuropsychiatric disorder risk loci. This atlas helps dissect transcriptional and chromatin dynamics in sEOAD, providing a key resource for AD research.

Single-nucleus multiomics reveals the disrupted regulatory programs in three brain regions of sporadic early-onset Alzheimer's disease.

Sporadic early-onset Alzheimer's disease (sEOAD) represents a significant but less-studied subtype of Alzheimer's disease (AD). Here, we generated a single-nucleus multiome atlas derived from the postmortem prefrontal cortex, entorhinal cortex, and hippocampus of nine individuals with or without sEOAD. Comprehensive analyses were conducted to delineate cell type-specific transcriptomic changes and linked candidate cis- regulatory elements (cCREs) across brain regions. We prioritized seven conservative transcription factors in glial cells in multiple brain regions, including RFX4 in astrocytes and IKZF1 in microglia, which are implicated in regulating sEOAD-associated genes. Moreover, we identified the top 25 altered intercellular signaling between glial cells and neurons, highlighting their regulatory potential on gene expression in receiver cells. We reported 38 cCREs linked to sEOAD-associated genes overlapped with late-onset AD risk loci, and sEOAD cCREs enriched in neuropsychiatric disorder risk loci. This atlas helps dissect transcriptional and chromatin dynamics in sEOAD, providing a key resource for AD research.

Investigating gene functions and single-cell expression profiles of de novo variants in orofacial clefts.

Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs' deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p < 0.0002), while only DNV genes were enriched in the endothelium (adjusted p = 0.032). Although we did not achieve significant findings using CPO gene sets, which was mainly due to the limited number of DNV genes, scRNA-seq analysis implicated various expression patterns among DNV and OMIM genes. Our results suggest that combinatory pathway and scRNA-seq data analyses are helpful for contextualizing genes in OFC development.

The HSP90-MYC-CDK9 network drives therapeutic resistance in mantle cell lymphoma.

Brexucabtagene autoleucel CAR-T therapy is highly efficacious in overcoming resistance to Bruton's tyrosine kinase inhibitors (BTKi) in mantle cell lymphoma. However, many patients relapse post CAR-T therapy with dismal outcomes. To dissect the underlying mechanisms of sequential resistance to BTKi and CAR-T therapy, we performed single-cell RNA sequencing analysis for 66 samples from 25 patients treated with BTKi and/or CAR-T therapy and conducted in-depth bioinformatics™ analysis. Our analysis revealed that MYC activity progressively increased with sequential resistance. HSP90AB1 (Heat shock protein 90 alpha family class B member 1), a MYC target, was identified as early driver of CAR-T resistance. CDK9 (Cyclin-dependent kinase 9), another MYC target, was significantly upregulated in Dual-R samples. Both HSP90AB1 and CDK9 expression were correlated with MYC activity levels. Pharmaceutical co-targeting of HSP90 and CDK9 synergistically diminished MYC activity, leading to potent anti-MCL activity. Collectively, our study revealed that HSP90-MYC-CDK9 network is the primary driving force of therapeutic resistance.