Fat metabolism is associated with telomere length in six population-based studies.

Telomeres are repetitive DNA sequences located at the end of chromosomes, which are associated to biological aging, cardiovascular disease, cancer and mortality. Lipid and fatty acid metabolism have been associated with telomere shortening. We have conducted an in-depth study investigating the association of metabolic biomarkers with telomere length (LTL). We performed an association analysis of 226 metabolic biomarkers with LTL using data from 11 775 individuals from six independent population-based cohorts (BBMRI-NL consortium). Metabolic biomarkers include lipoprotein lipids and subclasses, fatty acids, amino acids, glycolysis measures and ketone bodies. LTL was measured by quantitative polymerase chain reaction or FlowFISH. Linear regression analysis was performed adjusting for age, sex, lipid-lowering medication and cohort-specific covariates (model 1) and additionally for body mass index (BMI) and smoking (model 2), followed by inverse variance-weighted meta-analyses (significance threshold Pmeta = 6.5 × 10-4). We identified four metabolic biomarkers positively associated with LTL, including two cholesterol to lipid ratios in small VLDL (S-VLDL-C % and S-VLDL-CE %) and two omega-6 fatty acid ratios (FAw6/FA and LA/FA). After additionally adjusting for BMI and smoking, these metabolic biomarkers remained associated with LTL with similar effect estimates. In addition, cholesterol esters in very small VLDL (XS-VLDL-CE) became significantly associated with LTL (P = 3.6 × 10-4). We replicated the association of FAw6/FA with LTL in an independent dataset of 7845 individuals (P = 1.9 × 10-4). To conclude, we identified multiple metabolic biomarkers involved in lipid and fatty acid metabolism that may be involved in LTL biology. Longitudinal studies are needed to exclude reversed causation.

Genome-wide Association Analysis in Humans Links Nucleotide Metabolism to Leukocyte Telomere Length.

Leukocyte telomere length (LTL) is a heritable biomarker of genomic aging. In this study, we perform a genome-wide meta-analysis of LTL by pooling densely genotyped and imputed association results across large-scale European-descent studies including up to 78,592 individuals. We identify 49 genomic regions at a false dicovery rate (FDR) < 0.05 threshold and prioritize genes at 31, with five highlighting nucleotide metabolism as an important regulator of LTL. We report six genome-wide significant loci in or near SENP7, MOB1B, CARMIL1, PRRC2A, TERF2, and RFWD3, and our results support recently identified PARP1, POT1, ATM, and MPHOSPH6 loci. Phenome-wide analyses in >350,000 UK Biobank participants suggest that genetically shorter telomere length increases the risk of hypothyroidism and decreases the risk of thyroid cancer, lymphoma, and a range of proliferative conditions. Our results replicate previously reported associations with increased risk of coronary artery disease and lower risk for multiple cancer types. Our findings substantially expand current knowledge on genes that regulate LTL and their impact on human health and disease.

Somatic TARDBP variants as a cause of semantic dementia.

The aetiology of late-onset neurodegenerative diseases is largely unknown. Here we investigated whether de novo somatic variants for semantic dementia can be detected, thereby arguing for a more general role of somatic variants in neurodegenerative disease. Semantic dementia is characterized by a non-familial occurrence, early onset (<65 years), focal temporal atrophy and TDP-43 pathology. To test whether somatic variants in neural progenitor cells during brain development might lead to semantic dementia, we compared deep exome sequencing data of DNA derived from brain and blood of 16 semantic dementia cases. Somatic variants observed in brain tissue and absent in blood were validated using amplicon sequencing and digital PCR. We identified two variants in exon one of the TARDBP gene (L41F and R42H) at low level (1-3%) in cortical regions and in dentate gyrus in two semantic dementia brains, respectively. The pathogenicity of both variants is supported by demonstrating impaired splicing regulation of TDP-43 and by altered subcellular localization of the mutant TDP-43 protein. These findings indicate that somatic variants may cause semantic dementia as a non-hereditary neurodegenerative disease, which might be exemplary for other late-onset neurodegenerative disorders.

Reduced penetrance of pathogenic ACMG variants in a deeply phenotyped cohort study and evaluation of ClinVar classification over time.

