Automatic cell-type harmonization and integration across Human Cell Atlas datasets.

Harmonizing cell types across the single-cell community and assembling them into a common framework is central to building a standardized Human Cell Atlas. Here, we present CellHint, a predictive clustering tree-based tool to resolve cell-type differences in annotation resolution and technical biases across datasets. CellHint accurately quantifies cell-cell transcriptomic similarities and places cell types into a relationship graph that hierarchically defines shared and unique cell subtypes. Application to multiple immune datasets recapitulates expert-curated annotations. CellHint also reveals underexplored relationships between healthy and diseased lung cell states in eight diseases. Furthermore, we present a workflow for fast cross-dataset integration guided by harmonized cell types and cell hierarchy, which uncovers underappreciated cell types in adult human hippocampus. Finally, we apply CellHint to 12 tissues from 38 datasets, providing a deeply curated cross-tissue database with ∼3.7 million cells and various machine learning models for automatic cell annotation across human tissues.

Molecular topography of an entire nervous system.

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

Lactate Is a Natural Suppressor of RLR Signaling by Targeting MAVS.

RLR-mediated type I IFN production plays a pivotal role in elevating host immunity for viral clearance and cancer immune surveillance. Here, we report that glycolysis, which is inactivated during RLR activation, serves as a barrier to impede type I IFN production upon RLR activation. RLR-triggered MAVS-RIG-I recognition hijacks hexokinase binding to MAVS, leading to the impairment of hexokinase mitochondria localization and activation. Lactate serves as a key metabolite responsible for glycolysis-mediated RLR signaling inhibition by directly binding to MAVS transmembrane (TM) domain and preventing MAVS aggregation. Notably, lactate restoration reverses increased IFN production caused by lactate deficiency. Using pharmacological and genetic approaches, we show that lactate reduction by lactate dehydrogenase A (LDHA) inactivation heightens type I IFN production to protect mice from viral infection. Our study establishes a critical role of glycolysis-derived lactate in limiting RLR signaling and identifies MAVS as a direct sensor of lactate, which functions to connect energy metabolism and innate immunity.

Atad3a suppresses Pink1-dependent mitophagy to maintain homeostasis of hematopoietic progenitor cells.

Although deletion of certain autophagy-related genes has been associated with defects in hematopoiesis, it remains unclear whether hyperactivated mitophagy affects the maintenance and differentiation of hematopoietic stem cells (HSCs) and committed progenitor cells. Here we report that targeted deletion of the gene encoding the AAA+-ATPase Atad3a hyperactivated mitophagy in mouse hematopoietic cells. Affected mice showed reduced survival, severely decreased bone-marrow cellularity, erythroid anemia and B cell lymphopenia. Those phenotypes were associated with skewed differentiation of stem and progenitor cells and an enlarged HSC pool. Mechanistically, Atad3a interacted with the mitochondrial channel components Tom40 and Tim23 and served as a bridging factor to facilitate appropriate transportation and processing of the mitophagy protein Pink1. Loss of Atad3a caused accumulation of Pink1 and activated mitophagy. Notably, deletion of Pink1 in Atad3a-deficient mice significantly 'rescued' the mitophagy defect, which resulted in restoration of the progenitor and HSC pools. Our data indicate that Atad3a suppresses Pink1-dependent mitophagy and thereby serves a key role in hematopoietic homeostasis.

Immune microniches shape intestinal Treg function.

The intestinal immune system is highly adapted to maintaining tolerance to the commensal microbiota and self-antigens while defending against invading pathogens1,2. Recognizing how the diverse network of local cells establish homeostasis and maintains it in the complex immune environment of the gut is critical to understanding how tolerance can be re-established following dysfunction, such as in inflammatory disorders. Although cell and molecular interactions that control T regulatory (Treg) cell development and function have been identified3,4, less is known about the cellular neighbourhoods and spatial compartmentalization that shapes microorganism-reactive Treg cell function. Here we used in vivo live imaging, photo-activation-guided single-cell RNA sequencing5-7 and spatial transcriptomics to follow the natural history of T cells that are reactive towards Helicobacter hepaticus through space and time in the settings of tolerance and inflammation. Although antigen stimulation can occur anywhere in the tissue, the lamina propria-but not embedded lymphoid aggregates-is the key microniche that supports effector Treg (eTreg) cell function. eTreg cells are stable once their niche is established; however, unleashing inflammation breaks down compartmentalization, leading to dominance of CD103+SIRPα+ dendritic cells in the lamina propria. We identify and validate the putative tolerogenic interaction between CD206+ macrophages and eTreg cells in the lamina propria and identify receptor-ligand pairs that are likely to govern the interaction. Our results reveal a spatial mechanism of tolerance in the lamina propria and demonstrate how knowledge of local interactions may contribute to the next generation of tolerance-inducing therapies.

