Institute of Genetics and Cytology of the National Academy of Sciences of Belarus

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

Micronucleus (MN) expression in peripheral blood lymphocytes is well established as a standard method for monitoring chromosome damage in human populations. The first results of an analysis of pooled data from laboratories using the cytokinesis-block micronucleus (CBMN) assay and participating in the HUMN (HUman MicroNucleus project) international collaborative study are presented. The effects of laboratory protocol, scoring criteria, and host factors on baseline micronucleated binucleate cell (MNC) frequency are evaluated, and a reference range of "normal" values against which future studies may be compared is provided. Primary data from historical records were submitted by 25 laboratories distributed in 16 countries. This resulted in a database of nearly 7000 subjects. Potentially significant differences were present in the methods used by participating laboratories, such as in the type of culture medium, the concentration of cytochalasin-B, the percentage of fetal calf serum, and in the culture method. Differences in criteria for scoring micronuclei were also evident. The overall median MNC frequency in nonexposed (i.e., normal) subjects was 6.5 per thousand and the interquartile range was between 3 and 12 per thousand. An increase in MNC frequency with age was evident in all but two laboratories. The effect of gender, although not so evident in all databases, was also present, with females having a 19% higher level of MNC frequency (95% confidence interval: 14-24%). Statistical analyses were performed using random-effects models for correlated data. Our best model, which included exposure to genotoxic factors, host factors, methods, and scoring criteria, explained 75% of the total variance, with the largest contribution attributable to laboratory methods.

The Slavic branch of the Balto-Slavic sub-family of Indo-European languages underwent rapid divergence as a result of the spatial expansion of its speakers from Central-East Europe, in early medieval times. This expansion-mainly to East Europe and the northern Balkans-resulted in the incorporation of genetic components from numerous autochthonous populations into the Slavic gene pools. Here, we characterize genetic variation in all extant ethnic groups speaking Balto-Slavic languages by analyzing mitochondrial DNA (n = 6,876), Y-chromosomes (n = 6,079) and genome-wide SNP profiles (n = 296), within the context of other European populations. We also reassess the phylogeny of Slavic languages within the Balto-Slavic branch of Indo-European. We find that genetic distances among Balto-Slavic populations, based on autosomal and Y-chromosomal loci, show a high correlation (0.9) both with each other and with geography, but a slightly lower correlation (0.7) with mitochondrial DNA and linguistic affiliation. The data suggest that genetic diversity of the present-day Slavs was predominantly shaped in situ, and we detect two different substrata: 'central-east European' for West and East Slavs, and 'south-east European' for South Slavs. A pattern of distribution of segments identical by descent between groups of East-West and South Slavs suggests shared ancestry or a modest gene flow between those two groups, which might derive from the historic spread of Slavic people.

Phytosterols are natural sterols widely found in plants that have a variety of physiological functions, and their role in reducing cholesterol absorption has garnered much attention. Although the bioavailability of phytosterols is only 0.5-2%, they can still promote cholesterol balance in the body. A mechanism of phytosterols for lowering cholesterol has now been proposed. They not only reduce the uptake of cholesterol in the intestinal lumen and affect its transport, but also regulate the metabolism of cholesterol in the liver. In addition, phytosterols can significantly reduce the plasma concentration of total cholesterol, triglycerides, and low-density lipoprotein cholesterol (LDL-C), with a dose-response relationship. Ingestion of 3 g of phytosterols per day can reach the platform period, and this dose can reduce LDL-C by about 10.7%. On the other hand, phytosterols can also activate the liver X receptor α-CPY7A1 mediated bile acids excretion pathway and accelerate the transformation and metabolism of cholesterol. This article reviews the research progress of phytosterols as a molecular regulator of cholesterol and the mechanism of action for this pharmacological effect.