PURPOSE: We studied the penetrance of pathogenically classified variants in an elderly Dutch population from the Rotterdam Study, for which deep phenotyping is available. We screened the 59 actionable genes for which reporting of known pathogenic variants was recommended by the American College of Medical Genetics and Genomics (ACMG), and demonstrate that determining what constitutes a known pathogenic variant can be quite challenging. METHODS: We defined "known pathogenic" as classified pathogenic by both ClinVar and the Human Gene Mutation Database (HGMD). In 2628 individuals, we performed exome sequencing and identified known pathogenic variants. We investigated the clinical records of carriers and evaluated clinical events during 25 years of follow-up for evidence of variant pathogenicity. RESULTS: Of 3815 variants detected in the 59 ACMG genes, 17 variants were considered known pathogenic. For 14/17 variants the ClinVar classification had changed over time. Of 24 confirmed carriers of these variants, we observed at least one clinical event possibly caused by the variant in only three participants (13%). CONCLUSION: We show that the definition of "known pathogenic" is often unclear and should be approached carefully. Additionally variants marked as known pathogenic do not always have clinical impact on their carriers. Definition and classification of true (individual) expected pathogenic impact should be defined carefully.

Short-Term, Combined Fasting and Exercise Improves Body Composition in Healthy Males.

Fasting enhances the beneficial metabolic outcomes of exercise; however, it is unknown whether body composition is favorably modified on the short term. A baseline-follow-up study was carried out to assess the effect of an established protocol involving short-term combined exercise with fasting on body composition. One hundred seven recreationally exercising males underwent a 10-day intervention across 15 fitness centers in the Netherlands involving a 3-day gradual decrease of food intake, a 3-day period with extremely low caloric intake, and a gradual 4-day increase to initial caloric intake, with daily 30-min submaximal cycling. Using dual-energy X-ray absorptiometry analysis, all subjects substantially lost total body mass (-3.9 ± 1.9 kg; p < .001) and fat mass (-3.3 ± 1.3 kg; p < .001). Average lean mass was lost (-0.6 ± 1.5 kg; p < .001), but lean mass as a percentage of total body mass was not reduced. The authors observed a loss of -3.9 ± 1.9% android fat over total fat mass (p < .001), a loss of -2.2 ± 1.9% gynoid over total fat mass (p < .001), and reduced android/gynoid ratios (-0.05 ± 0.1; p < .001). Analyzing 15 preselected single-nucleotide polymorphisms in 13 metabolism-related genes revealed trending associations for thyroid state-related single-nucleotide polymorphisms rs225014 (deiodinase 2) and rs35767 (insulin-like growth factor1), and rs1053049 (PPARD). In conclusion, a short period of combined fasting and exercise leads to a substantial loss of body and fat mass without a loss of lean mass as a percentage of total mass.

Characteristics of de novo structural changes in the human genome.

Small insertions and deletions (indels) and large structural variations (SVs) are major contributors to human genetic diversity and disease. However, mutation rates and characteristics of de novo indels and SVs in the general population have remained largely unexplored. We report 332 validated de novo structural changes identified in whole genomes of 250 families, including complex indels, retrotransposon insertions, and interchromosomal events. These data indicate a mutation rate of 2.94 indels (1-20 bp) and 0.16 SVs (>20 bp) per generation. De novo structural changes affect on average 4.1 kbp of genomic sequence and 29 coding bases per generation, which is 91 and 52 times more nucleotides than de novo substitutions, respectively. This contrasts with the equal genomic footprint of inherited SVs and substitutions. An excess of structural changes originated on paternal haplotypes. Additionally, we observed a nonuniform distribution of de novo SVs across offspring. These results reveal the importance of different mutational mechanisms to changes in human genome structure across generations.

Vitamin D Receptor Polymorphisms Are Associated with Reduced Esophageal Vitamin D Receptor Expression and Reduced Esophageal Adenocarcinoma Risk.

