Population-specific genetic variation in large sequencing data sets: why more data is still better.

We have generated a next-generation whole-exome sequencing data set of 2628 participants of the population-based Rotterdam Study cohort, comprising 669 737 single-nucleotide variants and 24 019 short insertions and deletions. Because of broad and deep longitudinal phenotyping of the Rotterdam Study, this data set permits extensive interpretation of genetic variants on a range of clinically relevant outcomes, and is accessible as a control data set. We show that next-generation sequencing data sets yield a large degree of population-specific variants, which are not captured by other available large sequencing efforts, being ExAC, ESP, 1000G, UK10K, GoNL and DECODE.

The transcriptional landscape of age in human peripheral blood.

Disease incidences increase with age, but the molecular characteristics of ageing that lead to increased disease susceptibility remain inadequately understood. Here we perform a whole-blood gene expression meta-analysis in 14,983 individuals of European ancestry (including replication) and identify 1,497 genes that are differentially expressed with chronological age. The age-associated genes do not harbor more age-associated CpG-methylation sites than other genes, but are instead enriched for the presence of potentially functional CpG-methylation sites in enhancer and insulator regions that associate with both chronological age and gene expression levels. We further used the gene expression profiles to calculate the 'transcriptomic age' of an individual, and show that differences between transcriptomic age and chronological age are associated with biological features linked to ageing, such as blood pressure, cholesterol levels, fasting glucose, and body mass index. The transcriptomic prediction model adds biological relevance and complements existing epigenetic prediction models, and can be used by others to calculate transcriptomic age in external cohorts.

Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

A genome-wide association study identifies an osteoarthritis susceptibility locus on chromosome 7q22.

OBJECTIVE: To identify novel genes involved in osteoarthritis (OA), by means of a genome-wide association study. METHODS: We tested 500,510 single-nucleotide polymorphisms (SNPs) in 1,341 Dutch Caucasian OA cases and 3,496 Dutch Caucasian controls. SNPs associated with at least 2 OA phenotypes were analyzed in 14,938 OA cases and approximately 39,000 controls. Meta-analyses were performed using the program Comprehensive Meta-analysis, with P values <1 x 10(-7) considered genome-wide significant. RESULTS: The C allele of rs3815148 on chromosome 7q22 (minor allele frequency 23%; intron 12 of the COG5 gene) was associated with a 1.14-fold increased risk (95% confidence interval 1.09-1.19) of knee and/or hand OA (P = 8 x 10(-8)) and also with a 30% increased risk of knee OA progression (95% confidence interval 1.03-1.64) (P = 0.03). This SNP is in almost complete linkage disequilibrium with rs3757713 (68 kb upstream of GPR22), which is associated with GPR22 expression levels in lymphoblast cell lines (P = 4 x 10(-12)). Immunohistochemistry experiments revealed that G protein-coupled receptor protein 22 (GPR22) was absent in normal mouse articular cartilage or synovium. However, GPR22-positive chondrocytes were found in the upper layers of the articular cartilage of mouse knee joints that were challenged with in vivo papain treatment or methylated bovine serum albumin treatment. GPR22-positive chondrocyte-like cells were also found in osteophytes in instability-induced OA. CONCLUSION: Our findings identify a novel common variant on chromosome 7q22 that influences susceptibility to prevalence and progression of OA. Since the GPR22 gene encodes a G protein-coupled receptor, this is potentially an interesting therapeutic target.

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.

Low dietary riboflavin but not folate predicts increased fracture risk in postmenopausal women homozygous for the MTHFR 677 T allele.

