Genome-wide association analysis for chronic venous disease identifies EFEMP1 and KCNH8 as susceptibility loci

Omics Data Integration Uncovers mRNA-miRNA Interaction Regions in Genes Associated with Chronic Venous Insufficiency

Background/Objectives: Chronic venous insufficiency (CVI), a chronic vascular dysfunction, is a common health problem that causes serious complications such as painful varicose veins and even skin ulcers. Identifying the underlying genetic and epigenetic factors is important for improving the quality of life of individuals with CVI. In the literature, many genes, variants, and miRNAs associated with CVI have been identified through genomic and transcriptomic studies. Despite molecular pathogenesis studies, how the genes associated with CVI are regulated by miRNAs and the effect of variants in binding regions on expression levels are still not fully understood. In this study, previously identified genes, variants, and miRNAs associated with CVI, common variants in the mRNA-miRNA binding regions, were investigated using in silico analyses. Methods: For this purpose, miRNA research tools, MBS (miRNA binding site) database, genome browsers, and the eQTL Calculator in the GTEx portal were used. Results: We identified SNVs associated with CVI that may play a direct role in the miRNA-mediated regulation of the ZNF664, COL1A2, HFE, MDN, MTHFR, SRPX, TDRD5, TSPYL4, VEGFA, and APOE genes. In addition, when the common SNVs in the mRNA binding region of 75 unique CVI related-miRNAs in five candidate genes associated with CVI were examined, seven miRNAs associated with the expression profiles of ABCA1, PIEZO1, and CASZ1 genes were identified. Conclusions: In conclusion, the relationship between genetic markers identified in the literature that play a role in the pathogenesis of the CVI and the expression profiles was evaluated for the first time in the mRNA-miRNA interaction axis.

Biallelic and monoallelic variants in EFEMP1 can cause a severe and distinct subtype of heritable connective tissue disorder

Variants in EFEMP1, encoding Fibulin-3, were previously reported as a rare cause of heritable connective tissue disorder (HCTD) with recurrent hernias and joint hypermobility. We report three new cases with biallelic or monoallelic EFEMP1 variants and severe hernia phenotypes. Two male siblings of 10 and 13 years old presented with marfanoid habitus, recurrent inguinal and umbilical hernias, generalized joint hypermobility, and scoliosis. Parents and halfsiblings reported joint hypermobility and umbilical hernias. The eldest boy died at age 16 from incarcerated gastrointestinal herniation complicated by gastric and bowel necrosis with perforation. Autopsy revealed widespread intestinal diverticula. Immunohistochemistry of skin and fascia tissue did not reveal any abnormalities, including normal staining of elastic fibers. Both siblings harbored compound heterozygous likely pathogenic EFEMP1 variants (c.1320 + 2T > A, p.? and c.698G > A, p.Gly233Asp). An unrelated 58-year-old male had marfanoid features, high myopia, recurrent diaphragmatic and inguinal hernias, and chronic gastrointestinal dilatation with severe malabsorption. Both his dizygotic twin-brother and mother had recurrent hernias and high myopia. This man died at 59 years of age, and autopsy showed extensive diaphragmatic herniation, bowel diverticula, and pulmonary emphysema. A heterozygous EFEMP1 splice-variant (c.81 + 1G > A, p.?) was identified, causing exon skipping leading to a start-loss. Targeted genome reanalysis nor RNA-sequencing revealed a second variant at the other allele. The reported individuals expand the clinical and pathological phenotypes of EFEMP1-related disease, a distinct entity within the spectrum of HCTD. The severe and recurrent hernias, gastrointestinal dilatation, and diverticulosis result in an increased risk for life-threatening complications, demanding early recognition and close monitoring.

Glycocalyx disruption, endothelial dysfunction and vascular remodeling as underlying mechanisms and treatment targets of chronic venous disease

The glycocalyx is an essential structural and functional component of endothelial cells. Extensive hemodynamic changes cause endothelial glycocalyx disruption and vascular dysfunction, leading to multiple arterial and venous disorders. Chronic venous disease (CVD) is a common disorder of the lower extremities with major health and socio-economic implications, but complex pathophysiology. Genetic aberrations accentuated by environmental factors, behavioral tendencies, and hormonal disturbances promote venous reflux, valve incompetence, and venous blood stasis. Increased venous hydrostatic pressure and changes in shear-stress cause glycocalyx injury, endothelial dysfunction, secretion of adhesion molecules, leukocyte recruitment/activation, and release of cytokines, chemokines, and hypoxia-inducible factor, causing smooth muscle cell switch from contractile to synthetic proliferative phenotype, imbalance in matrix metalloproteinases (MMPs), degradation of collagen and elastin, and venous tissue remodeling, leading to venous dilation and varicose veins. In the advanced stages of CVD, leukocyte infiltration of the vein wall causes progressive inflammation, fibrosis, disruption of junctional proteins, accumulation of tissue metabolites and reactive oxygen and nitrogen species, and iron deposition, leading to skin changes and venous leg ulcer (VLU). CVD management includes compression stockings, venotonics, and surgical intervention. In addition to its antithrombotic and fibrinolytic properties, literature suggests sulodexide benefits in reducing inflammation, promoting VLU healing, improving endothelial function, exhibiting venotonic properties, and inhibiting MMP-9. Understanding the role of glycocalyx, endothelial dysfunction, and vascular remodeling should help delineate the underlying mechanisms and develop improved biomarkers and targeted therapy for CVD and VLU.

Retracted article: The role of second generation sequencing technology and nanomedicine in the monitoring and treatment of lower extremity deep vein thrombosis susceptibility genes

Rong Chang, Chunsheng Wang, Xiangqi Kong, Wenhui Li and Jinchun Wu. The role of second generation sequencing technology and nanomedicine in the monitoring and treatment of lower extremity deep vein thrombosis susceptibility genes. Bioengineered. 2021 Nov. doi: 10.1080/21655979.2021.2003926.Since publication, significant concerns have been raised about the compliance with ethical policies for human research and the integrity of the data reported in the article.When approached for an explanation, the authors provided some original data but were not able to provide all the necessary supporting information. As verifying the validity of published work is core to the scholarly record's integrity, we are retracting the article. All authors listed in this publication have been informed.We have been informed in our decision-making by our editorial policies and the COPE guidelines.The retracted article will remain online to maintain the scholarly record, but it will be digitally watermarked on each page as 'Retracted.'