BACKGROUND: The genetic origins of Uralic speakers from across a vast territory in the temperate zone of North Eurasia have remained elusive. Previous studies have shown contrasting proportions of Eastern and Western Eurasian ancestry in their mitochondrial and Y chromosomal gene pools. While the maternal lineages reflect by and large the geographic background of a given Uralic-speaking population, the frequency of Y chromosomes of Eastern Eurasian origin is distinctively high among European Uralic speakers. The autosomal variation of Uralic speakers, however, has not yet been studied comprehensively. RESULTS: Here, we present a genome-wide analysis of 15 Uralic-speaking populations which cover all main groups of the linguistic family. We show that contemporary Uralic speakers are genetically very similar to their local geographical neighbours. However, when studying relationships among geographically distant populations, we find that most of the Uralic speakers and some of their neighbours share a genetic component of possibly Siberian origin. Additionally, we show that most Uralic speakers share significantly more genomic segments identity-by-descent with each other than with geographically equidistant speakers of other languages. We find that correlated genome-wide genetic and lexical distances among Uralic speakers suggest co-dispersion of genes and languages. Yet, we do not find long-range genetic ties between Estonians and Hungarians with their linguistic sisters that would distinguish them from their non-Uralic-speaking neighbours. CONCLUSIONS: We show that most Uralic speakers share a distinct ancestry component of likely Siberian origin, which suggests that the spread of Uralic languages involved at least some demic component.

Abstract In order lo investigate the origin of two of the German 1RS. 1BL wheat‐rye translocations used world‐wide in breeding, a number of DNA probes were considered which (a) were critical for the short arm of the rye chromosome 1 R and (b) should show a specificity for the gene pool of Petkus rye. The DNA probe CDO580 was revealed as a specific one. (1) It clearly differentiated 1RS.1AL (‘Amigo’). 1RS.lBL (‘Salmon’) and 1RS.1DL (‘Gabo’) from the two German sources. (2) Both translocation wheats deriving from the Weihenstephan (Munich) and from the Salzmünde (Halle/S.) origin showed an identical DNA fragment which was typical for the gene pool of Petkus rye. It is supposed that both German sources have one progenitor in common.

The current method for scoring micronuclei as a measure of genetic damage rate in peripheral blood cells is to enumerate this end point in cytokinesis-blocked binucleated cultured lymphocytes. However, one can expect that, due to chronic exposure to genotoxins or inherent genetic instability, micronuclei may be expressed continually in vivo in dividing cell populations such as the progenitor cell lineages leading to mature lymphocytes or erythrocytes. Consequently, micronuclei may already be expressed in peripheral blood lymphocytes prior to culture. In view of these considerations, we have performed a study in children living in regions of Belarus that are contaminated by radionuclides from the Chernobyl disaster and compared their micronucleus frequency in erythrocytes, nondivided lymphocytes, and cultured cytokinesis-blocked binucleated lymphocytes to that of controls living in noncontaminated areas. Preliminary data presented in this paper indicate a significant two- to fourfold increase in micronucleus expression (P < 0.05) in exposed children relative to controls in erythrocytes or peripheral blood lymphocytes in blood smears as well as in mononuclear and cytokinesis-blocked binucleated lymphocytes in cultures. The measurement of micronuclei in nondivided mononuclear lymphocytes represents chromosomal damage expressed during in vivo divisions. The micronuclei in binucleated cultured cells represent micronuclei expressed ex-vivo and may include micronuclei already present in a cell prior to tissue culture. These preliminary data suggest that a different spectrum and level of damage may be observed in nondivided mononuclear lymphocytes, binucleated lymphocytes, and erythrocytes and that a combination of these approaches may provide a more comprehensive assessment of the extent of genetic damage induced by chronic exposure to radionuclides or other genotoxins in haematopoietic tissue.

OBJECTIVES: To investigate the relationship between body mass index (BMI) and 9-year mortality in older (> or = 65) Americans with and without disability. DESIGN: Cohort study. SETTING: The unique disability-focused National Long Term Care Survey (NLTCS) data that assessed the health and well-being of older individuals in 1994 were analyzed. PARTICIPANTS: Four thousand seven hundred ninety-one individuals in the 1994 survey. MEASUREMENTS: BMI (kg/m2) was calculated from self- or proxy reports of height and weight. The analysis was adjusted for 1-year change in BMI and demographic and health-related factors, as well as reports by proxies, and death occurring during the first 2 years after the interview. RESULTS: The relative risk of death as a function of BMI formed a nonsymmetric U-shaped pattern, with larger risks associated with lower BMI (< 22.0) and minimal risks for BMI of 25.0 to 34.9. (BMI 22.0-24.9 was the reference.) Adjustments for demographic and health-related factors had little effect on this pattern. Nondisabled individuals exhibited a similar U-shaped pattern but with lower risks associated with lower BMI. For disabled individuals, the mortality-risk pattern was higher for lower BMI (< 22.0) and flat for higher BMI, thus exhibiting an inverse J shape. BMI patterns were age sensitive, with disability status affecting sensitivity. CONCLUSION: Overweight or mild (grade 1) obesity was not a risk factor for 9-year mortality in older Americans participating in the 1994 NLTCS. A flatter BMI pattern of the relative risk of death for disabled than for nondisabled individuals suggests that optimal body weight can be sensitive to age and health and well-being.

