Genetic and Environmental Factors in Autoimmune Thyroid Disease: Exploring Associations with Selenium Levels and Novel Loci in a Latvian Cohort.

The interplay of genetic, immune and environmental factors strongly contributes to the development of autoimmune thyroid disease (AITD), which can be classified as Graves' disease (GD) or Hashimoto thyroiditis (HT). One of the most studied exogenous factors in the pathogenesis of AITD is selenium, which, in the form of selenoproteins, strengthens the antioxidative defence system of thyroid cells against superoxide production. Furthermore, it modulates inflammatory cytokine release and autoantibody production. The aim of this study was to assess the associations of genetic factors with selenium levels in a cohort of adults with HT and GD and healthy controls from Latvia. A total of 148 GD patients, 102 HT patients and 2442 control participants were included in the study. The genotypes were determined using genome-wide genotyping; imputation was carried out using the TOPMed r2 imputation panel; and association analysis was performed with PLINK v1.90b7. We found three loci associated with GD (LSAMP, HNRNPA3P5, and NTN1) and one locus associated with HT (VAT1L); furthermore, one locus was associated with a serum selenium concentration > 80 µg/L (LINC01544/RNF152/PIGN). The detected associations could be attributed to population-specific effects or unknown stratification in our cohort, and further assessment of these results is required to explain the relationships of genetic traits with AITD and other phenotypes.

Evaluating the Efficacy of Type 2 Diabetes Polygenic Risk Scores in an Independent European Population.

Numerous type 2 diabetes (T2D) polygenic risk scores (PGSs) have been developed to predict individuals' predisposition to the disease. An independent assessment and verification of the best-performing PGS are warranted to allow for a rapid application of developed models. To date, only 3% of T2D PGSs have been evaluated. In this study, we assessed all (n = 102) presently published T2D PGSs in an independent cohort of 3718 individuals, which has not been included in the construction or fine-tuning of any T2D PGS so far. We further chose the best-performing PGS, assessed its performance across major population principal component analysis (PCA) clusters, and compared it with newly developed population-specific T2D PGS. Our findings revealed that 88% of the published PGSs were significantly associated with T2D; however, their performance was lower than what had been previously reported. We found a positive association of PGS improvement over the years (p-value = 8.01 × 10-4 with PGS002771 currently showing the best discriminatory power (area under the receiver operating characteristic (AUROC) = 0.669) and PGS003443 exhibiting the strongest association PGS003443 (odds ratio (OR) = 1.899). Further investigation revealed no difference in PGS performance across major population PCA clusters and when compared with newly developed population-specific PGS. Our findings revealed a positive trend in T2D PGS performance, consistently identifying high-T2D-risk individuals in an independent European population.

Whole-Genome Sequencing of 502 Individuals from Latvia: The First Step towards a Population-Specific Reference of Genetic Variation.

Despite rapid improvements in the accessibility of whole-genome sequencing (WGS), understanding the extent of human genetic variation is limited by the scarce availability of genome sequences from underrepresented populations. Developing the population-scale reference database of Latvian genetic variation may fill the gap in European genomes and improve human genomics research. In this study, we analysed a high-coverage WGS dataset comprising 502 individuals selected from the Genome Database of the Latvian Population. An assessment of variant type, location in the genome, function, medical relevance, and novelty was performed, and a population-specific imputation reference panel (IRP) was developed. We identified more than 18.2 million variants in total, of which 3.3% so far are not represented in gnomAD and dbSNP databases. Moreover, we observed a notable though distinct clustering of the Latvian cohort within the European subpopulations. Finally, our findings demonstrate the improved performance of imputation of variants using the Latvian population-specific reference panel in the Latvian population compared to established IRPs. In summary, our study provides the first WGS data for a regional reference genome that will serve as a resource for the development of precision medicine and complement the global genome dataset, improving the understanding of human genetic variation.

Solitary fibrous tumor with IGF-II-induced non-islet cell tumor hypoglycemia: a case report and molecular characterization by next-generation sequencing.