Exploring the causal relationship between plasma lipids and varicose veins of lower extremity: A comprehensive two-sample Mendelian randomization study

Varicose veins of the lower extremities (VVs) is a common chronic vascular disease, with high prevalence rates in some countries; however, their pathogenesis remains unclear. Some studies have identified associations between changes in specific plasma lipid molecules, such as phosphatidylethanolamine (PE), phosphatidylcholine (PC), and sphingomyelin (SM), and the onset of VVs, but due to confounders and reverse causality, the causal relationship remains unclear. Meanwhile, studies on the potential link between other plasma lipids beyond PE, PC, and SM and the risk of VVs in the lower extremities are lacking. This study aimed to explore the potential causal relationship between VVs and plasma lipid levels to provide theoretical insights into the interrelation of plasma lipids and VVs in their occurrence and progression. We conducted a two-sample Mendelian randomization (MR) analysis to assess the potential connection between genetically predicted levels of individual plasma lipids and the risk of developing VVs. We utilized data from a large-scale genome-wide association study involving 7174 Finnish individuals for 179 plasma lipidomes along with VVs genome-wide association study data from 408,455 UK individuals. MR analysis employed methods, such as inverse-variance weighting, weighted median, Bayesian Weighted Mendelian Randomization, and MR-Egger regression. The inverse-variance weighting method was primarily used to assess causality. The validity of the results was demonstrated through sensitivity analysis. In total, 12 lipids were found to have their plasma levels associated with an increased risk of VVs. This includes 3 types of PE, 7 types of PC, and 2 types of phosphatidylinositol. However, no significant causal relationship was found between the plasma levels of 11 types of SM and VVs. These results support the existence of a potential causal relationship between specific types of lipid levels and the risk of VVs, which can provide clues for further studies on biological mechanisms and the exploration of potential therapeutic targets.

Exome sequencing identifies novel genetic variants associated with varicose veins

BACKGROUND

Varicose veins (VV) are one of the common human diseases, but the role of genetics in its development is not fully understood.

METHODS

We conducted an exome-wide association study of VV using whole-exome sequencing data from the UK Biobank, and focused on common and rare variants using single-variant association analysis and gene-level collapsing analysis.

FINDINGS

A total of 13,823,269 autosomal genetic variants were obtained after quality control. We identified 36 VV-related independent common variants mapping to 34 genes by single-variant analysis and three rare variant genes (PIEZO1, ECE1, FBLN7) by collapsing analysis, and most associations between genes and VV were replicated in FinnGen. PIEZO1 was the closest gene associated with VV (P = 5.05 × 10-31), and it was found to reach exome-wide significance in both single-variant and collapsing analyses. Two novel rare variant genes (ECE1 and METTL21A) associated with VV were identified, of which METTL21A was associated only with females. The pleiotropic effects of VV-related genes suggested that body size, inflammation, and pulmonary function are strongly associated with the development of VV.

CONCLUSIONS

Our findings highlight the importance of causal genes for VV and provide new directions for treatment.

Key Regulators of Angiogenesis and Inflammation Are Dysregulated in Patients with Varicose Veins

Varicose veins (VVs) are the most common manifestation of chronic venous disease (CVD) and appear as abnormally enlarged and tortuous superficial veins. VVs result from functional abnormalities in the venous circulation of the lower extremities, such as venous hypertension, venous valve incompetence, and venous reflux. Previous studies indicate that enhanced angiogenesis and inflammation contribute to the progression and onset of VVs; however, dysregulations in signaling pathways associated with these processes in VVs patients are poorly understood. Therefore, in our study, we aimed to identify key regulators of angiogenesis and inflammation that are dysregulated in patients with VVs. Expression levels of 18 genes were analyzed in peripheral blood mononuclear cells (PBMC) using real-time PCR, as well as plasma levels of 6 proteins were investigated using ELISA. Higher levels of CCL5, PDGFA, VEGFC, TGF-alpha, TGF-beta 1, and VEGF-A, as well as lower levels of VEGFB and VEGF-C, were found to be statistically significant in the VV group compared to the control subjects without VVs. None of the analyzed factors was associated with the venous localization of the varicosities. The presented study identified dysregulations in key angiogenesis- and inflammation-related factors in PBMC and plasma from VVs patients, providing new insight into molecular mechanisms that could contribute to the development of VVs and point out promising candidates for circulatory biomarkers of this disease.

Vicious Circle With Venous Hypertension, Irregular Flow, Pathological Venous Wall Remodeling, and Valve Destruction in Chronic Venous Disease: A Review

Substantial advances occurred in phlebological practice in the last two decades. With the use of modern diagnostic equipment, the patients' venous hemodynamics can be examined in detail in everyday practice. Application of venous segments for arterial bypasses motivated studies on the effect of hemodynamic load on the venous wall. New animal models have been developed to study hemodynamic effects on the venous system. In vivo and in vitro studies revealed cellular phase transitions of venous endothelial, smooth muscle, and fibroblastic cells and changes in connective tissue composition, under hemodynamic load and at different locations of the chronically diseased venous system. This review is an attempt to integrate our knowledge from epidemiology, paleoanthropology and anthropology, clinical and experimental hemodynamic studies, histology, cell physiology, cell pathology, and molecular biology on the complex pathomechanism of this frequent disease. Our conclusion is that the disease is initiated by limited genetic adaptation of mankind not to bipedalism but to bipedalism in the unmoving standing or sitting position. In the course of the disease several pathologic vicious circles emerge, sustained venous hypertension inducing cellular phase transitions, chronic wall inflammation, apoptosis of cells, pathologic dilation, and valvular damage which, in turn, further aggravate the venous hypertension.

Replication study identified EFEMP1 association with varicose vein predisposition among Indians

BACKGROUND

Varicose vein is a chronic condition that affects the lower extremities of the human body. Several factors have been implicated in the development of this disease, viz age, gender, weight, height and prolonged standing. Recently, genome-wide studies have identified genetic biomarkers that are associated with varicose veins in different ethnic groups. Such genetic studies are lacking in South Asians specifically in Indians where the prevalence of varicose veins is high, and it is important to replicate these variants in the stated population. The study aimed to replicate the association of genetic variants associated with varicose veins in this target population, which were found to be associated with the other ethnic groups.