Inositol serves as a natural inhibitor of mitochondrial fission by directly targeting AMPK.

Mitochondrial dynamics regulated by mitochondrial fusion and fission maintain mitochondrial functions, whose alterations underline various human diseases. Here, we show that inositol is a critical metabolite directly restricting AMPK-dependent mitochondrial fission independently of its classical mode as a precursor for phosphoinositide generation. Inositol decline by IMPA1/2 deficiency elicits AMPK activation and mitochondrial fission without affecting ATP level, whereas inositol accumulation prevents AMPK-dependent mitochondrial fission. Metabolic stress or mitochondrial damage causes inositol decline in cells and mice to elicit AMPK-dependent mitochondrial fission. Inositol directly binds to AMPKγ and competes with AMP for AMPKγ binding, leading to restriction of AMPK activation and mitochondrial fission. Our study suggests that the AMP/inositol ratio is a critical determinant for AMPK activation and establishes a model in which AMPK activation requires inositol decline to release AMPKγ for AMP binding. Hence, AMPK is an inositol sensor, whose inactivation by inositol serves as a mechanism to restrict mitochondrial fission.

Challenges and perspectives in computational deconvolution of genomics data.

Deciphering cell-type heterogeneity is crucial for systematically understanding tissue homeostasis and its dysregulation in diseases. Computational deconvolution is an efficient approach for estimating cell-type abundances from a variety of omics data. Despite substantial methodological progress in computational deconvolution in recent years, challenges are still outstanding. Here we enlist four important challenges related to computational deconvolution: the quality of the reference data, generation of ground truth data, limitations of computational methodologies, and benchmarking design and implementation. Finally, we make recommendations on reference data generation, new directions of computational methodologies, and strategies to promote rigorous benchmarking.

Tiam1 methylation by NSD2 promotes Rac1 signaling activation and colon cancer metastasis.

Metastasis is a major cause of cancer therapy failure and mortality. However, targeting metastatic seeding and colonization remains a significant challenge. In this study, we identified NSD2, a histone methyltransferase responsible for dimethylating histone 3 at lysine 36, as being overexpressed in metastatic tumors. Our findings suggest that NSD2 overexpression enhances tumor metastasis both in vitro and in vivo. Further analysis revealed that NSD2 promotes tumor metastasis by activating Rac1 signaling. Mechanistically, NSD2 combines with and activates Tiam1 (T lymphoma invasion and metastasis 1) and promotes Rac1 signaling by methylating Tiam1 at K724. In vivo and in vitro studies revealed that Tiam1 K724 methylation could be a predictive factor for cancer prognosis and a potential target for metastasis inhibition. Furthermore, we have developed inhibitory peptide which was proved to inhibit tumor metastasis through blocking the interaction between NSD2 and Tiam1. Our results demonstrate that NSD2-methylated Tiam1 promotes Rac1 signaling and cancer metastasis. These results provide insights into the inhibition of tumor metastasis.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

Bi-allelic variants in chromatoid body protein TDRD6 cause spermiogenesis defects and severe oligoasthenoteratozoospermia in humans.