BACKGROUND: (CYP11A1) encoded by the CYP11A1 gene catalyzes the conversion of cholesterol into pregnenolone - the general precursor of all the steroid hormones, starting with progesterone. This stage is missing in plants where mitochondrial cytochromes P450 (the mito CYP clan) have not been found. Generating transgenic plants with a mitochondrial type P450 from animals would offer an interesting option to verify whether plant mitochondria could serve as another site of P450 monooxygenase reaction for the steroid hormones biosynthesis. RESULTS: For a more detailed comparison of steroidogenic systems of Plantae and Animalia, we have created and studied transgenic tobacco and tomato plants efficiently expressing mammalian CYP11A1 cDNA. The detailed phenotypic characterization of plants obtained has shown that through four generations studied, the transgenic tobacco plants have reduced a period of vegetative development (early flowering and maturation of bolls), enlarged biomass and increased productivity (quantity and quality of seeds) as compared to the only empty-vector containing or wild type plants. Moreover, the CYP11A1 transgenic plants show resistance to such fungal pathogen as Botrytis cinerea. Similar valuable phenotypes (the accelerated course of ontogenesis and/or stress resistance) are also visible in two clearly distinct transgenic tomato lines expressing CYP11A1 cDNA: one line (No. 4) has an accelerated rate of vegetative development, while the other (No. 7) has enhanced immunity to abiotic and biotic stresses. The progesterone level in transgenic tobacco and tomato leaves is 3-5 times higher than in the control plants of the wild type. CONCLUSIONS: For the first time, we could show the compatibility in vivo of even the most specific components of the systems of biosynthesis of steroid hormones in Plantae and Animalia. The hypothesis is proposed and substantiated that the formation of the above-noted special phenotypes of transgenic plants expressing mammalian CYP11A1 cDNA is due to the increased biosynthesis of progesterone that can be considered as a very ancient bioregulator of plant cells and the first real hormone common to plants and animals.

Twenty different Pseudomonas strains utilizing m-toluate were isolated from oil-contaminated soil samples near Minsk, Belarus. Seventeen of these isolates carried plasmids ranging in size from 78 to about 200 kb (assigned pSVS plasmids) and encoding the meta cleavage pathway for toluene metabolism. Most plasmids were conjugative but of unknown incompatibility groups, except for one, which belonged to the IncP9 group. The organization of the genes for toluene catabolism was determined by restriction analysis and hybridization with xyl gene probes of pWW0. The majority of the plasmids carried xyl-type genes highly homologous to those of pWW53 and organized in a similar manner (M. T. Gallegos, P. A. Williams, and J. L. Ramos, J. Bacteriol. 179:5024-5029, 1997), with two distinguishable meta pathway operons, one upper pathway operon, and three xylS-homologous regions. All of these plasmids also possessed large areas of homologous DNA outside the catabolic genes, suggesting a common ancestry. Two other pSVS plasmids carried only one meta pathway operon, one upper pathway operon, and one copy each of xylS and xylR. The backbones of these two plasmids differed greatly from those of the others. Whereas these parts of the plasmids, carrying the xyl genes, were mostly conserved between plasmids of each group, the noncatabolic parts had undergone intensive DNA rearrangements. DNA sequencing of specific regions near and within the xylTE and xylA genes of the pSVS plasmids confirmed the strong homologies to the xyl genes of pWW53 and pWW0. However, several recombinations were discovered within the upper pathway operons of the pSVS plasmids and pWW0. The main genetic mechanisms which are thought to have resulted in the present-day configuration of the xyl operons are discussed in light of the diversity analysis carried out on the pSVS plasmids.

