Design and synthesis of large Stokes shift DNA dyes with reduced genotoxicity.

Despite intensive search over the past decades, only a few small-molecule DNA fluorescent dyes were found with large Stokes shifts. These molecules, however, are often too toxic for widespread usage. Here, we designed DNA-specific fluorescent dyes rooted in benzimidazole architectures with a hitherto unexplored molecular framework based on thiazole-benzimidazole scaffolding. We further incorporated a pyrazole ring with an extended sidechain to prevent cell penetration. These novel benzimidazole derivatives were predicted by quantum calculations and subsequently validated to have large Stokes shifts ranging from 135 to 143 nm, with their emission colors changed from capri blue for the Hoechst reference compound to iguana green. These readily-synthesized compounds, which displayed improved DNA staining intensity and detection limits along with a complete loss of capability for cellular membrane permeation and negligible mutagenic effects as designed, offer a safer alternative to the existing high-performance small-molecule DNA fluorescent dyes.

MARS and RNAcmap3: The Master Database of All Possible RNA Sequences Integrated with RNAcmap for RNA Homology Search.

Recent success of AlphaFold2 in protein structure prediction relied heavily on co-evolutionary information derived from homologous protein sequences found in the huge, integrated database of protein sequences (Big Fantastic Database). In contrast, the existing nucleotide databases were not consolidated to facilitate wider and deeper homology search. Here, we built a comprehensive database by incorporating the non-coding RNA (ncRNA) sequences from RNAcentral, the transcriptome assembly and metagenome assembly from metagenomics RAST (MG-RAST), the genomic sequences from Genome Warehouse (GWH), and the genomic sequences from MGnify, in addition to the nucleotide (nt) database and its subsets in National Center of Biotechnology Information (NCBI). The resulting Master database of All possible RNA sequences (MARS) is 20-fold larger than NCBI's nt database or 60-fold larger than RNAcentral. The new dataset along with a new split-search strategy allows a substantial improvement in homology search over existing state-of-the-art techniques. It also yields more accurate and more sensitive multiple sequence alignments (MSAs) than manually curated MSAs from Rfam for the majority of structured RNAs mapped to Rfam. The results indicate that MARS coupled with the fully automatic homology search tool RNAcmap will be useful for improved structural and functional inference of ncRNAs and RNA language models based on MSAs. MARS is accessible at https://ngdc.cncb.ac.cn/omix/release/OMIX003037, and RNAcmap3 is accessible at http://zhouyq-lab.szbl.ac.cn/download/.

Examining the correlates of cigarette smoking, e-cigarette use and dual use among Canadian post-secondary students.

Many Canadians use nicotine products such as cigarettes and e-cigarettes. A particular subpopulation of concern is post-secondary students given they have a higher prevalence of use. Many correlates of cigarette smoking and e-cigarette use have been identified. However, less focus has been on examining the correlates of cigarette smoking, e-cigarette use and dual use. This study explores the correlates of different nicotine modality use in post-secondary students. Using data from the Canadian Campus Wellbeing Survey (CCWS; n = 27,164), a multi-level nominal regression assessed the correlates of nicotine modality use. In comparison to individuals who were <20, individuals 20-24 (OR = .448, 95% CI .321, .625), 25-29 (OR = .140, 95% CI .093, .212), 30-34 (OR = .076, 95% CI .046, .125) and over 35 (OR = .041, 95% CI .024, .071) had lower odds of e-cigarette use compared to cigarette smoking. Identifying as a woman (OR = 1.553, 95% CI 1.202, 2.006), non-heterosexual (OR = .642, 95% CI = .485,0.851), current cannabis user (OR = 1.651, 95% CI 1.296, 2.104), and being an international student (OR = .350, 95% CI .251, .487) also impacted the odds of e-cigarette use vs only cigarette smoking. When considering dual use vs cigarette smoking, individuals aged 20-24 (OR = .491, 95% CI .337, .717), 25-29 (OR = .221, 95% CI .137, .357), 30-34 (OR = .163, 95% CI .091, .292) and over 35 (OR = .122, 95% CI .065, .230) had lower odds than individuals <20. Current cannabis use (OR = 1.680, 95% CI = 1.209, 2.138), binge drinking (OR = 1.885, 95% CI 1.384, 2.568), and international student status (OR = .689, 95% CI .476, .996) also impacted cigarette smoking vs dual-use. Overall, a minority of young adults (11.5%) at post-secondary institutions in our sample use nicotine products, and the higher prevalence of e-cigarette use warrants continued monitoring. Health promotion campaigns addressing e-cigarette use are required. Additionally, tailored intervention efforts could prioritize the treatment needs of international students studying in Canada.