BACKGROUND: The association between the TDRD6 variants and human infertility remains unclear, as only one homozygous missense variant of TDRD6 was found to be associated with oligoasthenoteratozoospermia (OAT). METHODS: Whole-exome sequencing and Sanger sequencing were employed to identify potential pathogenic variants of TDRD6 in infertile men. Histology, immunofluorescence, immunoblotting and ultrastructural analyses were conducted to clarify the structural and functional abnormalities of sperm in mutated patients. Tdrd6-knockout mice were generated using the CRISPR-Cas9 system. Total RNA-seq and single-cell RNA-seq (scRNA-seq) analyses were used to elucidate the underlying molecular mechanisms, followed by validation through quantitative RT-PCR and immunostaining. Intracytoplasmic sperm injection (ICSI) was also used to assess the efficacy of clinical treatment. RESULTS: Bi-allelic TDRD6 variants were identified in five unrelated Chinese individuals with OAT, including homozygous loss-of-function variants in two consanguineous families. Notably, besides reduced concentrations and impaired motility, a significant occurrence of acrosomal hypoplasia was detected in multiple spermatozoa among five patients. Using the Tdrd6-deficient mice, we further elucidate the pivotal role of TDRD6 in spermiogenesis and acrosome identified. In addition, the mislocalisation of crucial chromatoid body components DDX4 (MVH) and UPF1 was also observed in round spermatids from patients harbouring TDRD6 variants. ScRNA-seq analysis of germ cells from a patient with TDRD6 variants revealed that TDRD6 regulates mRNA metabolism processes involved in spermatid differentiation and cytoplasmic translation. CONCLUSION: Our findings strongly suggest that TDRD6 plays a conserved role in spermiogenesis and confirms the causal relationship between TDRD6 variants and human OAT. Additionally, this study highlights the unfavourable ICSI outcomes in individuals with bi-allelic TDRD6 variants, providing insights for potential clinical treatment strategies.

Posttranslational regulation of liver kinase B1 in human cancer.

Liver kinase B1 (LKB1) is a serine-threonine kinase that participates in multiple cellular and biological processes, including energy metabolism, cell polarity, cell proliferation, cell migration, and many others. LKB1 is initially identified as a germline-mutated causative gene in Peutz-Jeghers syndrome and is commonly regarded as a tumor suppressor due to frequent inactivation in a variety of cancers. LKB1 directly binds and activates its downstream kinases including the AMP-activated protein kinase (AMPK) and AMPK-related kinases by phosphorylation, which has been intensively investigated for the past decades. An increasing number of studies have uncovered the posttranslational modifications (PTMs) of LKB1 and consequent changes in its localization, activity, and interaction with substrates. The alteration in LKB1 function as a consequence of genetic mutations and aberrant upstream signaling regulation leads to tumor development and progression. Here, we review current knowledge about the mechanism of LKB1 in cancer and the contributions of PTMs, such as phosphorylation, ubiquitination, SUMOylation, acetylation, prenylation, and others, to the regulation of LKB1 function, offering new insights into the therapeutic strategies in cancer.

How well do neural signatures of resting-state EEG detect consciousness? A large-scale clinical study.

The assessment of consciousness states, especially distinguishing minimally conscious states (MCS) from unresponsive wakefulness states (UWS), constitutes a pivotal role in clinical therapies. Despite that numerous neural signatures of consciousness have been proposed, the effectiveness and reliability of such signatures for clinical consciousness assessment still remains an intense debate. Through a comprehensive review of the literature, inconsistent findings are observed about the effectiveness of diverse neural signatures. Notably, the majority of existing studies have evaluated neural signatures on a limited number of subjects (usually below 30), which may result in uncertain conclusions due to small data bias. This study presents a systematic evaluation of neural signatures with large-scale clinical resting-state electroencephalography (EEG) signals containing 99 UWS, 129 MCS, 36 emergence from the minimally conscious state, and 32 healthy subjects (296 total) collected over 3 years. A total of 380 EEG-based metrics for consciousness detection, including spectrum features, nonlinear measures, functional connectivity, and graph-based measures, are summarized and evaluated. To further mitigate the effect of data bias, the evaluation is performed with bootstrap sampling so that reliable measures can be obtained. The results of this study suggest that relative power in alpha and delta serve as dependable indicators of consciousness. With the MCS group, there is a notable increase in the phase lag index-related connectivity measures and enhanced functional connectivity between brain regions in comparison to the UWS group. A combination of features enables the development of an automatic detector of conscious states.

SIRPG expression positively associates with an inflamed tumor microenvironment and response to PD-1 blockade.