Background: Non-islet cell tumor-induced hypoglycemia (NICTH) is a rare, life-threatening medical condition caused by excessive insulin-like growth factor II (IGF-II) secretion from tumors of most commonly mesenchymal origin. Using next-generation sequencing, we have characterized the genome and transcriptome of the resected IGF-II-secreting solitary fibrous tumor from a patient with severe hypoglycemia accompanied by hypoglycemia unawareness. Case presentation: A 69-year-old male patient presenting with abdominal discomfort was examined using computer tomography, revealing a large lesion at the lesser pelvis extending above the umbilicus. As no bone and lymph node metastases were detected, the patient was scheduled for laparotomy. Before surgery, the patient presented with symptoms of severe hypoglycemia. Suppressed C-peptide levels and subsequent hypokalemia indicated a possible case of NICTH. The patient was treated with methylprednisolone (8 mg) to assess hypoglycemia. After the surgery, mild hypoglycemia was present for the postoperative period, and no radiological recurrences were observed 3 and 12 months after discharge. Histopathological examination results were consistent with the diagnosis of malignant solitary fibrous tumor (SFT). Overexpression of IGF-II was confirmed by both immunohistochemistry and RNA sequencing. Further NGS analysis revealed an SFT characteristic alteration-NAB2-STAT6 fusion. Additionally, three deleterious missense variants were detected in oncogenes BIRC6, KIT, and POLQ, and one homozygous in-frame deletion in the RBM10 tumor suppressor gene. Conclusion: While the NAB2-STAT6 fusions are well characterized, the mutational landscape of SFTs remains understudied. This study reports the importance of NGS to characterize SFTs as we detected four coding variants in genes (BIRC6, KIT, POLQ, and RBM10) associated with tumorigenesis that could potentially contribute to the overall pathogenesis of SFT.

Transcriptome of GH-producing pituitary neuroendocrine tumours and models are significantly affected by somatostatin analogues.

Pituitary neuroendocrine tumours (PitNETs) are neoplasms of the pituitary that overproduce hormones or cause unspecific symptoms due to mass effect. Growth hormone overproducing GH-producing PitNETs cause acromegaly leading to connective tissue, metabolic or oncologic disorders. The medical treatment of acromegaly is somatostatin analogues (SSA) in specific cases combined with dopamine agonists (DA), but almost half of patients display partial or full SSA resistance and potential causes of this are unknown. In this study we investigated transcriptomic landscape of GH-producing PitNETs on several levels and functional models-tumour tissue of patients with and without SSA preoperative treatment, tumour derived pituispheres and GH3 cell line incubated with SSA to study effect of medication on gene expression. MGI sequencing platform was used to sequence total RNA from PitNET tissue, pituispheres, mesenchymal stromal stem-like cells (MSC), and GH3 cell cultures, and data were analysed with Salmon-DeSeq2 pipeline. We observed that the GH-producing PitNETs have distinct changes in growth hormone related pathways related to its functional status alongside inner cell signalling, ion transport, cell adhesion and extracellular matrix characteristic patterns. In pituispheres model, treatment regimens (octreotide and cabergoline) affect specific cell proliferation (MKI67) and core functionality pathways (RYR2, COL8A2, HLA-G, ARFGAP1, TGFBR2). In GH3 cells we observed that medication did not have transcriptomic effects similar to preoperative treatment in PitNET tissue or pituisphere model. This study highlights the importance of correct model system selection for cell transcriptomic profiling and data interpretation that could be achieved in future by incorporating NGS methods and detailed cell omics profiling in PitNET model research.

Genome wide analysis of circulating miRNAs in growth hormone secreting pituitary neuroendocrine tumor patients' plasma.