Epitope delimitation: A new method for defining epitopes of human IgG-reactive antigenic peptides based on rabbit-recognized epitope motifs.

The use of precise epitope peptides as antigens is essential for accurate serological diagnosis of viral-infected individuals, but now it remains an unsolvable problem for mapping precise B cell epitopes (BCEs) recognized by human serum. To address this challenge, we propose a novel epitope delimitation (ED) method to uncover BCEs in the delineated human IgG-reactive (HR) antigenic peptides (APs). Specifically, the method based on the rationale of similarities in humoral immune responses between mammalian species consists of a pair of elements: experimentally delineated HR-AP and rabbit-recognized (RR) BCE motif and corresponding pair of sequence alignment analysis. As a result of using the ED approach, after decoding four RR-epitomes of human papillomavirus types 16/18-E6 and E7 proteins utilizing rabbit serum against each recombinant protein and sequence alignment analysis of HR-APs and RR-BCEs, 19 fine BCEs in 17 of 22 known HR-APs were defined based on each corresponding RR-BCE motifs, including the type-specificity of each delimited BCE in homologous proteins. The test with 22 known 16/20mer HR-APs demonstrated that the ED method is effective and efficient, indicating that it can be used as an alternative method to the conventional identification of fine BCEs using overlapping 8mer peptides.

Multiple sequence alignment-based RNA language model and its application to structural inference.

Compared with proteins, DNA and RNA are more difficult languages to interpret because four-letter coded DNA/RNA sequences have less information content than 20-letter coded protein sequences. While BERT (Bidirectional Encoder Representations from Transformers)-like language models have been developed for RNA, they are ineffective at capturing the evolutionary information from homologous sequences because unlike proteins, RNA sequences are less conserved. Here, we have developed an unsupervised multiple sequence alignment-based RNA language model (RNA-MSM) by utilizing homologous sequences from an automatic pipeline, RNAcmap, as it can provide significantly more homologous sequences than manually annotated Rfam. We demonstrate that the resulting unsupervised, two-dimensional attention maps and one-dimensional embeddings from RNA-MSM contain structural information. In fact, they can be directly mapped with high accuracy to 2D base pairing probabilities and 1D solvent accessibilities, respectively. Further fine-tuning led to significantly improved performance on these two downstream tasks compared with existing state-of-the-art techniques including SPOT-RNA2 and RNAsnap2. By comparison, RNA-FM, a BERT-based RNA language model, performs worse than one-hot encoding with its embedding in base pair and solvent-accessible surface area prediction. We anticipate that the pre-trained RNA-MSM model can be fine-tuned on many other tasks related to RNA structure and function.

IDR-targeting compounds suppress HPV genome replication via disruption of phospho-BRD4 association with DNA damage response factors.

Compounds binding to the bromodomains of bromodomain and extra-terminal (BET) family proteins, particularly BRD4, are promising anticancer agents. Nevertheless, side effects and drug resistance pose significant obstacles in BET-based therapeutics development. Using high-throughput screening of a 200,000-compound library, we identified small molecules targeting a phosphorylated intrinsically disordered region (IDR) of BRD4 that inhibit phospho-BRD4 (pBRD4)-dependent human papillomavirus (HPV) genome replication in HPV-containing keratinocytes. Proteomic profiling identified two DNA damage response factors-53BP1 and BARD1-crucial for differentiation-associated HPV genome amplification. pBRD4-mediated recruitment of 53BP1 and BARD1 to the HPV origin of replication occurs in a spatiotemporal and BRD4 long (BRD4-L) and short (BRD4-S) isoform-specific manner. This recruitment is disrupted by phospho-IDR-targeting compounds with little perturbation of the global transcriptome and BRD4 chromatin landscape. The discovery of these protein-protein interaction inhibitors (PPIi) not only demonstrates the feasibility of developing PPIi against phospho-IDRs but also uncovers antiviral agents targeting an epigenetic regulator essential for virus-host interaction and cancer development.

Diversity and genetic characterization of orthohantavirus from small mammals and humans during 2012-2022 in Hubei Province, Central China.

