Test of linkage and/or association between the estrogen receptor alpha gene with bone mineral density in Caucasian nuclear families

Intronic SNP in ESR1 encoding human estrogen receptor alpha is associated with brain ESR1 mRNA isoform expression and behavioral traits

Genetic variants of ESR1 have been implicated in multiple diseases, including behavioral disorders, but causative variants remain uncertain. We have searched for regulatory variants affecting ESR1 expression in human brain, measuring allelic ESR1 mRNA expression in human brain tissues with marker SNPs in exon4 representing ESR1-008 (or ESRα-36), and in the 3'UTR of ESR1-203, two main ESR1 isoforms in brain. In prefrontal cortex from subjects with bipolar disorder, schizophrenia, and controls (n = 35 each; Stanley Foundation brain bank), allelic ESR1 mRNA ratios deviated from unity up to tenfold at the exon4 marker SNP, with large allelic ratios observed primarily in bipolar and schizophrenic subjects. SNP scanning and targeted sequencing identified rs2144025, associated with large allelic mRNA ratios (p = 1.6E10-6). Moreover, rs2144025 was significantly associated with ESR1 mRNA levels in the Brain eQTL Almanac and in brain regions in the Genotype-Tissue Expression project. In four GWAS cohorts, rs2104425 was significantly associated with behavioral traits, including: hypomanic episodes in female bipolar disorder subjects (GAIN bipolar disorder study; p = 0.0004), comorbid psychological symptoms in both males and females with attention deficit hyperactivity disorder (GAIN ADHD, p = 0.00002), psychological diagnoses in female children (eMERGE study of childhood health, subject age ≥9, p = 0.0009), and traits in schizophrenia (e.g., grandiose delusions, GAIN schizophrenia, p = 0.0004). The first common ESR1 variant (MAF 12-33% across races) linked to regulatory functions, rs2144025 appears conditionally to affect ESR1 mRNA expression in the brain and modulate traits in behavioral disorders.

Bivariate genome-wide linkage analysis for traits BMD and AAM: effect of menopause on linkage signals

Osteoporosis is an age-related systemic skeletal disease, characterized by low bone mineral density (BMD). Low BMD is closely associated with late age at menarche (AAM). Our previous bivariate genome-wide linkage analyses (GWLAs) between BMD and AAM identified two shared genomic regions in 2584 Caucasian females including both pre- and post-menopausal females. However, menopause often causes dramatic bone loss in post-menopausal females; this may introduce some confounding effects on the bivariate GWLA for BMD and AAM. To address the effect of menopause on the identification of genetic factors shared by BMD and AAM, we segregated the previously studied population of 2584 females into two separate subgroups consisting of 1462 pre-menopause subjects and 1122 post-menopausal subjects, and performed further bivariate GWLAs. The BMD was measured by Hologic Dual-energy X-ray (DXA) scanners (Hologic Inc., Bedford, MA, USA). Based on the genome-wide thresholds corrected for multiple testing, we found more significant genomic regions in the pre-menopausal group than in total group (including pre- and post-menopausal women), e.g., we found 4, 1, and 2 shared by spine BMD and AAM, femoral neck (FNK) BMD and AAM and ultra distal (UD) BMD and AAM, respectively. We did not found any significant linkage signals in the post-menopausal group. Importantly, the linkage signals at all significant regions were much stronger in pre-menopausal group than in the other groups: post-menopausal females and total females. For example, the linkage LOD score for FNK BMD and AAM is as high as 4.88 in pre-menopausal females, but only 0.24 and 0.31 in post-menopausal and total females, respectively. These results suggest that menopause introduces some noise signals into GWLAs when estimating the shared genetic factors by BMD and AAM. Therefore, it is very important to classify female subjects properly according to their menopause stage when performing such studies.

A bivariate whole-genome linkage scan suggests several shared genomic regions for obesity and osteoporosis

CONTEXT