Epidemiological studies indicate that vitamin D exerts a protective effect on the development of various solid cancers. However, concerns have been raised regarding the potential deleterious role of high vitamin D levels in the development of esophageal adenocarcinoma (EAC). This study investigated genetic variation in the vitamin D receptor (VDR) in relation to its expression and risk of Barrett esophagus (BE) and EAC. VDR gene regulation was investigated by immunohistochemistry, reverse transcriptase-polymerase chain reaction (RT-PCR) and gel shift assays. Fifteen haplotype tagging single-nucleotide polymorphisms (SNPs) of the VDR gene were analyzed in 858 patients with reflux esophagitis (RE), BE or EAC and 202 healthy controls. VDR mRNA expression was higher in BE compared with squamous epithelium. VDR protein was located in the nucleus in BE. An rs1989969T/rs2238135G haplotype was identified in the 5' regulatory region of the VDR gene. It was associated with an approximately two-fold reduced risk of RE, BE and EAC. Analysis of a replication cohort was done for BE that confirmed this. The rs1989969T allele causes a GATA-1 transcription factor binding site to appear. The signaling of GATA-1, which is regarded as a negative transcriptional regulator, could explain the findings for rs1989969. The rs2238135G allele was associated with a significantly reduced VDR expression in BE; for the rs1989969T allele, a trend in reduced VDR expression was observed. We identified a VDR haplotype associated with reduced esophageal VDR expression and a reduced incidence of RE, BE and EAC. This VDR haplotype could be useful in identifying individuals who benefit most from vitamin D chemoprevention.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus.

BACKGROUND: Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk. AIM: To identify CNVs associated with osteoporotic bone fracture risk. METHOD: We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies. RESULTS: A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk. CONCLUSIONS: These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

A new variant in the ZCCHC8 gene: diverse clinical phenotypes and expression in the lung.

Introduction: Pulmonary fibrosis is a severe disease which can be familial. A genetic cause can only be found in ∼40% of families. Searching for shared novel genetic variants may aid the discovery of new genetic causes of disease. Methods: Whole-exome sequencing was performed in 152 unrelated patients with a suspected genetic cause of pulmonary fibrosis from the St Antonius interstitial lung disease biobank. Variants of interest were selected by filtering for novel, potentially deleterious variants that were present in at least three unrelated pulmonary fibrosis patients. Results: The novel c.586G>A p.(E196K) variant in the ZCCHC8 gene was observed in three unrelated patients: two familial patients and one sporadic patient, who was later genealogically linked to one of the families. The variant was identified in nine additional relatives with pulmonary fibrosis and other telomere-related phenotypes, such as pulmonary arterial venous malformations, emphysema, myelodysplastic syndrome, acute myeloid leukaemia and dyskeratosis congenita. One family showed incomplete segregation, with absence of the variant in one pulmonary fibrosis patient who carried a PARN variant. The majority of ZCCHC8 variant carriers showed short telomeres in blood. ZCCHC8 protein was located in different lung cell types, including alveolar type 2 (AT2) pneumocytes, the culprit cells in pulmonary fibrosis. AT2 cells showed telomere shortening and increased DNA damage, which was comparable to patients with sporadic pulmonary fibrosis and those with pulmonary fibrosis carrying a telomere-related gene variant, respectively. Discussion: The ZCCHC8 c.586G>A variant confirms the involvement of ZCCHC8 in pulmonary fibrosis and short-telomere syndromes and underlines the importance of including the ZCCHC8 gene in diagnostic gene panels for these diseases.

Disruption of TUFT1, a Desmosome-Associated Protein, Causes Skin Fragility, Woolly Hair, and Palmoplantar Keratoderma.

Desmosomes are dynamic complex protein structures involved in cellular adhesion. Disruption of these structures by loss-of-function variants in desmosomal genes leads to a variety of skin- and heart-related phenotypes. In this study, we report TUFT1 as a desmosome-associated protein, implicated in epidermal integrity. In two siblings with mild skin fragility, woolly hair, and mild palmoplantar keratoderma but without a cardiac phenotype, we identified a homozygous splice-site variant in the TUFT1 gene, leading to aberrant mRNA splicing and loss of TUFT1 protein. Patients' skin and keratinocytes showed acantholysis, perinuclear retraction of intermediate filaments, and reduced mechanical stress resistance. Immunolabeling and transfection studies showed that TUFT1 is positioned within the desmosome and that its location is dependent on the presence of the desmoplakin carboxy-terminal tail. A Tuft1-knockout mouse model mimicked the patients' phenotypes. Altogether, this study reveals TUFT1 as a desmosome-associated protein, whose absence causes skin fragility, woolly hair, and palmoplantar keratoderma.