Hemorrhagic fever with renal syndrome (HFRS) is a significant public health problem in Hubei Province, China, where a novel strain of orthohantavirus, HV004, was reported in 2012. However, no systematic study has investigated the prevalence and variation of orthohantavirus in rodents and humans. Herein, 2137 small mammals were collected from ten HFRS epidemic areas in Hubei Province from 2012 to 2022, and 143 serum samples from patients with suspected hemorrhagic fever were collected from two hospitals from 2017 to 2021. Orthohantavirus RNA was recovered from 134 lung tissue samples from five rodent species, with a 6.27 % prevalence, and orthohantavirus was detected in serum samples from 25 patients. Genetic analyses revealed that orthohantavirus hantanense (HTNV), orthohantavirus seoulense (SEOV), and orthohantavirus dabieshanense (DBSV) are co-circulating in rodents in Hubei, and HTNV and SEOV were identified in patient serum. Phylogenetic analysis showed that most of the HTNV sequences were clustered with HV004, indicating that HV004-like orthohantavirus was the main HNTV subtype in rodents. Two genetic reassortments and six recombination events were observed in Hubei orthohantaviruses. In summary, this study identified the diversity of orthohantaviruses circulating in Hubei over the past decade, with the HV004-like subtype being the main genotype in rodents and patients. These findings highlight the need for continued attention and focus on orthohantaviruses, especially concerning newly identified strains.

SPIN-CGNN: Improved fixed backbone protein design with contact map-based graph construction and contact graph neural network.

Recent advances in deep learning have significantly improved the ability to infer protein sequences directly from protein structures for the fix-backbone design. The methods have evolved from the early use of multi-layer perceptrons to convolutional neural networks, transformers, and graph neural networks (GNN). However, the conventional approach of constructing K-nearest-neighbors (KNN) graph for GNN has limited the utilization of edge information, which plays a critical role in network performance. Here we introduced SPIN-CGNN based on protein contact maps for nearest neighbors. Together with auxiliary edge updates and selective kernels, we found that SPIN-CGNN provided a comparable performance in refolding ability by AlphaFold2 to the current state-of-the-art techniques but a significant improvement over them in term of sequence recovery, perplexity, deviation from amino-acid compositions of native sequences, conservation of hydrophobic positions, and low complexity regions, according to the test by unseen structures, "hallucinated" structures and diffusion models. Results suggest that low complexity regions in the sequences designed by deep learning, for generated structures in particular, remain to be improved, when compared to the native sequences.

Serum sirtuin1: a potential blood biomarker for early diagnosis of Alzheimer's disease.

BACKGROUND: Sirtuin 1, a nicotinamide adenine dinucleotide-dependent deacetylase that is highly expressed in the hippocampus and anterior cortex tissues related to Alzheimer's Disease pathology, can cross the blood-brain barrier and is a promising biomarker. METHODS: A 1:1:1 case-control study was conducted and serum fasting blood glucose, triglyceride, total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, SIRT1, IL-6, Aß1-42, T-tau and P-tau-181 levels were evaluated in blood samples of 26 patients form the Alzheimer's Disease group, 26 patients form the mild cognitive impairment group, and 26 individuals form the normal control group. Receiver operator characteristic curves were used to evaluate the diagnostic significance. RESULTS: Serum SIRT1 level was significantly down-regulated in the mild cognitive impairment patients and Alzheimer's Disease patients compared with that in the normal control group (P<0.05). ROC curve analysis demonstrated that SIRT1 was a promising biomarker to distinguish Alzheimer's Disease patients from the mild cognitive impairment patients and the normal control group. In addition, SIRT1 was estimated to perform well in the diagnosis of Alzheimer's Disease ([AUC] = 0.742). CONCLUSIONS: In summary, the present study suggested that serum SIRT1 might be an early promising diagnostic biomarker for Alzheimer's Disease.

Metabolomics analysis of mycelial exudates provides insights into fungal antagonists of Armillaria.

The genus Armillaria has high edible and medical values, with zones of antagonism often occurring when different species are paired in culture on agar media, while the antagonism-induced metabolic alteration remains unclear. Here, the metabolome of mycelial exudates of two Chinese Armillaria biological species, C and G, co-cultured or cultured separately was analysed to discover the candidate biomarkers and the key metabolic pathways involved in Armillaria antagonists. A total of 2,377 metabolites were identified, mainly organic acids and derivatives, lipids and lipid-like molecules, and organoheterocyclic compounds. There were 248 and 142 differentially expressed metabolites between group C-G and C, C-G, and G, respectively, and fourteen common differentially expressed metabolites including malate, uracil, Leu-Gln-Arg, etc. Metabolic pathways like TCA cycle and pyrimidine metabolism were significantly affected by C-G co-culture. Additionally, 156 new metabolites (largely organic acids and derivatives) including 32 potential antifungal compounds, primarily enriched into biosynthesis of secondary metabolites pathways were identified in C-G co-culture mode. We concluded that malate and uracil could be used as the candidate biomarkers, and TCA cycle and pyrimidine metabolism were the key metabolic pathways involved in Armillaria antagonists. The metabolic changes revealed in this study provide insights into the mechanisms underlying fungal antagonists.