A genome-wide bivariate analysis was conducted for body fat mass (BFM) and bone mineral density (BMD) in a large Caucasian sample. We found some quantitative trait loci shared by BFM and BMD in the total sample and the gender-specific subgroups, and quantitative trait loci with potential pleiotropy were disclosed. BFM and BMD, as the respective measure for obesity and osteoporosis, are phenotypically and genetically correlated. However, specific genomic regions accounting for their genetic correlation are unknown.

OBJECTIVE

To identify systemically the shared genomic regions for BFM and BMD, we performed a bivariate whole-genome linkage scan in 4498 Caucasian individuals from 451 families for BFM and BMD at the hip, spine, and wrist, respectively. Linkage analyses were performed in the total sample and the male and female subgroups, respectively.

RESULTS

In the entire sample, suggestive linkages were detected at 7p22-p21 (LOD 2.69) for BFM and spine BMD, 6q27 (LOD 2.30) for BFM and hip BMD, and 11q13 (LOD 2.64) for BFM and wrist BMD. Male-specific suggestive linkages were found at 13q12 (LOD 3.23) for BFM and spine BMD and at 7q21 (LOD 2.59) for BFM and hip BMD. Female-specific suggestive LOD scores were 3.32 at 15q13 for BFM and spine BMD and 3.15 at 6p25-24 for BFM and wrist BMD.

CONCLUSIONS

Several shared genomic regions for BFM and BMD were identified here. Our data may benefit further positional and functional studies, aimed at eventually uncovering the complex mechanism underlying the shared genetic determination of obesity and osteoporosis.

Molecular genetic studies of gene identification for osteoporosis: a 2004 update

This review summarizes comprehensively the most important and representative molecular genetics studies of gene identification for osteoporosis published up to the end of December 2004. It is intended to constitute a sequential update of our previously published review covering the available data up to the end of 2002. Evidence from candidate gene association studies and genome-wide linkage studies in humans, as well as quantitative trait locus mapping animal models are reviewed separately. Studies of transgenic and knockout mice models relevant to osteoporosis are summarized. An important extension of this update is incorporation of functional genomic studies (including DNA microarrays and proteomics) on osteogenesis and osteoporosis, in light of the rapid advances and the promising prospects of the field. Comments are made on the most notable findings and representative studies for their potential influence and implications on our present understanding of genetics of osteoporosis. The format adopted by this review should be ideal for accommodating future new advances and studies.

Is a gene important for bone resorption a candidate for obesity? An association and linkage study on the RANK (receptor activator of nuclear factor-kappaB) gene in a large Caucasian sample

In light of findings that osteoporosis and obesity may share some common genetic determination and previous reports that RANK (receptor activator of nuclear factor-kappaB) is expressed in skeletal muscles which are important for energy metabolism, we hypothesize that RANK, a gene essential for osteoclastogenesis, is also important for obesity. In order to test the hypothesis with solid data we first performed a linkage analysis around the RANK gene in 4,102 Caucasian subjects from 434 pedigrees, then we genotyped 19 SNPs in or around the RANK gene. A family-based association test (FBAT) was performed with both a quantitative measure of obesity [fat mass, lean mass, body mass index (BMI), and percentage fat mass (PFM)] and a dichotomously defined obesity phenotype-OB (OB if BMI > or = 30 kg/m(2)). In the linkage analysis, an empirical P = 0.004 was achieved at the location of the RANK gene for BMI. Family-based association analysis revealed significant associations of eight SNPs with at least one obesity-related phenotype (P < 0.05). Evidence of association was obtained at SNP10 (P = 0.002) and SNP16 (P = 0.001) with OB; SNP1 with fat mass (P = 0.003); SNP1 (P = 0.003) and SNP7 (P = 0.003) with lean mass; SNP1 (P = 0.002) and SNP7 (P = 0.002) with BMI; SNP1 (P = 0.003), SNP4 (P = 0.007), and SNP7 (P = 0.002) with PFM. In order to deal with the complex multiple testing issues, we performed FBAT multi-marker test (FBAT-MM) to evaluate the association between all the 18 SNPs and each obesity phenotype. The P value is 0.126 for OB, 0.033 for fat mass, 0.021 for lean mass, 0.016 for BMI, and 0.006 for PFM. The haplotype data analyses provide further association evidence. In conclusion, for the first time, our results suggest that RANK is a novel candidate for determination of obesity.

