Quantitative and structural characteristics of mitochondrial DNA in varicose veins.

OBJECTIVE: We examined quantitative (in terms of mtDNA/nuclear DNA) and structural (in terms of common deletions in the MT-ND4 gene region) characteristics of mitochondrial DNA (mtDNA) in varicose veins (VVs) and venous wall layers by comparing mitochondrial genome parameters, as well as mitochondrial function (in terms of mitochondrial membrane potential (MtMP)), in varicose vein (VV) vs. non-varicose vein (NV) tissue samples. METHODS: We analyzed paired great saphenous vein samples (VV vs. NV segments from each patient left after venous surgery) harvested from patients with VVs. Relative mtDNA level and the proportion of no-deletion mtDNA were determined by a multiplex quantitative PCR (qPCR), confirming the latter with a more sensitive method - droplet digital PCR (ddPCR). Mitochondria's functional state in VVs was assessed using fluorescent (dependent on MtMP) live-staining of mitochondria in venous tissues. RESULTS: Total mtDNA level was lower in VV than in NV samples (predominantly in the t. media layer). ddPCR analysis showed lower proportion of no-deletion mtDNA in VVs. Because of the decrease in relative MtMP in VVs, our results suggest a possible reduction of mitochondrial function in VVs. CONCLUSION: Quantitative and structural changes (copy number and integrity) of mtDNA are plausibly involved in VV pathogenesis. Future clinical studies implementing the mitochondrial targeting may be eventually fostered after auxiliary mechanistic studies.

Mendelian randomization analysis of plasma levels of CD209 and MICB proteins and the risk of varicose veins of lower extremities.

Varicose veins of lower extremities (VVs) are a highly prevalent condition, the pathogenesis of which is still not fully elucidated. Mendelian randomization (MR) can provide useful preliminary information on the traits that are potentially causally related to the disease. The aim of the present study is to replicate the effects of the plasma levels of MHC class I polypeptide-related sequence B (MICB) and cluster of differentiation 209 (CD209) proteins reported in a previous hypothesis-free MR study. We conducted MR analysis using a fixed effects inverse-variance weighted meta-analysis of Wald ratios method. For MICB and CD209, we used data from a large-scale genome-wide association study (GWAS) for plasma protein levels (N = 3,301). For VVs, we used GWAS data obtained in the FinnGen project (N = 128,698), the eMERGE network (phase 3, N = 48,429), and the UK Biobank data available in the Gene ATLAS (N = 452,264). The data used in the study were obtained in individuals of European descent. The results for MICB did not pass criteria for statistical significance and replication. The results for CD209 passed all statistical significance thresholds, indicating that the genetically predicted increase in CD209 level is associated with increased risk of VVs (ßMR (SE) = 0.07 (0.01), OR (95% CI) = 1.08 (1.05-1.10), P-value = 5.9 ×10-11 in the meta-analysis of three cohorts). Our findings provide further support that CD209 can potentially be involved in VVs. In future studies, independent validation of our results using data from more powerful GWASs for CD209 measured by different methods would be beneficial.

In silico genome-wide gene-based association analysis reveals new genes predisposing to coronary artery disease.

Genome-wide association study (GWAS) have identified more than 300 single nucleotide polymorphisms at 163 independent loci associated with coronary artery disease (CAD). However, there is no full understanding about the causal genes for CAD and the mechanisms of their action. We aimed to perform a post GWAS analysis to identify genes whose polymorphism may influence the risk of CAD. Using the UK Biobank GWAS summary statistics, we performed a gene-based association analysis. We found 63 genes significantly associated with CAD due to their within-gene polymorphisms. Many of these genes are well known. Some known CAD genes such as FURIN and SORT1 did not show the gene-based association because their variants had low GWAS signals or gene-based association was inflated by the strong GWAS signal outside the gene. For several known CAD genes, we demonstrated that their effects could be explained not only or not at all by their own variants but by the variants within the neighboring genes controlling their expression. Using several bioinformatics techniques, we suggested potential mechanisms underlying gene-CAD associations. Three genes, CDK19, NCALD, and ARHGEF12 were not previously associated with CAD. The role of these genes should be clarified in further studies.