Glycosylated, Lipid-Binding, CDR-Like Domains of SARS-CoV-2 ORF8 Indicate Unique Sites of Immune Regulation.

The outbreak of the novel coronavirus SARS-CoV-2 has posed a significant threat to human health and the global economy since the end of 2019. Unfortunately, due to the virus's rapid evolution, preventingand controlling the epidemic remains challenging. The ORF8 protein is a unique accessory protein in SARS-CoV-2 that plays a crucial role in immune regulation, but its molecular details are still largely unknown. In this study, we successfully expressed SARS-CoV-2 ORF8 in mammalian cells and determined its structure using X-ray crystallography at a resolution of 2.3 Å. Our findings reveal several novel features of ORF8. We found that four pairs of disulfide bonds and glycosylation at residue N78 are essential for stabilizing ORF8's protein structure. Additionally, we identified a lipid-binding pocket and three functional loops that tend to form CDR-like domains that may interact with immune-related proteins to regulate the host immune system. On cellular experiments also demonstrated that glycosylation at N78 regulats of ORF8's ability to bind to monocytes cells. These novel features of ORF8 provide structural insights to into its immune-related function and may serve as new targets for developing ORF8-mediated immune regulation inhibitors. IMPORTANCE COVID-19, caused by the novel coronavirus SARS-CoV-2 virus, has triggered a global outbreak. The virus's continuous mutation increases its infectivity and may be directly related to the immune escape response of viral proteins. In this study, we used X-ray crystallography to determine the structure of SARS-CoV-2 ORF8 protein, a unique accessory protein expressed in mammalian cells, at a resolution of 2.3 Å. Our novel structure reveals important structure details that shed light on ORF8's involvement in immune regulation, including conservation disulfide bonds, a glycosylation site at N78, a lipid-binding pocket, and three functional loops that tend to form CDR-like domains that may interact with immune-related proteins to modulate the host immune system. We also conducted preliminary validation experiments on immune cells. These new insights into ORF8's structure and function provide potential targets for developing inhibitors to block the ORF8-mediated immune regulation between viral protein and host, ultimately contributing to the development of novel therapeutics for COVID-19.

HNF4G stimulates the development of pancreatic cancer by promoting IGF2BP2 transcription

Objective Pancreatic cancer is a devastating and lethal malignancy. Our study investigated the effective mechanism of HNF4G on pancreatic cancer cell functions through the IGF2BP2 transcription. Methods HNF4G and IGF2BP2 expressions in pancreatic cancer were examined. The relationship between HNF4G expression and pancreatic cancer patients’ clinicopathological characteristics was evaluated. After interfering with HNF4G expression in pancreatic cancer cells, the cell proliferative, migratory, and invasive capabilities were evaluated. Also, the expression of proliferation-related gene PCNA and migration and invasion-related gene MMP2 was determined. The binding relation between HNF4G and HNF4G promoter was forecasted and testified. A tumorigenesis assay in nude mice was performed to detect the HNF4G interference’s effect on the subcutaneous tumorigenic capacity of pancreatic cancer cells. Results HNF4G and IGF2BP2 expressions were up-regulated in pancreatic cancer. Specifically, interfering with HNF4G inhibited PANC-1 cell proliferative, invasive and migratory behaviors, and decreased PCNA and MMP2 expression. Mechanistically, HNF4G as a transcription factor could specifically bind to IGF2BP2 and promote its expression. Rescue assay findings showed that IGF2BP2 overexpression could reverse the inhibiting effect of HNF4G interference on pancreatic cancer cells. For the in vivo finding, interfering HNF4G expression retarded the subcutaneous tumorigenic ability of pancreatic cancer cells. Conclusion We summarize that HNF4G as a transcription factor regulates IGF2BP2 expression to promote pancreatic cancer cell proliferation and migration capacities (AU)

Consistency of the anesthesia consciousness index versus the bispectral index during laparoscopic gastrointestinal surgery with sevoflurane anesthesia: A prospective multi-center randomized controlled clinical study.