Ethnic Belarusians make up more than 80% of the nine and half million people inhabiting the Republic of Belarus. Belarusians together with Ukrainians and Russians represent the East Slavic linguistic group, largest both in numbers and territory, inhabiting East Europe alongside Baltic-, Finno-Permic- and Turkic-speaking people. Till date, only a limited number of low resolution genetic studies have been performed on this population. Therefore, with the phylogeographic analysis of 565 Y-chromosomes and 267 mitochondrial DNAs from six well covered geographic sub-regions of Belarus we strove to complement the existing genetic profile of eastern Europeans. Our results reveal that around 80% of the paternal Belarusian gene pool is composed of R1a, I2a and N1c Y-chromosome haplogroups - a profile which is very similar to the two other eastern European populations - Ukrainians and Russians. The maternal Belarusian gene pool encompasses a full range of West Eurasian haplogroups and agrees well with the genetic structure of central-east European populations. Our data attest that latitudinal gradients characterize the variation of the uniparentally transmitted gene pools of modern Belarusians. In particular, the Y-chromosome reflects movements of people in central-east Europe, starting probably as early as the beginning of the Holocene. Furthermore, the matrilineal legacy of Belarusians retains two rare mitochondrial DNA haplogroups, N1a3 and N3, whose phylogeographies were explored in detail after de novo sequencing of 20 and 13 complete mitogenomes, respectively, from all over Eurasia. Our phylogeographic analyses reveal that two mitochondrial DNA lineages, N3 and N1a3, both of Middle Eastern origin, might mark distinct events of matrilineal gene flow to Europe: during the mid-Holocene period and around the Pleistocene-Holocene transition, respectively.

BACKGROUND AND OBJECTIVE. Variation of osteoporosis in the population is the result of an interaction between the genotype and the environment, and the genetic causes of osteoporosis are being widely investigated. The aim of this study was to analyze the association between the polymorphisms of the vitamin D receptor (VDR), type I collagen (COL1A1), and lactase (LCT) genes and severe postmenopausal osteoporosis as well as bone mineral density (BMD). MATERIAL AND METHODS. A total of 54 women with severe postmenopausal osteoporosis and 77 controls (mean age, 58.3 years [SD, 6.2] and 56.7 years [SD, 7.42], respectively) were included into the study. The subjects were recruited at the City Center for Osteoporosis Prevention (Minsk, Belarus). Dual-energy x-ray absorptiometry was used to measure bone mineral density at the lumbar spine and the femoral neck. Severe osteoporosis was diagnosed in the women with the clinical diagnosis of postmenopausal osteoporosis and at least 1 fragility fracture. The control group included women without osteoporosis. Polymorphic sites in osteoporosis predisposition genes (ApaI, BsmI, TaqI, and Cdx2 of the VDR gene, G2046T of the COL1A1 gene, and T-13910C of the LCT gene) were determined using the polymerase chain reaction on the deoxyribonucleic acid isolated from dried bloodspots. RESULTS. The data showed that the ApaI and BsmI polymorphisms of the VDR gene and T- 13910C of the LCT gene were associated with severe postmenopausal osteoporosis in the analyzed Belarusian women (P<0.01). A statistically significant positive correlation between the VDR risk genotypes ApaI and TaqI and bone mineral density was found (P<0.05). CONCLUSIONS. The findings of this study suggest that at least the ApaI and BsmI polymorphisms of the VDR gene and T-13910C of the LCT gene are associated with the risk of postmenopausal osteoporosis in our sample of the Belarusian women.