Background: Circulating plasma miRNAs have been increasingly studied in the field of pituitary neuroendocrine tumor (PitNET) research. Our aim was to discover circulating plasma miRNAs species associated with growth hormone (GH) secreting PitNETs versus assess how the plasma levels of discovered miRNA candidates are impacted by SSA therapy and whether there is a difference in their levels between GH secreting PitNETs versus other PitNET types and healthy individuals. Design: We compared plasma miRNA content and levels before and after surgery focusing on GH secreting PitNET patients. Selected miRNA candidates from our data and literature were then tested in a longitudinal manner in somatostatin analogues (SSA) treatment group. Additionally, we validated selected targets in an independent GH secreting PitNET group. Methods: miRNA candidates were discovered using the whole miRNA sequencing approach and differential expression analysis. Selected miRNAs were then analyzed using real-time polymerase chain reaction (qPCR). Results: Whole miRNA sequencing discovered a total of 16 differentially expressed miRNAs (DEMs) in GH secreting PitNET patients' plasma 24 hours after surgery and 19 DEMs between GH secreting PitNET patients' plasma and non-functioning (NF) PitNET patients' plasma. Seven miRNAs were selected for further testing of which miR-625-5p, miR-503-5p miR-181a-2-3p and miR-130b-3p showed a significant downregulation in plasma after 1 month of SSA treatment. mir-625-5p was found to be significantly downregulated in plasma of GH secreting PitNET patients vs. NF PitNET patients. miR-625-5p alongside miR-130b-3p were also found to be downregulated in GH PitNETs compared to healthy individuals. Conclusions: Our study suggests that expression of plasma miRNAs miR-625-5p, miR-503-5p miR-181a-2-3p and miR-130b-3p in GH secreting PitNETs is affected by SSA treatment. Additionally, miR-625-5p can distinguish GH secreting PitNETs from other PitNET types and healthy controls warranting further research on these miRNAs for treatment efficacy.

Whole exome sequencing reveals novel risk genes of pituitary neuroendocrine tumors.

Somatic genetic alterations in pituitary neuroendocrine tumors (PitNET) tissues have been identified in several studies, but detection of overlapping somatic PitNET candidate genes is rare. We sequenced and by employing multiple data analysis methods studied the exomes of 15 PitNET patients to improve discovery of novel factors involved in PitNET development. PitNET patients were recruited to the study before PitNET removal surgery. For each patient, two samples for DNA extraction were acquired: venous blood and PitNET tissue. Exome sequencing was performed using Illumina NexSeq 500 sequencer and data analyzed using two separate workflows and variant calling algorithms: GATK and Strelka2. A combination of two data analysis pipelines discovered 144 PitNET specific somatic variants (mean = 9.6, range 0-19 per PitNET) of which all were SNVs. Also, we detected previously known GNAS PitNET mutation and identified somatic variants in 11 genes, which have contained somatic variants in previous WES and WGS studies of PitNETs. Noteworthy, this is the third study detecting somatic variants in gene RYR1 in the exomes of PitNETs. In conclusion, we have identified two novel PitNET candidate genes (AC002519.6 and AHNAK) with recurrent somatic variants in our PitNET cohort and found 13 genes overlapping from previous PitNET studies that contain somatic variants. Our study demonstrated that the use of multiple sequencing data analysis pipelines can provide more accurate identification of somatic variants in PitNETs.

Case Report: Micro-RNAs in Plasma From Bilateral Inferior Petrosal Sinus Sampling and Peripheral Blood From Corticotroph Pituitary Neuroendocrine Tumors.

Objective: Circulating miRNAs are found in bodily fluids including plasma and can serve as biomarkers for diseases. The aim of this study was to provide the first insight into the landscape of circulating miRNAs in close proximity to the adrenocorticotropic hormone (ACTH) secreting PitNET. To achieve this objective next-generation sequencing of miRNAs in plasma from bilateral inferior petrosal sinus sampling (BIPSS) - a gold standard in diagnosing ACTH-secreting PitNETs was carried out and selected miRNA candidates were further tested by RT-qPCR in independent patient cohorts. Methods: Sinistral (left) and dextral (right) BIPSS blood samples of the patient were collected in three time points: before the administration of corticotropin-releasing hormone, 5 and 15 minutes after stimulation. In differential expression analysis, sinistral plasma was compared with dextral. The selected miRNA candidates were tested in plasma by RT-qPCR in two patient groups: 1) in five ACTH secreting PitNET patients with plasma samples taken before and 24 hours after surgery, 2) in 12 ACTH secreting PitNET patients vs. 9 non-functioning PitNET patients. Results: BIPSS concluded that the highest amount of ACTH was released in the sinistral side at the 5th minute mark indicating a presence of a tumor. The highest amount of differentially expressed miRNAs was observed 5 minutes after stimulation (20 upregulated, 14 downregulated). At the 5th minute mark in sinistral plasma, two miRNAs were identified: hsa-miR-7-5p and hsa-miR-375-3p that were highly upregulated compared to other BIPSS samples and peripheral plasma samples. Further testing by qPCR revealed significant reduction of miR-7-5p in plasma 24 hours after surgery and upregulation in plasma of ACTH secreting PitNET patients compared to non-functioning PitNET patients (P =0.0013). Conclusions: By stimulating the ACTH secreting PitNET with CRH a rapid increase of two miRNAs (hsa-mir-7-5p, hsa-mir-375-3p) and ACTH can be observed in sinistral inferior petrosal (tumor side). A decrease of miR-7-5p in plasma after surgery and upregulation in plasma of ACTH secreting PitNET patients was discovered implying that further studies of this miRNA as diagnostic marker is needed.