Background: This study aimed to compare the consistency of anesthesia consciousness index (Ai) with that of bispectral index (BIS) in monitoring the depth of anesthesia (DOA) during sevoflurane anesthesia, to reveal the optimal cutoff values in different states of consciousness, and explore the stability of DOA monitoring during intraoperative injurious stimulation. Methods: We enrolled 145 patients (97 men and 48 women) from 10 medical centers. General anesthesia was induced using intravenous anesthetics and maintained with sevoflurane. Ai and BIS values were recorded. Results: The mean difference between the Ai and BIS was-0.1747 (95% confidence interval, -0.6660 to 0.3166; p = 0.4857). The regression equation of Ai and BIS from the Deming regression analysis was y = 5.6387 + 0.9067x (y is BIS, x is Ai), and the slope and intercept were statistically significant. Meanwhile, the receiver operating characteristic curve analysis of anesthesia-induced unconsciousness, loss of consciousness, and recovery of consciousness revealed that the accuracy of Ai and BIS were similar. In addition, the optimal cutoff values of the different states of consciousness were not sensitive to age, and both Ai and BIS had no correlation with hemodynamics. Conclusion: We conclude that Ai and BIS show no systematic deviation in readings with high consistency, similar accuracy, and good stability; these insights provide more data for clinical application.

Inhibiting the transcription and replication of Ebola viruses by disrupting the nucleoprotein and VP30 protein interaction with small molecules.

Ebola virus (EBOV) causes hemorrhagic fever in humans with high morbidity and fatality. Although over 45 years have passed since the first EBOV outbreak, small molecule drugs are not yet available. Ebola viral protein VP30 is a unique RNA synthesis cofactor, and the VP30/NP interaction plays a critical role in initiating the transcription and propagation of EBOV. Here, we designed a high-throughput screening technique based on a competitive binding assay to bind VP30 between an NP-derived peptide and a chemical compound. By screening a library of 8004 compounds, we obtained two lead compounds, Embelin and Kobe2602. The binding of these compounds to the VP30-NP interface was validated by dose-dependent competitive binding assay, surface plasmon resonance, and thermal shift assay. Moreover, the compounds were confirmed to inhibit the transcription and replication of the Ebola genome by a minigenome assay. Similar results were obtained for their two respective analogs (8-gingerol and Kobe0065). Interestingly, these two structurally different molecules exhibit synergistic binding to the VP30/NP interface. The antiviral efficacy (EC50) increased from 1 µM by Kobe0065 alone to 351 nM when Kobe0065 and Embelin were combined in a 4:1 ratio. The synergistic anti-EBOV effect provides a strong incentive for further developing these lead compounds in future studies.

HNF4G stimulates the development of pancreatic cancer by promoting IGF2BP2 transcription.

OBJECTIVE: Pancreatic cancer is a devastating and lethal malignancy. Our study investigated the effective mechanism of HNF4G on pancreatic cancer cell functions through the IGF2BP2 transcription. METHODS: HNF4G and IGF2BP2 expressions in pancreatic cancer were examined. The relationship between HNF4G expression and pancreatic cancer patients' clinicopathological characteristics was evaluated. After interfering with HNF4G expression in pancreatic cancer cells, the cell proliferative, migratory, and invasive capabilities were evaluated. Also, the expression of proliferation-related gene PCNA and migration and invasion-related gene MMP2 was determined. The binding relation between HNF4G and HNF4G promoter was forecasted and testified. A tumorigenesis assay in nude mice was performed to detect the HNF4G interference's effect on the subcutaneous tumorigenic capacity of pancreatic cancer cells. RESULTS: HNF4G and IGF2BP2 expressions were up-regulated in pancreatic cancer. Specifically, interfering with HNF4G inhibited PANC-1 cell proliferative, invasive and migratory behaviors, and decreased PCNA and MMP2 expression. Mechanistically, HNF4G as a transcription factor could specifically bind to IGF2BP2 and promote its expression. Rescue assay findings showed that IGF2BP2 overexpression could reverse the inhibiting effect of HNF4G interference on pancreatic cancer cells. For the in vivo finding, interfering HNF4G expression retarded the subcutaneous tumorigenic ability of pancreatic cancer cells. CONCLUSION: We summarize that HNF4G as a transcription factor regulates IGF2BP2 expression to promote pancreatic cancer cell proliferation and migration capacities.

High-throughput mapping of RNA solvent accessibility at the single-nucleotide resolution by RtcB ligation between a fixed 5'-OH-end linker and unique 3'-P-end fragments from hydroxyl radical cleavage.

