Identification of the TIFAB Gene as a Susceptibility Locus for Coronary Artery Aneurysm in Patients with Kawasaki Disease

Sex-Specific Susceptibility Loci Associated With Coronary Artery Aneurysms in Patients With Kawasaki Disease

BACKGROUND AND OBJECTIVES

Kawasaki disease (KD) is an acute vasculitis that primarily affects children under age 5 years. Approximately 20-25% of untreated children with KD and 3-5% of those treated with intravenous immunoglobulin therapy develop coronary artery aneurysms (CAAs). The prevalence of CAAs is much higher in male than in female patients with KD, but the underlying factors contributing to susceptibility to CAAs in patients with KD remain unclear. This study aimed to identify sex-specific susceptibility loci associated with CAAs in KD patients.

METHODS

A sex-stratified genome-wide association study (GWAS) was performed using previously obtained GWAS data from 296 KD patients and a new replication study in an independent set of 976 KD patients by comparing KD patients without CAA (controls) and KD patients with aneurysms (internal diameter ≥5 mm) (cases).

RESULTS

Six male-specific susceptibility loci, PDE1C, NOS3, DLG2, CPNE8, FUNDC1, and GABRQ (odds ratios [ORs], 2.25-9.98; p=0.00204-1.96×10-6), and 2 female-specific susceptibility loci, SMAD3 (OR, 4.59; p=0.00016) and IL1RAPL1 (OR, 4.35; p=0.00026), were significantly associated with CAAs in patients with KD. In addition, the numbers of CAA risk alleles additively contributed to the development of CAAs in patients with KD.

CONCLUSIONS

A sex-stratified GWAS identified 6 male-specific (PDE1C, NOS3, DLG2, CPNE8, FUNDC1, and GABRQ) and 2 female-specific (SMAD3 and IL1RAPL1) CAA susceptibility loci in patients with KD.

Site-specific genetic and functional signatures of aortic endothelial cells at aneurysm predilection sites in healthy and AngII ApoE-/- mice

Aortic aneurysm is characterized by a pathological dilation at specific predilection sites of the vessel and potentially results in life-threatening vascular rupture. Herein, we established a modified "Häutchen method" for the local isolation of endothelial cells (ECs) from mouse aorta to analyze their spatial heterogeneity and potential role in site-specific disease development. When we compared ECs from aneurysm predilection sites of healthy mice with adjacent control segments we found regulation of genes related to extracellular matrix remodeling, angiogenesis and inflammation, all pathways playing a critical role in aneurysm development. We also detected enhanced cortical stiffness of the endothelium at these sites. Gene expression of ECs from aneurysms of the AngII ApoE-/- model when compared to sham animals mimicked expression patterns from predilection sites of healthy animals. Thus, this work highlights a striking genetic and functional regional heterogeneity in aortic ECs of healthy mice, which defines the location of aortic aneurysm formation in disease.

Pharmacogenomics of coronary artery response to intravenous gamma globulin in kawasaki disease

Kawasaki disease (KD) is a multisystem inflammatory illness of infants and young children that can result in acute vasculitis. The mechanism of coronary artery aneurysms (CAA) in KD despite intravenous gamma globulin (IVIG) treatment is not known. We performed a Whole Genome Sequencing (WGS) association analysis in a racially diverse cohort of KD patients treated with IVIG, both using AHA guidelines. We defined coronary aneurysm (CAA) (N = 234) as coronary z ≥ 2.5 and large coronary aneurysm (CAA/L) (N = 92) as z ≥ 5.0. We conducted logistic regression models to examine the association of genetic variants with CAA/L during acute KD and with persistence >6 weeks using an additive model between cases and 238 controls with no CAA. We adjusted for age, gender and three principal components of genetic ancestry. The top significant variants associated with CAA/L were in the intergenic regions (rs62154092 p < 6.32E-08 most significant). Variants in SMAT4, LOC100127, PTPRD, TCAF2 and KLRC2 were the most significant non-intergenic SNPs. Functional mapping and annotation (FUMA) analysis identified 12 genomic risk loci with eQTL or chromatin interactions mapped to 48 genes. Of these NDUFA5 has been implicated in KD CAA and MICU and ZMAT4 has potential functional implications. Genetic risk score using these 12 genomic risk loci yielded an area under the receiver operating characteristic curve (AUC) of 0.86. This pharmacogenomics study provides insights into the pathogenesis of CAA/L in IVIG-treated KD and shows that genomics can help define the cause of CAA/L to guide management and improve risk stratification of KD patients.

Pharmacogenomics of Coronary Artery Response to Intravenous Gamma Globulin in Kawasaki Disease

Background

