Haplotype structure, LD blocks, and uneven recombination within the LRP5 gene

The role of Wnt pathway in obesity induced inflammation and diabetes: a review

Diabetes has become a major killer worldwide and at present, millions are affected by it. Being a chronic disease it increases the risk of other diseases ranging from pulmonary disorders to soft tissue infections. The loss of insulin-producing capacity of the pancreatic β-cells is the main reason for the development of the disease. Obesity is a major complication that can give rise to several other diseases such as cancer, diabetes, etc. Visceral adiposity is one of the major factors that play a role in the development of insulin resistance. Obesity causes a chronic low-grade inflammation in the tissues that further increases the chances of developing diabetes. Several pathways have been associated with the development of diabetes due to inflammation caused by obesity. The Wnt pathway is one such candidate pathway that is found to have a controlling effect on the development of insulin resistance. Moreover, the pathway has also been linked to obesity and inflammation. This review aims to find a connection between obesity, inflammation, and diabetes by taking the wnt pathway as the connecting link.

LDkit: a parallel computing toolkit for linkage disequilibrium analysis

BACKGROUND

Linkage disequilibrium (LD) analysis is broadly utilized in genetics to understand the evolutionary and demographic history and helps geneticists identify genes associated with interested inherited traits, such as diseases. There are some tools for linkage disequilibrium analysis either in a local or online way; however, there has been no such tool supporting both graphical user interface (GUI) and parallel computing.

RESULTS

We developed a GUI software called LDkit for LD analysis, which supports parallel computing. The LDkit supports both variant call format (VCF) and PLINK 'ped + map' format. At the same time, users could also just analyze a subset of individuals from the whole population. The LDkit reads the data by block and then paralleled the computation process by monitoring the usage of processes. Assessment on the Human 1000 genome data showed that when paralleled with 32 threads, the running time was reduced to less than 6 minutes from ~77 minutes using the chromosome 22 dataset with 1,103,547 SNPs and 2504 individuals.

CONCLUSIONS

The software LDkit can be effectively used to calculate and plot LD decay, LD block, and linkage disequilibrium analysis between a site and a given region. Most importantly, both graphical user interface (GUI) and stand-alone packages are available for users' convenience. LDkit was written in JAVA language under cross-platform support.

Risk of Wnt/β-catenin signalling pathway gene polymorphisms in primary Sjögren's syndrome

OBJECTIVE

To explore genetic polymorphisms of the Wnt/β-catenin signalling pathway in primary SS (PSS).

METHODS

We included 98 patients with PSS and 165 healthy volunteers. Genomic DNA was extracted from peripheral blood samples. Through an open-array platform of low density, we genotyped 25 polymorphisms from 14 genes (WISP1, DKK1, SOST, FRZB, LRP1, LRP4, LRP5, LRP6, GSKB, ADAMTS5, GDF5, FMN2, ADIPOQ and COL11A1) involved in the Wnt/β-catenin signalling pathway. We compared the allelic and genotypic frequencies with Fisher's exact test and logistic regression analysis adjusted by age, gender and individual admixture, as well as bootstrap-resampling analysis. We assessed the gene-gene interaction by the multifactor dimensionality reduction method.

RESULTS

We found a positive significant association with four polymorphisms: LRP5 rs606989, FRZB rs409238, GSK3B rs2037547 and ADIPOQ rs2241766. All of them conferred risk for PSS, being the highest among subjects carrying three to four risk alleles (P < 0.001). According to a multifactor dimensionality reduction analysis, the best models included the LRP5 (rs606989), FRZB (rs409238) and ADIPOQ (rs2241766) polymorphisms.

CONCLUSION

LRP5, FRZB and ADIPOQ genes related in the Wnt/β-catenin signalling pathway increased the risk of PSS. Further research is needed to establish their functional role in this clinical entity.

WNT Signaling in Disease

Developmental signaling pathways control a vast array of biological processes during embryogenesis and in adult life. The WNT pathway was discovered simultaneously in cancer and development. Recent advances have expanded the role of WNT to a wide range of pathologies in humans. Here, we discuss the WNT pathway and its role in human disease and some of the advances in WNT-related treatments.

Genome-Wide Association Studies for Pasmo Resistance in Flax (Linum usitatissimum L.)

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of QTL with positive-effect or favorite alleles (NPQTL) in accessions was significantly correlated with PR (R 2 = 0.55), suggesting that these QTL effects are mainly additive. NPQTL was also significantly associated with morphotype (R 2 = 0.52) and major QTL with positive effect alleles were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

The low-density lipoprotein receptor-related protein 5 (LRP5) 4037C>T polymorphism: candidate for susceptibility to type 1 diabetes mellitus

OBJECTIVE

The present study has investigated the association between low-density lipoprotein receptor-related protein 5 (LRP5) 4037C>T polymorphism and type 1 diabetes mellitus (T1DM) susceptibility in a Brazilian population.

SUBJECTS AND METHODS

A total number of 134 T1DM patients and 180 normoglycemic individuals (NG) aged 6-20 years were studied. Glycated hemoglobin and glucose levels were determined. Genotyping of LRP5 4037C>T (rs3736228) was performed.

RESULTS

T1DM patients showed poor glycemic control. Genotypes in the codominant (CT: OR = 2.99 [CI 95%: 1.71-5.24], p < 0.001; TT: OR = 5.34 [CI 95%: 1.05-2702], p < 0.001), dominant (CT + TT: OR = 3.16 [CI 95%: 1.84-5.43], p < 0.001) and log-additive (OR = 2.78 [CI 95%: 1.70-4.52], p < 0.001) models, and LRP5 4037T allele (OR = 2.88, [CI 95%: 1.78-4.77], p < 0.001) were associated with an increased risk of developing T1DM. LRP5 4037CT and CT+TT carriers in T1DM group showed higher concentrations of serum glucose and glycated hemoglobin when compared with CC carriers.

CONCLUSION

The LRP5 4037C>T may represent a candidate for T1DM susceptibility, as well as poor glycemic control.

A new haplotype block detection method for dense genome sequencing data based on interval graph modeling of clusters of highly correlated SNPs

Motivation

Linkage disequilibrium (LD) block construction is required for research in population genetics and genetic epidemiology, including specification of sets of single nucleotide polymorphisms (SNPs) for analysis of multi-SNP based association and identification of haplotype blocks in high density sequencing data. Existing methods based on a narrow sense definition do not allow intermediate regions of low LD between strongly associated SNP pairs and tend to split high density SNP data into small blocks having high between-block correlation.

Results

We present Big-LD, a block partition method based on interval graph modeling of LD bins which are clusters of strong pairwise LD SNPs, not necessarily physically consecutive. Big-LD uses an agglomerative approach that starts by identifying small communities of SNPs, i.e. the SNPs in each LD bin region, and proceeds by merging these communities. We determine the number of blocks using a method to find maximum-weight independent set. Big-LD produces larger LD blocks compared to existing methods such as MATILDE, Haploview, MIG ++, or S-MIG ++ and the LD blocks better agree with recombination hotspot locations determined by sperm-typing experiments. The observed average runtime of Big-LD for 13 288 240 non-monomorphic SNPs from 1000 Genomes Project autosome data (286 East Asians) is about 5.83 h, which is a significant improvement over the existing methods.

Availability and implementation

Source code and documentation are available for download at http://github.com/sunnyeesl/BigLD.

Contact