Toll-Like Receptor 1 Locus Re-examined in a Genome-Wide Association Study Update on Anti-Helicobacter pylori IgG Titers.

BACKGROUND & AIMS: A genome-wide significant association between anti-Helicobacter pylori (H pylori) IgG titers and Toll-like receptor (TLR1/6/10) locus on 4p14 was demonstrated for individuals of European ancestry, but not uniformly replicated. We re-investigated this association in an updated genome-wide association study (GWAS) meta-analysis for populations with low gastric cancer incidence, address potential causes of cohort heterogeneity, and explore functional implications of genetic variation at the TLR1/6/10 locus. METHODS: The dichotomous GWAS (25% individuals exhibiting highest anti-H pylori IgG titers vs remaining 75%) included discovery and replication sampls of, respectively, n = 15,685 and n = 9676, all of European ancestry. Longitudinal analysis of serologic data was performed on H pylori-eradicated subjects (n = 132) and patients under surveillance for premalignant gastric lesions (n = 107). TLR1/6/10 surface expression, TLR1 mRNA, and cytokine levels were measured in leukocyte subsets of healthy subjects (n = 26) genotyped for TLR1/6/10 variants. RESULTS: The association of the TLR1/6/10 locus with anti-H pylori IgG titers (rs12233670; ß = -0.267 ± SE 0.034; P = 4.42 × 10-15) presented with high heterogeneity and failed replication. Anti-H pylori IgG titers declined within 2-4 years after eradication treatment (P = 0.004), and decreased over time in patients with premalignant gastric lesions (P < 0.001). Variation at the TLR1/6/10 locus affected TLR1-mediated cytokine production and TLR1 surface expression on monocytes (P = 0.016) and neutrophils (P = 0.030), but not mRNA levels. CONCLUSIONS: The association between anti-H pylori IgG titers and TLR1/6/10 locus was not replicated across cohorts, possibly owing to dependency of anti-H pylori IgG titers on therapy, clearance, and antibody decay. H pylori-mediated immune cell activation is partly mediated via TLR1 signaling, which in turn is affected by genetic variation.

Large Scale Molecular Studies of Pituitary Neuroendocrine Tumors: Novel Markers, Mechanisms and Translational Perspectives.

Pituitary neuroendocrine tumors (PitNETs) are non-metastatic neoplasms of the pituitary, which overproduce hormones leading to systemic disorders, or tumor mass effects causing headaches, vertigo or visual impairment. Recently, PitNETs have been investigated in large scale (exome and genome) molecular analyses (transcriptome microarrays and sequencing), to uncover novel markers. We performed a literature analysis on these studies to summarize the research data and extrapolate overlapping gene candidates, biomarkers, and molecular mechanisms. We observed a tendency in samples with driver mutations (GNAS, USP8) to have a smaller overall mutational rate, suggesting driver-promoted tumorigenesis, potentially changing transcriptome profiles in tumors. However, direct links from drivers to signaling pathways altered in PitNETs (Notch, Wnt, TGF-ß, and cell cycle regulators) require further investigation. Modern technologies have also identified circulating nucleic acids, and pinpointed these as novel PitNET markers, i.e., miR-143-3p, miR-16-5p, miR-145-5p, and let-7g-5p, therefore these molecules must be investigated in the future translational studies. Overall, large-scale molecular studies have provided key insight into the molecular mechanisms behind PitNET pathogenesis, highlighting previously reported molecular markers, bringing new candidates into the research field, and reapplying traditional perspectives to newly discovered molecular mechanisms.