Multivariate genome-wide analysis of immunoglobulin G N-glycosylation identifies new loci pleiotropic with immune function.

The N-glycosylation of immunoglobulin G (IgG) affects its structure and function. It has been demonstrated that IgG N-glycosylation patterns are inherited as complex quantitative traits. Genome-wide association studies identified loci harboring genes encoding enzymes directly involved in protein glycosylation as well as loci likely to be involved in regulation of glycosylation biochemical pathways. Many of these loci could be linked to immune functions and risk of inflammatory and autoimmune diseases. The aim of the present study was to discover and replicate new loci associated with IgG N-glycosylation and to investigate possible pleiotropic effects of these loci onto immune function and the risk of inflammatory and autoimmune diseases. We conducted a multivariate genome-wide association analysis of 23 IgG N-glycosylation traits measured in 8090 individuals of European ancestry. The discovery stage was followed up by replication in 3147 people and in silico functional analysis. Our study increased the total number of replicated loci from 22 to 29. For the discovered loci, we suggest a number of genes potentially involved in the control of IgG N-glycosylation. Among the new loci, two (near RNF168 and TNFRSF13B) were previously implicated in rare immune deficiencies and were associated with levels of circulating immunoglobulins. For one new locus (near AP5B1/OVOL1), we demonstrated a potential pleiotropic effect on the risk of asthma. Our findings underline an important link between IgG N-glycosylation and immune function and provide new clues to understanding their interplay.

Replication of 15 loci involved in human plasma protein N-glycosylation in 4802 samples from four cohorts.

Human protein glycosylation is a complex process, and its in vivo regulation is poorly understood. Changes in glycosylation patterns are associated with many human diseases and conditions. Understanding the biological determinants of protein glycome provides a basis for future diagnostic and therapeutic applications. Genome-wide association studies (GWAS) allow to study biology via a hypothesis-free search of loci and genetic variants associated with a trait of interest. Sixteen loci were identified by three previous GWAS of human plasma proteome N-glycosylation. However, the possibility that some of these loci are false positives needs to be eliminated by replication studies, which have been limited so far. Here, we use the largest set of samples so far (4802 individuals) to replicate the previously identified loci. For all but one locus, the expected replication power exceeded 95%. Of the 16 loci reported previously, 15 were replicated in our study. For the remaining locus (near the KREMEN1 gene), the replication power was low, and hence, replication results were inconclusive. The very high replication rate highlights the general robustness of the GWAS findings as well as the high standards adopted by the community that studies genetic regulation of protein glycosylation. The 15 replicated loci present a good target for further functional studies. Among these, eight loci contain genes encoding glycosyltransferases: MGAT5, B3GAT1, FUT8, FUT6, ST6GAL1, B4GALT1, ST3GAL4 and MGAT3. The remaining seven loci offer starting points for further functional follow-up investigation into molecules and mechanisms that regulate human protein N-glycosylation in vivo.

Prioritization of causal genes for coronary artery disease based on cumulative evidence from experimental and in silico studies.

Genome-wide association studies have led to a significant progress in identification of genomic loci affecting coronary artery disease (CAD) risk. However, revealing the causal genes responsible for the observed associations is challenging. In the present study, we aimed to prioritize CAD-relevant genes based on cumulative evidence from the published studies and our own study of colocalization between eQTLs and loci associated with CAD using SMR/HEIDI approach. Prior knowledge of candidate genes was extracted from both experimental and in silico studies, employing different prioritization algorithms. Our review systematized information for a total of 51 CAD-associated loci. We pinpointed 37 genes in 36 loci. For 27 genes we infer they are causal for CAD, and for 10 further genes we judge them most likely causal. Colocalization analysis showed that for 18 out of these loci, association with CAD can be explained by changes in gene expression in one or more CAD-relevant tissues. Furthermore, for 8 out of 36 loci, existing evidence suggested additional CAD-associated genes. For the remaining 15 loci, we concluded that evidence for gene prioritization remains inconsistent, insufficient, or absent. Our results provide deeper insights into the genetic etiology of CAD and demonstrate knowledge gaps where further research is warranted.