METHODOLOGY

The studied cohort is of the Indian population comprising unrelated 104 varicose veins cases and 448 non-varicose vein controls. The samples were genotyped using the Illumina Global Screening Array. Using the genomic data from UK BioBank and 23andMe studied cohorts; eight genetic variants were selected to replicate in our dataset. The allelic association was performed to identify the effective allele and risk was estimated using odds ratio and p-value as level of significance. Multifactor Dimensionality Reduction was used to estimate the cumulative effect of variants in Indians.

RESULT

Variant rs3791679 of EFEMP1 was found to be associated with varicose veins in Indians. After observing the association of the EFEMP1 with varicose veins, we further ensued to identify all genetic variants within EFEMP1 to uncover the additional variants associated with this trait. Interestingly, we identified six new variants of EFEMP1 gene that have shown association. Moreover, the cumulative effect of all associated variations was estimated and the risk was 2.7 times higher in cases than controls whereas independently their effect ranges from 0.37-1.58.

CONCLUSION

This study identifies EFEMP1 as a potential gene related to the risk of varicose veins in Indians. It also highlights that evaluating the maximum number of variants of a gene rather than focusing solely on replicating single variations offers a more comprehensive and nuanced understanding of the genetic factors contributing to a complex trait like varicose veins.

SKAP2-A Molecule at the Crossroads for Integrin Signalling and Immune Cell Migration and Function

Src-kinase associated protein 2 (SKAP2) is an intracellular scaffolding protein that is broadly expressed in immune cells and is involved in various downstream signalling pathways, including, but not limited to, integrin signalling. SKAP2 has a wide range of binding partners and fine-tunes the rearrangement of the cytoskeleton, thereby regulating cell migration and immune cell function. Mutations in SKAP2 have been associated with several inflammatory disorders such as Type 1 Diabetes and Crohn's disease. Rodent studies showed that SKAP2 deficient immune cells have diminished pathogen clearance due to impaired ROS production and/or phagocytosis. However, there is currently no in-depth understanding of the functioning of SKAP2. Nevertheless, this review summarises the existing knowledge with a focus of its role in signalling cascades involved in cell migration, tissue infiltration and immune cell function.

Genome-wide association study of varicose veins identifies a protective missense variant in GJD3 enriched in the Finnish population

Varicose veins is the most common manifestation of chronic venous disease that displays female-biased incidence. To identify protein-inactivating variants that could guide identification of drug target genes for varicose veins and genetic evidence for the disease prevalence difference between the sexes, we conducted a genome-wide association study of varicose veins in Finns using the FinnGen dataset with 17,027 cases and 190,028 controls. We identified 50 associated genetic loci (P < 5.0 × 10-8) of which 29 were novel including one near ERG with female-specificity (rs2836405-G, OR[95% CI] = 1.09[1.05-1.13], P = 3.1 × 10-8). These also include two X-chromosomal (ARHGAP6 and SRPX) and two autosomal novel loci (TGFB2 and GJD3) with protein-coding lead variants enriched above 56-fold in Finns over non-Finnish non-Estonian Europeans. A low-frequency missense variant in GJD3 (p.Pro59Thr) is exclusively associated with a lower risk for varicose veins (OR = 0.62 [0.55-0.70], P = 1.0 × 10-14) in a phenome-wide scan of the FinnGen data. The absence of observed pleiotropy and its membership of the connexin gene family underlines GJD3 as a potential connexin-modulating therapeutic strategy for varicose veins. Our results provide insights into varicose veins etiopathology and highlight the power of isolated populations, including Finns, to discover genetic variants that inform therapeutic development.

Molecular Determinants of Chronic Venous Disease: A Comprehensive Review

Chronic Venous Disease (CVD) refers to several pathological and hemodynamic alterations of the veins of lower limbs causing a wide range of symptoms and signs with a high prevalence in the general population and with disabling consequences in the most severe forms. The etiology and pathophysiology of CVD is complex and multifactorial, involving genetic, proteomic, and cellular mechanisms that result in changes to the venous structure and functions. Expressions of several genes associated with angiogenesis, vascular development, and the regulation of veins are responsible for the susceptibility to CVD. Current evidence shows that several extracellular matrix alterations (ECM) could be identified and in some cases pharmacologically targeted. This review shows the most up to date information on molecular determinants of CVD in order to provide a complete overview of the current knowledge on this topic. In particular, the article explores the genetic influence, the hormonal influence, ECM imbalance, and histopathology of CVD and the role of endothelial dysfunction in CVD.

Genetics of varicose veins reveals polygenic architecture and genetic overlap with arterial and venous disease

Varicose veins represent a common cause of cardiovascular morbidity, with limited available medical therapies. Although varicose veins are heritable and epidemiologic studies have identified several candidate varicose vein risk factors, the molecular and genetic basis remains uncertain. Here we analyzed the contribution of common genetic variants to varicose veins using data from the Veterans Affairs Million Veteran Program and four other large biobanks. Among 49,765 individuals with varicose veins and 1,334,301 disease-free controls, we identified 139 risk loci. We identified genetic overlap between varicose veins, other vascular diseases and dozens of anthropometric factors. Using Mendelian randomization, we prioritized therapeutic targets via integration of proteomic and transcriptomic data. Finally, topological enrichment analyses confirmed the biologic roles of endothelial shear flow disruption, inflammation, vascular remodeling and angiogenesis. These findings may facilitate future efforts to develop nonsurgical therapies for varicose veins.

Fibulin-3 Deficiency Protects Against Myocardial Injury Following Ischaemia/ Reperfusion in in vitro Cardiac Spheroids