Meta-analysis of genome-wide scans for human adult stature identifies novel Loci and associations with measures of skeletal frame size.

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value = 6.1x10(-8) and rs910316 in TMED10, P-value = 1.4x10(-7)) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value = 3x10(-7) and rs849141 in JAZF1, P-value = 3.2x10(-11)). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value = 4x10(-5) and rs6817306 in LCORL, P-value = 4x10(-4)), hip axis length (including rs6830062 at LCORL, P-value = 4.8x10(-4) and rs4911494 at UQCC, P-value = 1.9x10(-4)), and femur length (including rs710841 at PRKG2, P-value = 2.4x10(-5) and rs10946808 at HIST1H1D, P-value = 6.4x10(-6)). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.

Impact of SNP microarray analysis of compromised DNA on kinship classification success in the context of investigative genetic genealogy.

Single nucleotide polymorphism (SNP) data generated with microarray technologies have been used to solve murder cases via investigative leads obtained from identifying relatives of the unknown perpetrator included in accessible genomic databases, an approach referred to as investigative genetic genealogy (IGG). However, SNP microarrays were developed for relatively high input DNA quantity and quality, while DNA typically obtainable from crime scene stains is of low DNA quantity and quality, and SNP microarray data obtained from compromised DNA are largely missing. By applying the Illumina Global Screening Array (GSA) to 264 DNA samples with systematically altered quantity and quality, we empirically tested the impact of SNP microarray analysis of compromised DNA on kinship classification success, as relevant in IGG. Reference data from manufacturer-recommended input DNA quality and quantity were used to estimate genotype accuracy in the compromised DNA samples and for simulating data of different degree relatives. Although stepwise decrease of input DNA amount from 200 ng to 6.25 pg led to decreased SNP call rates and increased genotyping errors, kinship classification success did not decrease down to 250 pg for siblings and 1st cousins, 1 ng for 2nd cousins, while at 25 pg and below kinship classification success was zero. Stepwise decrease of input DNA quality via increased DNA fragmentation resulted in the decrease of genotyping accuracy as well as kinship classification success, which went down to zero at the average DNA fragment size of 150 base pairs. Combining decreased DNA quantity and quality in mock casework and skeletal samples further highlighted possibilities and limitations. Overall, GSA analysis achieved maximal kinship classification success from 800 to 200 times lower input DNA quantities than manufacturer-recommended, although DNA quality plays a key role too, while compromised DNA produced false negative kinship classifications rather than false positive ones.

Three genome-wide association studies and a linkage analysis identify HERC2 as a human iris color gene.

Human iris color was one of the first traits for which Mendelian segregation was established. To date, the genetics of iris color is still not fully understood and is of interest, particularly in view of forensic applications. In three independent genome-wide association (GWA) studies of a total of 1406 persons and a genome-wide linkage study of 1292 relatives, all from the Netherlands, we found that the 15q13.1 region is the predominant region involved in human iris color. There were no other regions showing consistent genome-wide evidence for association and linkage to iris color. Single nucleotide polymorphisms (SNPs) in the HERC2 gene and, to a lesser extent, in the neighboring OCA2 gene were independently associated to iris color variation. OCA2 has been implicated in iris color previously. A replication study within two populations confirmed that the HERC2 gene is a new and significant determinant of human iris color variation, in addition to OCA2. Furthermore, HERC2 rs916977 showed a clinal allele distribution across 23 European populations, which was significantly correlated to iris color variation. We suggest that genetic variants regulating expression of the OCA2 gene exist in the HERC2 gene or, alternatively, within the 11.7 kb of sequence between OCA2 and HERC2, and that most iris color variation in Europeans is explained by those two genes. Testing markers in the HERC2-OCA2 region may be useful in forensic applications to predict eye color phenotypes of unknown persons of European genetic origin.

Recognition of scared faces and the serotonin transporter gene in young children: the Generation R Study.