Analysis of genetically independent phenotypes identifies shared genetic factors associated with chronic musculoskeletal pain conditions.

Chronic musculoskeletal pain affects all aspects of human life. However, mechanisms of its genetic control remain poorly understood. Genetic studies of pain are complicated by the high complexity and heterogeneity of pain phenotypes. Here, we apply principal component analysis to reduce phenotype heterogeneity of chronic musculoskeletal pain at four locations: the back, neck/shoulder, hip, and knee. Using matrices of genetic covariances, we constructed four genetically independent phenotypes (GIPs) with the leading GIP (GIP1) explaining 78.4% of the genetic variance of the analyzed conditions, and GIP2-4 explain progressively less. We identified and replicated five GIP1-associated loci and one GIP2-associated locus and prioritized the most likely causal genes. For GIP1, we showed enrichment with multiple nervous system-related terms and genetic correlations with anthropometric, sociodemographic, psychiatric/personality traits and osteoarthritis. We suggest that GIP1 represents a biopsychological component of chronic musculoskeletal pain, related to physiological and psychological aspects and reflecting pain perception and processing.

Absence of EGFR C797S Mutation in Tyrosine Kinase Inhibitor-Naïve Non-Small Cell Lung Cancer Tissues.

EGFR tyrosine-kinase inhibitors (TKIs) are used as targeted therapeutics for the treatment of advanced non-small cell lung cancer (NSCLC) with EGFR-activating mutations. EGFR C797S is common causes of acquired resistance to third-generation TKIs. There is wide-spread opinion that resistance-conferring mutation present even in a small proportion of cancer cells before the start of therapy could potentially predict poor response to a targeted drug. In our study, we tested whether C797S can be found in previously untreated NSCLCs. We analyzed DNA samples extracted from formalin-fixed paraffin-embedded (FFPE) tumor tissue sections of 470 lung adenocarcinoma patients, including 235 samples with activating EGFR mutations. Screening was performed using highly sensitive droplet digital PCR assay. No tumor samples with baseline C797S were identified. C797S does not occur in TKI-naïve NSCLCs and provide evidence that screening for this mutation before TKIs administration may not be necessary.

Defining the genetic control of human blood plasma N-glycome using genome-wide association study.

Glycosylation is a common post-translational modification of proteins. Glycosylation is associated with a number of human diseases. Defining genetic factors altering glycosylation may provide a basis for novel approaches to diagnostic and pharmaceutical applications. Here we report a genome-wide association study of the human blood plasma N-glycome composition in up to 3811 people measured by Ultra Performance Liquid Chromatography (UPLC) technology. Starting with the 36 original traits measured by UPLC, we computed an additional 77 derived traits leading to a total of 113 glycan traits. We studied associations between these traits and genetic polymorphisms located on human autosomes. We discovered and replicated 12 loci. This allowed us to demonstrate an overlap in genetic control between total plasma protein and IgG glycosylation. The majority of revealed loci contained genes that encode enzymes directly involved in glycosylation (FUT3/FUT6, FUT8, B3GAT1, ST6GAL1, B4GALT1, ST3GAL4, MGAT3 and MGAT5) and a known regulator of plasma protein fucosylation (HNF1A). However, we also found loci that could possibly reflect other more complex aspects of glycosylation process. Functional genomic annotation suggested the role of several genes including DERL3, CHCHD10, TMEM121, IGH and IKZF1. The hypotheses we generated may serve as a starting point for further functional studies in this research area.

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.

