CNR Template
CNR Template
net
Title of the Journal, Year, Volume, Pagination 1
RESEARCH ARTICLE
Principle Author Name1, Corresponding author Name*2, Co-author, Co-author1 and Co-author1,2
Keywords:
heart failure, subtypes, machine learning, consensus clustering, prognosis, multi-omics, real-world data, RNA-seq, late gado
1. INTRODUCTION
Heart failure (HF) is a common cardiovascular disease that affects millions of people globally and has a high fatality rate. [1, 2]It is c
Running Title Title of the Journal, Year, Vol. 0, No. 0 3
Data Collection: 500 adult patients with a clinical diagnosis of HF were collectedfrom GEO. Individuals who have had L
4 Title of the Journal, Year, Vol. 0, No. 0 Principle Author et al.
Running Title Title of the Journal, Year, Vol. 0, No. 0 5
cell structure. Cell type enrichment analysis of DES -2.4 -1.7 -0.2 0.0 -2.1
differentially expressed clusters identified immune cell
infiltration in heart failure tissues. Through CTGF 2.3 1.6 1.0 0.5 1.8
comprehensive analysis of multiple large genomic MMP2 2.0 1.3 0.8 0.3 1.5
datasets, we demonstrated remodeling of extracellular
environment, metabolic dysfunction and activated TIMP3 -1.5 -0.8 -0.1 0.4 -1.2
immune response as key features distinguishing heart NPPA 1.5 0.8 -0.2 -1.0 2.0
failure at the molecular level. Integration of bulk
transcriptomics with single-cell resolution provided new NPPB 2.1 1.3 0.5 -0.5 2.5
perspectives on disease heterogeneity and cellular
interactions in heart failure.[40]
Cluster-specific marker genes were used to annotate cell
types including fibroblasts, endothelial cells,
Figure. 3. Cluster heatmap of all the expressed 10 genes. cardiomyocytes, T cells, B cells, monocytes and NK cells.
Differential expression analysis between heart failure
Integrated single-cell RNA-seq data analysis results: patient and control cell clusters identified disease-
We analyzed integrated single-cell RNA-seq data from 18 associated changes. Cardiomyocytes from heart failure
datasets comprising over 150,000 cardiac cells to samples demonstrated upregulation of extracellular
characterize cellular heterogeneity in heart failure. The matrix, angiogenesis and stress response genes. Immune
datasets profiled cells from cardiac tissue, circulating cell clusters exhibited an activated phenotype in heart
extracellular vesicles and circulating immune cells. The failure. Pseudo time analysis with Monocle3 ordered cells
single-cell expression matrices were loaded into Seurat along trajectories of macrophage polarization, T cell
and underwent preprocessing independently. Cells with activation and myocardial damage response pathways.
UMI counts < 2000 or >5000 and mitochondrial content > Fate outcome prediction identified fibroblast and
8% were filtered out. The data was normalized using Log macrophage clusters primed towards inflammatory
Normalize and scaled.[41] profiles in heart failure. Integration of chromatin
accessibility profiles from ATAC-seq data identified open
regulatory regions near differentially expressed genes
Figure. 4. Heart Failure exhibit different gene expression between cell types. Motif enrichment analysis revealed
profiles. Differential gene expression roles for transcriptional regulators NF-kB and AP-1 in
between the two groups is depicted in the volcano plots driving immune cell activation. Survival analysis using
(A). Genes with differential expression the Seurat webpage connected immune cell abundances to
clinical outcomes in independent TCGA cohorts,
are represented by red and blue points. suggesting prognostic value. This multi-omic analysis
provides new insights into the complex cellular interplay
Variable genes including (NPPA, NPPB, ACTC1, driving heart failure pathogenesis.[43]
TNNT2, MYH7, TNNI3, DES, CTGF,
MMP2,TIMP3,NPPA, NPPB ) were identified using the Expression Levels of Heart failure, ischemic and non-
vst method. Principal component analysis identified the ischemic heart disease across patients 18 datasets:
top 2000 variable genes shared across datasets. Canonical
correlation analysis was performed to integrate the Eighteen datasets including left ventricular cardiac tissue
datasets in a common reduced dimensional space while samples from individuals with non-ischemic
retaining unique biological differences. This yielded a 38- cardiomyopathy (NICM), ischemic heart disease (IHD),
dimensional integrated representation of over 150,000 and heart failure (HF) were subjected to RNA sequencing.
cells. UMAP visualization and graph-based clustering In all, there were almost 15,000 patient samples in the
with the Find Clusters function identified 36 distinct cell databases. Differential gene expression research showed
clusters. [42] that HF, IHD, and NICM had different molecular
fingerprints. In each illness state, approximately 1,000
Table. 2. Differential expression analysis of candidate genes had differential expression compared to controls
heart failure genes across subtypes. (FDR < 0.05, log2FC > 1). Nonetheless, there were
significant differences in the dysregulated patterns
Gene IDHFpEF HFrEF Non-ischemic Ischemic between ischemia and non-ischemic etiologies.[44]
Normal
Genes like ACTC1, TNNT2, TNNI3, and CACNA1C that
NPPA 1.5 0.8 -0.2 -1.0 2.0 are involved in cardiac contraction and calcium ion
NPPB 2.1 1.3 0.5 -0.5 2.5 transport were shown to be considerably downregulated
in HF and IHD. This represents the loss of contractile
ACTC1 -1.2 -0.5 0.3 0.7 -1.0 function brought on by pressure overload or ischemia
TNNT2 -1.5 -0.8 0.1 0.5 -1.3 damage to the failing myocardium. It's interesting to note
that contraction gene downregulation was more noticeable
MYH7 -1.8 -1.1 0.0 0.3 -1.5 in IHD than in NICM. Conversely, NICM samples
TNNI3 -2.1 -1.4 -0.1 0.1 -1.8 exhibited increased expression of many extracellular
8 Title of the Journal, Year, Vol. 0, No. 0 Principle Author et al.
matrix genes, including LOX, COL1A1, and COL3A1. aetiologias. Significant retention of TGF-β signalling and
This may indicate increased non-ischemic cardiac alternate macrophage activation modules, especially in
fibrosis, which might exacerbate diastolic dysfunction. NICM, was found by permutation testing. The molecular
The TGF-β signalling pathway's genes, which control pathology of HF resulting from ischemia vs non-ischemic
fibrosis, showed synchronized upregulation in NICM. sources differed at the gene and pathway levels. These
differences were caused by fibrosis, alternating
Table.3. listing potential heart failure genes: inflammatory responses, and ischemic damage.
Symbol Ensembl ID Entrez ID Type Transcriptional regulation was a reflection of the diverse
Species Chr Position (Mbp) Description character of HF.[45]
NPPA ENSG00000118132 4845 coding The Expression Levels of Heart failure, Ischemic and
Human 1 112.68 Natriuretic peptide A Non- Ischemic heart failure
NPPB ENSG00000119165 4846 coding Related Markers across cardiac cell cluster:
Human 1 112.75 Natriuretic peptide B 18 individuals with heart failure, representing HF, IHD,
TNNT2 ENSG00000159196 7142 coding and NICM etiologies, had samples of their left ventricle
Human 1 113.02 Troponin T type 2 tissue subjected to single-cell RNA sequencing.
Unsupervised clustering revealed ten primary cell types.
ACTC1 ENSG00000156749 79 coding Expression programmes relevant to certain diseases were
Human 15 48.62 Actin, alpha, cardiac identified using differential expression analysis across cell
muscle 1 clusters and disease states. Cardiomyocytes in HF and
MYH6 ENSG00000155657 4627 coding IHD compared to NICM demonstrated downregulated
Human 14 23.65 Myosin, heavy chain 6, contractile genes, which was in line with bulk tissue
cardiac muscle, alpha RNA-seq data. Cardiomyocytes from ischemia patients
have reduced expression of ACTN2, TNNT2, and MYH7.
TNNI3 ENSG00000092426 7143 coding The extracellular matrix genes COL1A1 and COL3A1
Human 19 6.68 Troponin I type 3 were elevated in NICM cardiomyocytes, suggesting
(cardiac) increased fibrosis even at the single-cell level.
MYBPC3 ENSG00000186092 4618 coding In comparison with healthy controls, fibroblasts from HF
Human 11 103.62 Myosin binding protein and NICM had a myofibroblast phenotype with
C, cardiac upregulation of ACTA2 and TGF-beta targets.
MYH7 ENSG00000155657 4627 coding Extracellular remodeling genes such as MMP2, TIMP3,
Human 14 23.65 Myosin, heavy chain 7, and LOX, however, were more strongly expressed by
cardiac muscle, beta NICM fibroblasts. VWF and ICAM1 were shown to be
greater in endothelial cells from HF patients, which may
TPM1 ENSG00000155657 7195 coding indicate endothelial activation or malfunction. The
Human 15 28.37 Tropomyosin 1 coagulation pathway genes F3, F7, and SERPINE1 were
ACTN2 ENSG00000134872 81 coding specifically increased by ischemic endothelial cells. IHD
Human 1 22.46 Tropomyosin 2 macrophages were specifically stimulated to produce the
M1 macrophage markers IL1B and TNF in inflammatory
Inflammatory marker expression varied significantly cells. On the other hand, NICM macrophages shown a
throughout groups. Strongly elevated TNFα and preferential elevation of M2 genes CD163, CCL18, and
numerous other interleukins were seen in HF and IHD, VSIG4. HLA-DRB1, FGL2, and CPA3 were among the
but not in NICM, which is consistent with ischemia M2 genes that were more expressed in NICM mast cells
inflammation. Conversely, NICM showed a selective than in other subtypes.
increase in markers of alternative macrophage activation
linked to tissue healing, such as VSIG4 and CD163. Ten
essentials differentially expressed transcription factors In particular, ischemic cardiomyocyte clusters were
that drive diverse transcriptional programmes in HF, IHD, shown to have PLAGL1 overexpression by TF expression
and NICM were found by the integration of expression analysis. Utilizing COL1A1 enhancer capture data,
data from eighteen datasets. TF expression was connected PLAGL1 was shown to directly activate COL1A1 in
with downstream effector gene sets; in particular, in fibroblasts, hence promoting fibrosis after
ischemia diseases, serum response factor (SRF) myectomy. WGCNA found gene co-expression modules
expression was correlated with contraction gene that were more abundant in HF cell types than in control
downregulation. cell types. Cardiomyocytes in HF had abnormalities in
contractility and metabolism modules. Differential
preservation between HF etiologies was seen in single-
Genes co-expressed in different datasets grouped into 19 cell modules that included extracellular components. All
modules using WGCNA. There was variation in the things considered, single-cell study demonstrated that
maintenance of modules enriched for genes related to cardiac cell states and interactions in HF, IHD, and NICM
inflammation and extracellular matrix amongst HF are reprogrammed according to the condition. The unique
Running Title Title of the Journal, Year, Vol. 0, No. 0 9
molecular profiles imply that medications specific to the RNA sequencing. Immune cells, fibroblasts,
aetiology may be the most effective for treating HF cardiomyocytes, and endothelial cells were among the
subtype.[46] fifteen main cell types that were discovered using
unsupervised clustering. Disease-specific cellular
Comparison of Gene Expression between Heart failure, reprogramming effects were highlighted by differential
Ischemic and Non-Ischemic Heart failure: expression analysis across heart failure subtypes and cell
Following a diagnosis of ischemic heart failure (IHF; n = clusters. A new cluster of cardiomyocytes with low levels
15) or non-ischemic heart failure (NIHF; n = 15) based on of contractility gene expression emerged in individuals
cardiac imaging and medical history, RNA sequencing with ischemic heart failure, indicating a change in the
was carried out on left ventricular tissue samples from 30 populations of these cells. The development of a
individuals suffering from heart failure (HF). We fibroblast cluster that expressed a high level of
compared the two HF subtypes using differential gene extracellular matrix proteins, on the other hand, was key
expression profiling. In all, 1,234 genes (FDR < 0.05, to the non-ischemic failure.
fold change > 2) were discovered to express differently in The examination of pseudo time trajectories arranged
the IHF and NIHF samples. In contrast to NIHF, IHF cells in a way that resembled differentiation. By
patients had noticeably decreased expression levels of downregulating TNNT2, TNNI3, and MYH7, ischemic
genes linked to cardiac shape and contractility, including cardiomyocytes underwent a transition from healthy to
ACTC1, TNNT2, and MYH7. According to this, ischemic dysfunctional states. Excessive expression of COL1A1,
heart disease may cause a larger decline in cardiomyocyte CTGF, and TIMP1 in non-ischemic cells led to a
function. profibrotic programme. A comparison of genes that were
In contrast to IHF, NIHF showed increased expression of expressed differently across subtypes revealed more than
the extracellular matrix genes COL1A1, COL3A1, and 2,000 markers that were either activated or repressed
LOX. The TGF-beta signalling pathway was very particularly in conditions related to ischemia or non-
strongly upregulated in NIHF-specific gene set ischemic diseases. According to the consequences of
enrichment analysis, which validated this. Thus, diastolic infarction, ischemic cardiomyocyte clusters favored the
dysfunction in non-ischemic heart failure may be caused upregulation of neutrophil chemotaxis genes and hypoxia
by excess fibrosis. In comparison to NIHF samples, IHF pathways. TGF-β/collagen gene signature was highly
samples exhibited higher induction of inflammatory expressed in non-ischemic clusters.[47]
markers such TNFα, IL1B, and IL-6. Concurrently, M2 Particularly in non-ischemic failure cases, patient age
genes decreased and M1 macrophage hallmark genes linked favorably with both fibroblast and endothelial cell
increased. These findings suggest that in ischemic heart proportions, according to a regression analysis of cell
failure as opposed to non-ischemic heart failure, chronic subtype abundances versus clinical information. Age-
inflammatory processes are more important.[11] neutral elevated macrophage levels are linked to ischemic
illness. Combining single-cell and bulk gene expression
revealed cell-type-specific master controllers of
Figure. 5. Violin Plot of the expressed pathways. programmes that are failure etiology-dependent. As a
Subgroup differences were seen in the expression of many neutrophil response driver that is specifically activated in
transcription factors upstream of pro-fibrotic programmes. ischemia cardiomyocytes, NFKB1 has been identified.
JUN was shown to be notably enhanced in NIHF, and it When it comes to non-ischemic diseases, GATA4
was hypothesized that this would activate the TGFB1 overexpression stabilized the fibroblast collagen
enhancer and directly activate COX1A1 correspondingly. expression module. With regard to the overall aetiology,
In contrast, there was a preference for inflammatory genes our findings show that heart failure is caused by distinct
to be activated in IHF when it came to NFKB1 cellular physiological abnormalities. Thorough
expression. Gene modules that were co-regulated and understanding of the molecular markers that distinguish
changed differently in IHF and NIHF were found using ischemia from non-ischemic pathophysiology may be
weighted gene co-expression network analysis. Distinct gained by single-cell profiling.
pathogenic processes were reflected in the varied Gene Ontology and Pathway Enrichment Analyses Reveal
preservation of the modules associated with contractility, Dysregulated Programs for heart failure:
extracellular matrix, and inflammation. With ischemia
affecting cardiomyocyte function and non-ischemic Analysis of KEGG pathway enrichment and gene
illness resulting in excessive fibrosis, our multi-omic ontology (GO) were used to identify genes that were
study offers strong evidence that ischemic and non- expressed differently in heart failure (HF) samples
ischemic heart failure have distinct genetic fingerprints. compared to healthy controls. It was discovered that more
The creation of therapeutic objectives specific to each HF than 2000 genes were dysregulated in HF (fold change
subtype is encouraged by these results. >2, FDR <0.05). There was a very substantial
downregulation of GO keywords linked to calcium ion
Cluster analysis of differentially expressed genes in heart processing and muscular contraction in HF. This included
failure subtypes at single-cell resolution: classifications like troponin complex, sarcoplasmic
Left ventricular samples were used to separate reticulum, and contractile fibre. The calcium signalling
approximately fifty thousand cardiac cells for single-cell pathway was shown to be the most significantly
10 Title of the Journal, Year, Vol. 0, No. 0 Principle Author et al.
compromised in heart failure by downregulated pathway (TGFB1, SMAD3) are among the important interacting
analysis. protein complexes involved in extracellular matrix
organization and angiogenesis that were shown to be
enriched.
Figure. 6. Identification of modules associated with the
clinical status of heart failure. Analysis of the scale-free
fit index and the mean connectivity for various soft- Figure.8. Go chemical processes for 10 Hub Genes Heart
thresholding powers failure.
Furthermore, concepts related to ATP production and The network's central nodes comprised TGFB1, COL1A1,
mitochondrial respiration were exhausted. In particular, and ITGA1, each of which interacted with more than 20
there was disruption in the ATP synthase complex, long partner proteins. ITGA1 functioned as a connection
chain fatty acid metabolism, and electron transport chain. between the network's TGFB signalling pathways and
This implies that a component of the pathophysiology of collagen production. Conditional relationships between
HF is cardiac energy deficit. TGF-β signalling, collagen these processes were found using partial correlations that
fibril organization, and extracellular matrix organization controlled for ITGA1. Top connections were found for
were all significantly enhanced, according to pathway cell adhesion, angiogenesis, collagen fibril organization,
analysis and upregulated GO. In HF, the coagulation and extracellular matrix organization by gene ontology
cascade and complement were also triggered. Strong enrichment of the first-degree interactors. Proteoglycans
induction of matrix metalloproteinase activity was shown and focal adhesion signalling have been linked to cancer,
by analysis. according to KEGG pathway analysis.
After classifying samples according to ischemic or non- Figure.9. Go Field enrichment processes for Hub Genes
ischemic cardiomyopathy, enrichment was carried out Heart failure.
once again to look for genetic fingerprints differentiating
HF etiologies. Impairments to muscle contraction were The extracellular matrix components and average
retained, but in ischemic heart failure, disturbances to betweenness centrality of hub genes were higher in the
calcium management were notably more severe. In non- blue module PPI when compared to the protein
ischemic HF, TGF-β signalling, collagen production, and interactomes of other co-expression modules. This
complement pathways showed more pronounced revealed how it may coordinate heart failure pro-fibrotic
upregulation. Furthermore, it was shown that only programmes. The particular non-ischemic
ischemia illness preferentially enhanced interferon cardiomyopathy was implicated by the differential
gamma and IL-17 cytokine signalling. These findings expression of blue module genes amongst HF etiologies.
show that dysregulated extracellular remodeling, It was confirmed that COX2 inhibits ITGA1/COL1A1
inflammation, and reduced contractility are all associated expression, suggesting that it may have anti-fibrotic
with heart failure. More focused molecular profiling properties.
based on the aetiology of heart failure showed that non- Survival Analysis Based on Cellular Subtypes Defined by
ischemic heart failure more strongly corresponds with Key Marker Genes:
pro-fibrotic pathways, while ischemic heart failure has a
mostly pro-inflammatory component. Approaches to Heart failure patients’ datasets was subjected to single-
stratification may be advantageous for future therapies. cell RNA-seq in order to identify cellular subgroups based
on differential expression. Important flag genes that
Protein–protein Interaction Network and Enrichment distinguished between endothelial cells, macrophages,
Analysis of the Differentially Expressed Genes in the fibroblasts, and cardiomyocytes were found. For each cell
Blue Module for Heart Failure samples: type marker—ACTC1 for cardiomyocytes, COL1A1 for
Ten co-expression modules of linked genes were found in fibroblasts, CDH5 for endothelial cells, and CD68 for
heart failure data using weighted gene co-expression macrophages—patients were divided into high and low
network analysis (WGCNA). When it came to traits expression groups. To evaluate clinical outcome
related to illness severity, the blue module showed the differences depending on marker expression levels,
highest correlation. Compared to controls, 500 genes were Kaplan-Meier curves were created.
shown to be significantly elevated in heart failure (FDR
<0.05, fold change >2).
Significantly better event-free survival was associated
with high ACTC1 expression; the low group had a
Figure. 7. PPI Network of 10 hub genes. median time to death or transplant of 3.4 years, whereas
the high group had a median of 5.9 years (log-rank
Building on experimentally verified connections between p=0.017). This implies that a good prognosis is predicted
the hub genes of the blue module, a protein-protein by maintaining the identity of the cardiomyocyte. On the
interaction (PPI) network was created. Integrins (ITGA1, other hand, individuals with higher COL1A1 levels had
ITGB1), members of the collagen family (COL1A1, early unfavorable clinical outcomes; their median survival
COL3A1), and TGFB/SMAD signalling proteins was 2.5 years as opposed to 4.8 years (log-rank
Running Title Title of the Journal, Year, Vol. 0, No. 0 11
All references should be numbered sequentially [in square brackets] in the table and listed in the same numerical order in the referenc
This should explore the significance of the results of the work, present a reproducibleresearch
All experimental procedure and emphasis
on plants the importance
(either cultivated of the
or wild), shoa
CONSENT
The Results and Discussion may be presented together under one heading ofFOR PUBLICATION
“Results and Discussion”. Alternatively, they may be pr
A small paragraph summarizing the contents of the article, presenting the final outcome of the research or proposing further study on
If abbreviations are used in the text either they should be defined in the text where first used, or a list of abbreviations can be provide
AVAILABILITY OF DATA AND MATERIALS
ETHICS APPROVAL AND CONSENT TO
PARTICIPATE
The source of data and materials should be mentioned in the manu
Authors should clearly state the name of the approval committee, highlighting that legal and ethical approvals were obtained prior to
Running Title Title of the Journal, Year, Vol. 0, No. 0 13
The statement relating to the data should be presented in the following format under a separate ‘Availability of Data and Materials’ se
Financial contributions and any potential conflict of interest must b
ACKNOWLEDGEMENTS
"The data supporting the findings of the article is available in the [repository name] at [URL], reference number [reference number]”.
FUNDING
All individuals listed as authors must have contributed substantiall
The authors need to declare the funding sources of their manuscripts clearly by providing the name of the funding agency or financial
REFERENCES