UNLABELLED: The MTHFR C677T polymorphism is associated with mildly elevated homocysteine levels when folate and/or riboflavin status is low. Furthermore, a mildly elevated homocysteine level is a risk factor for osteoporotic fractures. We studied whether dietary intake of riboflavin and folate modifies the effects of the MTHFR C677T variant on fracture risk in 5,035 men and women from the Rotterdam Study. We found that the MTHFR C677T variant interacts with dietary riboflavin intake to influence fracture risk in women. INTRODUCTION: The MTHFR C677T polymorphism is associated with mildly elevated homocysteine (Hcy) levels in the presence of low folate and/or riboflavin status. A mildly elevated Hcy level was recently identified as a modifiable risk factor for osteoporotic fracture. We studied whether dietary intake of riboflavin and folate modifies the effects of the MTHFR C677T polymorphism on BMD and fracture risk. MATERIALS AND METHODS: We studied 5,035 individuals from the Rotterdam Study, >or=55 yr of age, who had data available on MTHFR, nutrient intake, and fracture risk. We performed analysis on Hcy levels in a total of 666 individuals, whereas BMD data were present for 4,646 individuals (2,692 women). RESULTS: In the total population, neither the MTHFR C677T polymorphism nor low riboflavin intake was associated with fracture risk and BMD. However, in the lowest quartile of riboflavin intake, female 677-T homozygotes had a 1.8 (95% CI: 1.1-2.9, p = 0.01) times higher risk for incident osteoporotic fractures and a 2.6 (95% CI: 1.3-5.1, p = 0.01) times higher risk for fragility fractures compared with the 677-CC genotype (interaction, p = 0.0002). This effect was not seen for baseline BMD in both men and women. No significant influence was found for dietary folate intake on the association between the MTHFR C677T genotype and fracture risk or BMD. In the lowest quartile of dietary riboflavin intake, T-homozygous individuals (men and women combined) had higher (22.5%) Hcy levels compared with C-homozygotes (mean difference = 3.44 microM, p = 0. 01; trend, p = 0.02). CONCLUSIONS: In this cohort of elderly whites, the MTHFR C677T variant interacts with dietary riboflavin intake to influence fracture risk in women.

The catechol-O-methyltransferase Met158 low-activity allele and association with nonvertebral fracture risk in elderly men.

CONTEXT: Because sex steroids play an important role in bone development, variants in genes encoding proteins involved in estrogen synthesis and metabolism could contribute to interindividual variation in bone parameters and fracture risk. An example is catechol-O-methyltransferase (COMT), an estrogen-degrading enzyme involved in inactivation of catechol-estrogens. Its gene contains a functional valine to methionine substitution at codon 158. OBJECTIVE: The aim of our study was to determine whether this polymorphism is associated with bone parameters and fracture risk in elderly subjects. METHODS: COMT genotypes were determined using TaqMan allelic discrimination in 2515 men and 3554 women from the Rotterdam Study, a population-based cohort study of individuals aged 55 and older. Associations with bone mineral density (BMD) and bone loss were analyzed using ANOVA or analysis of covariance, whereas fracture risk was analyzed using Cox's proportional hazard regression analysis. COMT mRNA expression in three osteoblastic cell lines (SaOS, MG63, and SVHFO) was analyzed by RT-PCR. RESULTS: Male carriers of the Met(158) allele had an increased risk for osteoporotic fractures (hazard ratio = 1.6; 95% confidence interval, 1.0-2.4) and for fragility fractures (hazard ratio = 2.7; 95% confidence interval, 1.3-5.9), with evidence for a dominant effect. Adjustments for age, height, weight, and BMD did not change the risk estimates. In women, this association was weaker and not significant. BMD was not significantly associated with the variant in either men or women. COMT mRNA was expressed in all three osteoblastic cell lines tested. CONCLUSION: The COMT Val158Met polymorphism is associated with fracture risk in elderly men, through a mechanism independent of BMD.

Proteomic analysis of proteins regulated by TRPS1 transcription factor in DU145 prostate cancer cells.

The aim of the present study was to identify proteins differentially regulated by TRPS1 in human prostate cancer cells in order to better understand the role of TRPS1 in prostate cancer development. The proteomes of androgen-independent DU145 prostate cancer cells, that do not express TRPS1 and of genetically engineered DU145 cells that stable and inducible express recombinant TRPS1 protein, were compared. Using two-dimensional electrophoresis followed by mass spectrometric analysis, 13 proteins that were differentially expressed between these two cell lines were identified. These proteins represent a dominant reduction of expression of antioxidant proteins, including superoxide dismutase, protein disulfide isomerase A3 precursor, endoplasmin precursor and annexin A2. Furthermore, regulation was observed for mitochondrion-associated proteins, glycolytic enzymes, a cytoskeleton-associated protein, a nuclear protein and proteins involved in apoptosis. Our data indicate that overexpression of TRPS1 protein is correlated with reduced protein expression of certain antioxidants. This suggests a possible involvement of TRPS1 in oxidative stress, and possibly in apoptosis in androgen-independent DU145 prostate cancer cells.