Myocardial infarction (MI, or heart attack) is a leading cause of death worldwide. Myocardial ischaemia reperfusion (I/R) injury typical of MI events is also associated with the development of cardiac fibrosis and heart failure in patients. Fibulin-3 is an extracellular matrix component that plays a role in regulating MI response in the heart. In this study, we generated and compared in vitro cardiac spheroids (CSs) from wild type (WT) and fibulin-3 knockout (Fib-3 KO) mice. These were then exposed to pathophysiological changes in oxygen (O₂) concentrations to mimic an MI event. We finally measured changes in contractile function, cell death, and mRNA expression levels of cardiovascular disease genes between WT and Fib-3 KO CSs. Our results demonstrated that there are significant differences in growth kinetics and endothelial network formation between WT and Fib-3 KO CSs, however, they respond similarly to changes in O₂ concentrations. Fib-3 deficiency resulted in an increase in viability of cells and improvement in contraction frequency and fractional shortening compared to WT I/R CSs. Gene expression analyses demonstrated that Fib-3 deficiency inhibits I/R injury and cardiac fibrosis and promotes angiogenesis in CSs. Altogether, our findings suggest that Fib-3 deficiency makes CSs resistant to I/R injury and associated cardiac fibrosis and helps to improve the vascular network in CSs.

Genome-wide association analysis and replication in 810,625 individuals with varicose veins

Varicose veins affect one-third of Western society, with a significant subset of patients developing venous ulceration, costing $14.9 billion annually in the USA. Current management consists of either compression stockings, or surgical ablation for more advanced disease. Most varicose veins patients report a positive family history, and heritability is ~17%. We describe the largest two-stage genome-wide association study of varicose veins in 401,656 individuals from UK Biobank, and replication in 408,969 individuals from 23andMe (total 135,514 cases and 675,111 controls). Forty-nine signals at 46 susceptibility loci were discovered. We map 237 genes to these loci, several of which are biologically plausible and tractable to therapeutic targeting. Pathway analysis identified enrichment in extracellular matrix biology, inflammation, (lymph)angiogenesis, vascular smooth muscle cell migration, and apoptosis. Using a polygenic risk score (PRS) derived in an independent cohort, we demonstrate its predictive utility and correlation with varicose veins surgery.

A genome-wide association study for varicose veins

Background

The aim was to compare the genetic information of varicose vein patients with that of a healthy population attempting to identify certain significant genetic associations.

Method

Patients’ clinical characteristics and demographics were collected, and their genetic samples were examined. The results were compared to the genetic information of one thousand sex-matched healthy controls from Taiwan Biobank database. The Clinical-Etiology-Anatomy-Pathophysiology classification was applied for further subgroup analysis.

Results

After comparison of genetic information of ninety-six patients to that of healthy controls, two significant single nucleotide polymorphisms (SNPs) were identified. One was in DPYSL2 gene, and the other was in VSTM2L gene. A further comparison between C2-3 patient subgroup and C4-6 subgroup identified another four significant SNPs, which were located in ZNF664-FAM101A, PHF2, ACOT11, and TOM1L1 genes.

Conclusion

Our preliminary result identified six significant SNPs located in six different genes. All of them and their genetic products may warrant further investigations.

Single-Nucleus Transcriptional Profiling of Chronic Kidney Disease after Cisplatin Nephrotoxicity

Cisplatin induces both acute and chronic nephrotoxicity during chemotherapy in patients with cancer. Presented here is the first study of single-nucleus RNA sequencing (snRNA-seq) of cisplatin-induced nephrotoxicity. Repeated low-dose cisplatin treatment (RLDC) led to decreases in renal function and kidney weight in mice at 9 weeks. The kidneys of these mice showed tubular degeneration and dilation. snRNA-seq identified 16 cell types and 17 cell clusters in these kidneys. Cluster-by-cluster comparison demonstrated cell type-specific changes in gene expression and identified a unique proximal tubule (PT) injury/repair cluster that co-expressed the injury marker kidney injury molecule-1 (Kim1) and the proliferation marker Ki-67. Compared with control, post-RLDC kidneys had 424 differentially expressed genes in PT cells, including tubular transporters and cytochrome P450 enzymes involved in lipid metabolism. snRNA-seq also revealed transcriptional changes in potential PT injury markers (Krt222, Eda2r, Ltbp2, and Masp1) and repair marker (Bex4). RLDC induced inflammation and proinflammatory cytokines (RelB, TNF-α, Il7, Ccl2, and Cxcl2) and the expression of fibrosis markers (fibronectin, collagen I, connective tissue growth factor, vimentin, and α-smooth muscle actin). Together, these results provide new insights into RLDC-induced transcriptional changes at the single-cell level that may contribute to the development of chronic kidney problems in patients with cancer after cisplatin chemotherapy.

Integrative analysis prioritizes the relevant genes and risk factors for chronic venous disease

OBJECTIVES

Chronic venous disease (CVD) refers to a range of symptoms resulting from long-term morphological and functional abnormalities of the venous system. However, the mechanism of CVD development remains largely unknown. Here we aim to provide more information on CVD pathogenesis, prevention strategies and therapy development through the integrative analysis of large-scale genetic data.

METHODS

Genetic data were obtained from publicly accessible databases. We utilized different approaches, including FUMA, DEPICT, Sherlock, SMR/HEIDIS, DEPICT and NetWAS to identify possible causal genes for CVD. Candidate genes were prioritized to further literature-based review. The differential expression of prioritized genes was validated by microarray from the Gene Expression Omnibus (GEO), a public genomics data repository" and Real-time quantitative PCR (qPCR) of varicose veins (VVs) specimens. The causal relationships between risk factors and CVD were assessed using the Two-sample Mendelian randomization (MR) approach.

RESULTS

We identified 46 lead single-nucleotide polymorphisms (SNPs) and 26 plausible causal genes for CVD. Microarray data indicated differential expression of possible causal genes in CVD when compared to controls. The expression levels of WDR92, RSPO3, LIMA, ABCB10, DNAJC7, C1S, CXCL1 were significantly down-regulated (P<0.05). PHLDA1 and SERPINE1 were significantly upregulated (P<0.05). Dysregulated expression of WDR92, RSPO3 and CASZ1 was also found in varicose vein specimens by qPCR. Two-sample MR suggested causative effects of BMI (OR, 1.008, 95%CI, 1.005-1.010), standing height (OR, 1.009, 95%CI, 1.007-1.011), college degree (OR, 0.983, 95%CI, 0.991-0.976), insulin (OR, 0.858, 95%CI, 0.794-0.928) and metformin (OR, 0.944, 95%CI, 0.904-0.985) on CVD.

CONCLUSIONS