Given the challenges for the experimental determination of RNA tertiary structures, probing solvent accessibility has become increasingly important to gain functional insights. Among various chemical probes developed, backbone-cleaving hydroxyl radical is the only one that can provide unbiased detection of all accessible nucleotides. However, the readouts have been based on reverse transcription (RT) stop at the cleaving sites, which are prone to false positives due to PCR amplification bias, early drop-off of reverse transcriptase, and the use of random primers in RT reaction. Here, we introduced a fixed-primer method called RL-Seq by performing RtcB Ligation (RL) between a fixed 5'-OH-end linker and unique 3'-P-end fragments from hydroxyl radical cleavage prior to high-throughput sequencing. The application of this method to E. coli ribosomes confirmed its ability to accurately probe solvent accessibility with high sensitivity (low required sequencing depth) and accuracy (strong correlation to structure-derived values) at the single-nucleotide resolution. Moreover, a near-perfect correlation was found between the experiments with and without using unique molecular identifiers, indicating negligible PCR biases in RL-Seq. Further improvement of RL-Seq and its potential transcriptome-wide applications are discussed.

Identification of the Immune Status of Hypertrophic Cardiomyopathy by Integrated Analysis of Bulk- and Single-Cell RNA Sequencing Data.

Objectives: Hypertrophic cardiomyopathy (HCM) is the most common hereditary cardiomyopathy and immune infiltration is considered an indispensable factor involved in its pathogenesis. In this study, we attempted to combine bulk sequencing and single-cell sequencing to map the immune infiltration-related genes in hypertrophic cardiomyopathy. Methods: The GSE36961, GSE160997, and GSE122930 datasets were obtained from the Gene Expression Omnibus database. The compositional patterns of the 18 types of immune cell fraction and pathway enrichment score in control and HCM patients were estimated based on the GSE36961 cohort using xCell algorithm. The Weighted Gene Coexpression Network Analysis (WGCNA) was performed to identify genes associated with immune infiltration for hypertrophic cardiomyopathy. The area under the curve (AUC) value was obtained and used to evaluate the discriminatory ability of common immune-related DEGs. "NetworkAnalyst" platform was used to identify TF-gene and TF-miRNA interaction with identified common genes. Heat map was used to determine the association between common DEGs and various immune cells. Results: Immune infiltration analysis by the xCell algorithm showed a higher level of CD8+ naive T cells, CD8+ T cells, as well as a lower level of activated dendritic cells (aDC), dendritic cells (DC), immature dendritic cells (iDC), conventional dendritic cells (cDC), macrophages, M1 macrophages, monocytes, and NKT cell in HCM compared with the control group in GSE36961 dataset. aDC, macrophages, and M1 macrophages were the top three discriminators between HCM and control groups with the area under the curve (AUC) of 0.907, 0.867, and 0.941. WGCNA analysis showed that 1258 immune-related genes were included in four different modules. Of these modules, the turquoise module showed a pivotal correlation with HCM. 13 common immune-related DEGs were found by intersecting common DEGs in GSE36961 and GSE160997 datasets with genes from the genes in turquoise module. 5 hub immune-related genes (S100A9, TYROBP, FCER1G, CD14, and S100A8) were identified by protein interaction network. Through analysis of single-cell sequencing data, S100a9, TYROBP, FCER1G, and S100a8 were mainly expressed by infiltrated M1 proinflammatory cells, especially Ccr2-M1 proinflammatory macrophage cells in the heart immune microenvironment while Cd14 was expressed by infiltrated M1 proinflammatory macrophage cells and M2 macrophages in transverse aortic constriction (TAC) mice at 1 week. Higher M2 macrophage and M1 proinflammatory macrophage infiltration as well as lower Ccr2-M1 proinflammatory macrophage and dendritic cells were shown in TAC 1week mice compared with sham mice. Conclusions: There was a difference in immune infiltration between HCM patients and normal groups. aDC, macrophages, and M1 macrophages were the top three discriminator immune cell subsets between HCM and control groups. S100A9, TYROBP, FCER1G, CD14, and S100A8 were identified as potential biomarkers to discriminate HCM from the control group. S100a9, TYROBP, FCER1G, and S100a8 were mainly expressed by infiltrated M1 proinflammatory cells, especially Ccr2-M1 proinflammatory cells in the heart immune microenvironment while Cd14 was expressed by M2 macrophages in transverse aortic constriction (TAC) mice at 1 week.