Oleic acid availability impacts thymocyte preprogramming and subsequent peripheral Treg cell differentiation.

The nature of activation signals is essential in determining T cell subset differentiation; however, the features that determine T cell subset preference acquired during intrathymic development remain elusive. Here we show that naive CD4+ T cells generated in the mouse thymic microenvironment lacking Scd1, encoding the enzyme catalyzing oleic acid (OA) production, exhibit enhanced regulatory T (Treg) cell differentiation and attenuated development of experimental autoimmune encephalomyelitis. Scd1 deletion in K14+ thymic epithelia recapitulated the enhanced Treg cell differentiation phenotype of Scd1-deficient mice. The dearth of OA permitted DOT1L to increase H3K79me2 levels at the Atp2a2 locus of thymocytes at the DN2-DN3 transition stage. Such epigenetic modification persisted in naive CD4+ T cells and facilitated Atp2a2 expression. Upon T cell receptor activation, ATP2A2 enhanced the activity of the calcium-NFAT1-Foxp3 axis to promote naive CD4+ T cells to differentiate into Treg cells. Therefore, OA availability is critical for preprogramming thymocytes with Treg cell differentiation propensities in the periphery.

Proteogenomic characterization of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

Lung mesenchymal cells elicit lipid storage in neutrophils that fuel breast cancer lung metastasis.

Acquisition of a lipid-laden phenotype by immune cells has been defined in infectious diseases and atherosclerosis but remains largely uncharacterized in cancer. Here, in breast cancer models, we found that neutrophils are induced to accumulate neutral lipids upon interaction with resident mesenchymal cells in the premetastatic lung. Lung mesenchymal cells elicit this process through repressing the adipose triglyceride lipase (ATGL) activity in neutrophils in prostaglandin E2-dependent and -independent manners. In vivo, neutrophil-specific deletion of genes encoding ATGL or ATGL inhibitory factors altered neutrophil lipid profiles and breast tumor lung metastasis in mice. Mechanistically, lipids stored in lung neutrophils are transported to metastatic tumor cells through a macropinocytosis-lysosome pathway, endowing tumor cells with augmented survival and proliferative capacities. Pharmacological inhibition of macropinocytosis significantly reduced metastatic colonization by breast tumor cells in vivo. Collectively, our work reveals that neutrophils serve as an energy reservoir to fuel breast cancer lung metastasis.

KRAS Allelic Imbalance Enhances Fitness and Modulates MAP Kinase Dependence in Cancer.

Investigating therapeutic "outliers" that show exceptional responses to anti-cancer treatment can uncover biomarkers of drug sensitivity. We performed preclinical trials investigating primary murine acute myeloid leukemias (AMLs) generated by retroviral insertional mutagenesis in KrasG12D "knockin" mice with the MEK inhibitor PD0325901 (PD901). One outlier AML responded and exhibited intrinsic drug resistance at relapse. Loss of wild-type (WT) Kras enhanced the fitness of the dominant clone and rendered it sensitive to MEK inhibition. Similarly, human colorectal cancer cell lines with increased KRAS mutant allele frequency were more sensitive to MAP kinase inhibition, and CRISPR-Cas9-mediated replacement of WT KRAS with a mutant allele sensitized heterozygous mutant HCT116 cells to treatment. In a prospectively characterized cohort of patients with advanced cancer, 642 of 1,168 (55%) with KRAS mutations exhibited allelic imbalance. These studies demonstrate that serial genetic changes at the Kras/KRAS locus are frequent in cancer and modulate competitive fitness and MEK dependency.

Lung fibroblasts facilitate pre-metastatic niche formation by remodeling the local immune microenvironment.

