High degree of polyclonality hinders somatic mutation calling in lung brush samples of COPD cases and controls.

Chronic obstructive pulmonary disease (COPD) is induced by cigarette smoking and characterized by inflammation of airway tissue. Since smokers with COPD have a higher risk of developing lung cancer than those without, we hypothesized that they carry more mutations in affected tissue. We called somatic mutations in airway brush samples from medium-coverage whole genome sequencing data from healthy never and ex-smokers (n = 8), as well as from ex-smokers with variable degrees of COPD (n = 4). Owing to the limited concordance of resulting calls between the applied tools we built a consensus, a strategy that was validated with high accuracy for cancer data. However, consensus calls showed little promise of representing true positives due to low mappability of corresponding sequence reads and high overlap with positions harbouring known genetic polymorphisms. A targeted re-sequencing approach suggested that only few mutations would survive stringent verification testing and that our data did not allow the inference of any difference in the mutational load of bronchial brush samples between former smoking COPD cases and controls. High polyclonality in airway brush samples renders medium-depth sequencing insufficient to provide the resolution to detect somatic mutations. Deep sequencing data of airway biopsies are needed to tackle the question.

Whole-genome fingerprint of the DNA methylome during human B cell differentiation.

We analyzed the DNA methylome of ten subpopulations spanning the entire B cell differentiation program by whole-genome bisulfite sequencing and high-density microarrays. We observed that non-CpG methylation disappeared upon B cell commitment, whereas CpG methylation changed extensively during B cell maturation, showing an accumulative pattern and affecting around 30% of all measured CpG sites. Early differentiation stages mainly displayed enhancer demethylation, which was associated with upregulation of key B cell transcription factors and affected multiple genes involved in B cell biology. Late differentiation stages, in contrast, showed extensive demethylation of heterochromatin and methylation gain at Polycomb-repressed areas, and genes with apparent functional impact in B cells were not affected. This signature, which has previously been linked to aging and cancer, was particularly widespread in mature cells with an extended lifespan. Comparing B cell neoplasms with their normal counterparts, we determined that they frequently acquire methylation changes in regions already undergoing dynamic methylation during normal B cell differentiation.

Whole-epigenome analysis in multiple myeloma reveals DNA hypermethylation of B cell-specific enhancers.

While analyzing the DNA methylome of multiple myeloma (MM), a plasma cell neoplasm, by whole-genome bisulfite sequencing and high-density arrays, we observed a highly heterogeneous pattern globally characterized by regional DNA hypermethylation embedded in extensive hypomethylation. In contrast to the widely reported DNA hypermethylation of promoter-associated CpG islands (CGIs) in cancer, hypermethylated sites in MM, as opposed to normal plasma cells, were located outside CpG islands and were unexpectedly associated with intronic enhancer regions defined in normal B cells and plasma cells. Both RNA-seq and in vitro reporter assays indicated that enhancer hypermethylation is globally associated with down-regulation of its host genes. ChIP-seq and DNase-seq further revealed that DNA hypermethylation in these regions is related to enhancer decommissioning. Hypermethylated enhancer regions overlapped with binding sites of B cell-specific transcription factors (TFs) and the degree of enhancer methylation inversely correlated with expression levels of these TFs in MM. Furthermore, hypermethylated regions in MM were methylated in stem cells and gradually became demethylated during normal B-cell differentiation, suggesting that MM cells either reacquire epigenetic features of undifferentiated cells or maintain an epigenetic signature of a putative myeloma stem cell progenitor. Overall, we have identified DNA hypermethylation of developmentally regulated enhancers as a new type of epigenetic modification associated with the pathogenesis of MM.

Functional relevance of the BMD-associated polymorphism rs312009: novel involvement of RUNX2 in LRP5 transcriptional regulation.

LRP5 is an osteoporosis susceptibility gene. Association analyses reveal that individual single-nucleotide polymorphisms (SNPs) determine variation in bone mineral density (BMD) among individuals as well as fracture risk. In a previous study, we identified a lumbar spine BMD-associated SNP, rs312009, located in the LRP5 5' region. A RUNX2 binding site was identified in this region by gel-shift experiments. Here we test the functionality of this SNP and examine whether RUNX2 is indeed a regulator of LRP5 expression. Gene reporter assays were used to test rs312009 functionality. Bioinformatic predictive tools and gel-shift and gene reporter assays were used to identify and characterize additional RUNX2 binding elements in the 3.3-kb region upstream of LRP5. Allelic differences in the transcriptional activity of rs312009 were observed in two osteoblastic cell lines, the T allele being a better transcriber than the C allele. RUNX2 cotransfection in HeLa cells revealed that the LRP5 5' region responded to RUNX2 in a dose-dependent manner and that the previously identified RUNX2 binding site participated in this response. Also, RUNX2 inhibition by RNAi led to nearly 60% reduction of endogenous LRP5 mRNA in U-2 OS cells. Four other RUNX2 binding sites were identified in the 5' region of LRP5. Luciferase experiments revealed the involvement of each of them in the RUNX2 response. The allelic differences observed point to the involvement of rs312009 as a functional SNP in the observed association. To our knowledge, this is the first time that the direct action of RUNX2 on LRP5 has been described. This adds evidence to previously described links between two important bone-regulating systems: the RUNX2 transcription-factor cascade and the Wnt signaling pathway.

