Molecular and transcriptional basis of bidirectional CD4+ T cell exhaustion in oropharyngeal squamous cell carcinoma.

Tumor-infiltrating CD4+ T cells orchestrate the adaptive immune response through remarkable plasticity, and the expression patterns of exhaustion-related inhibitory receptors in these cells differ significantly from those of CD8+ T cells. Thus, a better understanding of the molecular basis of CD4+ T cell exhaustion and their responses to immune checkpoint blockade (ICB) is required. Here, we integrated multiomics approaches to define the phenotypic and molecular profiles of exhausted CD4+ T cells in oropharyngeal squamous cell carcinoma (OPSCC). Two distinct immune-promoting (Module 1) and immunosuppressive (Module 2) functional modules in tumor-infiltrating CD4+ T cells were identified, and both the immune-promoting function of Module 1 cells and immunosuppressive function of Module 2 cells were positively associated with their corresponding exhaustion states. Furthermore, the application of ICBs targeting effector CD4+ T cells in Module 1 (αPD-1) and Treg cells in Module 2 (αCTLA-4) in mouse models could help reinvigorate the effector function of Module 1-exhausted CD4+ T cells and reduce the immunosuppressive function of Module 2-exhausted CD4+ T cells, ultimately promoting OPSCC tumor regression. Taken together, our study provides a crucial cellular basis for the selection of optimal ICB in treating OPSCC.

Bacterial interactome disturbance in chronic obstructive pulmonary disease clinical stability and exacerbations.

RATIONALE: Our understanding of airway dysbiosis in chronic obstructive pulmonary disease (COPD) remains incomplete, which may be improved by unraveling the complexity in microbial interactome. OBJECTIVES: To characterize reproducible features of airway bacterial interactome in COPD at clinical stability and during exacerbation, and evaluate their associations with disease phenotypes. METHODS: We performed weighted ensemble-based co-occurrence network analysis of 1742 sputum microbiomes from published and new microbiome datasets, comprising two case-control studies of stable COPD versus healthy control, two studies of COPD stability versus exacerbation, and one study with exacerbation-recovery time series data. RESULTS: Patients with COPD had reproducibly lower degree of negative bacterial interactions, i.e. total number of negative interactions as a proportion of total interactions, in their airway microbiome compared with healthy controls. Evaluation of the Haemophilus interactome showed that the antagonistic interaction networks of this established pathogen rather than its abundance consistently changed in COPD. Interactome dynamic analysis revealed reproducibly reduced antagonistic interactions but not diversity loss during COPD exacerbation, which recovered after treatment. In phenotypic analysis, unsupervised network clustering showed that loss of antagonistic interactions was associated with worse clinical symptoms (dyspnea), poorer lung function, exaggerated neutrophilic inflammation, and higher exacerbation risk. Furthermore, the frequent exacerbators (≥ 2 exacerbations per year) had significantly reduced antagonistic bacterial interactions while exhibiting subtle compositional changes in their airway microbiota. CONCLUSIONS: Bacterial interactome disturbance characterized by reduced antagonistic interactions, rather than change in pathogen abundance or diversity, is a reproducible feature of airway dysbiosis in COPD clinical stability and exacerbations, which suggests that we may target interactome rather than pathogen alone for disease treatment.

A fine-scale Arabidopsis chromatin landscape reveals chromatin conformation-associated transcriptional dynamics.

Plants, as sessile organisms, deploy transcriptional dynamics for adapting to extreme growth conditions such as cold stress. Emerging evidence suggests that chromatin architecture contributes to transcriptional regulation. However, the relationship between chromatin architectural dynamics and transcriptional reprogramming in response to cold stress remains unclear. Here, we apply a chemical-crosslinking assisted proximity capture (CAP-C) method to elucidate the fine-scale chromatin landscape, revealing chromatin interactions within gene bodies closely associated with RNA polymerase II (Pol II) densities across initiation, pausing, and termination sites. We observe dynamic changes in chromatin interactions alongside Pol II activity alterations during cold stress, suggesting local chromatin dynamics may regulate Pol II activity. Notably, cold stress does not affect large-scale chromatin conformations. We further identify a comprehensive promoter-promoter interaction (PPI) network across the genome, potentially facilitating co-regulation of gene expression in response to cold stress. Our study deepens the understanding of chromatin conformation-associated gene regulation in plant response to cold.