Primary tumors are drivers of pre-metastatic niche formation, but the coordination by the secondary organ toward metastatic dissemination is underappreciated. Here, by single-cell RNA sequencing and immunofluorescence, we identified a population of cyclooxygenase 2 (COX-2)-expressing adventitial fibroblasts that remodeled the lung immune microenvironment. At steady state, fibroblasts in the lungs produced prostaglandin E2 (PGE2), which drove dysfunctional dendritic cells (DCs) and suppressive monocytes. This lung-intrinsic stromal program was propagated by tumor-associated inflammation, particularly the pro-inflammatory cytokine interleukin-1ß, supporting a pre-metastatic niche. Genetic ablation of Ptgs2 (encoding COX-2) in fibroblasts was sufficient to reverse the immune-suppressive phenotypes of lung-resident myeloid cells, resulting in heightened immune activation and diminished lung metastasis in multiple breast cancer models. Moreover, the anti-metastatic activity of DC-based therapy and PD-1 blockade was improved by fibroblast-specific Ptgs2 deletion or dual inhibition of PGE2 receptors EP2 and EP4. Collectively, lung-resident fibroblasts reshape the local immune landscape to facilitate breast cancer metastasis.

Transport of sugars.

Soluble sugars serve five main purposes in multicellular organisms: as sources of carbon skeletons, osmolytes, signals, and transient energy storage and as transport molecules. Most sugars are derived from photosynthetic organisms, particularly plants. In multicellular organisms, some cells specialize in providing sugars to other cells (e.g., intestinal and liver cells in animals, photosynthetic cells in plants), whereas others depend completely on an external supply (e.g., brain cells, roots and seeds). This cellular exchange of sugars requires transport proteins to mediate uptake or release from cells or subcellular compartments. Thus, not surprisingly, sugar transport is critical for plants, animals, and humans. At present, three classes of eukaryotic sugar transporters have been characterized, namely the glucose transporters (GLUTs), sodium-glucose symporters (SGLTs), and SWEETs. This review presents the history and state of the art of sugar transporter research, covering genetics, biochemistry, and physiology-from their identification and characterization to their structure, function, and physiology. In humans, understanding sugar transport has therapeutic importance (e.g., addressing diabetes or limiting access of cancer cells to sugars), and in plants, these transporters are critical for crop yield and pathogen susceptibility.

CARD9S12N facilitates the production of IL-5 by alveolar macrophages for the induction of type 2 immune responses.

The adaptor CARD9 functions downstream of C-type lectin receptors (CLRs) for the sensing of microbial infection, which leads to responses by the TH1 and TH17 subsets of helper T cells. The single-nucleotide polymorphism rs4077515 at CARD9 in the human genome, which results in the substitution S12N (CARD9S12N), is associated with several autoimmune diseases. However, the function of CARD9S12N has remained unknown. Here we generated CARD9S12N knock-in mice and found that CARD9S12N facilitated the induction of type 2 immune responses after engagement of CLRs. Mechanistically, CARD9S12N mediated CLR-induced activation of the non-canonical transcription factor NF-κB subunit RelB, which initiated production of the cytokine IL-5 in alveolar macrophages for the recruitment of eosinophils to drive TH2 cell-mediated allergic responses. We identified the homozygous CARD9 mutation encoding S12N in patients with allergic bronchopulmonary aspergillosis and revealed activation of RelB and production of IL-5 in peripheral blood mononuclear cells from these patients. Our study provides genetic and functional evidence demonstrating that CARD9S12N can turn alveolar macrophages into IL-5-producing cells and facilitates TH2 cell-mediated pathologic responses.

Parental histone transfer caught at the replication fork.

In eukaryotes, DNA compacts into chromatin through nucleosomes1,2. Replication of the eukaryotic genome must be coupled to the transmission of the epigenome encoded in the chromatin3,4. Here we report cryo-electron microscopy structures of yeast (Saccharomyces cerevisiae) replisomes associated with the FACT (facilitates chromatin transactions) complex (comprising Spt16 and Pob3) and an evicted histone hexamer. In these structures, FACT is positioned at the front end of the replisome by engaging with the parental DNA duplex to capture the histones through the middle domain and the acidic carboxyl-terminal domain of Spt16. The H2A-H2B dimer chaperoned by the carboxyl-terminal domain of Spt16 is stably tethered to the H3-H4 tetramer, while the vacant H2A-H2B site is occupied by the histone-binding domain of Mcm2. The Mcm2 histone-binding domain wraps around the DNA-binding surface of one H3-H4 dimer and extends across the tetramerization interface of the H3-H4 tetramer to the binding site of Spt16 middle domain before becoming disordered. This arrangement leaves the remaining DNA-binding surface of the other H3-H4 dimer exposed to additional interactions for further processing. The Mcm2 histone-binding domain and its downstream linker region are nested on top of Tof1, relocating the parental histones to the replisome front for transfer to the newly synthesized lagging-strand DNA. Our findings offer crucial structural insights into the mechanism of replication-coupled histone recycling for maintaining epigenetic inheritance.