Novel susceptibility loci identified in a genome-wide association study of type 2 diabetes complications in population of Latvia.

BACKGROUND: Type 2 diabetes complications cause a serious emotional and economical burden to patients and healthcare systems globally. Management of both acute and chronic complications of diabetes, which dramatically impair the quality of patients' life, is still an unsolved issue in diabetes care, suggesting a need for early identification of individuals with high risk for developing diabetes complications. METHODS: We performed a genome-wide association study in 601 type 2 diabetes patients after stratifying them according to the presence or absence of four types of diabetes complications: diabetic neuropathy, diabetic nephropathy, macrovascular complications, and ophthalmic complications. RESULTS: The analysis revealed ten novel associations showing genome-wide significance, including rs1132787 (GYPA, OR = 2.71; 95% CI = 2.02-3.64) and diabetic neuropathy, rs2477088 (PDE4DIP, OR = 2.50; 95% CI = 1.87-3.34), rs4852954 (NAT8, OR = 2.27; 95% CI = 2.71-3.01), rs6032 (F5, OR = 2.12; 95% CI = 1.63-2.77), rs6935464 (RPS6KA2, OR = 2.25; 95% CI = 6.69-3.01) and macrovascular complications, rs3095447 (CCDC146, OR = 2.18; 95% CI = 1.66-2.87) and ophthalmic complications. By applying the targeted approach of previously reported susceptibility loci we managed to replicate three associations: MAPK14 (rs3761980, rs80028505) and diabetic neuropathy, APOL1 (rs136161) and diabetic nephropathy. CONCLUSIONS: Together these results provide further evidence for the implication of genetic factors in the development of type 2 diabetes complications and highlight several potential key loci, able to modify the risk of developing these conditions. Moreover, the candidate variant approach proves a strong and consistent effect for multiple variants across different populations.

Pituispheres Contain Genetic Variants Characteristic to Pituitary Adenoma Tumor Tissue.

The most common type of pituitary neoplasms is benign pituitary adenoma (PA). Clinically significant PAs affect around 0.1% of the population. Currently, there is no established human PA cell culture available and when PA tumor cells are cultured they form two distinct types depending on culturing conditions either free-floating aggregates also known as pituispheres or cells adhering to the surface of cell plates and displaying mesenchymal stem-like properties. The aim of this study was to trace the origin of sphere-forming and adherent pituitary cell cultures and characterize the potential use of these surgery derived cell lines as PA model. We carried out a paired-end exome sequencing of patients' tumor and germline DNA using Illumina NextSeq followed by characterization of corresponding PA cell cultures. Variation analysis revealed a low amount of somatic mutations (mean = 5.2, range 3-7) in exomes of PAs. Somatic mutations of the primary surgery material can be detected in the exomes of respective pituispheres, but not in exomes of respective mesenchymal stem-like cells. For the first time, we show that the genome of pituispheres represents genome of PA while mesenchymal stem cells derived from the PA tissue do not contain mutations characteristic to PA in their genome, therefore, most likely representing normal cells of pituitary or surrounding tissues. This finding indicates that pituispheres can be used as a human model of PA cells, but combination of cell culturing techniques and NGS needs to be employed to adjust for disability to propagate spheres in culturing conditions.

Case report: recurrent pituitary adenoma has increased load of somatic variants.