Kawasaki disease (KD) is a multisystem inflammatory illness of infants and young children that can result in acute vasculitis. The pathological walls of afflicted coronary arteries show propensity for forming thrombosis and aneurysms. The mechanism of coronary artery aneurysms (CAA) despite intravenous gamma globulin (IVIG) treatment is not known.

Methods

We performed a Whole Genome Sequencing (WGS) association analysis in a racially diverse cohort of KD patients treated with IVIG, both using AHA guidelines. We defined coronary aneurysm (CAA) (N = 234) as coronary z>2.5 and large coronary aneurysm (CAA/L) (N = 92) as z>5.0. We conducted logistic regression models to examine the association of genetic variants with CAA/L during acute KD and with persistence >6 weeks using an additive model between cases and 238 controls with no CAA. We adjusted for age, gender and three principal components of genetic ancestry. We performed functional mapping and annotation (FUMA) analysis and further assessed the predictive risk score of genomic risk loci using the area under the receiver operating characteristic curve (AUC).

Results

The top significant variants associated with CAA/L were in the intergenic regions (rs62154092 p<6.32E-08 most significant). Variants in SMAT4, LOC100127 , PTPRD, TCAF2 and KLRC2 were the most significant non-intergenic SNPs. FUMA identified 12 genomic risk loci with eQTL or chromatin interactions mapped to 48 genes. Of these NDUFA5 has been implicated in KD CAA and MICU and ZMAT4 has potential functional implications. Genetic risk score using these 12 genomic risk loci yielded an AUC of 0.86.

Conclusions

This pharmacogenomics study provides insights into the pathogenesis of CAA/L in IVIG-treated KD patients. We have identified multiple novel SNPs associated with CAA/L and related genes with potential functional implications. The study shows that genomics can help define the cause of CAA/L to guide management and improve risk stratification of KD patients.

Role of lncRNAs in Remodeling of the Coronary Artery Plaques in Patients with Atherosclerosis

INTRODUCTION

Cardiovascular diseases (CVDs) are the leading cause of death worldwide according to World Health Organization (WHO) data. Atherosclerosis is considered as a chronic inflammatory disease that develops in response to damage to the vascular intima-media layer in most cases. In recent years, epigenetic events have emerged as important players in the development and progression of CVDs. Since noncoding RNA (ncRNAs) are important regulators in the organization of the pathophysiological processes of the cardiovascular system, they have the potential to be used as therapeutic targets, diagnostic and prognostic biomarkers. In this study long noncoding RNA (lncRNA) and mRNA gene expression were compared between coronary atherosclerotic plaques (CAP) and the internal mammary artery (IMA)  which has the same genetic makeup and is exposed to the same environmental stress conditions with CAP in the same individual.

METHODS

lncRNA and mRNA gene expressions were determined using the microarray in the samples. Microarray results were validated by RT-qPCR. Differentially expressed genes (DEGs; lncRNAs and mRNAs) were determined by GeneSpring (Ver 3.0) [p values < 0.05 and fold change (FC) > 2]. DAVID bioinformatics program was used for Gene Ontology (GO) annotation and enrichment analyses of statistically significant genes between CAP and IMA tissue.

RESULTS AND CONCLUSIONS

In our study, 345 DEGs were found to be statistically significant (p < 0.05; FC > 2) between CAP and IMA. Of these, 65 were lncRNA and 280 were mRNA. Thirty-three lncRNAs were upregulated, while 32 lncRNAs were downregulated. Some of the important mRNAs are SPP1, CYP4B1, CHRDL1, MYOC, and ALKAL2, while some of the lncRNAs are LOC105377123, LINC01857, DIO3OS, LOC101928134, and KCNA3 between CAP and IMA tissue. We also identified genes that correlated with statistically significant lncRNAs. The results of this study are expected to be an important source of data in the development of new genetically based drugs to prevent atherosclerotic plaque. In addition, the data obtained may contribute to the explanation of the epigenetic mechanisms that play a role in the pathological basis of the process that protects the IMA from atherosclerosis.

Identification of rare coding variants associated with Kawasaki disease by whole exome sequencing

Kawasaki disease (KD) is an acute pediatric vasculitis that affects genetically susceptible infants and children. To identify coding variants that influence susceptibility to KD, we conducted whole exome sequencing of 159 patients with KD and 902 controls, and performed a replication study in an independent 586 cases and 732 controls. We identified five rare coding variants in five genes (FCRLA, PTGER4, IL17F, CARD11, and SIGLEC10) associated with KD (odds ratio [OR], 1.18 to 4.41; p = 0.0027–0.031). We also performed association analysis in 26 KD patients with coronary artery aneurysms (CAAs; diameter > 5 mm) and 124 patients without CAAs (diameter < 3 mm), and identified another five rare coding variants in five genes (FGFR4, IL31RA, FNDC1, MMP8, and FOXN1), which may be associated with CAA (OR, 3.89 to 37.3; p = 0.0058–0.0261). These results provide insights into new candidate genes and genetic variants potentially involved in the development of KD and CAA.