Absence of near-ambient superconductivity in LuH2±xNy.

A recent study demonstrated near-ambient superconductivity in nitrogen-doped lutetium hydride1. This stimulated a worldwide interest in exploring room-temperature superconductivity at low pressures. Here, by using a high-pressure and high-temperature synthesis technique, we have obtained nitrogen-doped lutetium hydride (LuH2±xNy), which has a dark-blue colour and a structure with the space group [Formula: see text] as evidenced by X-ray diffraction. This structure is the same as that reported in ref. 1, with a slight difference in lattice constant. Raman spectroscopy of our samples also showed patterns similar to those observed in ref. 1. Energy-dispersive X-ray spectroscopy confirmed the presence of nitrogen in the samples. We observed a metallic behaviour from 350 K to 2 K at ambient pressure. On applying pressures from 2.1 GPa to 41 GPa, we observed a gradual colour change from dark blue to violet to pink-red. By measuring the resistance at pressures ranging from 0.4 GPa to 40.1 GPa, we observed a progressively improved metallic behaviour; however, superconductivity was not observed above 2 K. Temperature dependence of magnetization at high pressure shows a very weak positive signal between 100 K and 320 K, and the magnetization increases with an increase in magnetic field at 100 K. All of these are not expected for superconductivity above 100 K. Thus, we conclude the absence of near-ambient superconductivity in this nitrogen-doped lutetium hydride at pressures below 40.1 GPa.

The ZATT-TOP2A-PICH Axis Drives Extensive Replication Fork Reversal to Promote Genome Stability.

Replication fork reversal is a global response to replication stress in mammalian cells, but precisely how it occurs remains poorly understood. Here, we show that, upon replication stress, DNA topoisomerase IIalpha (TOP2A) is recruited to stalled forks in a manner dependent on the SNF2-family DNA translocases HLTF, ZRANB3, and SMARCAL1. This is accompanied by an increase in TOP2A SUMOylation mediated by the SUMO E3 ligase ZATT and followed by recruitment of a SUMO-targeted DNA translocase, PICH. Disruption of the ZATT-TOP2A-PICH axis results in accumulation of partially reversed forks and enhanced genome instability. These results suggest that fork reversal occurs via a sequential two-step process. First, HLTF, ZRANB3, and SMARCAL1 initiate limited fork reversal, creating superhelical strain in the newly replicated sister chromatids. Second, TOP2A drives extensive fork reversal by resolving the resulting topological barriers and via its role in recruiting PICH to stalled forks.

A statistical method for image-mediated association studies discovers genes and pathways associated with four brain disorders.

Brain imaging and genomics are critical tools enabling characterization of the genetic basis of brain disorders. However, imaging large cohorts is expensive and may be unavailable for legacy datasets used for genome-wide association studies (GWASs). Using an integrated feature selection/aggregation model, we developed an image-mediated association study (IMAS), which utilizes borrowed imaging/genomics data to conduct association mapping in legacy GWAS cohorts. By leveraging the UK Biobank image-derived phenotypes (IDPs), the IMAS discovered genetic bases underlying four neuropsychiatric disorders and verified them by analyzing annotations, pathways, and expression quantitative trait loci (eQTLs). A cerebellar-mediated mechanism was identified to be common to the four disorders. Simulations show that, if the goal is identifying genetic risk, our IMAS is more powerful than a hypothetical protocol in which the imaging results were available in the GWAS dataset. This implies the feasibility of reanalyzing legacy GWAS datasets without conducting additional imaging, yielding cost savings for integrated analysis of genetics and imaging.

The development of brain pericytes requires expression of the transcription factor nkx3.1 in intermediate precursors.

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation.

Comprehensive single-cell RNA-seq analysis using deep interpretable generative modeling guided by biological hierarchy knowledge.