Assessment of linkage and association of 13 genetic loci with bone mineral density

Bone mineral density (BMD), an important risk factor for osteoporosis, is a complex trait likely affected by multiple genes. The linkage and/or association of 13 polymorphic loci of seven candidate genes (estrogen receptor alpha [ERalpha] and beta [ERbeta], calcium-sensing receptor, vitamin D receptor, collagen type 1alpha1, low-density lipoprotein [LDL] receptor-related protein 5 [LRPS], and transforming growth factor beta1) were evaluated in 177 southern Chinese pedigrees of 674 subjects, with each pedigree identified through a proband having a BMD Z score of -1.28 or less at the hip or spine. A suggestive linkage was detected between the IVS1-351A/G polymorphism of ERalpha and spine BMD, and between the 1082G/A, 1730G/A, and D14S1026 polymorphisms of ERbeta and BMD at both spine and hip. The quantitative transmission disequilibrium test (QTDT) detected total family association between 1730G/A of ERbeta and BMD at spine and hip; between D14S1026 of ERbeta and hip BMD; and between the 266A/G and 2220C/T polymorphisms of LRP5 and hip BMD. Similar total family associations were detected when only the females were analyzed. In addition, the IVS1-397T/C polymorphism of ERalpha was associated with spine BMD, and the 266A/G and 2220C/T polymorphisms of LRP5 were associated with femoral neck BMD in the females. A within-family association was detected with the IVS1-397T/C polymorphism of ERalpha, and the 266A/G and 2220C/T polymorphisms of LRP5 in the females. The effect of each polymorphism on BMD variance ranged from 1% to 4%. In conclusion, ERalpha, ERbeta and LRP5 are important candidate genes determining BMD variation, especially in females.

Association analysis of estrogen receptor alpha gene polymorphisms with cross-sectional geometry of the femoral neck in Caucasian nuclear families

Bone geometry is a key factor in bone strength, which is the ultimate intrinsic determinant of fracture risk. Though the heritability of bone geometry is high, little effort has been spent on searching for the underlying genes. In this study, employing a sample of 1,873 subjects from 405 Caucasian nuclear families, we studied seven single nucleotide polymorphisms (SNPs) and their haplotypes of the ER-alpha gene for association with six hip geometric variables, namely, cross-sectional area (CSA), cortical thickness (CT), endocortical diameter (ED), subperiosteal width (W), sectional modulus (Z) and buckling ratio (BR). The major method used was the quantitative transmission disequilibrium test (QTDT). Our major findings were summarized below. The within-family association between SNP4 (rs1801132) in exon 4 with endocortical diameter and subperiosteal width was detected in single locus analyses (P=0.008 and 0.021, respectively) and verified in haplotype analyses (P=0.034 and 0.058, respectively). The total association of SNP4 with these two diameters was also observed in both single locus and haplotype analyses (P=0.005 and 0.031 for ED, plus P=0.003 and 0.070 for W). In addition, the total association between SNP5 (rs932477) in intron 4 with cortical thickness and buckling ratio was detected (single locus analyses: P=0.035 and 0.041, respectively). Haplotype analyses further supported the above association (P=0.010 and 0.004, respectively). Similar patterns of associations with the studied SNPs and their haplotypes were present in subsamples stratified by sex, too. However, after permutation tests, the empirical significance level was set as P<0.011, which renders most associations insignificant. Therefore, we concluded that polymorphisms in the ER-alpha gene were nominally associated with femoral neck (FN) geometry variables estimated from DXA. Such genetic effects on hip geometry were not sex specific.