Effect of IL-1beta, PGE(2), and TGF-beta1 on the expression of OPG and RANKL in normal and osteoporotic primary human osteoblasts.

The RANKL/RANK/OPG pathway is essential for bone remodeling regulation. Many hormones and cytokines are involved in regulating gene expression in most of the pathway components. Moreover, any deregulation of this pathway can alter bone metabolism, resulting in loss or gain of bone mass. Whether osteoblasts from osteoporotic and nonosteoporotic patients respond differently to cytokines is unknown. The aim of this study was to compare the effect of interleukin (IL)-1beta, proftaglandin E(2) (PGE(2)), and transforming growth factor-beta1 (TGF-beta1) treatments on OPG and RANKL gene expression in normal (n = 11) and osteoporotic (n = 8) primary osteoblasts. OPG and RANKL mRNA levels of primary human osteoblastic (hOB) cell cultures were assessed by real-time PCR. In all cultures, OPG mRNA increased significantly in response to IL-1beta treatment and decreased in response to TGF-beta1 whereas PGE(2) treatment had no effect. RANKL mRNA levels were significantly increased by all treatments. Differences in OPG and RANKL responses were observed between osteoporotic and nonosteoporotic hOB: in osteoporotic hOB, the OPG response to IL-1beta treatment was up to three times lower (P = 0.009), whereas that of RANKL response to TGF-beta1 was five times higher (P = 0.002) after 8 h of treatment, as compared with those in nonosteoporotic hOBs. In conclusion, osteoporotic hOB cells showed an anomalous response under cytokine stimulation, consistent with an enhanced osteoclastogenesis resulting in high levels of bone resorption.

A haplotype-based analysis of the LRP5 gene in relation to osteoporosis phenotypes in Spanish postmenopausal women.

LRP5 encodes the low-density lipoprotein receptor-related protein 5, a transmembrane protein involved in Wnt signaling. LRP5 is an important regulator of osteoblast growth and differentiation, affecting bone mass in vertebrates. Whether common variations in LRP5 are associated with normal BMD variation or osteoporotic phenotypes is of great relevance. We used a haplotype-based approach to search for common disease-associated variants in LRP5 in a cohort of 964 Spanish postmenopausal women. Twenty-four SNPs were selected, covering the LRP5 region, including the missense changes p.V667M and p.A1330V. The SNPs were genotyped and evaluated for association with BMD at the lumbar spine (LS) or femoral neck (FN) and with osteoporotic fracture, at single SNP and haplotype levels, by regression methods. Association with LS BMD was found for SNP 1, rs312009, located in the 5'-flanking region (p = 0.011, recessive model). SNP 6, rs2508836, in intron 1, was also associated with BMD, both at LS (p = 0.025, additive model) and FN (p = 0.031, recessive model). Two polymorphisms were associated with fracture: SNP 11, rs729635, in intron 1, and SNP 15, rs643892, in intron 5 (p = 0.007 additive model and p = 0.019 recessive model, respectively). Haplotype analyses did not provide additional information, except for haplotype "GC" of the block located at the 3'end of the gene. This haplotype spans intron 22 and the 3' untranslated region and was associated with FN BMD (p = 0.029, one copy of the haplotype versus none). In silico analyses showed that SNP 1 (rs312009) lies in a putative RUNX2 binding site. Electro-mobility shift assays confirmed RUNX2 binding to this site.

Promoter 2 -1025 T/C polymorphism in the RUNX2 gene is associated with femoral neck bmd in Spanish postmenopausal women.

Stimulation of bone formation is a key therapeutic target in osteoporosis. Runx2 is a runt domain transcription factor essential to osteoblast differentiation, bone remodeling, and fracture healing. Runx2 knockout mice exhibit a complete lack of ossification, while overexpression of this gene in transgenic mice results in an osteoporotic phenotype. Thus, RUNX2 is a good candidate for the genetic determination of osteoporosis. In this association study, the effects of the -330 G/T polymorphism in promoter 1 and the -1025 T/C polymorphism (rs7771980) in promoter 2 of RUNX2 were tested in relation to lumbar spine (LS) and femoral neck (FN) bone mineral density (BMD) in a cohort of 821 Spanish postmenopausal women. The minor allele frequencies for the two polymorphisms were 0.15 and 0.07, respectively. The two polymorphisms, located more than 90 kb apart, were not in linkage disequilibrium (D' = 0.27, r (2) = 0.028). In an ANCOVA test adjusting by weight, height, age, and years since menopause, the -330 G/T polymorphism was not associated with any of the phenotypes analyzed, while we found the -1025 T/C polymorphism to be associated with FN BMD (p = 0.001). In particular, individuals carrying the TC genotype had higher mean adjusted FN BMD values than those bearing the TT genotype. Our results highlight the importance of this RUNX2 promoter 2 polymorphism in FN BMD determination.