yyoo@snu.ac.kr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Effects of Single Nucleotide Polymorphism Marker Density on Haplotype Block Partition

Many researchers have found that one of the most important characteristics of the structure of linkage disequilibrium is that the human genome can be divided into non-overlapping block partitions in which only a small number of haplotypes are observed. The location and distribution of haplotype blocks can be seen as a population property influenced by population genetic events such as selection, mutation, recombination and population structure. In this study, we investigate the effects of the density of markers relative to the full set of all polymorphisms in the region on the results of haplotype partitioning for five popular haplotype block partition methods: three methods in Haploview (confidence interval, four gamete test, and solid spine), MIG++ implemented in PLINK 1.9 and S-MIG++. We used several experimental datasets obtained by sampling subsets of single nucleotide polymorphism (SNP) markers of chromosome 22 region in the 1000 Genomes Project data and also the HapMap phase 3 data to compare the results of haplotype block partitions by five methods. With decreasing sampling ratio down to 20% of the original SNP markers, the total number of haplotype blocks decreases and the length of haplotype blocks increases for all algorithms. When we examined the marker-independence of the haplotype block locations constructed from the datasets of different density, the results using below 50% of the entire SNP markers were very different from the results using the entire SNP markers. We conclude that the haplotype block construction results should be used and interpreted carefully depending on the selection of markers and the purpose of the study.

Adipose- and muscle-derived Wnts trigger pancreatic β-cell adaptation to systemic insulin resistance

Wnt signaling molecules are associated with obesity, hyperlipidemia, and type 2 diabetes (T2D). Here, we show that two Wnt proteins, WNT3a and WNT4, are specifically secreted by skeletal muscle and adipose tissue during the development of insulin resistance and play an important role in cross-talk between insulin-resistant tissues and pancreatic beta cells. The activation of Frizzled receptor and Wnt signaling in pancreatic islets via circulating WNT3a in blood resulted in higher insulin secretion and an increase in beta cell proliferation, thus leading to islet adaptation in a pre-diabetic state. Interestingly, in fully developed T2D, the expression profiles of Wnt3a and Wnt4 in adipose tissue and muscle cells and blood plasma levels of these proteins were opposite to the pre-diabetic state, thus favoring the downregulation of Wnt signaling in beta cells and resulting in dysfunctional pancreatic islets. These results demonstrate that alterations in the secretion profile of a canonical Wnt activator (WNT3a) and inhibitor (WNT4) from insulin-resistant tissues during the development of T2D are responsible for triggering progression from a pre-diabetic to a diabetic state. We also show here that WNT3a and WNT4 are potent myokines, and their expression and secretion are regulated in response to nutritional and metabolic changes.

Genetic analysis of the relationship between bone mineral density and low-density lipoprotein receptor-related protein 5 gene polymorphisms

Background

A number of studies have examined the association between the polymorphisms of the low-density lipoprotein receptor-related protein 5 gene (LRP5), but previous results have been inconclusive. Thus we performed a meta-analysis of studies on the association between the LRP5 polymorphisms and bone mineral density (BMD) to assess their pooled effects.

Methods

Published literature from PubMed, EMBASE and ISI web of science were searched for eligible publications. Weighted mean difference (WMD) and 95% confidence interval (CI) was calculated using fixed- or random-effects model.

Results

A total of 19 studies with 25773 subjects were considered in this meta-analysis. Of them, 17 examined the association between the A1330V polymorphism and BMD, 8 were focused on the V667M polymorphism, and 2 analyzed the Q89R polymorphism. Individuals with the A1330V AA genotype showed significantly higher BMD than those with the AV/VV genotypes [at lumbar spine (LS): WMD = 0.02g/cm², 95% CI = 0.01-0.03, P < 10^-4; at femur neck (FN): WMD = 0.01g/cm², 95% CI = 0.00-0.02, P = 0.01] or VV genotype (at LS: WMD = 0.02g/cm², 95% CI = 0.01-0.04, P = 0.01). Significant associations were also detected in the analysis for V667M (VV vs. VM/MM: WMD at LS = 0.02g/cm², 95% CI = 0.02-0.03, P < 10^-5; WMD at FN = 0.01g/cm², 95% CI = 0.01-0.02, P = 0.0002). As for Q89R, subjects with the QQ genotype tended to have higher BMD than those with the QR/RR genotypes at FN (WMD = 0.03g/cm², 95% CI = 0.01-0.05, P = 0.005).

Conclusion

This meta-analysis demonstrated that the LRP5 polymorphisms may be modestly associated with BMD of LS and FN.

Haplotype kernel association test as a powerful method to identify chromosomal regions harboring uncommon causal variants

For most complex diseases, the fraction of heritability that can be explained by the variants discovered from genome-wide association studies is minor. Although the so-called "rare variants" (minor allele frequency [MAF] < 1%) have attracted increasing attention, they are unlikely to account for much of the "missing heritability" because very few people may carry these rare variants. The genetic variants that are likely to fill in the "missing heritability" include uncommon causal variants (MAF < 5%), which are generally untyped in association studies using tagging single-nucleotide polymorphisms (SNPs) or commercial SNP arrays. Developing powerful statistical methods can help to identify chromosomal regions harboring uncommon causal variants, while bypassing the genome-wide or exome-wide next-generation sequencing. In this work, we propose a haplotype kernel association test (HKAT) that is equivalent to testing the variance component of random effects for distinct haplotypes. With an appropriate weighting scheme given to haplotypes, we can further enhance the ability of HKAT to detect uncommon causal variants. With scenarios simulated according to the population genetics theory, HKAT is shown to be a powerful method for detecting chromosomal regions harboring uncommon causal variants.

The involvement of the wnt signaling pathway and TCF7L2 in diabetes mellitus: The current understanding, dispute, and perspective

The Wnt signaling pathway was initially discovered for its role in tumorigenesis and the development of Drosophila and other eukaryotic organisms. The key effector of this pathway, the bipartite transcription factor β-cat/TCF, is formed by free β-catenin (β-cat) and a TCF protein, including TCF7L2. Extensive recent investigations have highlighted the role of the Wnt signaling pathway in metabolic homeostasis and its implication in diabetes and other metabolic diseases. Genome-wide association studies have shown that several key components of the Wnt signaling pathway are implicated in metabolic homeostasis and the development of type 2 diabetes (T2D). Despite controversial observations regarding the role of Wnt signaling in the development and function of pancreatic islets, the discovery of the association between certain single nucleotide polymorphisms of TCF7L2 and T2D susceptibility has fueled great efforts to explore the role of Wnt signaling in the function of pancreatic β-cells and glucose homeostasis. Here we have introduced our basic understanding of the canonical Wnt signaling pathway, summarized our current knowledge on its implication in metabolic homeostasis and T2D, discussed the work on TCF7L2 as a T2D susceptibility gene, and presented the controversial role of Wnt signaling and TCF7L2 in pancreatic islets as well as their potential metabolic function in other organs. We then expanded our view into the crosstalk among Wnt, insulin and FOXO signaling cascades, which further illustrates the complexity of the Wnt signaling pathway in metabolic homeostasis. Finally, we have presented our perspectives.

Similarity-based multimarker association tests for continuous traits

