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.

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.

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.

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.

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.

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.

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.

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.

Association of polymorphisms near the FOXC2 gene with the risk of varicose veins in ethnic Russians.

OBJECTIVE: To investigate the association of polymorphisms located near the FOXC2 gene with the risk of varicose veins in ethnic Russians. METHODS: Allele, genotype, and haplotype frequencies were determined in the sample of 474 patients with primary varicose veins and in the control group of 478 individuals without a history of chronic venous disease. RESULTS: Polymorphisms rs7189489, rs4633732, and rs1035550 showed the association with the increased risk of varicose veins, but none of the observed associations remained significant after correction for multiple testing. Haplotype analysis revealed the association of haplotype rs7189489 C-rs4633732 T-rs34221221 C-rs1035550 C-rs34152738 T-rs12711457 G with the increased risk of varicose veins (OR = 2.67, P = 0.01). CONCLUSIONS: Our results provide evidence that the studied polymorphisms do not play a major role in susceptibility to varicose veins development in the Russian population.

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.

TERT polymorphisms rs2853669 and rs7726159 influence on prostate cancer risk in Russian population.

Telomere length and telomerase activity have been hypothesized to play a role in cancer development. The aim of our study was to investigate the association of allelic variants of three functional polymorphisms rs2853669, rs2736100, and rs7726159 in the telomerase reverse transcriptase (TERT) gene with the risk of the breast cancer and prostate cancer in Russian population. Six hundred sixty women with breast cancer, 372 men with prostate cancer, and corresponding control groups of 523 women and 363 men were included in the present case-control study. We observed an association of allele rs2853669 C with increased risk of prostate cancer (co-dominant model TC vs. TT OR = 1.65, P = 0.002; additive model OR = 1.42, P = 0.005; dominant model: OR = 1.64, P = 0.001) and allele rs7726159 A with reduced risk of this malignancy (сo-dominant model: AA vs. CC OR = 0.42, P = 0.002; additive model: OR = 0.69, P = 0.002; dominant model: OR = 0.67, P = 0.01; recessive model: OR = 0.48, P = 0.005). None of the studied polymorphisms showed an association with the risk of breast cancer. Our results provide evidence that the TERT gene variability modulate prostate cancer predisposition in ethnical Russians.