BACKGROUND: Pituitary adenomas (PA) have an increased potential for relapse in one to 5 years after resection. In this study, we investigated the genetic differences in genomic DNA of primary and rapidly recurrent tumours in the same patient to explain the causality mechanisms of PA recurrence. CASE PRESENTATION: The patient was a 69-year-old female with non-functional pituitary macroadenoma with extension into the left cavernous sinus (Knosp grade 2) who underwent craniotomy and partial resection in August 2010. Two years later, the patient had prolonged tumour growth with an essential suprasellar extension (Knosp grade 2), and a second craniotomy with partial tumour resection was performed in September 2012. In both tumours, the KI-67 level was below 1.5%. Exome sequencing via semiconductor sequencing of patient germline DNA and somatic DNA from both tumours was performed. Tmap alignment and Platypus variant calling were performed followed by variant filtering and manual review with IGV software. We observed an increased load of missense variants in the recurrent PA tumour when compared to the original tumour. The number of detected variants increased from ten to 26 and potential clonal expansion of four variants was observed. Additionally, targeted SNP analysis revealed five rare missense SNPs with a potential impact on the function of the encoded proteins. CONCLUSIONS: In this case study, an SNP located in HRAS is the most likely candidate inducing rapid PA progression. The relapsed PA tumour had a higher variation load and fast tumour recurrence in this patient could be caused by clonal expansion of the leftover tumour tissue.

Medication for Acromegaly Reduces Expression of MUC16, MACC1 and GRHL2 in Pituitary Neuroendocrine Tumour Tissue.

Acromegaly is a disease mainly caused by pituitary neuroendocrine tumor (PitNET) overproducing growth hormone. First-line medication for this condition is the use of somatostatin analogs (SSAs), that decrease tumor mass and induce antiproliferative effects on PitNET cells. Dopamine agonists (DAs) can also be used if SSA treatment is not effective. This study aimed to determine differences in transcriptome signatures induced by SSA/DA therapy in PitNET tissue. We selected tumor tissue from twelve patients with somatotropinomas, with half of the patients receiving SSA/DA treatment before surgery and the other half treatment naive. Transcriptome sequencing was then carried out to identify differentially expressed genes (DEGs) and their protein-protein interactions, using pathway analyses. We found 34 upregulated and six downregulated DEGs in patients with SSA/DA treatment. Three tumor development promoting factors MUC16, MACC1, and GRHL2, were significantly downregulated in therapy administered PitNET tissue; this finding was supported by functional studies in GH3 cells. Protein-protein interactions and pathway analyses revealed extracellular matrix involvement in the antiproliferative effects of this type of the drug treatment, with pronounced alterations in collagen regulation. Here, we have demonstrated that somatotropinomas can be distinguished based on their transcriptional profiles following SSA/DA therapy, and SSA/DA treatment does indeed cause changes in gene expression. Treatment with SSA/DA significantly downregulated several factors involved in tumorigenesis, including MUC16, MACC1, and GRHL2. Genes that were upregulated, however, did not have a direct influence on antiproliferative function in the PitNET cells. These findings suggested that SSA/DA treatment acted in a tumor suppressive manner and furthermore, collagen related interactions and pathways were enriched, implicating extracellular matrix involvement in this anti-tumor effect of drug treatment.

Evaluation of the Possibility to Detect Circulating Tumor DNA From Pituitary Adenoma.

Objective: Circulating free DNA (cfDNA) in general and circulating tumor DNA (ctDNA) in particular is becoming an increasingly used form of liquid biopsy biomarkers. In this study, we are investigating the ability to detect ctDNA from the plasma of pituitary adenoma (PA) patients. Design: Tumor tissue samples were obtained from planed PA resections, before which blood plasma samples were taken. Somatic variants found in PA tissue samples were evaluated in related cfDNA, isolated from plasma samples. Methods: Sanger sequencing, as well as previously obtained whole-exome sequencing data, were used to evaluate somatic variants composition in tumor tissue samples. cfDNA was isolated from the same PA patients and competitive allele-specific TaqMan PCR and amplicon-based next-generation sequencing (NGS) approach were used for targeted detection of variants found in corresponding tumor tissue samples. Results: Using NGS-based analysis, we detected five out of 17 somatic variants in 40 to 60% of total reads, three variants in 0.50-5.00% of total read count, including GNAS c.601C>T, which was detected using ultra-deep NGS (1.78 million X) in 0.77% of amplicons reads. Nine variants were not detected. We also detected We were not able to detect variant found in PA tissue in cfDNA using cast-PCR, indicating that the portion of variant-containing ctDNA in total isolated cfDNA is too small to be detected with this method. Conclusions: For the first time, we demonstrate the possibility to detect somatic variants of PA in cfDNA isolated from patients' blood plasma. Whether the source of variant detected in cfDNA is PA should be further tested.