Mycoplasma hyorhinis reduces sensitivity of human lung carcinoma cells to Nutlin-3 and promotes their malignant phenotype.

PURPOSE: MDM2 inhibitors are promising anticancer agents that induce cell cycle arrest and tumor cells death via p53 reactivation. We examined the influence of Mycoplasma hyorhinis infection on sensitivity of human lung carcinoma cells NCI-H292 to MDM2 inhibitor Nutlin-3. In order to unveil possible mechanisms underlying the revealed effect, we investigated gene expression changes and signal transduction networks activated in NCI-H292 cells in response to mycoplasma infection. METHODS: Sensitivity of NCI-Н292 cells to Nutlin-3 was estimated by resazurin-based cell viability assay. Genome-wide transcriptional profiles of NCI-H292 and NCI-Н292Myc.h cell lines were determined using Illumina Human HT-12 v3 Expression BeadChip. Search for key transcription factors and key node molecules was performed using the geneXplain platform. Ability for anchorage-independent growth was tested by soft agar colony formation assay. RESULTS: NCI-Н292Myc.h cells were shown to be 1.5- and 5.2-fold more resistant to killing by Nutlin-3 at concentrations of 15 and 30 µM than uninfected NCI-Н292 cells (P < 0.05 and P < 0.001, respectively). Transcriptome analysis revealed differential expression of multiple genes involved in cancer progression and metastasis as well as epithelial-mesenchymal transition (EMT). Moreover, we have shown experimentally that NCI-Н292Myc.h cells were more capable of growing and dividing without binding to a substrate. The most likely mechanism explaining the observed changes was found to be TLR4- and IL-1b-mediated activation of NF-κB pathway. CONCLUSIONS: Our results provide evidence that mycoplasma infection is an important factor modulating the effect of MDM2 inhibitors on cancer cells and is able to induce EMT-related changes.

Genome-wide association study in ethnic Russians suggests an association of the MHC class III genomic region with the risk of primary varicose veins.

Heredity is a well-known risk factor for varicose veins, but genetic basis of this condition remains poorly studied. Our aim was to conduct a large-scale genetic association study for primary varicose veins (PVVs) in the population of ethnic Russians. An initial scan using Illumina HumanExome-12 v1.0 BeadChip was performed for 273 patients with PVVs and 250 controls without a history of chronic venous disease and other venous disorders. After quality control and removal of monomorphic markers, 25,424 common and 48,232 rare variants were included in the analysis. 42 single nucleotide polymorphisms (SNPs) were genotyped in the independent replication cohort of 447 PVVs patients and 443 controls. Association of common variants with PVVs was investigated by logistic regression, and the impact of rare variants was analyzed using sequence kernel association test. No effect of low frequency alleles has been revealed in our study. Common variant analysis identified a promising signal at chromosome 6 within classical major histocompatibility complex (MHC) class III subregion. The most strongly associated SNP in a combined analysis that reached a suggestive significance level of 3.2e-05 was polymorphism rs4151657 in the complement factor B gene. Testing for potential pleiotropy with other traits indicated that the same causal variant in this region increases the risk of rheumatoid arthritis and has a negative impact on human height. Our results provide suggestive evidence for the involvement of the MHC class III genes in the pathogenesis of PVVs. Further independent studies are needed to confirm our pilot findings.

Allele rs2010963 C of the VEGFA gene is associated with the decreased risk of primary varicose veins in ethnic Russians.