The diversity of inhibitory receptor co-expression patterns of exhausted CD8+ T cells in oropharyngeal carcinoma.

Exhausted CD8+ T cells (Texs) are characterized by the expression of various inhibitory receptors (IRs), whereas the functional attributes of these co-expressed IRs remain limited. Here, we systematically characterized the diversity of IR co-expression patterns in Texs from both human oropharyngeal squamous cell carcinoma (OPSCC) tissues and syngeneic OPSCC model. Nearly 60% of the Texs population co-expressed two or more IRs, and the number of co-expressed IRs was positively associated with superior exhaustion and cytotoxicity phenotypes. In OPSCC patients, programmed cell death-1 (PD-1) blockade significantly enhanced PDCD1-based co-expression with other IR genes, whereas dual blockades of PD-1 and cytotoxic T lymphocyte-associated protein 4 (CTLA-4) significantly upregulated CTLA4-based co-expression with other IR genes. Collectively, our findings demonstrate that highly diverse IR co-expression is a leading feature of Texs and represents their functional states, which might provide essential clues for the rational selection of immune checkpoint inhibitors in treating OPSCC.

iM-Seeker: a webserver for DNA i-motifs prediction and scoring via automated machine learning.

DNA, beyond its canonical B-form double helix, adopts various alternative conformations, among which the i-motif, emerging in cytosine-rich sequences under acidic conditions, holds significant biological implications in transcription modulation and telomere biology. Despite recognizing the crucial role of i-motifs, predictive software for i-motif forming sequences has been limited. Addressing this gap, we introduce 'iM-Seeker', an innovative computational platform designed for the prediction and evaluation of i-motifs. iM-Seeker exhibits the capability to identify potential i-motifs within DNA segments or entire genomes, calculating stability scores for each predicted i-motif based on parameters such as the cytosine tracts number, loop lengths, and sequence composition. Furthermore, the webserver leverages automated machine learning (AutoML) to effortlessly fine-tune the optimal i-motif scoring model, incorporating user-supplied experimental data and customised features. As an advanced, versatile approach, 'iM-Seeker' promises to advance genomic research, highlighting the potential of i-motifs in cell biology and therapeutic applications. The webserver is freely available at https://im-seeker.org.

The PTM profiling of CTCF reveals the regulation of 3D chromatin structure by O-GlcNAcylation.

CCCTC-binding factor (CTCF), a ubiquitously expressed and highly conserved protein, is known to play a critical role in chromatin structure. Post-translational modifications (PTMs) diversify the functions of protein to regulate numerous cellular processes. However, the effects of PTMs on the genome-wide binding of CTCF and the organization of three-dimensional (3D) chromatin structure have not been fully understood. In this study, we uncovered the PTM profiling of CTCF and demonstrated that CTCF can be O-GlcNAcylated and arginine methylated. Functionally, we demonstrated that O-GlcNAcylation inhibits CTCF binding to chromatin. Meanwhile, deficiency of CTCF O-GlcNAcylation results in the disruption of loop domains and the alteration of chromatin loops associated with cellular development. Furthermore, the deficiency of CTCF O-GlcNAcylation increases the expression of developmental genes and negatively regulates maintenance and establishment of stem cell pluripotency. In conclusion, these results provide key insights into the role of PTMs for the 3D chromatin structure.

Precise prediction of phase-separation key residues by machine learning.

Understanding intracellular phase separation is crucial for deciphering transcriptional control, cell fate transitions, and disease mechanisms. However, the key residues, which impact phase separation the most for protein phase separation function have remained elusive. We develop PSPHunter, which can precisely predict these key residues based on machine learning scheme. In vivo and in vitro validations demonstrate that truncating just 6 key residues in GATA3 disrupts phase separation, enhancing tumor cell migration and inhibiting growth. Glycine and its motifs are enriched in spacer and key residues, as revealed by our comprehensive analysis. PSPHunter identifies nearly 80% of disease-associated phase-separating proteins, with frequent mutated pathological residues like glycine and proline often residing in these key residues. PSPHunter thus emerges as a crucial tool to uncover key residues, facilitating insights into phase separation mechanisms governing transcriptional control, cell fate transitions, and disease development.

Prediction of DNA i-motifs via machine learning.

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine.

DASES: a database of alternative splicing for esophageal squamous cell carcinoma.