Our study integrates genetic and gene expression data to make an effective risk gene prediction and etiological inferences for CVD. Prioritized candidate genes provide more insights into CVD pathogenesis, and the causative effects of risk factors on CVD that deserve further investigation.

Varicose Veins and Chronic Venous Disease

Chronic venous disease is a worldwide problem associated with significant morbidity and is expected to increase in prevalence as the current population ages. This is a comprehensive review of the anatomy, pathophysiology, genomics, clinical classification, and treatment modalities of chronic venous disease.

Understanding Chronic Venous Disease: A Critical Overview of Its Pathophysiology and Medical Management

Chronic venous disease (CVD) is a multifactorial condition affecting an important percentage of the global population. It ranges from mild clinical signs, such as telangiectasias or reticular veins, to severe manifestations, such as venous ulcerations. However, varicose veins (VVs) are the most common manifestation of CVD. The explicit mechanisms of the disease are not well-understood. It seems that genetics and a plethora of environmental agents play an important role in the development and progression of CVD. The exposure to these factors leads to altered hemodynamics of the venous system, described as ambulatory venous hypertension, therefore promoting microcirculatory changes, inflammatory responses, hypoxia, venous wall remodeling, and epigenetic variations, even with important systemic implications. Thus, a proper clinical management of patients with CVD is essential to prevent potential harms of the disease, which also entails a significant loss of the quality of life in these individuals. Hence, the aim of the present review is to collect the current knowledge of CVD, including its epidemiology, etiology, and risk factors, but emphasizing the pathophysiology and medical care of these patients, including clinical manifestations, diagnosis, and treatments. Furthermore, future directions will also be covered in this work in order to provide potential fields to explore in the context of CVD.

Epidemiology and Genetics of Venous Thromboembolism and Chronic Venous Disease

Venous disease is a term that broadly covers both venous thromboembolic disease and chronic venous disease. The basic pathophysiology of venous thromboembolism and chronic venous disease differ as venous thromboembolism results from an imbalance of hemostasis and thrombosis while chronic venous disease occurs in the setting of tissue damage because of prolonged venous hypertension. Both diseases are common and account for significant mortality and morbidity, respectively, and collectively make up a large health care burden. Despite both diseases having well-characterized environmental components, it has been known for decades that family history is an important risk factor, implicating a genetic element to a patient's risk. Our understanding of the pathogenesis of these diseases has greatly benefited from an expansion of population genetic studies from pioneering familial studies to large genome-wide association studies; we now have multiple risk loci for each venous disease. In this review, we will highlight the current state of knowledge on the epidemiology and genetics of venous thromboembolism and chronic venous disease and directions for future research.

Genetic Polymorphisms Related to VO2max Adaptation Are Associated With Elite Rugby Union Status and Competitive Marathon Performance

PURPOSE

Genetic polymorphisms have been associated with the adaptation to training in maximal oxygen uptake (V˙O2max). However, the genotype distribution of selected polymorphisms in athletic cohorts is unknown, with their influence on performance characteristics also undetermined. This study investigated whether the genotype distributions of 3 polymorphisms previously associated with V˙O2max training adaptation are associated with elite athlete status and performance characteristics in runners and rugby athletes, competitors for whom aerobic metabolism is important.

METHODS

Genomic DNA was collected from 732 men including 165 long-distance runners, 212 elite rugby union athletes, and 355 nonathletes. Genotype and allele frequencies of PRDM1 rs10499043 C/T, GRIN3A rs1535628 G/A, and KCNH8 rs4973706 T/C were compared between athletes and nonathletes. Personal-best marathon times in runners, as well as in-game performance variables and playing position, of rugby athletes were analyzed according to genotype.

RESULTS

Runners with PRDM1 T alleles recorded marathon times ∼3 minutes faster than CC homozygotes (02:27:55 [00:07:32] h vs 02:31:03 [00:08:24] h, P = .023). Rugby athletes had 1.57 times greater odds of possessing the KCNH8 TT genotype than nonathletes (65.5% vs 54.7%, χ2 = 6.494, P = .013). No other associations were identified.

CONCLUSIONS

This study is the first to demonstrate that polymorphisms previously associated with V˙O2max training adaptations in nonathletes are also associated with marathon performance (PRDM1) and elite rugby union status (KCNH8). The genotypes and alleles previously associated with superior endurance-training adaptation appear to be advantageous in long-distance running and achieving elite status in rugby union.

Identification of Key Pathways and Genes in Obesity Using Bioinformatics Analysis and Molecular Docking Studies

Obesity is an excess accumulation of body fat. Its progression rate has remained high in recent years. Therefore, the aim of this study was to diagnose important differentially expressed genes (DEGs) associated in its development, which may be used as novel biomarkers or potential therapeutic targets for obesity. The gene expression profile of E-MTAB-6728 was downloaded from the database. After screening DEGs in each ArrayExpress dataset, we further used the robust rank aggregation method to diagnose 876 significant DEGs including 438 up regulated and 438 down regulated genes. Functional enrichment analysis was performed. These DEGs were shown to be significantly enriched in different obesity related pathways and GO functions. Then protein-protein interaction network, target genes - miRNA regulatory network and target genes - TF regulatory network were constructed and analyzed. The module analysis was performed based on the whole PPI network. We finally filtered out STAT3, CORO1C, SERPINH1, MVP, ITGB5, PCM1, SIRT1, EEF1G, PTEN and RPS2 hub genes. Hub genes were validated by ICH analysis, receiver operating curve (ROC) analysis and RT-PCR. Finally a molecular docking study was performed to find small drug molecules. The robust DEGs linked with the development of obesity were screened through the expression profile, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for obesity.

Investigation of association between genetic polymorphisms of MMP2, MMP8, MMP9 and TIMP2 and development of varicose veins in the Slovak Population - pilot study