Objective To study the association of polymorphisms rs699947, rs2010963, rs3025039 in the VEGFA gene region and rs1870377, rs2305949, rs2071559 in the VEGFR2 gene region with the risk of primary varicose veins in ethnic Russians. Methods Genotypes were determined by real-time PCR allelic discrimination. The case group consisted of 448 patients with primary varicose veins and the control group comprised 609 individuals without a history of chronic venous disease. Association was studied by logistic regression analysis. Results Allele rs2010963 C was associated with the decreased risk of varicose veins (additive model of inheritance: odds ratio = 0.73, 95% confidence interval = 0.59-0.91, P = 0.004). Conclusions Our results provide evidence that polymorphism rs2010963 located in the 5' untranslated region of the VEGFA gene can influence genetic susceptibility to primary varicose veins in Russians. Otherwise, it can be in linkage disequilibrium with another functional single nucleotide polymorphism that can alter the level of vascular endothelial growth factor A protein.

Polymorphisms in inflammation-related genes and the risk of primary varicose veins in ethnic Russians.

Inflammation was shown to be activated in varicose veins, although its role in the development of vein wall transformation remains inconclusive. We aimed to investigate the influence of 13 inflammation-related single nucleotide polymorphisms (SNPs) TNF rs1800629 and rs3093661, IL1A rs1800587, IL1RN rs4251961, IL6 rs1800795 and rs1800796, IFNG rs2430561, IL10 rs1800896, TGFB1 rs1800469, HIF1A rs11549465, NFKB1 rs28362491, and rs4648068 on the risk of primary varicose veins (PVVs) in ethnic Russians. We genotyped 709 patients with PVVs and 278 individuals without a history of chronic venous disease and performed a single SNP and a haplotype analysis. Several associations with P < 0.05 were revealed in our study. Variant allele HIF1A rs11549465 T, TNF rs3093661 A, and NFKB1 rs28362491 ATTG deletion showed the reverse association with PVV risk, and allele IL6 rs1800795 C was associated with the increased risk of the studied pathology. Haplotype analysis revealed associations of TNF haplotypes rs3093661 A-rs1800629 G and IL6 rs1800795 C-rs1800796 G with the decreased and the increased risk of PVVs, correspondingly. However, all the observed associations failed to reach statistical significance after the correction for multiple testing, which was set at a level of 10-3 due to many tests performed. Our study therefore provides evidence that investigated polymorphisms do not play a major role in susceptibility to PVVs.

Functional polymorphism rs1024611 in the MCP1 gene is associated with the risk of varicose veins of lower extremities.

OBJECTIVE: Monocyte chemoattractant protein 1 (MCP-1) is a chemokine responsible for monocyte, basophil, and T-lymphocyte attraction. Polymorphism rs1024611 located in the regulatory region of the MCP1 gene has previously been shown to be associated with increased MCP-1 production. In our study, we aimed to examine the association of rs1024611 with the risk of primary varicose veins (PVVs) of lower extremities. METHODS: The case group comprised 470 patients with PVVs, and the control group included 269 individuals without a history of chronic venous disease. All cases and controls were ethnic Russians. Genotypes were determined by real-time polymerase chain reaction allelic discrimination. Association was studied by logistic regression analysis. RESULTS: We revealed the association of genotype G/G with the increased risk of PVVs (odds ratio [OR], 1.87; 95% confidence interval [CI], 1.02-3.44; P = .04). In the subgroup analysis, association was revealed only in patients with C2 Clinical, Etiology, Anatomy, and Pathophysiology class (allele G: OR, 1.62 [95% CI, 1.13-2.33; P = .008]; genotype G/G: OR, 3.22 [95% CI, 1.43-7.27; P = .005]), in patients with age at onset of PVVs before 30 years (allele G: OR, 1.41 [95% CI, 1.08-1.85; P = .01]; genotype G/G: OR, 2.35 [95% CI, 1.22-4.55; P = .01]), and in patients who declared no family history (allele G: OR, 1.46 [95% CI, 1.02-2.09; P = .04]; genotype G/G: OR, 2.50 [95% CI, 1.11-5.63; P = .03]). CONCLUSIONS: Our results provide evidence for MCP-1 involvement in the development of PVVs and indicate that inflammation could be implicated in the pathogenesis of this condition.