BACKGROUND: Previous research highlights the significance of a functional polymorphism located in the promoter region (5-HTTLPR) of the serotonin transporter gene in emotional behaviour. This study examined the effect of the 5-HTTLPR polymorphism on emotion processing in a large number of healthy preschoolers. METHODS: The 5-HTTLPR genotype was classified in 605 children as homozygous for the short allele (SS), homozygous for the long allele (LL), or heterozygous (LS). Emotion-processing was assessed using age-appropriate computer tasks where children matched happy, sad, angry, and fearful facial expressions preceded by a shape-matching task to assess basic matching ability. RESULTS: We found that young children could differentiate between emotion categories (F = 12.1, p < .001). The effect of 5-HTTLPR genotype depended on the emotion category presented (F = 2.3, p = .031). This effect was explained by the finding that SS children were less accurate at recognising fearful faces than LL or LS children (F = 5.3, p = .005). We did not find any significant differences as a result of 5-HTTLPR genotype for happy, sad or angry expressions (p > .05). CONCLUSIONS: Results indicate that 5-HTTLPR allele status selectively impacts the processing of fearful but not other facial expressions. This pattern is already apparent in very young typically developing children. Results may signal an early vulnerability for affective problems before disorders emerge.

Telomere Length and the Risk of Alzheimer's Disease: The Rotterdam Study.

There is a wide interest in biomarkers that capture the burden of detrimental factors as these accumulate with the passage of time, i.e., increasing age. Telomere length has received considerable attention as such a marker, because it is easily quantified and it may aid in disentangling the etiology of dementia or serve as predictive marker. We determined the association of telomere length with risk of Alzheimer's disease and all-cause dementia in a population-based setting. Within the Rotterdam Study, we performed quantitative PCR to measure mean leukocyte telomere length in blood. We determined the association of telomere length with risk of Alzheimer's disease until 2016, using Cox regression models. Of 1,961 participants (mean age 71.4±9.3 years, 57.1% women) with a median follow-up of 8.3 years, 237 individuals were diagnosed with Alzheimer's disease. We found a U-shaped association between telomere length and risk of Alzheimer's disease: compared to the middle tertile the adjusted hazard ratio was 1.59 (95% confidence interval (CI), 1.13-2.23) for the lowest tertile and 1.47 (1.03-2.10) for the highest tertile. Results were similarly U-shaped but slightly attenuated for all-cause dementia. In conclusion, shorter and longer telomere length are both associated with an increased risk of Alzheimer's disease in the general population.

Exome Sequencing Analysis Identifies Rare Variants in ATM and RPL8 That Are Associated With Shorter Telomere Length.

Telomeres are important for maintaining genomic stability. Telomere length has been associated with aging, disease, and mortality and is highly heritable (∼82%). In this study, we aimed to identify rare genetic variants associated with telomere length using whole-exome sequence data. We studied 1,303 participants of the Erasmus Rucphen Family (ERF) study, 1,259 of the Rotterdam Study (RS), and 674 of the British Heart Foundation Family Heart Study (BHF-FHS). We conducted two analyses, first we analyzed the family-based ERF study and used the RS and BHF-FHS for replication. Second, we combined the summary data of the three studies in a meta-analysis. Telomere length was measured by quantitative polymerase chain reaction in blood. We identified nine rare variants significantly associated with telomere length (p-value < 1.42 × 10-7, minor allele frequency of 0.2-0.5%) in the ERF study. Eight of these variants (in C11orf65, ACAT1, NPAT, ATM, KDELC2, and EXPH5) were located on chromosome 11q22.3 that contains ATM, a gene involved in telomere maintenance. Although we were unable to replicate the variants in the RS and BHF-FHS (p-value ≥ 0.21), segregation analysis showed that all variants segregate with shorter telomere length in a family. In the meta-analysis of all studies, a nominally significant association with LTL was observed with a rare variant in RPL8 (p-value = 1.48 × 10-6), which has previously been associated with age. Additionally, a novel rare variant in the known RTEL1 locus showed suggestive evidence for association (p-value = 1.18 × 10-4) with LTL. To conclude, we identified novel rare variants associated with telomere length. Larger samples size are needed to confirm these findings and to identify additional variants.

Vitamin D binding protein genotype and osteoporosis.