Testing multiple markers simultaneously not only can capture the linkage disequilibrium patterns but also can decrease the number of tests and thus alleviate the multiple-testing penalty. If a gene is associated with a phenotype, subjects with similar genotypes in this gene should also have similar phenotypes. Based on this concept, we have developed a general framework that is applicable to continuous traits. Two similarity-based tests (namely, SIMc and SIMp tests) were derived as special cases of the general framework. In our simulation study, we compared the power of the two tests with that of the single-marker analysis, a standard haplotype regression, and a popular and powerful kernel machine regression. Our SIMc test outperforms other tests when the average R(2) (a measure of linkage disequilibrium) between the causal variant and the surrounding markers is larger than 0.3 or when the causal allele is common (say, frequency = 0.3). Our SIMp test outperforms other tests when the causal variant was introduced at common haplotypes (the maximum frequency of risk haplotypes >0.4). We also applied our two tests to an adiposity data set to show their utility.

Low‐density lipoprotein receptor‐related protein 5 variant Q89R is associated with hypertension in Japanese females

We examined the influence of the Q89R variant in low-density lipoprotein receptor-related protein 5 on hypertension in a large Japanese cohort. We used multivariate logistic regression analysis to adjust for the effects of other factors known to influence hypertension such as age, body mass index, blood chemistry and lifestyle. Our cohort consisted of 1444 males and 1161 females selected from 3834 male and 2591 female workers in a single company. Hypertension was defined as systolic blood pressure 140 mmHg and/or diastolic blood pressure 90 mmHg or the use of antihypertensive medication. Because the RR genotype was so rare (three normotensive males and three normotensive females), they were excluded from analyses. Genotype distributions for Q89R in hypertensive females (QR = 14, QQ = 60) were significantly different (p = 0.033) from normotensive females (QR = 113, QQ = 971). In contrast, those in hypertensive males (QR = 26, QQ = 230) and in normotensive males (QR = 145, QQ = 1040) were similar. Allele distributions were not significantly different in either gender. In females, multivariate logistic regression showed that the QR genotype was associated with hypertension with odds ratio of 2.1 compared to the QQ genotype. This study indicates that the Q89R polymorphism is an independent factor for hypertension in Japanese females.

Analysis of meiotic recombination products from human sperm

Traditional methods for surveying meiotic recombination in humans are limited to pedigree and linkage disequilibrium analyses. We have developed assays that allow the direct detection of crossover and gene conversion molecules in batches of sperm DNA. To date, we have characterized 26 recombination hotspots by allele-specific PCR and selectively amplified recombinant DNA molecules from these regions. These analyses have revealed that meiotic crossover hotspots in humans are highly localized and flanked by DNA segments where recombination is suppressed. The centers of crossover hotspots are also active in noncrossover recombination, displaying short conversion tracts.

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.

The WNT signalling pathway and diabetes mellitus

The WNT signalling pathway is involved in many physiological and pathophysiological activities. WNT ligands bind to Frizzled receptors and co-receptors (LDL receptor-related protein 5/6), triggering a cascade of signalling events. The major effector of the canonical WNT signalling pathway is the bipartite transcription factor beta-catenin/T cell transcription factor (beta-cat/TCF), formed by free beta-cat and one of the four TCFs. The WNT pathway is involved in lipid metabolism and glucose homeostasis, and mutations in LRP5 may lead to the development of diabetes and obesity. beta-Cat/TCF is also involved in the production of the incretin hormone glucagon-like peptide-1 in the intestinal endocrine L cells. More recently, genome-wide association studies have identified TCF7L2 as a diabetes susceptibility gene, and individuals carrying certain TCF7L2 single nucleotide polymorphisms could be more susceptible to the development of type 2 diabetes. Furthermore, beta-cat is able to interact with forkhead box transcription factor subgroup O (FOXO) proteins. Since FOXO and TCF proteins compete for a limited pool of beta-cat, enhanced FOXO activity during ageing and oxidative stress may attenuate WNT-mediated activities. These observations shed new light on the pathogenesis of type 2 diabetes as an age-dependent disease.

Large-scale association study between two coding LRP5 gene polymorphisms and bone phenotypes and fractures in men

Herein we investigated the association between polymorphisms in the LRP5 gene and bone phenotypes and fractures in three large male cohorts based on the rationale that mutations in LRP5 cause severe bone phenotypes. Results showed an association of the Val667Met SNP with spine BMD in 3,800 young and elderly men.

INTRODUCTION

The low-density lipoprotein receptor-related protein 5 (LRP5)-Wnt signalling system is of importance for regulating osteoblastic activity, which became clear after findings that inactivating mutations in LRP5 cause osteoporosis. The overall aim of this study was to investigate the association between polymorphisms in the LRP5 gene and bone mineral density (BMD) in three large cohorts of young and elderly men.

METHODS

The cohorts used were MrOS Sweden (n = 3014, aged 69-81 years) and MrOs Hong Kong (n = 2000, aged > 65 years) and the Swedish GOOD study (n = 1068, aged 18-20 years). The polymorphisms Val667Met and Ala1330Val were genotyped using a TaqMan assay.

RESULTS

When combining the data from the Swedish cohorts in a meta-analysis (n = 3,800), men carrying the 667Met-allele had 3% lower BMD at lumbar spine compared with non-carriers (p < 0.05). The Val667Met SNP was not polymorphic in the Hong Kong population and thus were not included. There were no associations between the Ala1330Val SNP and bone phenotypes in the study populations. No associations between the LRP5 polymorphisms and self-reported fractures were seen in MrOs Sweden.

CONCLUSIONS

Results from these three large cohorts indicate that the Val667Met polymorphism but not the Ala1330Val contributes to the observed variability in BMD in the Swedish populations.

Mammalian meiotic recombination hot spots

Our understanding of the details of mammalian meiotic recombination has recently advanced significantly. Sperm typing technologies, linkage studies, and computational inferences from population genetic data have together provided information in unprecedented detail about the location and activity of the sites of crossing-over in mice and humans. The results show that the vast majority of meiotic recombination events are localized to narrow DNA regions (hot spots) that constitute only a small fraction of the genome. The data also suggest that the molecular basis of hot spot activity is unlikely to be strictly determined by specific DNA sequence motifs in cis. Further molecular studies are needed to understand how hot spots originate, function and evolve.

Association of low-density lipoprotein receptor-related protein 5 (LRP5) promoter SNP with peak bone mineral density in Chinese women

OBJECTIVE

Low-density lipoprotein receptor-related protein 5 (LRP5) is important for osteoblast differentiation and mutations of the gene are associated with both low and high bone mass syndromes. Our study aimed to evaluate the importance of LRP5 in the determination of peak bone mass acquisition in Chinese females in the general population.

METHODS

A total of 286 young southern Chinese females (aged 22-44 years) with low bone mineral density (BMD) (defined by a BMD Z score < or =-1.28 at either the hip or spine) or high BMD (Z score > or =+1) were studied. The LRP5 gene was sequenced for single nucleotide polymorphisms (SNPs) and 4 SNPs were tagged from 8 genotyped SNPs for this study.

RESULTS

Single locus allele association tests revealed significant associations of rs682429 and rs686921 with BMD variation (p < 0.05). Omnibus test (likelihood ratio test) revealed overall significant association between LRP5 gene locus and total hip BMD, with rs682429 being most predictive. rs682429 is located in 5'UTR, 2 bases adjacent to a consensus recognition site for the Elk-1 binding element.

CONCLUSION

Common variations of the LRP5 promoter are associated with BMD in young women. These significant associations appear to be driven by rs682429. Functional studies are necessary to elucidate the role of this SNP on the effect of Elk-1 binding element transcriptional activity of LRP5 gene.