Recent advances in microfluidics and sequencing technologies allow researchers to explore cellular heterogeneity at single-cell resolution. In recent years, deep learning frameworks, such as generative models, have brought great changes to the analysis of transcriptomic data. Nevertheless, relying on the potential space of these generative models alone is insufficient to generate biological explanations. In addition, most of the previous work based on generative models is limited to shallow neural networks with one to three layers of latent variables, which may limit the capabilities of the models. Here, we propose a deep interpretable generative model called d-scIGM for single-cell data analysis. d-scIGM combines sawtooth connectivity techniques and residual networks, thereby constructing a deep generative framework. In addition, d-scIGM incorporates hierarchical prior knowledge of biological domains to enhance the interpretability of the model. We show that d-scIGM achieves excellent performance in a variety of fundamental tasks, including clustering, visualization, and pseudo-temporal inference. Through topic pathway studies, we found that d-scIGM-learned topics are better enriched for biologically meaningful pathways compared to the baseline models. Furthermore, the analysis of drug response data shows that d-scIGM can capture drug response patterns in large-scale experiments, which provides a promising way to elucidate the underlying biological mechanisms. Lastly, in the melanoma dataset, d-scIGM accurately identified different cell types and revealed multiple melanin-related driver genes and key pathways, which are critical for understanding disease mechanisms and drug development.

Imaging viscous flow of the Dirac fluid in graphene.

The electron-hole plasma in charge-neutral graphene is predicted to realize a quantum critical system in which electrical transport features a universal hydrodynamic description, even at room temperature1,2. This quantum critical 'Dirac fluid' is expected to have a shear viscosity close to a minimum bound3,4, with an interparticle scattering rate saturating1 at the Planckian time, the shortest possible timescale for particles to relax. Although electrical transport measurements at finite carrier density are consistent with hydrodynamic electron flow in graphene5-8, a clear demonstration of viscous flow at the charge-neutrality point remains elusive. Here we directly image viscous Dirac fluid flow in graphene at room temperature by measuring the associated stray magnetic field. Nanoscale magnetic imaging is performed using quantum spin magnetometers realized with nitrogen vacancy centres in diamond. Scanning single-spin and wide-field magnetometry reveal a parabolic Poiseuille profile for electron flow in a high-mobility graphene channel near the charge-neutrality point, establishing the viscous transport of the Dirac fluid. This measurement is in contrast to the conventional uniform flow profile imaged in a metallic conductor and also in a low-mobility graphene channel. Via combined imaging and transport measurements, we obtain viscosity and scattering rates, and observe that these quantities are comparable to the universal values expected at quantum criticality. This finding establishes a nearly ideal electron fluid in charge-neutral, high-mobility graphene at room temperature4. Our results will enable the study of hydrodynamic transport in quantum critical fluids relevant to strongly correlated electrons in high-temperature superconductors9. This work also highlights the capability of quantum spin magnetometers to probe correlated electronic phenomena at the nanoscale.

MAE-seq refines regulatory elements across the genome.

Proper cell fate determination relies on precise spatial and temporal genome-wide cooperation between regulatory elements (REs) and their targeted genes. However, the lengths of REs defined using different methods vary, which indicates that there is sequence redundancy and that the context of the genome may be unintelligible. We developed a method called MAE-seq (Massive Active Enhancers by Sequencing) to experimentally identify functional REs at a 25-bp scale. In this study, MAE-seq was used to identify 626879, 541617 and 554826 25-bp enhancers in mouse embryonic stem cells (mESCs), C2C12 and HEK 293T, respectively. Using ∼1.6 trillion 25 bp DNA fragments and screening 12 billion cells, we identified 626879 as active enhancers in mESCs as an example. Comparative analysis revealed that most of the histone modification datasets were annotated by MAE-Seq loci. Furthermore, 33.85% (212195) of the identified enhancers were identified as de novo ones with no epigenetic modification. Intriguingly, distinct chromatin states dictate the requirement for dissimilar cofactors in governing novel and known enhancers. Validation results show that these 25-bp sequences could act as a functional unit, which shows identical or similar expression patterns as the previously defined larger elements, Enhanced resolution facilitated the identification of numerous cell-specific enhancers and their accurate annotation as super enhancers. Moreover, we characterized novel elements capable of augmenting gene activity. By integrating with high-resolution Hi-C data, over 55.64% of novel elements may have a distal association with different targeted genes. For example, we found that the Cdh1 gene interacts with one novel and two known REs in mESCs. The biological effects of these interactions were investigated using CRISPR-Cas9, revealing their role in coordinating Cdh1 gene expression and mESC proliferation. Our study presents an experimental approach to refine the REs at 25-bp resolution, advancing the precision of genome annotation and unveiling the underlying genome context. This novel approach not only advances our understanding of gene regulation but also opens avenues for comprehensive exploration of the genomic landscape.