BACKGROUND: This study aimed to investigate the relationship between signal regulatory protein gamma (SIRPG) and tumor immune microenvironment phenotypes or T cell mediated-adaptive antitumor immunity, and its predictive value for response to PD-1 blockade in cancers. METHODS: Pan-cancer analysis of SIRPG expression and immune deconvolution was performed using transcriptomic data across 33 tumor types. Transcriptomic and clinical data from 157 patients with non-small-cell lung cancer (NSCLC) and melanoma received PD-1 blockade were analyzed. Expression characteristics of SIRPG were investigated using single-cell RNA sequencing (scRNA-seq) data of 103,599 cells. The effect of SIRPG expression was evaluated via SIRPG knockdown or overexpression in Jurkat T cells. RESULTS: The results showed that most cancers with high SIRPG expression had significantly higher abundance of T cells, B cells, NK cells, M1 macrophages and cytotoxic lymphocytes and increased expression level of immunomodulatory factors regulating immune cell recruitment, antigen presentation, T cell activation and cytotoxicity, but markedly lower abundance of neutrophils, M2 macrophages, and myeloid-derived suppressor cells. High SIRPG expression was associated with favorable response to PD-1 blockade in both NSCLC and melanoma. scRNA-seq data suggested SIRPG was mainly expressed in CD8+ exhausted T and CD4+ regulatory T cells, and positively associated with immune checkpoint expression including PDCD1 and CTLA4. In vitro test showed SIRPG expression in T cells could facilitate expression of PDCD1 and CTLA4. CONCLUSION: High SIRPG expression is associated with an inflamed immune phenotype in cancers and favorable response to PD-1 blockade, suggesting it would be a promising predictive biomarker for PD-1 blockade and novel immunotherapeutic target.

Superior radiation tolerance via reversible disordering-ordering transition of coherent superlattices.

Materials capable of sustaining high radiation doses at a high temperature are required for next-generation fission and future fusion energy. To date, however, even the most promising structural materials cannot withstand the demanded radiation environment due to irreversible radiation-driven microstructure degradation. Here we report a counterintuitive strategy to achieve exceptionally high radiation tolerance at high temperatures by enabling reversible local disordering-ordering transition of the introduced superlattice nanoprecipitates in metallic materials. As particularly demonstrated in martensitic steel containing a high density of B2-ordered superlattices, no void swelling was detected even after ultrahigh-dose radiation damage at 400-600 °C. The reordering process of the low-misfit superlattices in highly supersaturated matrices occurs through the short-range reshuffling of radiation-induced point defects and excess solutes right after rapid, ballistic disordering. This dynamic process stabilizes the microstructure, continuously promotes in situ defect recombination and efficiently prevents the capillary-driven long-range diffusion process. The strategy can be readily applied into other materials and pave the pathway for developing materials with high radiation tolerance.

Novel MEIOB pathogenic variants including a homozygous non-canonical splicing variant, cause meiotic arrest and human non-obstructive azoospermia.

Non-obstructive azoospermia (NOA) is the most severe form of human male infertility, and the genetic causes of NOA with meiotic arrest remain largely unclear. In this study, we identified novel compound heterozygous MEIOB variants (c.814C > T: p.R272X and c.976G > A: p.A326T) and a previously undescribed homozygous non-canonical splicing variant of MEIOB (c.528 + 3A > C) in two NOA-affected individuals from two irrelevant Chinese families. MEIOB missense variant (p.A326T) significantly reduced protein abundance and nonsense variant (p.R272X) produced a truncated protein. Both of two variants impaired the MEIOB-SPATA22 interaction. The MEIOB non-canonical splicing variant resulted in whole Exon 6 skipping by minigene assay, which was predicted to produce a frameshift truncated protein (p.S111Rfs*32). Histological and immunostaining analysis indicated that both patients exhibited a similar phenotype as we previously reported in Meiob mutant mice, that is, absence of spermatids in seminiferous tubules and meiotic arrest. Our study identified three novel pathogenic variants of MEIOB in NOA patients, extending the mutation spectrum of the MEIOB and highlighting the contribution of meiotic recombination related genes in human fertility.

Coriaceumins A-D, the First Nitrogen-Containing Crenulide Diterpenoids from the Brown Alga Dictyota coriacea Collected in the East China Sea.