Identification of novel locus associated with coronary artery aneurysms and validation of loci for susceptibility to Kawasaki disease

Kawasaki disease (KD) is a paediatric vasculitis associated with coronary artery aneurysms (CAA). Genetic variants influencing susceptibility to KD have been previously identified, but no risk alleles have been validated that influence CAA formation. We conducted a genome-wide association study (GWAS) for CAA in KD patients of European descent with 200 cases and 276 controls. A second GWAS for susceptibility pooled KD cases with healthy paediatric controls from vaccine trials in the UK (n = 1609). Logistic regression mixed models were used for both GWASs. The susceptibility GWAS was meta-analysed with 400 KD cases and 6101 controls from a previous European GWAS, these results were further meta-analysed with Japanese GWASs at two putative loci. The CAA GWAS identified an intergenic region of chromosome 20q13 with multiple SNVs showing genome-wide significance. The risk allele of the most associated SNV (rs6017006) was present in 13% of cases and 4% of controls; in East Asian 1000 Genomes data, the allele was absent or rare. Susceptibility GWAS with meta-analysis with previously published European data identified two previously associated loci (ITPKC and FCGR2A). Further meta-analysis with Japanese GWAS summary data from the CASP3 and FAM167A genomic regions validated these loci in Europeans showing consistent effects of the top SNVs in both populations. We identified a novel locus for CAA in KD patients of European descent. The results suggest that different genes determine susceptibility to KD and development of CAA and future work should focus on the function of the intergenic region on chromosome 20q13.

Multisystem inflammatory syndrome in children and Kawasaki disease: a critical comparison

Children and adolescents infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are predominantly asymptomatic or have mild symptoms compared with the more severe coronavirus disease 2019 (COVID-19) described in adults. However, SARS-CoV-2 is also associated with a widely reported but poorly understood paediatric systemic vasculitis. This multisystem inflammatory syndrome in children (MIS-C) has features that overlap with myocarditis, toxic-shock syndrome and Kawasaki disease. Current evidence indicates that MIS-C is the result of an exaggerated innate and adaptive immune response, characterized by a cytokine storm, and that it is triggered by prior SARS-CoV-2 exposure. Epidemiological, clinical and immunological differences classify MIS-C as being distinct from Kawasaki disease. Differences include the age range, and the geographical and ethnic distribution of patients. MIS-C is associated with prominent gastrointestinal and cardiovascular system involvement, admission to intensive care unit, neutrophilia, lymphopenia, high levels of IFNγ and low counts of naive CD4⁺ T cells, with a high proportion of activated memory T cells. Further investigation of MIS-C will continue to enhance our understanding of similar conditions associated with a cytokine storm.

Aberrant Whole Blood Gene Expression in the Lumen of Human Intracranial Aneurysms

The rupture of an intracranial aneurysm (IA) causes devastating hemorrhagic strokes. Yet, most IAs remain asymptomatic and undetected until they rupture. In the search for circulating biomarkers of unruptured IAs, we previously performed transcriptome profiling on whole blood and identified an IA-associated panel of 18 genes. In this study, we seek to determine if these genes are also differentially expressed within the IA lumen, which could provide a mechanistic link between the disease and the observed circulating gene expression patterns. To this end, we collected blood from the lumen of 37 IAs and their proximal parent vessels in 31 patients. The expression levels of 18 genes in the lumen and proximal vessel were then measured by quantitative polymerase chain reaction. This analysis revealed that the expression of 6/18 genes (CBWD6, MT2A, MZT2B, PIM3, SLC37A3, and TNFRSF4) was significantly higher in intraluminal blood, while the expression of 3/18 genes (ST6GALNAC1, TCN2, and UFSP1) was significantly lower. There was a significant, positive correlation between intraluminal and proximal expression of CXCL10, MT2A, and MZT2B, suggesting local increases of these genes is reflected in the periphery. Expression of ST6GALNAC1 and TIFAB was significantly positively correlated with IA size, while expression of CCDC85B was significantly positively correlated with IA enhancement on post-contrast MRI, a metric of IA instability and risk. In conclusion, intraluminal expression differences in half of the IA-associated genes observed in this study provide evidence for IA tissue-mediated transcriptional changes in whole blood. Additionally, some genes may be informative in assessing IA risk, as their intraluminal expression was correlated to IA size and aneurysmal wall enhancement.

TIFA and TIFAB: FHA-domain proteins involved in inflammation, hematopoiesis, and disease