Multiple less common genetic variants explain the association of the cholesteryl ester transfer protein gene with coronary artery disease

OBJECTIVES

The objective of this study was to identify associations of the cholesteryl ester transfer protein (CETP) gene with coronary artery disease (CAD) with tagging (t) single nucleotide polymorphisms (SNPs) chosen to optimally account for intra-genic variation.

BACKGROUND

The CETP gene plays a critical role in lipoprotein metabolism, but the common and well-studied TaqIB variant is inconsistently predictive of CAD.

METHODS

From a deoxyribonucleic acid bank of 10,020 individuals, nondiabetic nonsmoking patients (n = 4,811) with angiographically defined, clinically significant CAD (> or =70% stenosis) or normal coronaries were genotyped for 11 CETP tSNPs. Myocardial infarction (MI) and lipid levels were evaluated as secondary end points.

RESULTS

Analysis of single tSNPs, corrected for multiple comparisons (p < 0.00485), identified allele +1086A to be associated with CAD (p = 0.0034). Suggestive allelic and significant genotypic associations were found for -631AA (odds ratio [OR] = 3.95, p = 0.004 vs. CC) and +2389GA (OR = 1.21, p = 0.003 vs. GG). Haplotype analysis by linkage disequilibrium (LD) group revealed a CAD association for LD group B (p = 0.0025 across T+1086A, C+878T, C+408T) and near significance for LD group A (p = 0.013 across C-631A, MspI, G+2389A). A weak protective trend for TaqIB was eliminated by adjustment for other tSNPs, and haplotype analyses suggested that TaqIB was simply a marker for other tSNPs or haplotypes. No tSNP or haplotype associations with MI were found.

CONCLUSIONS