A phytochemical investigation of an East China Sea collection of the brown alga Dictyota coriacea has led to the isolation of four novel nitrogen-containing crenulide diterpenoids, named coriaceumins A-D (1-4), two rare nitrogenous xenicane diterpenoids, dictyolactams C (5) and D (6), and one known crenulide diterpenoid, hydroxycrenulide (7). The structures of the new compounds were elucidated by detailed spectroscopic data analyses, including HRESIMS and 1D/2D NMR. The absolute configurations were determined by a comparison of the experimental ECD spectra with the spectra computed by DFT-based quantum chemical calculations. Coriaceumins A-D (1-4) represent the first examples of nitrogen-containing crenulide diterpenoids. In a bioassay, compounds 2, 3, 5, and 7 were found to exhibit different levels of inhibitory effects against protein tyrosine phosphatase 1B (PTP1B) with IC50 values ranging from 7.3 to 19 µM. In addition, the primary structure-activity relationships of all the isolates were summarized.

Statistical learning in patients in the minimally conscious state.

When listening to speech, cortical activity can track mentally constructed linguistic units such as words, phrases, and sentences. Recent studies have also shown that the neural responses to mentally constructed linguistic units can predict the outcome of patients with disorders of consciousness (DoC). In healthy individuals, cortical tracking of linguistic units can be driven by both long-term linguistic knowledge and online learning of the transitional probability between syllables. Here, we investigated whether statistical learning could occur in patients in the minimally conscious state (MCS) and patients emerged from the MCS (EMCS) using electroencephalography (EEG). In Experiment 1, we presented to participants an isochronous sequence of syllables, which were composed of either 4 real disyllabic words or 4 reversed disyllabic words. An inter-trial phase coherence analysis revealed that the patient groups showed similar word tracking responses to real and reversed words. In Experiment 2, we presented trisyllabic artificial words that were defined by the transitional probability between words, and a significant word-rate EEG response was observed for MCS patients. These results suggested that statistical learning can occur with a minimal conscious level. The residual statistical learning ability in MCS patients could potentially be harnessed to induce neural plasticity.

TRAF6 Restricts p53 Mitochondrial Translocation, Apoptosis, and Tumor Suppression.

Mitochondrial p53 is involved in apoptosis and tumor suppression. However, its regulation is not well studied. Here, we show that TRAF6 E3 ligase is a crucial factor to restrict mitochondrial translocation of p53 and spontaneous apoptosis by promoting K63-linked ubiquitination of p53 at K24 in cytosol, and such ubiquitination limits the interaction between p53 and MCL-1/BAK. Genotoxic stress reduces this ubiquitination in cytosol by S13/T330 phosphorylation-dependent translocation of TRAF6 from cytosol to nucleus, where TRAF6 also facilitates the K63-linked ubiquitination of nuclear p53 and its transactivation by recruiting p300 for p53 acetylation. Functionally, K63-linked ubiquitination of p53 compromised p53-mediated apoptosis and tumor suppression. Colorectal cancer samples with WT p53 reveal that TRAF6 overexpression negatively correlates with apoptosis and predicts poor response to chemotherapy and radiotherapy. Together, our study identifies TRAF6 as a critical gatekeeper to restrict p53 mitochondrial translocation, and such mechanism may contribute to tumor development and drug resistance.

The effects of obstructive sleep apnea on blood pressure variability and load in patients with hypertension.

BACKGROUND: Hypertension frequently coexists with obstructive sleep apnea (OSA), and their interplay substantially impacts the prognosis of affected individuals. Investigating the influence of OSA on blood pressure variability (BPV) and blood pressure load (BPL) in hypertensive patients has become a focal point of clinical research. METHODS: This cross-sectional study recruited hypertensive patients (n = 265) without discrimination and classified them into four groups based on their apnea-hypopnea index (AHI): control group (n = 40), AHI < 5; mild group (n = 74), 5 ≤ AHI ≤ 15; moderate group (n = 68), 15 < AHI ≤ 30; severe group (n = 83), AHI > 30. All participants underwent comprehensive assessments, including polysomnography (PSG) monitoring, 24-h ambulatory blood pressure (ABP) monitoring, cardiac Doppler ultrasound, and additional examinations when indicated. RESULTS: BPV and BPL exhibited significant elevations in the moderate and severe OSA groups compared to the control and mild OSA groups (P < 0.05). Moreover, interventricular septum thickness and left ventricular end-diastolic volume (LVEDV) demonstrated higher values in the moderate and severe OSA groups (P < 0.05). Multiple stepwise regression analysis identified noteworthy risk factors for elevated BPV in hypertensive patients with OSA, including AHI, maximum apnea time, total times of oxygen reduction, and mean time of apnea. CONCLUSION: Hypertensive patients with moderate to severe OSA exhibited substantially increased BPV and BPL. Moreover, BPV was correlated with AHI, maximum apnea time, total times of oxygen reduction, and mean time of apnea in hypertensive patients with OSA.