Matrix metalloproteinases (MMPs) are a family of zinc-dependent metalloendopeptidases that degrades extracellular matrix (ECM) components. MMPs are associated with venous wall remodelling, proliferation, migration, phenotypic and functional transformation of vascular smooth muscle cells and ECM organization under the physiological and pathophysiological conditions. We investigated possible association of genetic promoter polymorphisms of MMP2 (rs243866), MMP8 (rs11225395), MMP9 (rs3918242) and TIMP2 (rs8179090) to varicose veins development in the Slovak population. Genomic DNA from 276 Slovak individuals (138 cases, 138 controls) was genotyped for selected SNPs (rs243866, rs11225395, rs3918242 and rs8179090) using the PCR-RFLP analysis. The data were analysed by chi-squared (chi2) test, logistic regression, and Mann-Whitney test. The risk of varicose veins development was evaluated in dominant, codominant and recessive genetic models. The statistical evaluation of selected polymorphisms in patients in all three genetic models has not shown a significant risk of varicose veins development. Our study has not shown the association between selected polymorphisms and increased risk of varicose veins development in Slovak population. More evidence with broaden sample size is needed.

Fibulin-3 affects vascular endothelial function and is regulated by angiotensin II

OBJECTIVE

The aim of the present study was to investigate the effect of fibulin-3 on vascular endothelial function, and to explore the relevant underlying mechanism with regard to the involvement of angiotensin II (AngII).

METHODS

One hundred and eight patients with essential hypertension (EH) and 31 controls were included to measure the flow-mediated dilatation (FMD). Serum fibulin-3 and AngII were examined using enzyme-linked immunosorbent assay (ELISA) and radioimmunoassay. Stable transfection of fibulin-3 was conducted on human umbilical vein endothelial cells (HUVECs) and SV40T-transformed HUVECs (PUMC-HUVEC-T1 cells). Cell counting kit-8 assay, cell cycle assay, wound healing assay, Transwell assay, apoptosis assay, and tube formation assay were subsequently performed. The expression of angiogenesis-associated genes [endothelial nitric oxide synthase (eNOS) and vascular endothelial growth factor A (VEGFA)] were measured by western blot analysis. HUVECs and PUMC-HUVEC-T1 cells were treated with AngII, and with or without an inhibitor of nuclear factor κB (NF-κB), BAY 11-7082. Pro-inflammatory cytokines [interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α)] were detected by ELISA. The expression levels of fibulin-3 and p65 were then measured by western blotting.

RESULTS

Lower levels of serum fibulin-3 were accompanied by poorer FMD and higher levels of serum AngII in patients with EH. Fibulin-3 overexpression promoted cell proliferation, migration, and angiogenesis, but led to an inhibition of apoptosis. By contrast, fibulin-3 downregulation inhibited cell proliferation, migration and angiogenesis, but promoted apoptosis. AngII induced inflammation and inhibited the expression of fibulin-3. BAY 11-7082 eliminated the inhibitory effect of AngII on fibulin-3.

CONCLUSIONS

Taken together, the results of the present study have shown that serum fibulin-3 may be a predictor of vascular endothelial function in patients with EH. Fibulin-3 gene may also have a beneficial role in repairing the vascular endothelium. Furthermore, the results also suggested that fibulin-3 gene was suppressed by AngII via the NF-κB signaling pathway.

The Pathophysiological Significance of Fibulin-3

Fibulin-3 (also known as EGF-containing fibulin extracellular matrix protein 1 (EFEMP1)) is a secreted extracellular matrix glycoprotein, encoded by the EFEMP1 gene that belongs to the eight-membered fibulin protein family. It has emerged as a functionally unique member of this family, with a diverse array of pathophysiological associations predominantly centered on its role as a modulator of extracellular matrix (ECM) biology. Fibulin-3 is widely expressed in the human body, especially in elastic-fibre-rich tissues and ocular structures, and interacts with enzymatic ECM regulators, including tissue inhibitor of metalloproteinase-3 (TIMP-3). A point mutation in EFEMP1 causes an inherited early-onset form of macular degeneration called Malattia Leventinese/Doyne honeycomb retinal dystrophy (ML/DHRD). EFEMP1 genetic variants have also been associated in genome-wide association studies with numerous complex inherited phenotypes, both physiological (namely, developmental anthropometric traits) and pathological (many of which involve abnormalities of connective tissue function). Furthermore, EFEMP1 expression changes are implicated in the progression of numerous types of cancer, an area in which fibulin-3 has putative significance as a therapeutic target. Here we discuss the potential mechanistic roles of fibulin-3 in these pathologies and highlight how it may contribute to the development, structural integrity, and emergent functionality of the ECM and connective tissues across a range of anatomical locations. Its myriad of aetiological roles positions fibulin-3 as a molecule of interest across numerous research fields and may inform our future understanding and therapeutic approach to many human diseases in clinical settings.

CEAP clinical classes C0S-C4: differences, similarities and role of Ruscus + HMC + vitamin C in patients with chronic venous disease

Since the publication of the CEAP classification, new research has enriched our knowledge; notably on the heritability of CVD and the genetic and environmental factors involved in this condition, as well as the symptoms apparent within the spectrum of the CEAP clinical classes and the benefits of medical treatment. Using the CEAP classification as a special theme, a symposium with the same title as the present paper was held at the annual meeting of the 2019 European Venous Forum. The lectures presented much valuable information, from which some key points can be extracted. The influence of environmental factors was demonstrated, and the fact that a large amount of information can be obtained from comprehensive history taking. There is robust evidence for heritability. Many candidate genes/loci have been identified, potentially offering new targets for treatment. More research is needed, notably using genome-wide association studies and also on microbiota, which may play a role in CVD through the inflammation pathway. Ruscus + HMC + vitamin C acts by increasing venous and lymphatic tone, protecting microcirculation, and reducing inflammation. It improves quality of life in C0S to C3 CVD patients, while a review of clinical studies and a meta-analysis have confirmed its clinical efficacy across a wide spectrum of CVD clinical classes: C0S, C1S, C2, C3 and C4. It has been awarded a Grade 1A recommendation by the international guidelines.

A variant of the castor zinc finger 1 (CASZ1) gene is differentially associated with the clinical classification of chronic venous disease