Multiple, less common SNPs and haplotype variants underlie CETP-related CAD risk, for which the common TaqIB variant is simply a poor marker. The occurrence of risk-related variants on separate haplotypes suggests genetic-risk complexity and allelic heterogeneity. (Database Registry of the Intermountain Heart Collaborative Study; http://clinicaltrials.gov/ct/show/NCT00406185?order=1; NCT00406185).

A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle

Analysis of data on 1000 Holstein-Friesian bulls genotyped for 15,036 single-nucleotide polymorphisms (SNPs) has enabled genomewide identification of haplotype blocks and tag SNPs. A final subset of 9195 SNPs in Hardy-Weinberg equilibrium and mapped on autosomes on the bovine sequence assembly (release Btau 3.1) was used in this study. The average intermarker spacing was 251.8 kb. The average minor allele frequency (MAF) was 0.29 (0.05-0.5). Following recent precedents in human HapMap studies, a haplotype block was defined where 95% of combinations of SNPs within a region are in very high linkage disequilibrium. A total of 727 haplotype blocks consisting of > or =3 SNPs were identified. The average block length was 69.7 +/- 7.7 kb, which is approximately 5-10 times larger than in humans. These blocks comprised a total of 2964 SNPs and covered 50,638 kb of the sequence map, which constitutes 2.18% of the length of all autosomes. A set of tag SNPs, which will be useful for further fine-mapping studies, has been identified. Overall, the results suggest that as many as 75,000-100,000 tag SNPs would be needed to track all important haplotype blocks in the bovine genome. This would require approximately 250,000 SNPs in the discovery phase.

Low-density lipoprotein receptor-related protein 5 (LRP5) gene polymorphisms are associated with bone mass in both Chinese and whites

In this study, the associations of novel LRP5 variants with BMD variation were detected and some replicated in the two ethnic groups of Chinese and white origins, respectively. These data support the concept that LRP5 variation can contribute to minor and major variation in bone structure.

INTRODUCTION

Mutations in the low-density lipoprotein receptor-related protein 5 (LRP5) gene have been shown to cause both high and low bone mass. However, it is still controversial whether LRP5 is associated with normal BMD variation. This study explored the association of LRP5 with BMD phenotypes at three clinically important skeletal sites-the spine, hip, and ultradistal radius (UD)-in two independent populations of Chinese and white ethnicities, respectively.

MATERIALS AND METHODS

The Chinese sample consisted of 733 unrelated subjects. The white sample was made up of 1873 subjects from 405 nuclear families. High-density single nucleotide polymorphisms (SNPs) across the whole LRP5 gene were genotyped and analyzed in both samples.

RESULTS

Linkage disequilibrium (LD) analyses showed that the haplotype structures of LRP5 between Chinese and whites were in good agreement. Association tests showed that polymorphisms in block 5 spanning intron 7 to intron 19 of LRP5 significantly associated with spine BMD variation in both samples. Particularly, the significant association of SNP rs491347 in intron 7 with spine BMD in the Chinese sample (p=0.002) was replicated in whites, even after adjusting for multiple testing (p=0.005). Its strongly associated SNP rs1784235 could cause the loss of an estrogen receptor alpha (ERalpha) binding site in LRP5, which could partially explain the above replicated association. However, we did not observe any significant replication with BMD variation at the hip and UD. After accounting for multiple testing, associations with BMD variation at these two sites were mainly found in Chinese. Sex-stratified analyses further revealed that the LRP5 associations with BMD in Chinese and whites were driven by male and female subjects, respectively.

CONCLUSIONS

Our work supported LRP5 genetic variants as possible susceptibility factors for osteoporosis and fractures in humans. Especially, the SNP rs491347 and its strongly associated SNPs (e.g., rs1784235) could be important to human osteoporosis phenotypes.

Genetic variation at the low-density lipoprotein receptor-related protein 5 (LRP5) locus modulates Wnt signaling and the relationship of physical activity with bone mineral density in men

Polymorphisms in the LRP5 gene have been associated with bone mineral density (BMD) in men and/or women. However, the functional basis for this association remains obscure. We hypothesized that LRP5 alleles could modulate Wnt signaling and the relationship between physical activity and BMD. This genetic association study was performed in the population-based Framingham Study Offspring Cohort, and included a subset of 1797 unrelated individuals who provided blood samples for DNA and who had BMD measurements of the hip and spine. Ten single-nucleotide polymorphisms (SNPs) spanning the LRP5 gene were genotyped and used for association and interaction analyses with BMD by regression methods. LRP5 haplotypes were transiently co-expressed with Wnt3a, MesD and Dkk1 in HEK293 cells and their activity evaluated by the TCF-Lef reporter assay. Six out of ten SNPs in LRP5 were associated with one or more of the femur or spine BMDs in men or women after adjustment for covariates, and these associations differed between genders. In men< or =age 60 years, 3 SNPs were significantly associated with BMD: rs2306862 on Exon 10 with femoral neck BMD (p=0.01) and Ward's BMD (p=0.01); rs4988321/p. V667M with Ward's BMD (p=0.02); and intronic rs901825 with trochanter BMD (p=0.03). In women, 3 SNPs in intron 2 were significantly associated with BMD: rs4988330 for trochanter (p=0.01) and spine BMD (p=0.003); rs312778 with femoral neck BMD (p=0.05); and rs4988331 with spine BMD (p=0.04). For each additional rare allele, BMD changed by 3-5% in males and 2-4% in females. Moreover, there was a significant interaction between physical activity and rs2306862 in exon 10 (p for interaction=0.02) and rs3736228/p. A1330V in exon 18 (p for interaction=0.05) on spine BMD in men. In both cases, the TT genotype was associated with lower BMD in men with higher physical activity scores, conversely with higher BMD in men with lower physical activity scores. In vitro, TCF-Lef activity in presence of Wnt3a was significantly reduced in cells expressing LRP5 haplotypes carrying the T allele of exon 10 and 18 compared to the wild-type allele, whereas co-expression of Dkk1 completely inhibited Wnt3a response through all LRP5 haplotypes. In summary, genetic variation in exons 10 and 18 of the LRP5 gene modulates Wnt signaling and the relationship between physical activity and BMD in men. These observations suggest that Wnt-LRP5 may play a role in the adaptation of bone to mechanical load in humans, and may explain some gender-related differences in bone mass.

Stepwise haplotype analysis: are LD patterns repeatable?

A variety of techniques exist to describe and depict patterns of pairwise linkage disequilibrium (LD). In the current paper, a new log-linear framework is proposed for the summarisation of local interactions among single nucleotide polymorphisms (SNPs). Our approach provides a straightforward means of capturing the diversity of higher-order LD relationships for small numbers of loci by investigating inter-marker interactions. Our method was applied to a dataset of 76 SNP markers spanning a genomic interval of length 2.8 megabases. The analysis of three short sub-regions is described in detail here. Model and graphical representations of contiguous markers in medium to high LD are presented. In the regions studied, evidence for sub-structure was detected, supporting the view that the genomic reality is complex. Interestingly, a critical evaluation of the method by bootstrapping showed that while some LD relationships were captured in a highly repeatable fashion, the majority were not. Large numbers of small interactions, both direct and indirect, mean that many models can adequately summarise the data at hand. Our results suggest that repeatability should be further investigated in the application of LD-based approaches.

G-protein beta 3 subunit polymorphism C1429T and low-density lipoprotein receptor-related protein 5 polymorphism A1330V are risk factors for hypercholesterolemia in Japanese males--a prospective study over 5 years

We examined the relationship between the C825T, C1429T, and A-350G variants in the G-protein beta 3 subunit (GNB3) gene, the A1330V and Q89R variants in the low-density lipoprotein receptor-related protein 5 (LRP5) gene, and the risk of hypercholesterolemia in a prospective study in Japanese workers. This study included observations over a 5-year period from 1997 to 2002 on 936 males and 662 females who were not hypercholesterolemic on entry. Hypercholesterolemia was defined as a serum total cholesterol level of 240 mg/dL or higher. Pooled logistic regression analyses were performed using either of the gene variants with age, body mass index, smoking, alcohol consumption, and habitual exercise as the covariates. The risk of the development of hypercholesterolemia was 2.27 times higher in males with the TT genotype of GNB3/C1429T than in males with the CC genotype (95% confidence interval, 1.04-4.94), after adjustment for the effects of other potential covariates. Simultaneously, the risk was 1.49 times higher in males with the AV genotype of LRP5/A1330V than in males with the AA genotype (95% confidence interval, 1.04-2.12) after adjustment for the effects of other potential covariates. This study indicates the GNB3/C1429T and LRP5/A1330V are independent risk factors for hypercholesterolemia in Japanese males and suggests that targeting these polymorphisms may be beneficial when attempting to prevent hypercholesterolemia in the general Japanese male population.

The bovine PPARGC1A gene: molecular characterization and association of an SNP with variation of milk fat synthesis

Several studies in a variety of breeds have reported at least two QTL for milk production traits, including milk fat synthesis on bovine chromosome 6 (BTA6), comprising a region that comparatively has been mapped to equivalent syntenic chromosome intervals in human, pig, and mouse harboring loci associated with type II diabetes and obesity-related traits. We identified the bovine peroxysome proliferator-activated receptor-gamma coactivator-1alpha gene (PPARGC1A) as a plausible positional and functional candidate gene for a previously described QTL for milk fat yield on BTA6 because of its chromosomal position and its key role in energy, fat, and glucose metabolism. To analyze the role of the bovine PPARGC1A gene in regulation of milk fat synthesis in dairy cattle, we determined its cDNA sequence, genomic organization, chromosomal localization, and expression pattern. The bovine PPARGC1A gene is organized in 13 exons comprising 6,261 bp and is expressed at different levels in a large number of tissues. Bovine PPARGC1A cDNA and protein sequences showed substantial similarity (92-95%) to its respective orthologs from human, rat, and mouse. Screening for polymorphisms in the coding sequence, exon/intron boundaries, 5'- and 3'-untranslated regions, and promoter region of the PPARGC1A gene in sires with a different genotype at the QTL for milk fat yield as well as in a multibreed panel revealed a total of 11 polymorphic loci. A significant association between an SNP in intron 9 of the PPARGC1A gene and milk fat yield was observed in a major dairy cattle population, indicating that the PPARGC1A gene could be involved in genetic variation underlying the QTL for milk fat synthesis on BTA6.

Low-density lipoprotein receptor-related protein 5 variant A1330V is a determinant of blood pressure in Japanese males

We examined the influence of the A1330V variant in the low-density lipoprotein receptor-related protein 5 gene on blood pressure in a large cohort of Japanese workers. This study used analysis of covariance in a multivariate general linear model to adjust for other potential factors such as age, body mass index, blood chemistry and lifestyle. The target subjects were 1440 males and 1169 females selected from 3834 male and 2591 female workers in a single company. Hypertension was defined as systolic blood pressure >or=140 mm Hg and/or diastolic blood pressure >or=90 mm Hg or the use of antihypertensive medications. Genotype distributions for A1330V in hypertensive males (AA=139(54.5%), AV=101(39.6%), VV=15(5.9%)) and females (AA=48(63.2%), AV=24(31.6%), VV=4(5.3%)) were not significantly different from normotensive males (AA=594(50.1%), AV=488(41.2%), VV=103(8.7%)) and females (AA=568(52.0%), AV=441(40.3%), VV=84(7.7%)). Allele distributions were not significantly different in either sex. In males, analysis of covariance showed that the VV genotype was associated with a 2.5 mm Hg lower diastolic blood pressure and a 2.3 mm Hg lower mean blood pressure than the AA genotype. This study indicates that the 1330V allele is an independent factor for lower diastolic and mean blood pressure in Japanese males.

Strong correlation between meiotic crossovers and haplotype structure in a 2.5-Mb region on the long arm of chromosome 21

Although the haplotype structure of the human genome has been studied in great detail, very little is known about the mechanisms underlying its formation. To investigate the role of meiotic recombination on haplotype block formation, single nucleotide polymorphisms were selected at a high density from a 2.5-Mb region of human chromosome 21. Direct analysis of meiotic recombination by high-throughput multiplex genotyping of 662 single sperm identifies 41 recombinants. The crossovers were nonrandomly distributed within 16 small areas. All, except one, of these crossovers fall in areas where the haplotype structure exhibits breakdown, displaying a strong statistically positive association between crossovers and haplotype block breaks. The data also indicate a particular clustered distribution of recombination hotspots within the region. This finding supports the hypothesis that meiotic recombination makes a primary contribution to haplotype block formation in the human genome.

Modeling haplotype block variation using Markov chains

Models of background variation in genomic regions form the basis of linkage disequilibrium mapping methods. In this work we analyze a background model that groups SNPs into haplotype blocks and represents the dependencies between blocks by a Markov chain. We develop an error measure to compare the performance of this model against the common model that assumes that blocks are independent. By examining data from the International Haplotype Mapping project, we show how the Markov model over haplotype blocks is most accurate when representing blocks in strong linkage disequilibrium. This contrasts with the independent model, which is rendered less accurate by linkage disequilibrium. We provide a theoretical explanation for this surprising property of the Markov model and relate its behavior to allele diversity.

Missense mutations in LRP5 are not a common cause of idiopathic osteoporosis in adult men

We studied whether the LRP5 gene contributes to the clinical phenotype of IO in men. Mutation analysis in 66 IO men revealed a range of sequence variants, of which two missense variants were shown to be of functional relevance.

INTRODUCTION

Mutations in the LDL receptor-related protein 5 (LRP5) gene have been associated with extreme bone phenotypes, which makes LRP5 a plausible candidate gene for idiopathic osteoporosis (IO).

MATERIALS AND METHODS

In 66 men with IO, all 23 exons and exon-intron boundaries of the LRP5 gene were screened for mutations, and functional analyses were performed for those that were putatively involved in the phenotype.

RESULTS

Mutation analysis in the IO probands revealed five missense mutations, of which 1067C>T (S356L), 1364C>T (S455L), and 4609G>A (A1537T) were of potential functional significance because they were located in highly conserved regions of LRP5 and not found in a control panel. Segregation analysis in the respective families could not exclude their possible causality for IO. Furthermore, functional analyses clearly showed an inhibitory effect of mutations 1067C>T and 1364C>T on Wnt signal transduction. These effects are most likely caused by impaired LRP5 synthesis in the case of 1067C>T and failure of protein trafficking to the cell surface for 1364C>T.

CONCLUSIONS

For 2 of 66 IO probands, a mutation in the LRP5 gene with proven functionality was found. The findings indicate that carrying an LRP5 mutation is a risk factor for IO, but that overall, IO in men is infrequently underlied by such a mutation.

Single-nucleotide polymorphisms and haplotype LD analysis of the 29-kb IGF2 region on chromosome 11p15.5 in the Korean population

OBJECTIVE

We investigated sequence variations of the 29-kb insulin-like growth factor 2 (IGF2) region in human chromosome region 11p15.5 in the Korean population. This region consists of IGF2, insulin-like growth factor 2 antisense (IGF2AS), and the insulin gene, all important candidate genes for various diseases, including cancer, obesity, diabetes, and coronary disease. While single nucleotide polymorphisms (SNPs) have been identified for this region and used in association studies, ethnic differences in genetic variation at this site have not been addressed. To date, SNPs for the entire 29-kb region in the Korean population have not been reported.

METHODS

We surveyed a population of 108 Koreans for SNPs in the 29-kb IGF2 region.

RESULTS

We identified 62 SNPs, consisting of 6 SNPs in the promoter region, 17 in the untranslated region, 19 in introns, and 20 in the intergenic region. We also analyzed linkage disequilibrium (LD) patterns and haplotypes using 36 high-frequency (> 5%)SNPs and found a well-defined LD block spanning about 13 kb that includes 8 kb of the IGF2AS gene, with two hot-spot regions flanking the LD block.

CONCLUSION

These SNPs may be useful as genetic markers in disease association studies in the Korean population.

Association of polymorphisms in low-density lipoprotein receptor-related protein 5 gene with bone mineral density in postmenopausal Chinese women

AIM

To investigate the possible association of Q89R, N740N and A1330V polymorphisms in low-density lipoprotein receptor-related protein 5 (LRP5) gene with bone mineral density (BMD) in postmenopausal Chinese women.

METHODS

Q89R, N740N and A1330V genotypes were determined by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) in 647 unrelated healthy postmenopausal Han Chinese women aged 43-76 years in Shanghai. BMD at lumbar spine 1-4 and the left proximal femur including the femoral neck, trochanter and Ward's triangle were measured by dual-energy X-ray absorptionmetry in all subjects.