Mitochondrial (mt) and chloroplast (ct) genome inheritance was studied in barley-wheat hybrids, as were their progenies obtained from backcrosses with different common wheat cultivars, by monitoring the composition of 4 mtDNA (coxI, a 5'-flanking region of cob, nad3-orf156, and 5'-upstream region of 18S/5S) and 2 ctDNA (simple-sequence repeat locus downstream of trnS and a 3'-flanking region of rbcL) loci. In male sterile F1 and BC1 plants, maternal barley mtDNA fragments were mainly detected and very low levels of paternal wheat fragments were occasionally detected by PCR in coxI, a 5'-flanking region of cob and nad3-orf156, whereas a 5'-upstream region of 18S/5S showed clear heteroplasmy, containing both maternal and paternal copies, with maternal copies prevailing. Plants showing such heteroplasmic mtDNA composition remained either semisterile or became completely sterile in the later backcross generations. Only maternal ctDNA copies were detected in these plants. In 3 stable, self-fertile, and vigourous lines obtained in the advanced backcross generations and possessing recombinant wheat nuclear genome, however, only mt- and ctDNA copies of wheat parents were detected; thus, the original alloplasmic condition appeared to be lost. Our results suggest that transmission followed by selective replication of the paternal wheat organellar DNA leads to a paternally oriented shift of the organellar DNA composition in barley-wheat hybrids, which correlates with the restoration of fertility and plant vigour. These 2 processes seem to be related to nucleocytoplasmic compatibility and to be under the control of the nuclear genome composition.

Abstract Two major genes ( eam8 and eam10 ) and two quantitative trait loci (QTL) determining flowering time in barley were associated with restriction fragment length polymorphism markers. The loci eam8 and eam10 were found to map in regions of chromosomes 1HL and 3HL, respectively, already estimated from previous classical linkage analyses. While investigating doubled haploid lines of a spring habit barley mapping population, two QTL for flowering time were detected on chromosomes 1HL and 7HS, respectively, when the material was grown under long photoperiod conditions. When growing the same lines under short photoperiod, no QTL were discernible. Allelic and homoeologous relationships with flowering time loci described earlier in barley and other Triticeae species are discussed.

Human mitochondrial DNA haplogroup U is among the initial maternal founders in Southwest Asia and Europe and one that best indicates matrilineal genetic continuity between late Pleistocene hunter-gatherer groups and present-day populations of Europe. While most haplogroup U subclades are older than 30 thousand years, the comparatively recent coalescence time of the extant variation of haplogroup U7 (~16-19 thousand years ago) suggests that its current distribution is the consequence of more recent dispersal events, despite its wide geographical range across Europe, the Near East and South Asia. Here we report 267 new U7 mitogenomes that - analysed alongside 100 published ones - enable us to discern at least two distinct temporal phases of dispersal, both of which most likely emanated from the Near East. The earlier one began prior to the Holocene (~11.5 thousand years ago) towards South Asia, while the later dispersal took place more recently towards Mediterranean Europe during the Neolithic (~8 thousand years ago). These findings imply that the carriers of haplogroup U7 spread to South Asia and Europe before the suggested Bronze Age expansion of Indo-European languages from the Pontic-Caspian Steppe region.

Caloric restriction (CR) has positive effects on health and longevity. CR in mammals implements time-restricted (TR) feeding, a short period of feeding followed by prolonged fasting. Periodic fasting, in the form of TR or mealtime, improves metabolism without reduction in caloric intake. In order to understand the relative contribution of reduced food intake and periodic fasting to the health benefits of CR, we compared physiological and metabolic changes induced by CR and TR (without reduced food intake) in mice. CR significantly reduced blood glucose and insulin around the clock, improved glucose tolerance, and increased insulin sensitivity (IS). TR reduced blood insulin and increased insulin sensitivity, but in contrast to CR, TR did not improve glucose homeostasis. Liver expression of circadian clock genes was affected by both diets while the mRNA expression of glucose metabolism genes was significantly induced by CR, and not by TR, which is in agreement with the minor effect of TR on glucose metabolism. Thus, periodic fasting contributes to some metabolic benefits of CR, but TR is metabolically different from CR. This difference might contribute to differential effects of CR and TR on longevity.

Based on the cross (Triticum aestivum L. x Secale cereale L.) x T. aestivum L., wheat-rye substitution lines (2n = 42) were produced with karyotypes containing, instead of a pair of homologous wheat chromosomes, a homeologous pair of rye chromosomes. The chromosome composition of these lines was described by GISH and C-banding methods, and SSR analysis. The results of genomic in situ hybridization demonstrated that karyotype of these lines included one pair of rye chromosomes each and lacked wheat--rye translocations. C-banding and SSR markers were used to identify rye chromosomes and determine the wheat chromosomes at which the substitution occurred. The lines were designated 1R(1D), 2R(2D)2, 2R(2D)3, 3R(3B), 6R(6A)2. The chromosome composition of lines IR(1A), 2R(W)1, 5R(W), 5R(5A), and 6R(W)1, which were earlier obtained according to the same scheme for crossing, was characterized using methods of telocentric analysis, GISH, C-banding, and SSR analysis. These lines were identified as 1R(1A), 2R(2D)1, 5R(5D), 5R(5A), and 6R(6A)1, C-banding of chromosomes belonging to line 1R(1A) revealed the presence of two translocated chromosomes (3DS.3DL-del. and 4AL.W) during simultaneous amplification of SSR markers located on 3DL and 4AS arms. The "combined" long arm of the newly derived chromosome 4A is assumed to be formed from the long arm of chromosome 4AS itself and a deleted segment 3DL. All examined lines are cytologically stable, except for 3R(3B), which does not affect the stability of rye 3R chromosome transfer. Chromosome identification and classification of the lines will permit them to be models for genetic studies that can be used thereafter as promising "secondary gene pools" for the purpose of plant breeding.

Plant fibers find wide application in various fields that demand specific parameters of fiber quality. To develop approaches for the improvement of yield and quality of bast fibers, the knowledge of the fiber developmental stages and of the key molecular players that are responsible for a certain parameter, are vitally important. In the present review the key stages of fiber development, such as initiation, intrusive growth, and formation of thickened cell wall layers (secondary and tertiary cell walls) are considered, as well as the impact of each of these stages on the final parameters of fiber yield and quality. The problems and perspectives of crop quality regulation are discussed.

Vitamin D plays an important role in bone metabolism and is important for the prevention of multifactorial pathologies, including osteoporosis (OP). The biological action of vitamin is realized through its receptor, which is coded by the VDR gene. VDR gene polymorphism can influence individual predisposition to OP and response to vitamin D supplementation. The aim of this work was to reveal the effects of VDR gene ApaI rs7975232, BsmI rs1544410, TaqI rs731236, FokI rs2228570, and Cdx2 rs11568820 variants on bone mineral density (BMD), 25-hydroxyvitamin D level, and OP risk in Belarusian women. Methods. The case group included 355 women with postmenopausal OP, and the control group comprised 247 women who met the inclusion criteria. TaqMan genotyping assay was used to determine VDR gene variants. Results. Rs7975232 A/A, rs1544410 T/T, and rs731236 G/G single variants and their A-T-G haplotype showed a significant association with increased OP risk (for A-T-G, OR = 1.8, p = 0.0001) and decreased BMD (A-T-G, −0.09 g/cm2, p = 0.0001). The rs11568820 A-allele showed a protective effect on BMD (+0.22 g/cm2, p = 0.027). A significant dose effect with 25(OH)D was found for rs1544410, rs731236, and rs11568820 genotypes. Rs731236 A/A was associated with the 25(OH)D deficiency state. Conclusion. Our novel data on the relationship between VDR gene variants and BMD, 25(OH)D level, and OP risk highlights the importance of genetic markers for personalized medicine strategy.

Vitamin D receptor (VDR) is one of the main mediators of vitamin D biological activity. VDR dysfunction might substantially contribute to development of postmenopausal osteoporosis (PMO). Numerous studies have revealed the effects of several VDR gene variants on osteoporosis risk, although significant variation in different ethnicities has been suggested. The main purpose of this work was to assess the frequency of distribution of VDR genetic variants with established effect and evaluate their haplotype association with the risk of PMO in a cohort of Belarusian and Lithuanian women. Case group included women with postmenopausal osteoporosis (n=149), the control group comprised women with normal bone mineral density (BMD) and without previous fragility fractures (n=172). Both groups were matched for age, height, sex, and BMI – no statistically significant differences observed. VDR gene polymorphic variants (ApaI rs7975232, BsmI rs1544410, TaqI rs731236, and Cdx2 rs11568820) were determined using polymerase chain reaction and restriction fragment length polymorphism. The lumbar spine (L1-L4) and femoral neck BMD was measured using dual-energy X-ray absorptiometry. Association between each VDR variant and PMO risk was assessed using multiple logistic regression. The genotyping revealed statistically significant difference in the rs7975232 genotype frequencies between the patients and the controls (homozygous C/C genotype was overrepresented in patients, P=0.008). Patients with osteoporosis were also three times more likely to carry the rs1544410 G/G genotype, when compared to controls. We found that rs7975232, rs1544410 and rs731236 variants were in a strong direct linkage disequilibrium (P<0.0001), suggesting that risk alleles of these markers are preferably inherited jointly. For the bearers of C-G-C haplotype (consisting of rs7975232, rs1544410 and rs731236 unfavorable alleles), the risk of PMO was significantly higher (OR=4.7, 95% CI 2.8 to 8.1, p<0.0001) compared to controls. This haplotype was significantly over-represented in PMO group compared to all other haplotypes. Our findings highlight the importance of identified haplotypes of VDR gene variants. Complex screening of these genetic markers can be used to implement personalized clinical approach for prevention, treatment and rehabilitation programs.