First report of pepper chlorosis-associated virus infecting tobacco (Nicotiana tabacum) plants in Sichuan Province in China.

Tobacco (Nicotiana tabacum) is an economically important crop in China, and more than 30 viruses have been reported to infect tobacco (Yin et al. 2022). In July 2022, we observed interveinal necrosis on tobacco leaves in fields in Sichuan Province (N 27.9172, E 105.6662) (Fig. 1). Total RNA was isolated from multiple leaves of one plant using an RNAprep Pure Polysaccharide Polyphenol Plant Total RNA Extraction Kit (TIANGEN, Beijing, China). Total RNAs were pooled, and a TruSeq Stranded Total RNA with RiboZero Gold Kit (Illumina, San Diego, CA, USA) was used to eliminate ribosomal RNA. An RNA-Seq library was constructed using VAHTS Universal V6 RNA-seq Library Prep (Nanjing Vazyme, China). High-throughput sequencing was performed on the Illumina DNBseq platform (BGI-ShenZhen, China), which yielded 20,102,087 reads with an average length of 150 nt (total size >6 Gb). Unaligned reads were assembled de novo using SPAdes (Bankevich et al. 2012). Contigs with length ≥200 nt were subjected to local BLASTn and BLASTx analyses against the GenBank nt and nr databases, respectively (Wang et al. 2022). A total of 23 contigs were identified through BLASTx (e-value cut-off = 10 -3), ranging from 631 to 1555 bp long, with 82% to 96% coverage to partial genomic sequences of pepper chlorosis-associated virus (PepCaV-Higashitsuno_2021; Accessions: LC719619 to LC719621) and one contig (6459 bp) with 99% similarity to tobacco mosaic virus (Accession: OP525281) isolate DSMZ PV-0109 from Germany. The complete genome sequence of PepCaV was obtained using primers based on the assembled contigs. The 5'- and 3'-terminal regions of the RNA genome were obtained by 5'- and 3'-rapid amplification of cDNA ends. These amplicons were cloned using the pEASY-Blunt Zero Cloning Kit (TRANSGEN, Nangjing, China) and sequenced by Sanger sequencing. Complete genome sequences of tripartite PepCaV from tobacco samples were 7697, 1808, and 1557 nucleotides long (Accession: OR451987 to OR451989) and showed genome organization typical of the genus Ophiovirus in the family Aspiviridae. The complete sequences of RNA1, RNA2 and RNA3 genome segments shared 92.36%, 90.43%, and 95.24%, nucleotide sequence identities, respectively, with the isolate PepCaV-Higashitsuno_2021 pepper isolate (Accession: LC719619 to LC719621) (Shimomoto et al. 2023), but PepCaV has not been reported to infect N. tabacum. In June 2023, 10 plants collected from each place of Macheng (N 27.9094, E 105.6740), Xiangyang (N 28.0936, E 105.6249) and Moni (N 27.8899, E 105.5936) showing interveinal necrosis symptoms were tested using RT-PCR using PepCaV-MP610-F (5'-TGTTCTCTGCTATGCGGTTG -3') and PepCaV-MP610-R (5'-AGCAATCTCGCACCTGAAGT-3') to product 610bp amplicon. Twenty-five tobacco plants were positive for PepCaV. Single sequence from each location was submitted to GenBank (Accession: PP728631 to PP728633). Sap extracts from the original field leaf samples collected from Sichuan Province were used to mechanically inoculate tobacco plants (10 plants) at the four-leaf stage. After 7 days, leaf samples were tested using RT-PCR assay specific to PepCaV and TMV while samples were positive only for TMV but failed to transmit PepCaV by mechanical inoculation. According to previous literature, ophioviruses may be transmitted though soilborne fungus (Jeong et al. 2014). Further research is needed to understand the transmission, epidemiology, and pathological properties of the PepCaV. To our knowledge, this is the first report documenting natural PepCaV infection of tobacco plants in China, providing a scientific basis for PepCaV infection control in tobacco plantations.