RESULTS

The distribution of the Q89R, N740N and A1330V genotypes in this population was as follows: QQ 80.5%, QR 18.7%, and RR 0.8%; TT 66.9%, TC 31.1%, and CC 2.0%; AA 68.0%, AV 29.7%, and VV 2.3%. The frequencies of the Q89R, N740N and A1330V genotypes and alleles did not deviate from the Hardy-Weinberg equilibrium. We found that the Q89R and A1330V polymorphisms were in linkage disequilibrium in our population (kappa2=13.50, P<0.01). Both before and after adjusting for age, years since menopause, height, and weight, the Q89R or N740N genotypes were significantly associated with BMD at the femoral neck (P<0.05). Subjects with the Q89R QQ genotype or the N740N TT genotype had a significantly higher BMD at the femoral neck, compared with those with the QR/RR or TC/CC genotypes, respectively. No significant association was found between A1330V polymorphism and BMD at any site.

CONCLUSION

Our findings suggest that the LRP5 gene is a candidate for the genetic determination of BMD in postmenopausal Chinese women.

Haplotypes of the low-density lipoprotein receptor-related protein 5 (LRP5) gene: are they a risk factor in osteoarthritis?

OBJECTIVE

Several genome-wide scans have revealed an osteoarthritis (OA)-susceptibility locus on chromosome 11q in close proximity to the low-density lipoprotein receptor-related protein 5 (LRP5) gene. The regulation of bone mass is under the control of LRP5 and since increased bone mass is thought to play a role in the pathology of OA we examined LRP5 polymorphisms and haplotypes to determine if variants of this locus may predispose to OA.

METHODS

A UK control population of 187 individuals was examined for five commonly occurring polymorphisms against a cohort of 158 DNAs from patients with knee OA. An additional UK cohort was also examined to confirm the findings of the first study; this second group consisted of 110 knee OA patients. Haplotype analysis was also performed on patient and control DNAs.

RESULTS

A study of individual polymorphisms revealed no association with disease. However, haplotype analysis of the initial two populations revealed a common haplotype (C-G-C-C-A) that provided a 1.6-fold increased risk of OA (P(c)=0.021). The data obtained from the second cohort confirmed the initial findings, with a 1.6-fold increased risk observed within this cohort for the risk haplotype (P=0.012).

CONCLUSIONS

A closer investigation of LRP5 and associated Wnt signalling molecules in OA will help determine disease aetiology and the development of novel treatment strategies that specifically target the bone compartment.

Pathogenic mutations and polymorphisms in the lipoprotein receptor-related protein 5 reveal a new biological pathway for the control of bone mass

PURPOSE OF REVIEW

This review summarizes recent findings concerning the genomic variations of the lipoprotein receptor-related protein 5 (LPR5) in relation to bone biology.

RECENT FINDINGS

Mutations in the LRP5 gene causing high bone mass (HBM) and osteoporosis-pseudoglioma (OPPG) underscored the role of the Wnt-LRP5 canonical signaling on bone formation. Additional LRP5 activating mutations have been identified in a variety of sclerosing bone dysplasias, improving the diagnostic classification of these disorders. Association of polymorphisms in LRP5 with bone mineral density indicate that LRP5 genetic variation contribute to the risk of osteoporosis. Transgenic mice carrying the LRP5 HBM mutation have improved bone biomechanical properties, and the molecular mechanisms by which this mutation exerts its effects have been clarified. A number of KO mice have shown the complex effects of the Wnt-LRP5 pathway on bone mass and skeletal morphology. In vitro studies indicate that osteoblasts produce a variety of Wnts, the LRP5 co-receptor frizzled (Fzd), as well as LRP5 and Wnt inhibitors, i.e. dickkopf (Dkk1) and frizzled-related proteins (Sfrps), respectively, and delineate the role of these molecules in regulating the commitment of mesenchymal stem cells along the osteoblastic lineage.