Recent reports have suggested a reproducible association between the rs11121615 SNP, located within an intron of the castor zinc finger 1 (CASZ1) gene, and varicose veins. This study aimed to determine if this variant is also differentially associated with the various clinical classifications of chronic venous disease (CVD). The rs11121615 SNP was genotyped in two independent cohorts from New Zealand (n = 1876 controls /1606 CVD cases) and the Netherlands (n = 1626/2966). Participants were clinically assessed using well-established CVD criteria. The association between the rs11121615 C-allele and varicose veins was validated in both cohorts. This was strongest in those with higher clinical severity classes and was not significant in those with non-varicose vein CVD. Functional analysis of the rs11121615 variant demonstrated that the risk allele was associated with increased enhancer activity. This study demonstrates that the CASZ1 gene associated C-allele of rs11121615 has a significant, reproducible, association with CVD (CEAP C ≥ 2 meta-odds ratio 1.31, 95% CI 1.27–1.34, P = 1 × 10⁻⁹⁸, PHet = 0.25), but not with non-varicose vein (CEAP C1, telangiectasia or reticular veins) forms of venous disease. The effect size of this association therefore appears to be susceptible to influence by phenotypic heterogeneity, particularly if a cohort includes a large number of cases with lower severity CVD.

Next-Generation Sequencing Profiles of the Methylome and Transcriptome in Peripheral Blood Mononuclear Cells of Rheumatoid Arthritis

Using next-generation sequencing to decipher methylome and transcriptome and underlying molecular mechanisms contributing to rheumatoid arthritis (RA) for improving future therapies, we performed methyl-seq and RNA-seq on peripheral blood mononuclear cells (PBMCs) from RA subjects and normal donors. Principal component analysis and hierarchical clustering revealed distinct methylation signatures in RA with methylation aberrations noted across chromosomes. Methylation alterations varied with CpG features and genic characteristics. Typically, CpG islands and CpG shores were hypermethylated and displayed the greatest methylation variance. Promoters were hypermethylated and enhancers/gene bodies were hypomethylated, with methylation variance associated with expression variance. RA genetically associated genes preferentially displayed differential methylation and differential expression or interacted with differentially methylated and differentially expressed genes. These differentially methylated and differentially expressed genes were enriched with several signaling pathways and disease categories. 10 genes (CD86, RAB20, XAF1, FOLR3, LTBR, KCNH8, DOK7, PDGFA, PITPNM2, CELSR1) with concomitantly differential methylation in enhancers/promoters/gene bodies and differential expression in B cells were validated. This integrated analysis of methylome and transcriptome identified novel epigenetic signatures associated with RA and highlighted the interaction between genetics and epigenetics in RA. These findings help our understanding of the pathogenesis of RA and advance epigenetic studies in regards to the disease.

Genome-wide association analysis of diverticular disease points towards neuromuscular, connective tissue and epithelial pathomechanisms

OBJECTIVE

Diverticular disease is a common complex disorder characterised by mucosal outpouchings of the colonic wall that manifests through complications such as diverticulitis, perforation and bleeding. We report the to date largest genome-wide association study (GWAS) to identify genetic risk factors for diverticular disease.

DESIGN

Discovery GWAS analysis was performed on UK Biobank imputed genotypes using 31 964 cases and 419 135 controls of European descent. Associations were replicated in a European sample of 3893 cases and 2829 diverticula-free controls and evaluated for risk contribution to diverticulitis and uncomplicated diverticulosis. Transcripts at top 20 replicating loci were analysed by real-time quatitative PCR in preparations of the mucosal, submucosal and muscular layer of colon. The localisation of expressed protein at selected loci was investigated by immunohistochemistry.

RESULTS

We discovered 48 risk loci, of which 12 are novel, with genome-wide significance and consistent OR in the replication sample. Nominal replication (p<0.05) was observed for 27 loci, and additional 8 in meta-analysis with a population-based cohort. The most significant novel risk variant rs9960286 is located near CTAGE1 with a p value of 2.3×10-10 and 0.002 (ORallelic=1.14 (95% CI 1.05 to 1.24)) in the replication analysis. Four loci showed stronger effects for diverticulitis, PHGR1 (OR 1.32, 95% CI 1.12 to 1.56), FAM155A-2 (OR 1.21, 95% CI 1.04 to 1.42), CALCB (OR 1.17, 95% CI 1.03 to 1.33) and S100A10 (OR 1.17, 95% CI 1.03 to 1.33).

CONCLUSION

In silico analyses point to diverticulosis primarily as a disorder of intestinal neuromuscular function and of impaired connective fibre support, while an additional diverticulitis risk might be conferred by epithelial dysfunction.

Varicose veins of lower extremities: Insights from the first large-scale genetic study

Varicose veins of lower extremities (VVs) are a common multifactorial vascular disease. Genetic factors underlying VVs development remain largely unknown. Here we report the first large-scale study of VVs performed on a freely available genetic data of 408,455 European-ancestry individuals. We identified the 12 reliably associated loci that explain 13% of the SNP-based heritability, and prioritized the most likely causal genes CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9. VVs-associated variants within these loci exhibited pleiotropic effects on several phenotypes including blood pressure/hypertension and blood cell traits. Gene set enrichment analysis revealed gene categories related to abnormal vasculogenesis. Genetic correlation analysis confirmed known epidemiological associations between VVs and deep venous thrombosis, weight, rough labor, and standing job, and found a genetic overlap with multiple traits that have not been previously suspected to share common genetic background with VVs. These traits included educational attainment, fluid intelligence and prospective memory scores, walking pace (negative correlation with VVs), smoking, height, number of operations, pain, and gonarthrosis (positive correlation with VVs). Finally, Mendelian randomization analysis provided evidence for causal effects of plasma levels of MICB and CD209 proteins, and anthropometric traits such as waist and hip circumference, height, weight, and both fat and fat-free mass. Our results provide novel insight into both VVs genetics and etiology. The revealed genes and proteins can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

Varicose veins of lower extremities (VVs) affect about 30% of adults in developed countries and cause both cosmetic and health problems. A strong body of evidence indicates that heredity plays an important role in the etiology of this condition. However, genetic basis of VVs remains poorly understood. Here, we present the results of the first large-scale genetic study for VVs. We identified genes which are the most likely involved in VVs pathogenesis. We show that VVs are correlated at a genetic level with numerous traits and phenotypes, including those already known from prior epidemiological studies (deep venous thrombosis, body mass index, standing job, etc.) as well as with those that have not been suspected to share common genetic background with VVs (fluid intelligence and prospective memory scores, smoking, walking pace, pain all over the body, and other traits). Finally, using genetic variants as instruments, we demonstrate direct causal effects of the traits related to anthropometry, such as height and weight, and plasma levels of immune-related proteins MICB and CD209. Our study provides novel insight into both VVs genetics and etiology. The revealed genes (CASZ1, PIEZO1, PPP3R1, EBF1, STIM2, HFE, GATA2, NFATC2, and SOX9) and proteins (MICB and CD209) can be considered as good candidates for follow-up functional studies and might be of interest as potential drug targets.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