Forkhead-associated (FHA) domain-containing proteins are widely expressed across eubacteria and in eukaryotes. FHA domains contain phosphopeptide recognition motifs, which operate in a variety of phosphorylation-dependent and -independent biological processes, including the DNA damage response, signal transduction, and regulation of the cell cycle. More recently, two FHA domain-containing proteins were discovered in mammalian cells as tumor necrosis factor receptor-associated factor (TRAF)-interacting proteins: TIFA and TIFAB. TIFA and TIFAB are important modifiers of the innate immune signaling through their regulation of TRAF proteins. Recent studies have also revealed distinct roles for TIFA and TIFAB in the context of immune cell function, chronic inflammation, hematopoiesis, and hematologic disorders. Collectively, these studies indicate the important role of TIFA- and TIFAB-dependent signaling in hematopoietic cells and their dysregulation in several human diseases. In this review, we summarize the molecular mechanisms and biological role of these FHA-domain homologues, placing them into the context of human disease.

A framework for understanding Kawasaki disease pathogenesis

Kawasaki disease (KD) is a common vasculitis of childhood, typically affecting children under the age of five. Despite many aspects of its presentation that bear resemblence to acute infection, no causative infectious agent has been identified despite years of intense scrutiny. Unlike most infections, however, there are significant differences in racial predilection that suggest a strong genetic influence. The inflammatory response in KD specifically targets the coronary arteries, also unusual for an infectious condition. In this review, we discuss recent hypotheses on KD pathogenesis as well as new insights into the innate immune response and mechanisms behind vascular damage. The pathogenesis is complex, however, and remains inadequately understood.

Current State of Precision Medicine in Primary Systemic Vasculitides

Precision medicine (PM) is an emerging data-driven health care approach that integrates phenotypic, genomic, epigenetic, and environmental factors unique to an individual. The goal of PM is to facilitate diagnosis, predict effective therapy, and avoid adverse reactions specific for each patient. The forefront of PM is in oncology; nonetheless, it is developing in other fields of medicine, including rheumatology. Recent studies on elucidating the genetic architecture of polygenic and monogenic rheumatological diseases have made PM possible by enabling physicians to customize medical treatment through the incorporation of clinical features and genetic data. For complex inflammatory disorders, the prevailing paradigm is that disease susceptibility is due to additive effects of common reduced-penetrance gene variants and environmental factors. Efforts have been made to calculate cumulative genetic risk score (GRS) and to relate specific susceptibility alleles for use of target therapies. The discovery of rare patients with single-gene high-penetrance mutations informed our understanding of pathways driving systemic inflammation. Here, we review the advances in practicing PM in patients with primary systemic vasculitides (PSVs). We summarize recent genetic studies and discuss current knowledge on the contribution of epigenetic factors and extracellular vesicles (EVs) in disease progression and treatment response. Implementation of PM in PSVs is a developing field that will require analysis of a large cohort of patients to validate data from genomics, transcriptomics, metabolomics, proteomics, and epigenomics studies for accurate disease profiling. This multi-omics approach to study disease pathogeneses should ultimately provide a powerful tool for stratification of patients to receive tailored optimal therapies and for monitoring their disease activity.

Identifying Patients With Kawasaki Disease Safe for Early Discharge: Development of a Risk Prediction Model at a US Children's Hospital

OBJECTIVES

To develop a model to predict risk of intravenous immunoglobulin (IVIg) nonresponse in patients with Kawasaki disease (KD) to assist in early discharge decision-making.

METHODS

Retrospective cohort study of 430 patients 0 to 18 years old discharged from a US children's hospital January 1, 2010, through July 31, 2017 with a diagnosis of KD. IVIg nonresponse was defined as at least 1 of the following: temperature ≥38.0°C between 36 hours and 7 days after initial IVIg dose, receipt of a second IVIg dose after a temperature ≥38.0°C at least 20 hours after initial IVIg dose, or readmission within 7 days with administration of a second IVIg dose. Backward stepwise logistic regression was used to select a predictive model.

RESULTS

IVIg nonresponse occurred in 19% (81 of 430) of patients. We identified a multivariate model (which included white blood cell count, hemoglobin level, platelet count, aspartate aminotransferase level, sodium level, albumin level, temperature within 6 hours of first IVIg dose, and incomplete KD) with good predictive ability (optimism-adjusted concordance index: 0.700) for IVIg nonresponse. Stratifying into 2 groups by a predictive probability cutoff of 0.10, we identified 26% of patients at low risk for IVIg nonresponse, with a sensitivity and specificity of 90% and 30%, respectively, and a negative predictive value of 93%.

CONCLUSIONS

We developed a model with good predictive value for identifying risk of IVIg nonresponse in patients with KD at a US children's hospital. Patients at lower risk may be considered for early discharge by using shared decision-making. Our model may be used to inform implementation of electronic health record tools and future risk prediction research.