SUMMARY

Identification of pathogenic mutations and allelic variations in LRP5 has improved our understanding of the physiology of bone mass acquisition and the pathophysiology of several bone diseases, including osteoporosis. Understanding how complex interactions between agonistic and inhibitory factors in the Wnt-LRP5 canonical pathway influence osteoblast functions has the potential of providing new anabolic treatments for osteoporosis.

Conservation of theRB1gene in human and primates

Mutations in the RB1 gene are associated with retinoblastoma, which has served as an important model for understanding hereditary predisposition to cancer. Despite the great scrutiny that RB1 has enjoyed as the prototypical tumor suppressor gene, it has never been the object of a comprehensive survey of sequence variation in diverse human populations and primates. Therefore, we analyzed the coding (2,787 bp) and adjacent intronic and untranslated (7,313 bp) sequences of RB1 in 137 individuals from a wide range of ethnicities, including 19 Asian Indian hereditary retinoblastoma cases, and five primate species. Aside from nine apparently disease-associated mutations, 52 variants were identified. They included six singleton, coding variants that comprised five amino acid replacements and one silent site. Nucleotide diversity of the coding region (pi=0.0763+/-1.35 x 10(-4)) was 52 times lower than that of the noncoding regions (pi=3.93+/-5.26 x 10(-4)), indicative of significant sequence conservation. The occurrence of purifying selection was corroborated by phylogeny-based maximum likelihood analysis of the RB1 sequences of human and five primates, which yielded an estimated ratio of replacement to silent substitutions (omega) of 0.095 across all lineages. RB1 displayed extensive linkage disequilibrium over 174 kb, and only four unique recombination events, two in Africa and one each in Europe and Southwest Asia, were observed. Using a parsimony approach, 15 haplotypes could be inferred. Ten were found in Africa, though only 12.4% of the 274 chromosomes screened were of African origin. In non-Africans, a single haplotype accounted for from 63 to 84% of all chromosomes, most likely the consequence of natural selection and a significant bottleneck in effective population size during the colonization of the non-African continents.

Genome-wide association studies: theoretical and practical concerns

To fully understand the allelic variation that underlies common diseases, complete genome sequencing for many individuals with and without disease is required. This is still not technically feasible. However, recently it has become possible to carry out partial surveys of the genome by genotyping large numbers of common SNPs in genome-wide association studies. Here, we outline the main factors - including models of the allelic architecture of common diseases, sample size, map density and sample-collection biases - that need to be taken into account in order to optimize the cost efficiency of identifying genuine disease-susceptibility loci.

Association analysis of signal transducer and activator of transcription 4 (STAT4) polymorphisms with asthma

The signal transducer and activator of transcription 4 (STAT4) on chromosome 2q32.2-q32.3 is known to be essential for mediating responses to interleukin 12 in lymphocytes and regulating the differentiation of T helper cells. In an effort to discover additional polymorphism(s) in genes in which variant(s) have been implicated in asthma, we investigated the genetic polymorphisms in STAT4 to evaluate it as a potential candidate gene for a host genetic study of asthma. By direct DNA sequencing in 24 individuals, we identified 12 sequence variants within introns and their flanking regions, including the 1.5 kb promoter region of STAT4. Among them, seven common polymorphic sites were selected for genotyping in our asthma cohort (502 asthmatic patients, 164 normal controls). Using logistic regression analysis for association with the risk of asthma, while controlling for age, gender, and smoking status as covariates, no significant associations with the risk of asthma were detected. However, one single nucleotide polymorphism (SNP) in intron 11 (+90089T--> C) and two haplotypes showed positive association (P= 0.03, 0.03 and 0.03, respectively) with production of specific IgE to Dermatophagoides farinae (D.f.) or Dermatophagoides pteronyssinus (D.p.) among asthmatic patients. The minor allele frequencies of +90089T--> C and BLOCK2-ht1 were higher (0.54 and 0.47, respectively) among individuals who produced specific IgE to D.f. or D.p. than frequencies (0.47 and 0.39, respectively) among individuals who did not produce specific IgE (OR=1.38 and 1.40, respectively). Our findings suggest that polymorphisms in STAT4 might be one of the genetic factors for the risk of production of specific IgE to mite allergens.

Selecting tagging SNPs for association studies using power calculations from genotype data

Recent studies have indicated that linkage disequilibrium (LD) between single nucleotide polymorphism (SNP) markers can be used to derive a reduced set of tagging SNPs (tSNPs) for genetic association studies. Previous strategies for identifying tSNPs have focused on LD measures or haplotype diversity, but the statistical power to detect disease-associated variants using tSNPs in genetic studies has not been fully characterized. We propose a new approach of selecting tSNPs based on determining the set of SNPs with the highest power to detect association. Two-locus genotype frequencies are used in the power calculations. To show utility, we applied this power method to a large number of SNPs that had been genotyped in Caucasian samples. We demonstrate that a significant reduction in genotyping efforts can be achieved although the reduction depends on genotypic relative risk, inheritance mode and the prevalence of disease in the human population. The tSNP sets identified by our method are remarkably robust to changes in the disease model when small relative risk and additive mode of inheritance are employed. We have also evaluated the ability of the method to detect unidentified SNPs. Our findings have important implications in applying tSNPs from different data sources in association studies.

Sequence analysis of the mannose-binding lectin (MBL2) gene reveals a high degree of heterozygosity with evidence of selection

Human mannose-binding protein (MBL) is a component of innate immunity. To capture the common genetic variants of MBL2, we resequenced a 10.0 kb region that includes MBL2 in 102 individuals representing four major US ethnic groups. In all, 87 polymorphic sites were observed, indicating a high level of heterozygosity (total pi=18.3 x 10(-4)). Estimates of linkage disequilibrium across MBL2 indicate that it is divided into two blocks, with a probable recombination hot spot in the 3' end. Three non-synonymous SNPs in exon 1 of the encoding MBL2 gene and three upstream SNPs form common 'secretor haplotypes' that can predict circulating levels. Common variants have been associated with increased susceptibility to infection and autoimmune diseases. The high frequencies of B, C and D alleles in certain populations suggest a possible selective advantage for heterozygosity. There is limited diversity of haplotype structure; the 'secretor haplotypes' lie on a restricted number of extended haplotypes, which could include additional linked SNPs, which might also have possible functional implications. There is evidence for gene conversion in the region between the two blocks, in the last exon. Our data should form the basis for conducting MBL2 candidate gene association studies using a locus-wide approach.

Contribution of the LRP5 gene to normal variation in peak BMD in women

The role of the LRP5 gene in rare BMD-related traits has recently been shown. We tested whether variation in this gene might play a role in normal variation in peak BMD. Association between SNPs in LRP5 and hip and spine BMD was measured in 1301 premenopausal women. Only a small proportion of the BMD variation was attributable to LRP5 in our sample.

INTRODUCTION

Mutations in the low-density lipoprotein receptor-related protein 5 (LRP5) gene have been implicated as the cause of multiple distinct BMD-related rare Mendelian phenotypes. We sought to examine whether the LRP5 gene contributes to the observed variation in peak BMD in the normal population.

MATERIALS AND METHODS