Common genetic variation of the low-density lipoprotein receptor-related protein 5 and 6 genes determines fracture risk in elderly white men.

UNLABELLED: Both LRP5 and LRP6 genes have been implicated to play a role in bone metabolism. In a large population-based study, we related common variation in both genes to bone parameters and fractures. LRP5 variation was associated to both BMD and frame size, whereas both LRP5 and 6 variations were associated with an increased fracture risk in males. INTRODUCTION: The low-density lipoprotein receptor-related protein 5 (LRP5) gene has a clear role in rare BMD traits and also in normal variation in peak BMD. We examined whether common variation in LRP5 and its close homolog, LRP6, plays a role in BMD in old age and fractures, the main clinical endpoint of osteoporosis. MATERIALS AND METHODS: We analyzed four variants of LRP5 and one amino acid variant of the LRP6 gene in a large prospective population-based cohort study of elderly subjects. RESULTS AND CONCLUSIONS: In men, the LRP5 1330-valine variant was associated with decreased BMD at the lumbar spine and the femoral neck with evidence for an allele-dose effect (p = 0.001 and 0.01, respectively). The Val allele was also associated with decreased vertebral body size and femoral neck width. Haplotype analysis of studied polymorphisms did not improve the association found and suggested that the 1330 variant was driving the association. We observed a borderline significant association of the LRP6 Ile1062Val polymorphism with height and vertebral body size in males. Male carriers of the LRP5 1330-valine variant had a 60% increased risk for fragility fractures, and the LRP6 1062-valine allele also conferred a 60% higher risk. Carriers of both the risk alleles of LRP5 and 6 had a 140% (p = 0.004) higher risk compared with noncarriers of both risk alleles and accounted for 10% of the fractures in males. The fracture risks were independent of age, height, weight, and BMD. In women, all of these associations were weaker and less consistent compared with men. The polymorphisms that were found associated were both situated in potentially important domains of the receptor and show considerable evolutionary conservation, which is evidence for functional importance of these residues.

Promoter and 3'-untranslated-region haplotypes in the vitamin d receptor gene predispose to osteoporotic fracture: the rotterdam study.

Polymorphisms of the vitamin D receptor gene (VDR) have been shown to be associated with several complex diseases, including osteoporosis, but the mechanisms are unknown and study results have been inconsistent. We therefore determined sequence variation across the major relevant parts of VDR, including construction of linkage disequilibrium blocks and identification of haplotype alleles. We analyzed 15 haplotype-tagging SNPs in relation to 937 clinical fractures recorded in 6,148 elderly whites over a follow-up period of 7.4 years. Haplotype alleles of the 5' 1a/1e, 1b promoter region and of the 3' untranslated region (UTR) were strongly associated with increased fracture risk. For the 16% of subjects who had risk genotypes at both regions, their risk increased 48% for clinical fractures (P = .0002), independent of age, sex, height, weight, and bone mineral density. The population-attributable risk varied between 1% and 12% for each block and was 4% for the combined VDR risk genotypes. Functional analysis of the variants demonstrated 53% lower expression of a reporter construct with the 1e/1a promoter risk haplotype (P = 5 x 10(-7)) in two cell lines and 15% lower mRNA level of VDR expression constructs carrying 3'-UTR risk haplotype 1 in five cell lines (P = 2 x 10(-6)). In a further analysis, we showed 30% increased mRNA decay in an osteoblast cell line for the construct carrying the 3'-UTR risk haplotype (P = .02). This comprehensive candidate-gene analysis demonstrates that the risk allele of multiple VDR polymorphisms results in lower VDR mRNA levels. This could impact the vitamin D signaling efficiency and might contribute to the increased fracture risk we observed for these risk haplotype alleles.

The TRPS1 transcription factor: androgenic regulation in prostate cancer and high expression in breast cancer.

