Next Article in Journal
Exploring Self-Care Management Practices among Patients Diagnosed with Type 2 Diabetes Mellitus at District General Hospital in Chilaw, Sri Lanka
Previous Article in Journal
Foliar Application with Plant-Derived Extracts Enhances Growth, Physiological Parameters, and Yield of Potatoes (Solanum tuberosum L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

Associations between Vitamin D Receptor (VDR) Polymorphisms and Gut Microbiota in a Spanish Population Sample †

1
Department of Biosciences, European University of Madrid, 28670 Madrid, Spain
2
Department of Biomedicine, Universidad Francisco de Vitoria, 28223 Madrid, Spain
3
Pre-Clinical Dentistry Department, European University of Madrid, 28670 Madrid, Spain
4
Nutrition and Food Sciences Department, Faculty of Pharmacy, University Complutense of Madrid, 28040 Madrid, Spain
5
Preventive Medicine and Public Health and Microbiology Center, Faculty of Medicine, University Autónoma of Madrid, 28049 Madrid, Spain
*
Author to whom correspondence should be addressed.
Presented at the 3rd International Electronic Conference on Nutrients, 1–15 November 2023; Available online: https://iecn2023.sciforum.net/.
Biol. Life Sci. Forum 2023, 29(1), 4; https://doi.org/10.3390/IECN2023-15530
Published: 31 October 2023
(This article belongs to the Proceedings of The 3rd International Electronic Conference on Nutrients)

Abstract

:
This study examined the association between genetic variations in the Vitamin D Receptor (VDR) gene and the composition of gut microbiota in a Spanish population. The VDR gene plays a role in mediating the effects of vitamin D on gut microbiota composition. We analyzed 87 healthy participants from the Spanish Caucasian population and found that the genetic variations were in Hardy–Weinberg equilibrium. Through metataxonomic sequencing, we identified that the GG genotype from the VDR polymorphism rs731236 seems to be associated with the opportunistic genus Solobacterium, indicating that genetic variations may influence the composition of gut microbiota.

1. Introduction

1,25-Dihydroxyvitamin D (1,25(OH)2D3), the active form of vitamin D, activates the widely expressed vitamin D receptor (VDR). The VDR gene encodes the vitamin D receptor, which plays a crucial role in mediating the biological effects of vitamin D. The vitamin D3 receptor belongs to a group of ligand-inducible transcription factors known as the nuclear hormone receptor superfamily [1]. The receptor not only functions as a receptor for vitamin D3 but also acts as a receptor for a secondary bile acid called lithocholic acid [2]. The main targets regulated by the vitamin D3 receptor are primarily involved in mineral metabolism, but it also plays a role in various other metabolic pathways such as immune response and cancer [3]. The vitamin D receptor (VDR) seems a key genetic factor for shaping the host microbiome as conditional VDR deletion severely changes metabolites specifically produced from carbohydrate, protein, lipid, and bile acid metabolism. In one study conducted in mice, gender differences were also observed regarding the effect of VDR deletion [4].
Nevertheless, low levels of vitamin D or inactivating polymorphisms in VDR have been associated with inflammatory and metabolic disorders. Mutations in this gene have been associated with a condition called type II vitamin D-resistant rickets [5]. Different transcript variants of the gene have been identified. A recent study has also provided evidence of translational readthrough in this gene, resulting in the expression of an extended isoform with a different C-terminus, utilizing an alternative translation termination codon within the gene’s coding sequence [6].
The rs731236 polymorphism, also known as the TaqI restriction site polymorphism (T/C), is one of the common genetic variations in the VDR gene. There are different studies suggesting a role of this polymorphism in health status. A nonsignificant increase in mortality of 55% was observed among patients homozygous for the rare allele C [7]. On the contrary, an analysis conducted among 111 Swedish breast cancer survivors below 37 years of age found a higher survival among estrogen receptor-positive tamoxifen-treated patients homozygous for the C allele [8]. VDR has also been suggested to be a potential key player in the pathogenetic mechanism of obesity. Vasilopoulos et al. found we found a strong association between the VDR TaqI (rs731236) ‘T’ allele and obesity when analyzing a Greek population [9].
A metanalysis from Wang et al. suggested a role of VDR in the composition of the gut microbiome and showed that Parabacteroides are the most significant taxa correlated with the VDR gene [10]. The lack of VDR leads to dysbiosis, indicating a critical role of VDR in the conformation of gut microbiota communities [11].
In this work, we have explored the connections between the gut microbiome and the VDR gene variations in a sample of the Caucasian Spanish population. By analyzing genetic variations in the VDR gene and their potential impact on the gut microbial composition, we aim to shed light on the interplay between these factors.

2. Materials and Methods

2.1. Study Population

This study included a total of 87 healthy participants, consisting of 57 men and 30 women, with ages ranging from 18 to 48 years. Strict exclusion criteria were applied, excluding individuals with any form of pathology within six months prior to the study, those who had undergone previous gastrointestinal surgery, those who had taken antibiotics within three months prior to the study, smokers, individuals using prebiotics, probiotics, or nutritional supplements, vegetarians, or vegans, and pregnant or lactating individuals. All participants in this study were of Caucasian ethnicity.

2.2. Sample Collection

Participants were provided with the Fe-Col® Fecal Sample Collection Kit (Alpha Laboratories, Hampshire, UK), along with insulated bags and ice blocks to maintain sample integrity during transportation. Stool samples were preserved at −80 °C until further processing.

2.3. DNA Extraction

Human and bacterial DNA extraction was conducted on 100 mg of stool sample utilizing the commercial E.Z.N.A.® Stool DNA Kit (Omega Biotek, Norcross, GA, USA) and a bead-beating homogenizer (Bullet Blender Storm, Next Advance, New York, NY, USA). DNA concentration and purity were assessed using the Quant-iT PicoGreen dsDNA Assay Kit (ThermoFisher Scientific, Waltham, MA, USA) and an FP-8300 spectrofluorimeter (Jasco, Tokyo, Japan). Bacterial DNA was employed for microbiota analysis, while human DNA was used for VDR genotyping.

2.4. VDR Genotyping

Applied Biosystems TaqMan® SNP Genotyping Assays (Applied Biosystems, Foster City, CA, USA; assay ID: C_2404008_10) were employed for allelic discrimination analysis of the TaqI VDR (rs731236) gene polymorphism. The StepOnePlus Real-Time PCR system (ThermoFisher Scientific) was utilized, with a protocol involving denaturation at 95 °C for 10 min, followed by 50 cycles of denaturation at 92 °C for 15 s, annealing/extension at 60 °C for 1 min, and a final extension step of 30 s at 60 °C. Fluorescence analysis was conducted with Allelic Discrimination 7500 software v.2.0.2. Following genotyping, participants were classified as follows based on the VDR genotype: TT as VDR-1, TC as VDR-2, and CC as VDR-3 for subsequent analysis.

2.5. Sequencing and Bioinformatics

The hypervariable V3 and V4 regions were amplified using the primer pair 5′-TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG-3′ and 5′ GTCTCGTGGGCTCGG AGATGTGTATAAGAGACAG-3′. The amplicon of 459 bp was visualized in a 0.8% agarose gel stained with ethidium bromide, and bands were cut and cleaned using the MinElute Gel Extraction Kit (Qiagen, Hilden, Germany). DNA amplicons were sequenced on a MiSeq Illumina platform (Illumina, San Diego, CA, USA). Sequence outputs were analyzed using the Quantitative Insights into Microbial Ecology (QIIME2) program, v2019.10. The 16s paired reads were imported into QIIME2 and processed with the DADA2 plugin], adjusting the maximum expected error threshold to 2.0 (both forward and reverse). The taxonomy assignments were performed with the classify-sklearn method and an in-house customized classifier based on the SILVA reference database [12,13,14]. To build the customized reference database, sequences, according to our primers (forward primer sequence: CCTACGGGNGGCWGCAG, reverse primer sequence: GACTACHVGGGTA TCTAATCC), were extracted from the SILVA 132 database clustered at 99% identity. The classifier was trained using our tailored reference reads and SILVA 7-levels for reference taxonomy, including the species probability (weights) likely to be observed for human stool (downloaded from https://github.com/BenKaehler/readytowear, accessed on 1 June 2023) [12,13].

2.6. Statistical Analysis

The allele and genotypes frequencies were calculated using the SNPStat program [15]. Statistical analysis was carried out using QIIME2 v2019.10, SPSS software v26.0 (SPSS, Chicago, IL, USA) and the R statistical package v4.1.1. Linear discriminant analysis coupled with effect size (LEfSe v1.0) was performed to identify bacterial associated pathways differentially represented between groups with default settings. Significance was set at p < 0.05.

3. Results

3.1. VDR Genotyping

Descriptive statistics were used to show the baseline characteristics of the participants in this study. Deviation from the Hardy–Weinberg equilibrium (HWE) was tested using SNPStat Software [15]. Genotyping results for the rs731236 (TaqI) polymorphism were obtained and are presented in Table 1. In comparison to the 1000 Genomes Project reference, the calculated allele frequencies (T = 0.61 and C = 0.39) were as expected for a European population sample. When examining genotypes, TT homozygotes were found at a frequency of 0.43, while TC heterozygotes and CT heterozygotes had frequencies of 0.38 and 0.20, respectively, in accordance with Hardy–Weinberg equilibrium (p = 0.071) (Table 1).

3.2. Bacterial Composition According to Genotypes

No differences were found in either body composition parameters or diet.
LEfSe (Linear Discriminant Analysis Effect Size) statistical analysis shows that the GG genotype is enriched in the genus Solobacterium, whilst the AA genotype is enriched in the genera Victivallaceae, Moryella, Enterococcus (Fam. Enterococcaceae) and the families Desulfovibrionaceae, Selenomonadaceae, and Ruminococcaceae (phylum Desulfobacterota; class Desulfovibrionia, order Desulfovibrionales) (Figure 1).

4. Discussion

The LEfSe analysis revealed a strong association between the Solobacterium genus and the presence of the rare GG genotype for the VDR gene. Some species of the Solobacterium genus are opportunistic pathogens and are part of the oral and intestinal microbiotas [16]. Solobacterium seems to be increased in obesity in some cases, and in a recent study, the intake of lactic acid appeared to decrease its presence [17]. However, in our study, we did not find associations between body composition and the genotypes, so further studies and analysis will be required to find consistent results showing a role of the VDR nuclear receptor in the gut microbiota composition and whether they lead to specific changes in the microbiota, such as the presence of opportunistic pathogens as it seems.

Author Contributions

Conceptualization, R.G.S.; methodology, C.B. and M.B.; software, R.G.S., C.B. and M.B.; validation, M.B. and R.G.S., experimental work, M.T. and R.G.S.; data curation, C.B. and M.B.; writing—original draft preparation, R.G.S.; writing—review and editing, R.G.S., M.L., S.C. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universidad Europea de Madrid, grant number 2018/UEM04.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of COMITÉ ÉTICO DE INVESTIGACIÓN CON MEDICAMENTOS REGIONAL DE LA COMUNIDAD DE MADRID (protocol number 47/560280.9/18 approved 10 September 2018).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Data Availability Statement

Data will be provided upon a collaborative research request.

Acknowledgments

The authors acknowledge Goyi Montalvo for the technical help and implication during the experimental work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pike, J.W.; Meyer, M.B. The vitamin D receptor: New paradigms for the regulation of gene expression by 1,25-dihydroxyvitamin D3. Endocrinol. Metab. Clin. N. Am. 2010, 39, 255–269. [Google Scholar] [CrossRef] [PubMed]
  2. Ishizawa, M.; Akagi, D.; Makishima, M. Lithocholic Acid Is a Vitamin D Receptor Ligand That Acts Preferentially in the Ileum. Int. J. Mol. Sci. 2018, 19, 1975. [Google Scholar] [CrossRef] [PubMed]
  3. Khammissa, R.A.G.; Fourie, J.; Motswaledi, M.H.; Ballyram, R.; Lemmer, J.; Feller, L. The Biological Activities of Vitamin D and Its Receptor in Relation to Calcium and Bone Homeostasis, Cancer, Immune and Cardiovascular Systems, Skin Biology, and Oral Health. BioMed Res. Int. 2018, 2018, 9276380. [Google Scholar] [CrossRef] [PubMed]
  4. Feldman, D.; Malloy, P.J. Mutations in the vitamin D receptor and hereditary vitamin D-resistant rickets. Bonekey Rep. 2014, 3, 510. [Google Scholar] [CrossRef] [PubMed]
  5. Chatterjee, I.; Lu, R.; Zhang, Y.; Zhang, J.; Dai, Y.; Xia, Y.; Sun, J. Vitamin D receptor promotes healthy microbial metabolites and microbiome. Sci. Rep. 2020, 10, 7340. [Google Scholar] [CrossRef] [PubMed]
  6. Loughran, G.; Jungreis, I.; Tzani, I.; Power, M.; Dmitriev, R.I.; Ivanov, I.P.; Kellis, M.; Atkins, J.F. Stop codon readthrough generates a C-terminally extended variant of the human vitamin D receptor with reduced calcitriol response. J. Biol. Chem. 2018, 293, 4434–4444. [Google Scholar] [CrossRef] [PubMed]
  7. Goode, E.L.; Dunning, A.M.; Kuschel, B.; Healey, C.S.; E Day, N.; Ponder, B.A.J.; Easton, D.F.; Pharoah, P.P.D. Effect of germ-line genetic variation on breast cancer survival in a population-based study. Cancer Res. 2002, 62, 3052–3057. [Google Scholar] [PubMed]
  8. Lundin, A.C.; Söderkvist, P.; Eriksson, B.; Bergman-Jungeström, M.; Wingren, S.; South-East Sweden Breast Cancer Group. Association of breast cancer progression with a vitamin D receptor gene polymorphism. Cancer Res. 1999, 59, 2332–2334. [Google Scholar] [PubMed]
  9. Vasilopoulos, Y.; Sarafidou, T.; Kotsa, K.; Papadimitriou, M.; Goutzelas, Y.; Stamatis, C.; Bagiatis, V.; Tsekmekidou, X.; Yovos, J.G.; Mamuris, Z. VDR TaqI is associated with obesity in the Greek population. Gene 2013, 512, 237–239. [Google Scholar] [CrossRef] [PubMed]
  10. Wang, J.; Thingholm, L.B.; Skiecevičienė, J.; Rausch, P.; Kummen, M.; Hov, J.R.; Degenhardt, F.; Heinsen, F.-A.; Rühlemann, M.C.; Szymczak, S.; et al. Genome-wide association analysis identifies variation in vitamin D receptor and other host factors influencing the gut microbiota. Nat. Genet. 2016, 48, 1396–1406. [Google Scholar] [CrossRef] [PubMed]
  11. Wu, S.; Zhang, Y.-G.; Lu, R.; Xia, Y.; Zhou, D.; O Petrof, E.; Claud, E.C.; Chen, D.; Chang, E.B.; Carmeliet, G.; et al. Intestinal epithelial vitamin D receptor deletion leads to defective autophagy in colitis. Gut 2015, 64, 1082–1094. [Google Scholar] [CrossRef] [PubMed]
  12. Kaehler, B.D.; Bokulich, N.A.; McDonald, D.; Knight, R.; Caporaso, J.G.; Huttley, G.A. Species-level microbial sequence classification is improved by source-environment information. Nat. Commun. 2019, 10, 4643. [Google Scholar] [CrossRef] [PubMed]
  13. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Caporaso, J.G. Optimizing taxonomic classification of marker gene sequences. Microbiome 2018, 6, 90. [Google Scholar] [CrossRef] [PubMed]
  14. Bailén, M.; Tabone, M.; Bressa, C.; Lominchar, M.G.M.; Larrosa, M.; González-Soltero, R. Unraveling Gut Microbiota Signatures Associated with PPARD and PARGC1A Genetic Polymorphisms in a Healthy Population. Genes 2022, 13, 289. [Google Scholar] [CrossRef] [PubMed]
  15. Sole, X.; Guino, E.; Valls, J.; Iniesta, R.; Moreno, V. SNPStats: A web tool for the analysis of association studies. Bioinformatics 2006, 22, 1928–1929. [Google Scholar] [CrossRef] [PubMed]
  16. Alauzet, C.; Aujoulat, F.; Lozniewski, A.; Brahim, S.B.; Domenjod, C.; Enault, C.; Lavigne, J.-P.; Marchandin, H. A new look at the genus solobacterium: A retrospective analysis of twenty-seven cases of infection involving s. moorei and a review of sequence databases and the literature. Microorganisms 2021, 9, 1229. [Google Scholar] [CrossRef] [PubMed]
  17. Wu, Y.; Chi, X.; Zhang, Q.; Chen, F.; Deng, X. Characterization of the salivary microbiome in people with obesity. PeerJ. 2018, 6, e4458. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Lefse analysis showing bacterial composition according to genotype.
Figure 1. Lefse analysis showing bacterial composition according to genotype.
Blsf 29 00004 g001
Table 1. Genotyping results showing the allele and genotype frequencies for the VDR rs731236 (TaqI) polymorphism (test for Hardy–Weinberg equilibrium; p = 0.071).
Table 1. Genotyping results showing the allele and genotype frequencies for the VDR rs731236 (TaqI) polymorphism (test for Hardy–Weinberg equilibrium; p = 0.071).
rs731236 (TaqI) (C/T)SNP Allele Frequencies (n = 87)GenotypesFrequencies
(H-W)
T0.39TT0.43
C0.61TC0.38
CC0.20
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

González Soltero, R.; Bressa, C.; Tabone, M.; Clemente, S.; Larrosa, M.; Bailén, M. Associations between Vitamin D Receptor (VDR) Polymorphisms and Gut Microbiota in a Spanish Population Sample. Biol. Life Sci. Forum 2023, 29, 4. https://doi.org/10.3390/IECN2023-15530

AMA Style

González Soltero R, Bressa C, Tabone M, Clemente S, Larrosa M, Bailén M. Associations between Vitamin D Receptor (VDR) Polymorphisms and Gut Microbiota in a Spanish Population Sample. Biology and Life Sciences Forum. 2023; 29(1):4. https://doi.org/10.3390/IECN2023-15530

Chicago/Turabian Style

González Soltero, Rocío, Carlo Bressa, Mariangela Tabone, Sara Clemente, Mar Larrosa, and María Bailén. 2023. "Associations between Vitamin D Receptor (VDR) Polymorphisms and Gut Microbiota in a Spanish Population Sample" Biology and Life Sciences Forum 29, no. 1: 4. https://doi.org/10.3390/IECN2023-15530

Article Metrics

Back to TopTop