We genotyped 12 single nucleotide polymorphisms (SNPs) in LRP5 using allele-specific PCR and mass spectrometry methods. Linkage disequilibrium between the genotyped LRP5 SNPs was measured. We tested for association between these SNPs and both hip and spine BMD (adjusted for age and body weight) in 1301 healthy premenopausal women who took part in a sibling pair study aimed at identifying the genes underlying peak bone mass. Our study used both population-based (ANOVA) and family-based (quantitative transmission disequilibrium test) association methodology.

RESULTS AND CONCLUSIONS

The linkage disequilibrium pattern and haplotype block structure within the LRP5 gene were consistent with that observed in other studies. Although significant evidence of association was found between LRP5 SNPs and both hip and spine BMD, only a small proportion of the total variation in these phenotypes was accounted for. The genotyped SNPs accounted for approximately 0.8% of the variation in femoral neck BMD and 1.1% of the variation in spine BMD. Results from our sample suggest that natural variation in and around LRP5 is not a major contributor to the observed variability in peak BMD at either the femoral neck or lumbar spine in white women.

Characterization of linkage disequilibrium structure, mutation history, and tagging SNPs, and their use in association analyses:ELAC2 and familial early-onset prostate cancer

In association analyses, it is critical that informative single-nucleotide polymorphisms (SNPs) be selected for study and utilized appropriately. We sequenced 38 kb, including exons of ELAC2, promoter region and conserved upstream intergenic sequences. A comprehensive characterization of linkage disequilibrium (LD) structure and mutation history was performed using our principal components analysis (PCA) method and a phylogenetic analysis. We identified a complex pattern of LD structure consistent with the occurrence of both recombination and mutation events within ELAC2. Four overlapping and noncontiguous LD groups were defined. Eight tagging SNPs (tSNPs) were identified, accounting for over 90% of the genetic variation of the 19 total variants. We tested associations between familial early-onset prostate cancer (PRCA) and each variant independently and in haplotypes. We performed these tests using all 19 variants and the 8 tSNPs; the results using tSNP haplotypes accurately represent the association evidence for the full haplotypes. We observed increased evidence for association when SNPs were analyzed in haplotypes. The phylogenetic analysis indicated three haplotypes, clustered farthest from the root-node, all of which were found more often in cases than controls. These three haplotypes together showed the best evidence of association with familial, early-onset PRCA (P=0.0024; odds ratio=2.23; 95% CI, 1.33-3.74), indicating possible allelic heterogeneity. Our results suggest that 8 tSNPs are required to comprehensively assess associations in ELAC2, and that haplotypes should be considered for analysis, and that a knowledge of mutation history may be helpful in parsing allelic heterogeneity and suggesting combinations of haplotypes to be tested.

Influence of LRP5 polymorphisms on normal variation in BMD

Genetic studies based on cohorts with rare and extreme bone phenotypes have shown that the LRP5 gene is an important genetic modulator of BMD. Using family-based and case-control approaches, this study examines the role of the LRP5 gene in determining normal population variation of BMD and describes significant association and suggestive linkage between LRP5 gene polymorphisms and BMD in >900 individuals with a broad range of BMD.

INTRODUCTION

Osteoporosis is a common, highly heritable condition determined by complex interactions of genetic and environmental etiologies. Genetic factors alone can account for 50-80% of the interindividual variation in BMD. Mutations in the LRP5 gene on chromosome 11q12-13 have been associated with rare syndromes characterized by extremely low or high BMD, but little is known about the contribution of this gene to the development of osteoporosis and determination of BMD in a normal population.

MATERIALS AND METHODS

To examine the entire spectrum of low to high BMD, 152 osteoporotic probands, their families (597 individuals), and 160 women with elevated BMD (T score > 2.5) were recruited. BMD at the lumbar spine, femoral neck, and hip were measured in each subject using DXA.

RESULTS

PAGE sequencing of the LRP5 gene revealed 10 single nucleotide polymorphisms (SNPs), 8 of which had allele frequencies of >5%, in exons 8, 9, 10, 15, and 18 and in introns 6, 7, and 21. Within families, a strong association was observed between an SNP at nucleotide C171346A in intron 21 and total hip BMD (p < 1 x 10(-5) in men only, p = 0.0019 in both men and women). This association was also observed in comparisons of osteoporotic probands and unrelated elevated BMD in women (p = 0.03), along with associations with markers in exons 8 (C135242T, p = 0.007) and 9 (C141759T, p = 0.02). Haplotypes composed of two to three of the SNPs G121513A, C135242T, G138351A, and C141759T were strongly associated with BMD when comparing osteoporotic probands and high BMD cases (p < 0.003). An SNP at nucleotide C165215T in exon 18 was linked to BMD at the lumbar spine, femoral neck, and total hip (parametric LOD scores = 2.8, 2.5, and 2.2 and nonparametric LOD scores = 0.3, 1.1, and 2.2, respectively) but was not genetically associated with BMD variation.

CONCLUSION

These results show that common LRP5 polymorphisms contribute to the determination of BMD in the general population.

Association between bone mineral density and LDL receptor-related protein 5 gene polymorphisms in young Korean men

Recently, It has been reported that the LDL receptor-related protein 5 (LRP5) regulates bone formation, and that mutations of the gene cause osteoporosis-pseudoglioma syndrome or high bone mass phenotypes. However, the mutations cannot explain a genetic trait for osteoporosis in the general population because of their rarity. From 219 Korean men aged 20-34 yr, we looked for six known polymorphisms causing amino acid changes in the LRP5 coding region, and investigated their association with bone mineral density (BMD) at the following anatomical sites: lumbar spine (L2-L4) and the left proximal femur (femoral neck, Ward's triangle, trochanter and shaft). We found that the Q89R polymorphism was significantly associated with BMD at the femoral neck and Ward's triangle (p=0.004 and <0.001, respectively). However, after adjusting for age, weight and height, a statistically significant association only occurred at the Ward's triangle (p=0.043), and a marginal association was observed at the femoral neck (p=0.098). No A400V, V667M, R1036Q and A1525V polymorphisms were found, and no statistically significant association was found between the A1330V polymorphism and BMD at any sites. Although we failed to demonstrate a clear association between the LRP5 polymorphism and peak bone mass in young men, the present study suggests that larger-scale studies on the Q89R polymorphism need to be performed.

LRP5, low-density-lipoprotein-receptor-related protein 5, is a determinant for bone mineral density

Osteoporosis is a multifactorial trait with low bone mineral density (BMD). We report results of an association study between BMD and nine candidate genes ( TGFB1, TGFBR2, SMAD2, SMAD3, SMAD4, IFNB1, IFNAR1, FOS and LRP5), as well as of a case-control study of osteoporosis. Samples for the former association study included 481 general Japanese women. Among the nine candidate genes examined, only LRP5 showed a significant association with BMD. We identified a strong linkage disequilibrium (LD) block within LRP5. Of five LPR5 single nucleotide polymorphisms (SNPs) that are located in the LD block, three gave relatively significant results: Women with the C/C genotype at the c.2220C>T SNP site had higher adjusted BMD (AdjBMD) value compared to those with C/T and T/T (p=0.022); and likewise, G/G at IVS17-30G>A and C/C women at c.3989C>T showed higher AdjBMD than those with G/A or A/A (p=0.039) and with C/T or T/T ( p=0.053), respectively. The case-control study in another series of samples consisting of 126 osteoporotic patients and 131 normal controls also gave a significant difference in allele frequency at c.2220C>T (kappa2=6.737, p=0.009). These results suggest that LRP5 is a BMD determinant and also contributes to a risk of osteoporosis.