TRPS1 mRNA is more highly expressed in androgen-dependent lymph node carcinoma of prostate-fast growing colony (LNCaP-FGC) compared with androgen-independent lymph node carcinoma of prostate-lymph node original (LNCaP-LNO) prostate cancer cell lines. Furthermore, TRPS1 mRNA expression is down-regulated by androgens in LNCaP-FGC cells, a process mediated by the androgen receptor (AR). Here, we present TRPS1 protein expression in human prostate cancer material derived from a panel of six androgen-dependent and eight androgen-independent human prostate cancer xenografts. TRPS1 protein is expressed in all androgen-dependent xenografts, which also express AR and prostate-specific antigen (PSA). Androgen withdrawal by castration resulted in an increase in TRPS1 protein in two androgen-dependent xenografts, indicating relieved repression by action of AR. TRPS1 protein is expressed in four androgen-independent xenografts and is low or absent in the other four androgen-independent xenografts. Androgen withdrawal by castration demonstrates that TRPS1 protein levels remain the same in 1 androgen-independent xenograft, most likely due to the lack of AR expression. These data show that TRPS1 protein expression is regulated by androgens via the AR in human prostate cancer xenografts. Analysis of TRPS1 mRNA expression in normal and tumour tissue of the prostate and 18 other human tissues, showed that TRPS1 had the highest mRNA expression levels in normal and tumour tissues of breast. In addition, high TRPS1 mRNA and protein expression levels were observed in four out of five human breast cancer cell lines. In conclusion, TRPS1 protein expression is down-regulated by androgens in human prostate cancer, and analysis of TRPS1 mRNA expression levels in several human tissues showed that the highest levels were observed in normal and tumour breast tissue.

The atypical GATA protein TRPS1 represses androgen-induced prostate-specific antigen expression in LNCaP prostate cancer cells.

Prostate-specific antigen (PSA) is considered as an important marker for prostate cancer. Regulation of PSA gene expression is mediated by androgens bound to androgen receptors via androgen response elements (AREs) in its promoter and far upstream enhancer regions. In addition, GATA proteins contribute to PSA gene transcription by interacting with GATA motifs present in the PSA enhancer sequence. The TRPS1 gene contains a single GATA zinc finger domain and not only binds to forward consensus GATA motifs but also to an inverse GATA motif overlapping the ARE III in the far upstream enhancer of the PSA gene. Overexpression of TRPS1 in androgen-dependent human LNCaP prostate cancer cells inhibited the transcription of a transiently transfected PSA enhancer/promoter-driven luciferase reporter construct. Furthermore, overexpression of TRPS1 reduced the androgen-induced endogenous PSA levels secreted in culture medium of LNCaP cells. Our results suggest a role of TRPS1 in androgen regulation of PSA gene expression.

A central dinucleotide within vitamin D response elements modulates DNA binding and transactivation by the vitamin D receptor in cellular response to natural and synthetic ligands.

There is considerable divergence in the sequences of steroid receptor response elements, including the vitamin D response elements (VDREs). Two major VDRE-containing and thus 1,25-dihydroxyvitamin D(3) (1,25-(OH)(2)D(3))-regulated genes are the two non-collagenous, osteoblast-derived bone matrix proteins osteocalcin and osteopontin. We observed a stronger induction of osteopontin than osteocalcin mRNA expression by 1,25-(OH)(2)D(3). Subsequently, we have shown that vitamin D receptor/retinoid X receptor alpha (VDR/RXRalpha) heterodimers bind more tightly to the osteopontin VDRE than to the osteocalcin VDRE. Studies using point mutants revealed that the internal dinucleotide at positions 3 and 4 of the proximal steroid half-element are most important for modulating the strength of receptor binding. In addition, studies with VDRE-driven luciferase reporter gene constructs revealed that the central dinucleotide influences the transactivation potential of VDR/RXRalpha with the same order of magnitude as that observed in the DNA binding studies. The synthetic vitamin D analog KH1060 is a more potent stimulator of transcription and inducer of VDRE binding of VDR/RXR in the presence of nuclear factors isolated from ROS 17/2.8 osteoblast-like cells than the natural ligand 1,25-(OH)(2)D(3). Interestingly, however, KH1060 is comparable or even less potent than 1,25-(OH)(2)D(3) in stimulating VDRE binding of in vitro synthesized VDR/RXRalpha. Thus, the extent of 1,25-(OH)(2)D(3)- and KH1060-dependent binding of VDR/RXRalpha is specified by a central dinucleotide in the VDRE, and the ligand-induced effects on DNA binding are in part controlled by the cellular context of nuclear proteins.