A deregulated m6A writer complex axis driven by BRD4 confers an epitranscriptomic vulnerability in combined DNA repair-targeted therapy.

Aberrant transcripts expression of the m6A methyltransferase complex (MTC) is widely found across human cancers, suggesting a dysregulated signaling cascade which integrates m6A epitranscriptome to drive tumorigenesis. However, the responsible transcriptional machinery directing the expression of distinct MTC subunits remains unclear. Here, we identified an unappreciated interplay between the histone acetyl-lysine reader BRD4 and the m6A writer complex across human cancers. BRD4 directly stimulates transcripts expression of seven MTC subunits, allowing the maintenance of the nuclear writer complex integrity. Upon BET inhibition, this BRD4-MTC signaling cascade accounts for global m6A reduction and the subsequent dynamic alteration of BRD4-dependent transcriptome, resulting in impaired DNA damage response that involves activation of homologous recombination (HR) repair and repression of apoptosis. We further demonstrated that the combined synergy upon BET/PARP inhibition largely relies on disrupted m6A modification of HR and apoptotic genes, counteracting PARP inhibitor (PARPi) resistance in patient-derived xenograft models. Our study revealed a widespread active cross-talk between BRD4-dependent epigenetic and MTC-mediated epitranscriptomic networks, which provides a unique therapeutic vulnerability that can be leveraged in combined DNA repair-targeted therapy.

cLD: Rare-variant linkage disequilibrium between genomic regions identifies novel genomic interactions.

Linkage disequilibrium (LD) is a fundamental concept in genetics; critical for studying genetic associations and molecular evolution. However, LD measurements are only reliable for common genetic variants, leaving low-frequency variants unanalyzed. In this work, we introduce cumulative LD (cLD), a stable statistic that captures the rare-variant LD between genetic regions, which reflects more biological interactions between variants, in addition to lack of recombination. We derived the theoretical variance of cLD using delta methods to demonstrate its higher stability than LD for rare variants. This property is also verified by bootstrapped simulations using real data. In application, we find cLD reveals an increased genetic association between genes in 3D chromatin interactions, a phenomenon recently reported negatively by calculating standard LD between common variants. Additionally, we show that cLD is higher between gene pairs reported in interaction databases, identifies unreported protein-protein interactions, and reveals interacting genes distinguishing case/control samples in association studies.

Nuclear VCP drives colorectal cancer progression by promoting fatty acid oxidation.

Fatty acid oxidation (FAO) fuels many cancers. However, knowledge of pathways that drive FAO in cancer remains unclear. Here, we revealed that valosin-containing protein (VCP) upregulates FAO to promote colorectal cancer growth. Mechanistically, nuclear VCP binds to histone deacetylase 1 (HDAC1) and facilitates its degradation, thus promoting the transcription of FAO genes, including the rate-limiting enzyme carnitine palmitoyltransferase 1A (CPT1A). FAO is an alternative fuel for cancer cells in environments exhibiting limited glucose availability. We observed that a VCP inhibitor blocked the upregulation of FAO activity and CPT1A expression triggered by metformin in colorectal cancer (CRC) cells. Combined VCP inhibitor and metformin prove more effective than either agent alone in culture and in vivo. Our study illustrates the molecular mechanism underlying the regulation of FAO by nuclear VCP and demonstrates the potential therapeutic utility of VCP inhibitor and metformin combination treatment for colorectal cancer.