Significantly altered peripheral blood cell DNA methylation profile as a result of immediate effect of metformin use in healthy individuals.

BACKGROUND: Metformin is a widely prescribed antihyperglycemic agent that has been also associated with multiple therapeutic effects in various diseases, including several types of malignancies. There is growing evidence regarding the contribution of the epigenetic mechanisms in reaching metformin's therapeutic goals; however, the effect of metformin on human cells in vivo is not comprehensively studied. The aim of our study was to examine metformin-induced alterations of DNA methylation profiles in white blood cells of healthy volunteers, employing a longitudinal study design. RESULTS: Twelve healthy metformin-naïve individuals where enrolled in the study. Genome-wide DNA methylation pattern was estimated at baseline, 10 h and 7 days after the start of metformin administration. The whole-genome DNA methylation analysis in total revealed 125 differentially methylated CpGs, of which 11 CpGs and their associated genes with the most consistent changes in the DNA methylation profile were selected: POFUT2, CAMKK1, EML3, KIAA1614, UPF1, MUC4, LOC727982, SIX3, ADAM8, SNORD12B, VPS8, and several differentially methylated regions as novel potential epigenetic targets of metformin. The main functions of the majority of top-ranked differentially methylated loci and their representative cell signaling pathways were linked to the well-known metformin therapy targets: regulatory processes of energy homeostasis, inflammatory responses, tumorigenesis, and neurodegenerative diseases. CONCLUSIONS: Here we demonstrate for the first time the immediate effect of short-term metformin administration at therapeutic doses on epigenetic regulation in human white blood cells. These findings suggest the DNA methylation process as one of the mechanisms involved in the action of metformin, thereby revealing novel targets and directions of the molecular mechanisms underlying the various beneficial effects of metformin. TRIAL REGISTRATION: EU Clinical Trials Register, 2016-001092-74. Registered 23 March 2017, https://www.clinicaltrialsregister.eu/ctr-search/trial/2016-001092-74/LV .

Single nucleotide polymorphisms in the intergenic region between metformin transporter OCT2 and OCT3 coding genes are associated with short-term response to metformin monotherapy in type 2 diabetes mellitus patients.

OBJECTIVES: High variability in clinical response to metformin is often observed in type 2 diabetes (T2D) patients, and it highlights the need for identification of genetic components affecting the efficiency of metformin therapy. Aim of this observational study is to evaluate the role of tagSNPs (tagging single nucleotide polymorphisms) from genomic regions coding for six metformin transporter genes with respect to the short-term efficiency. DESIGN: 102 tagSNPs in 6 genes coding for metformin transporters were genotyped in the group of 102 T2D patients treated with metformin for 3 months. METHODS: Most significant hits were analyzed in the group of 131 T2D patients from Slovakia. Pharmacokinetic study in 25 healthy nondiabetic volunteers was conducted to investigate the effects of identified polymorphisms. RESULTS: In the discovery group of 102 patients, minor alleles of rs3119309, rs7757336 and rs2481030 were significantly nominally associated with metformin inefficiency (P = 1.9 × 10-6 to 8.1 × 10-6). Effects of rs2481030 and rs7757336 did not replicate in the group of 131 T2DM patients from Slovakia alone, whereas rs7757336 was significantly associated with a reduced metformin response in combined group. In pharmacokinetic study, group of individuals harboring risk alleles of rs7757336 and rs2481030 displayed significantly reduced AUC∞ of metformin in plasma. CONCLUSIONS: For the first time, we have identified an association between the lack of metformin response and SNPs rs3119309 and rs7757336 located in the 5' flanking region of the genes coding for Organic cation transporter 2 and rs2481030 located in the 5' flanking region of Organic cation transporter 3 that was supported by the results of a pharmacokinetic study on 25 healthy volunteers.