Osteoporosis is a bone disease leading to an increased fracture risk. It is considered a complex multifactorial genetic disorder with interaction of environmental and genetic factors. As a candidate gene for osteoporosis, we studied vitamin D binding protein (DBP, or group-specific component, Gc), which binds to and transports vitamin D to target tissues to maintain calcium homeostasis through the vitamin D endocrine system. DBP can also be converted to DBP-macrophage activating factor (DBP-MAF), which mediates bone resorption by directly activating osteoclasts. We summarized the genetic linkage structure of the DBP gene. We genotyped two single-nucleotide polymorphisms (SNPs, rs7041 = Glu416Asp and rs4588 = Thr420Lys) in 6,181 elderly Caucasians and investigated interactions of the DBP genotype with vitamin D receptor (VDR) genotype and dietary calcium intake in relation to fracture risk. Haplotypes of the DBP SNPs correspond to protein variations referred to as Gc1s (haplotype 1), Gc2 (haplotype 2), and Gc1f (haplotype3). In a subgroup of 1,312 subjects, DBP genotype was found to be associated with increased and decreased serum 25-(OH)D(3) for haplotype 1 (P = 3 x 10(-4)) and haplotype 2 (P = 3 x 10(-6)), respectively. Similar associations were observed for 1,25-(OH)(2)D(3). The DBP genotype was not significantly associated with fracture risk in the entire study population. Yet, we observed interaction between DBP and VDR haplotypes in determining fracture risk. In the DBP haplotype 1-carrier group, subjects of homozygous VDR block 5-haplotype 1 had 33% increased fracture risk compared to noncarriers (P = 0.005). In a subgroup with dietary calcium intake <1.09 g/day, the hazard ratio (95% confidence interval) for fracture risk of DBP hap1-homozygote versus noncarrier was 1.47 (1.06-2.05). All associations were independent of age and gender. Our study demonstrated that the genetic effect of the DBP gene on fracture risk appears only in combination with other genetic and environmental risk factors for bone metabolism.

Intestinal microbiome composition and its relation to joint pain and inflammation.

Macrophage-mediated inflammation is thought to have a causal role in osteoarthritis-related pain and severity, and has been suggested to be triggered by endotoxins produced by the gastrointestinal microbiome. Here we investigate the relationship between joint pain and the gastrointestinal microbiome composition, and osteoarthritis-related knee pain in the Rotterdam Study; a large population based cohort study. We show that abundance of Streptococcus species is associated with increased knee pain, which we validate by absolute quantification of Streptococcus species. In addition, we replicate these results in 867 Caucasian adults of the Lifelines-DEEP study. Finally we show evidence that this association is driven by local inflammation in the knee joint. Our results indicate the microbiome is a possible therapeutic target for osteoarthritis-related knee pain.

The RIZ Pro704 insertion-deletion polymorphism, bone mineral density and fracture risk: the Rotterdam study.

Estrogens play a major role in the maintenance of bone and bone strength, and they exert their effects via estrogen receptors. Recently, an estrogen receptor alpha (ESR1) specific co-activator, retinoblastoma-interacting zinc-finger protein (RIZ1, 1p36), was shown to strongly enhance ESR1 function in vitro. The same study showed that a Proline insertion-deletion polymorphism at amino acid position 704 (Pro704 ins/del) in the RIZ1 gene was associated with heel BMD in young Swedish women. We tested the relation between the RIZ1 Pro704 ins/del polymorphism and BMD and fracture risk in Caucasian elderly men and women of the Rotterdam study. We also examined whether estradiol levels (measured in a subset) or genetic variation in ESR1 influenced this relation. In 2424 men and 3517 women from the Rotterdam study, RIZ1 genotypes were determined and associations with BMD (lumbar spine and femoral neck) and fracture risk were analysed. We recorded 374 vertebral fractures at baseline and during 6.4+/-0.4 (SD) years of follow-up, and 1219 incident non-vertebral fractures during 7.4+/-3.3 (SD) years of follow-up. The allele frequency of the Pro704 insertion was 41%, the genotype distribution was in Hardy-Weinberg Equilibrium (P=0.94). We found no association of this polymorphism with BMD or fracture risk. Stratification for gender, estradiol levels or interaction with ESR1 risk haplotype did not change these results. In conclusion, in this large study we observed no association of the RIZ1 Pro704 insertion-deletion polymorphism with BMD or fracture risk. This suggests this polymorphism to play a minor role, if any, as a genetic determinant of osteoporosis in elderly subjects.