Allele and genotype frequencies of polymorphisms in cytokine genes in ethnic Russian individuals from Moscow, Russia.

Two hundred and twenty eight ethnic Russian individuals from Moscow, Russia, were genotyped at 14 single nucleotide polymorphisms CCL2 A-2578G; VEGFA C-2578A, G-634C, and C+936T; TNF G+419A and G-308A; IL1A G-889A; IL1RN T+1018C; IL6G-174C and G-572C; IFNG T+874A; IL1B C-511T; IL10 A+1082G; TGFB1 C-509T. Genotypes were determined using real-time polymerase chain reaction with TaqMan probes and polymerase chain reaction followed by melting analysis of dual-labeled probe. Genotype distribution was in accordance with Hardy-Weinberg equilibrium for all studied polymorphisms. Genotype data are available in the Allele Frequencies Net Database under identifier AFND 3367 and the population name "Russia Moscow Cytokine".

Polymorphisms in the MTHFR and MTR genes and the risk of varicose veins in ethnical Russians.

OBJECTIVES: The objective of this study was to study the association of polymorphisms MTHFR C677T (rs1801133) and MTR A2756G (rs1805087) with the risk of varicose veins in ethnical Russians. METHODS: We genotyped 475 patients with varicose veins, 168 individual without chronic venous disease, and the population-based group of 896 subjects. Association was studied using logistic regression analysis adopting co-dominant, additive, recessive, and dominant models of inheritance. RESULTS: None of the polymorphisms showed a statistically significant association with the risk of varicose veins. CONCLUSIONS: Our results provide evidence that the studied polymorphisms do not contribute to genetic susceptibility to varicose veins in ethnical Russians.

Polymorphisms in DNA repair genes and breast cancer risk in Russian population: a case-control study.

Genetic variation in DNA repair genes can alter an individual's capacity to repair damaged DNA and influence the risk of cancer. We tested seven polymorphisms in DNA repair genes XRCC1, ERCC2, XRCC3, XRCC2, EXOI and TP53 for a possible association with breast cancer risk in a sample of 672 case and 672 control Russian women. An association was observed for allele A of the polymorphism XRCC1 (R399Q) rs25487 (co-dominant model AA vs. GG: OR 1.76, P = 0.003; additive model OR 1.28, P = 0.005; dominant model: OR 1.29, P = 0.03; recessive model OR 1.63, P = 0.008). Allele T of the polymorphism ERCC2 (D312N) rs1799793 was also associated with breast cancer risk (co-dominant model TT vs. CC: OR 1.43, P = 0.04; additive model OR 1.21, P = 0.02; dominant model: OR 1.30, P = 0.02), but the association became insignificant after applying Bonferroni correction. No association with breast cancer was found for the remaining SNPs. In summary, our study provides evidence that polymorphisms in DNA repair genes may play a role in susceptibility to breast cancer in the population of ethnical Russians.

HFE p.C282Y gene variant is associated with varicose veins in Russian population.

Recently, the association of polymorphism rs1800562 (p.C282Y) in the hemochromatosis (HFE) gene with the increased risk of venous ulceration was shown. We hypothesized that HFE gene polymorphism might be involved not only in ulceration process, but also in susceptibility to primary varicose veins. We genotyped HFE p.C282Y (rs1800562) and p.H63D (rs1799945) variants in patients with primary varicose veins (n = 463) and in the control group (n = 754). In our study, p.282Y variant (rs1800562 A allele) was significantly associated with the risk of varicose veins (OR 1.79, 95 % CI = 1.11-2.89, P = 0.02). A borderline significant reverse association of p.63D variant (rs1799945 G allele) with venous leg ulcer development was revealed in Russians (OR 0.25, 95 % CI = 0.06-1.00, P = 0.05), but not in the meta-analysis (P = 0.56). We conclude that the HFE gene polymorphism can affect the risk of developing primary varicose veins.