Esophageal carcinoma ranks as the sixth leading cause of cancer-related mortality globally, with esophageal squamous cell carcinoma (ESCC) being particularly prevalent among Asian populations. Alternative splicing (AS) plays a pivotal role in ESCC development and progression by generating diverse transcript isoforms. However, the current landscape lacks a specialized database focusing on alternative splicing events (ASEs) derived from a large number of ESCC cases. Additionally, most existing AS databases overlook the contribution of long non-coding RNAs (lncRNAs) in ESCC molecular mechanisms, predominantly focusing on mRNA-based ASE identification. To address these limitations, we deployed DASES (http://www.hxdsjzx.cn/DASES). Employing a combination of publicly available and in-house ESCC RNA-seq datasets, our extensive analysis of 346 samples, with 93% being paired tumor and adjacent non-tumor tissues, led to the identification of 257 novel lncRNAs in esophageal squamous cell carcinoma. Leveraging a paired comparison of tumor and adjacent normal tissues, DASES identified 59,094 ASEs that may be associated with ESCC. DASES fills a critical gap by providing comprehensive insights into ASEs in ESCC, encompassing lncRNAs and mRNA, thus facilitating a deeper understanding of ESCC molecular mechanisms and serving as a valuable resource for ESCC research communities.

Genetic characteristics of common variable immunodeficiency patients with autoimmunity.

Background: The pathogenesis of common variable immunodeficiency disorder (CVID) is complex, especially when combined with autoimmunity. Genetic factors may be potential explanations for this complex situation, and whole genome sequencing (WGS) provide the basis for this potential. Methods: Genetic information of patients with CVID with autoimmunity, together with their first-degree relatives, was collected through WGS. The association between genetic factors and clinical phenotypes was studied using genetic analysis strategies such as sporadic and pedigree. Results: We collected 42 blood samples for WGS (16 CVID patients and 26 first-degree relatives of healthy controls). Through pedigree, sporadic screening strategies and low-frequency deleterious screening of rare diseases, we obtained 9,148 mutation sites, including 8,171 single-nucleotide variants (SNVs) and 977 Insertion-deletions (InDels). Finally, we obtained a total of 28 candidate genes (32 loci), of which the most common mutant was LRBA. The most common autoimmunity in the 16 patients was systematic lupus erythematosis. Through KEGG pathway enrichment, we identified the top ten signaling pathways, including "primary immunodeficiency", "JAK-STAT signaling pathway", and "T-cell receptor signaling pathway". We used PyMOL to predict and analyse the three-dimensional protein structures of the NFKB1, RAG1, TIRAP, NCF2, and MYB genes. In addition, we constructed a PPI network by combining candidate mutants with genes associated with CVID in the OMIM database via the STRING database. Conclusion: The genetic background of CVID includes not only monogenic origins but also oligogenic effects. Our study showed that immunodeficiency and autoimmunity may overlap in genetic backgrounds. Clinical Trial Registration: identifier ChiCTR2100044035.

[Role of Liquid-Liquid Phase Separation in Cell Fate Transition and Diseases].

Liquid-liquid phase separation (LLPS), a novel mechanism of the organization and formation of cellular structures, plays a vital role in regulating cell fate transitions and disease pathogenesis and is gaining widespread attention. LLPS may lead to the assemblage of cellular structures with liquid-like fluidity, such as germ granules, stress granules, and nucleoli, which are classic membraneless organelles. These structures are typically formed through the high-concentration liquid aggregation of biomacromolecules driven by weak multivalent interactions. LLPS is involved in regulating various intracellular life activities and its dysregulation may cause the disruption of cellular functions, thereby contributing to the pathogenesis and development of neurodegenerative diseases, infectious diseases, cancers, etc. Herein, we summarized published findings on the LLPS dynamics of membraneless organelles in physiological and pathological cell fate transition, revealing their crucial roles in cell differentiation, development, and various pathogenic processes. This paper provides a fresh theoretical framework and potential therapeutic targets for LLPS-related studies, opening new avenues for future research.

Prevotella copri alleviates sarcopenia via attenuating muscle mass loss and function decline.

BACKGROUND: The gut microbiome and fecal metabolites have been found to influence sarcopenia, but whether there are potential bacteria that can alleviate sarcopenia has been under-investigated, and the molecular mechanism remains unclear. METHODS: To investigate the relationships between the gut microbiome, fecal metabolites and sarcopenia, subjects were selected from observational multi-ethnic study conducted in Western China. Sarcopenia was diagnosed according to the criteria of the Asian Working Group for Sarcopenia 2014. The gut microbiome was profiled by shotgun metagenomic sequencing. Untargeted metabolomic analysis was performed to analyse the differences in fecal metabolites. We investigated bacterium with the greatest relative abundance difference between healthy individuals and sarcopenia patients, and the differences in metabolites associated with the bacteria, to verify its effects on muscle mass and function in a mouse model. RESULTS: The study included 283 participants (68.90% females, mean age: 66.66 years old) with and without sarcopenia (141 and 142 participants, respectively) and from the Han (98 participants), Zang (88 participants) and Qiang (97 participants) ethnic groups. This showed an overall reduction (15.03% vs. 20.77%, P = 0.01) of Prevotella copri between the sarcopenia and non-sarcopenia subjects across the three ethnic groups. Functional characterization of the differential bacteria showed enrichment (odds ratio = 15.97, P = 0.0068) in branched chain amino acid (BCAA) metabolism in non-sarcopenia group. A total of 13 BCAA and their derivatives have relatively low levels in sarcopenia. In the in vivo experiment, we found that the blood BCAA level was higher in the mice gavaged with live P. copri (LPC) (P < 0.001). The LPC mice had significantly longer wire and grid hanging time (P < 0.02), longer time on rotor (P = 0.0001) and larger grip strength (P < 0.0001), indicating better muscle function. The weight of gastrocnemius mass and rectus femoris mass (P < 0.05) was higher in LPC mice. The micro-computed tomography showed a larger leg area (P = 0.0031), and a small animal analyser showed a higher lean mass ratio in LPC mice (P = 0.0157), indicating higher muscle mass. CONCLUSIONS: The results indicated that there were lower levels of both P. copri and BCAA in sarcopenia individuals. In vivo experiments, gavage with LPC could attenuate muscle mass and function decline, indicating alleviating sarcopenia. This suggested that P. copri may play a therapeutic potential role in the management of sarcopenia.

Simultaneous profiling of chromatin architecture and transcription in single cells.

The three-dimensional structure of chromatin plays a crucial role in development and disease, both of which are associated with transcriptional changes. However, given the heterogeneity in single-cell chromatin architecture and transcription, the regulatory relationship between the three-dimensional chromatin structure and gene expression is difficult to explain based on bulk cell populations. Here we develop a single-cell, multimodal, omics method allowing the simultaneous detection of chromatin architecture and messenger RNA expression by sequencing (single-cell transcriptome sequencing (scCARE-seq)). Applying scCARE-seq to examine chromatin architecture and transcription from 2i to serum single mouse embryonic stem cells, we observe improved separation of cell clusters compared with single-cell chromatin conformation capture. In addition, after defining the cell-cycle phase of each cell through chromatin architecture extracted by scCARE-seq, we find that periodic changes in chromatin architecture occur in parallel with transcription during the cell cycle. These findings highlight the potential of scCARE-seq to facilitate comprehensive analyses that may boost our understanding of chromatin architecture and transcription in the same single cell.

Proliferative exhausted CD8+ T cells exacerbate long-lasting anti-tumor effects in human papillomavirus-positive head and neck squamous cell carcinoma.

The survival prognosis of human papillomavirus (HPV)-positive and HPV-negative head and neck squamous cell carcinoma (HNSCC) is largely different, and little is known about the anti-tumor mechanism of tumor-infiltrated exhausted CD8+ T cells (Tex) in HNSCC. We performed cell-level multi-omics sequencing on human HNSCC samples to decipher the multi-dimensional characteristics of Tex cells. A proliferative exhausted CD8+ T cell cluster (P-Tex) which was beneficial to survival outcomes of patients with HPV-positive HNSCC was identified. Interestingly, P-Tex cells expressed CDK4 genes as high as cancer cells, which could be simultaneously inhibited by CDK4 inhibitors and might be a potential reason for the ineffectiveness of CDK4 inhibitors in treating HPV-positive HNSCC. P-Tex cells could aggregate in the antigen-presenting cell niches and activate certain signaling pathways. Together, our findings suggest a promising role for P-Tex cells in the prognosis of patients with HPV-positive HNSCC by providing modest but persistent anti-tumor effects.

The comprehension of relative clauses in Mandarin Children with suspected specific language impairment.

This paper investigates the comprehension of Relative Clauses (RCs) in 15 Mandarin children with suspected Specific Language Impairment (SLI) (aged between 4; 5 and 6; 0) and 29 typically developing (TD) controls. Results from a Character Picture Matching Task indicate that (i) the subject RC was better understood than the object RC in children with SLI, but there was no asymmetry in the comprehension of the two RCs in TD children; (ii) the performance of children with SLI was significantly worse than that of their TD peers; (iii) children with SLI were prone to committing thematic role reversal errors and middle errors. In order to overcome the shortcomings of previous accounts, we therefore put forward the Edge Feature Underspecification Hypothesis, which can not only explain the asymmetry of comprehension seen in children with SLI but also shed light on the nature of errors committed by them in the task.

G4Atlas: a comprehensive transcriptome-wide G-quadruplex database.

RNA G-quadruplex (rG4) is a vital RNA tertiary structure motif that involves the base pairs on both Hoogsteen and Watson-Crick faces of guanines. rG4 is of great importance in the post-transcriptional regulation of gene expression. Experimental technologies have advanced to identify in vitro and in vivo rG4s across diverse transcriptomes. Building on these recent advances, here we present G4Atlas, the first transcriptome-wide G-quadruplex database, in which we have collated, classified, and visualized transcriptome rG4 experimental data, generated from rG4-seq, chemical profiling and ligand-binding methods. Our comprehensive database includes transcriptome-wide rG4s generated from 82 experimental treatments and 238 samples across ten species. In addition, we have also included RNA secondary structure prediction information across both experimentally identified and unidentified rG4s to enable users to display any potential competitive folding between rG4 and RNA secondary structures. As such, G4Atlas will enable users to explore the general functions of rG4s in diverse biological processes. In addition, G4Atlas lays the foundation for further data-driven deep learning algorithms to examine rG4 structural features.

RNA G-quadruplex structure contributes to cold adaptation in plants.

Nucleotide composition is suggested to infer gene functionality and ecological adaptation of species to distinct environments. However, the underlying biological function of nucleotide composition dictating environmental adaptations is largely unknown. Here, we systematically analyze the nucleotide composition of transcriptomes across 1000 plants (1KP) and their corresponding habitats. Intriguingly, we find that plants growing in cold climates have guanine (G)-enriched transcriptomes, which are prone to forming RNA G-quadruplex structures. Both immunofluorescence detection and in vivo structure profiling reveal that RNA G-quadruplex formation in plants is globally enhanced in response to cold. Cold-responsive RNA G-quadruplexes strongly enhanced mRNA stability, rather than affecting translation. Disruption of individual RNA G-quadruplex promotes mRNA decay in the cold, leading to impaired plant cold response. Therefore, we propose that plants adopted RNA G-quadruplex structure as a molecular signature to facilitate their adaptation to the cold during evolution.

Guide RNA engineering enables efficient CRISPR editing with a miniature Syntrophomonas palmitatica Cas12f1 nuclease.

Gene therapy is limited by inefficient delivery of large clustered regularly interspaced short palindromic repeat (CRISPR) effectors, such as Cas9 and Cas12a nucleases. Cas12f nucleases are currently one of the most compact CRISPR genome editors. However, the available toolkit of efficient Cas12f editors is limited. Here, we report the characterization and engineering of a miniature CRISPR-Cas12f system from Syntrophomonas palmitatica (SpaCas12f1, 497 amino acids). We show that CRISPR-SpaCas12f1 cleaves double-stranded DNA (dsDNA) with 5' T-rich PAM specificity and is naturally active for genome editing in bacteria. We identify that CRISPR-SpaCas12f1 trans-activating CRISPR RNA (tracrRNA) harbors a unique head-to-toe hairpin structure, and the natural hairpin structure is a key factor in restricting genome editing by SpaCas12f1 in human cells. Systematical engineering of SpaCas12f1 guide RNA transforms CRISPR-SpaCas12f1 into an efficient genome editor comparable to Francisella novicida CRISPR-Cas12a. Our findings expand the mini CRISPR toolbox, paving the way for therapeutic applications of CRISPR-SpaCas12f1 and engineering compact genome manipulation technologies.