Novel biomarkers for cardiovascular risk

Cardiovascular disease refers to different diseases involving the heart and/or the arteries and/or the veins. Cardiovascular disease, overall considered, is a notable source of morbidity and mortality worldwide. Therefore, several research studies are dedicated to explore, by means of biomarkers, the possiblity to calculate the cardiovascular risk both for the onset and for the complications of the related clinical manifestations such as coronary artery disease, carotid artery stenosis, peripheral artery disease, arterial aneurysm, chronic venous disease and venous thromboembolism. This review discusses the most updated information in the area of the novel biomarkers related to omics, imaging techniques and clinical data, that may help physicians in order to improve the knowledge and the management of the cardiovascular risk.

A genome-wide association study identifies a susceptibility locus for biliary atresia on 2p16.1 within the gene EFEMP1

Biliary atresia (BA) is a rare pediatric cholangiopathy characterized by fibrosclerosing obliteration of the extrahepatic bile ducts, leading to cholestasis, fibrosis, cirrhosis, and eventual liver failure. The etiology of BA remains unknown, although environmental, inflammatory, infectious, and genetic risk factors have been proposed. We performed a genome-wide association study (GWAS) in a European-American cohort of 343 isolated BA patients and 1716 controls to identify genetic loci associated with BA. A second GWAS was performed in an independent European-American cohort of 156 patients with BA and other extrahepatic anomalies and 212 controls to confirm the identified candidate BA-associated SNPs. Meta-analysis revealed three genome-wide significant BA-associated SNPs on 2p16.1 (rs10865291, rs6761893, and rs727878; P < 5 ×10-8), located within the fifth intron of the EFEMP1 gene, which encodes a secreted extracellular protein implicated in extracellular matrix remodeling, cell proliferation, and organogenesis. RNA expression analysis showed an increase in EFEMP1 transcripts from human liver specimens isolated from patients with either BA or other cholestatic diseases when compared to normal control liver samples. Immunohistochemistry demonstrated that EFEMP1 is expressed in cholangiocytes and vascular smooth muscle cells in liver specimens from patients with BA and other cholestatic diseases, but it is absent from cholangiocytes in normal control liver samples. Efemp1 transcripts had higher expression in cholangiocytes and portal fibroblasts as compared with other cell types in normal rat liver. The identification of a novel BA-associated locus, and implication of EFEMP1 as a new BA candidate susceptibility gene, could provide new insights to understanding the mechanisms underlying this severe pediatric disorder.

Polymorphisms of genes involved in inflammation and blood vessel development influence the risk of varicose veins

Heredity plays an important role in the etiology of varicose veins (VVs). However, the genetic basis underlying this condition remains poorly understood. Our aim was to replicate top association signals from genome-wide association studies (GWASs) for VVs of lower extremities using 2 independent datasets-our sample of ethnic Russian individuals (709 cases and 278 controls) and a large cohort of British residents from UK Biobank (10 861 cases and 397 594 controls). Associations of polymorphisms rs11121615, rs6712038, rs507666, rs966562, rs7111987, rs6062618, and rs6905288 were validated in the UK Biobank individuals at a Bonferroni-corrected significance level. In Russian cohort, only rs11121615 reached a nominal significance level of P < .05. Results of original GWAS and replication studies were combined by a meta-analysis, and polymorphisms listed above as well as rs111434909 and rs4463578 passed a genome-wide significant threshold. Notably, the majority of these polymorphisms were located within or near genes involved in vascular development and remodeling, and regulation of inflammatory response. Our results confirm the role of these polymorphisms in genetic susceptibility to VVs and indicate the revealed genomic regions as good candidates for further fine-mapping studies and functional analysis. Moreover, our findings implicate inflammation and abnormal vascular architecture in VVs pathogenesis.

Identification of novel hyper- or hypomethylated CpG sites and genes associated with atherosclerotic plaque using an epigenome-wide association study

DNA methylation is an important epigenetic modification that has been implicated in the pathogenesis of atherosclerosis. Although previous studies have identified various CpG sites and genes whose methylation is associated with atherosclerosis in populations with European or Mexican ancestry, the genome‑wide pattern of DNA methylation in the atherosclerotic human aorta is yet to be elucidated in Japanese individuals. In the present study, a genome‑wide analysis of DNA methylation at ~853,000 CpG sites was performed using 128 postmortem aortic intima specimens obtained from 64 Japanese patients. To avoid the effects of interindividual variation, intraindividual paired comparisons were performed between atheromatous plaque lesions and corresponding plaque‑free tissue for each patient. Bisulfite‑modified genomic DNA was analyzed using a specific microarray for DNA methylation. DNA methylation at each CpG site was calculated as the β value, where β = (intensity of the methylated allele)/(intensity of the methylated allele + intensity of the unmethylated allele + 100). Bonferroni's correction for statistical significance of association was applied to compensate for multiple comparisons. The methylation of 2,679 CpG sites differed significantly (P<5.86x10‑8) between atheromatous plaque lesions and the corresponding plaque‑free intima, with 2,272 and 407 CpG sites in atheromatous plaques being hyper‑ or hypomethylated, respectively. A total of 5 hypermethylated CpG sites in atheromatous plaques were demonstrated to have a difference in β value of >0.15 (plaque lesion‑plaque‑free intima) and 11 had a β ratio of >1.50 (plaque/plaque‑free intima). A further 15 and 17 hypomethylated CpG sites in atheromatous plaques were observed to have a difference in β value of <‑0.15 or a β ratio of <0.67, respectively. According to these limits, a total of 16 novel genes that were significantly hyper‑ or hypomethylated in atheromatous plaque lesions compared with the plaque‑free intima were identified in the present study. The results of the present study suggest that the methylation of these genes may contribute to the pathogenesis of atherosclerosis in the Japanese population.