Meta-Analysis of Genomewide Association Studies Reveals Genetic Variants for Hip Bone Geometry.

Hip geometry is an important predictor of fracture. We performed a meta-analysis of GWAS studies in adults to identify genetic variants that are associated with proximal femur geometry phenotypes. We analyzed four phenotypes: (i) femoral neck length; (ii) neck-shaft angle; (iii) femoral neck width, and (iv) femoral neck section modulus, estimated from DXA scans using algorithms of hip structure analysis. In the Discovery stage, 10 cohort studies were included in the fixed-effect meta-analysis, with up to 18,719 men and women ages 16 to 93 years. Association analyses were performed with ∼2.5 million polymorphisms under an additive model adjusted for age, body mass index, and height. Replication analyses of meta-GWAS significant loci (at adjusted genomewide significance [GWS], threshold p ≤ 2.6 × 10-8 ) were performed in seven additional cohorts in silico. We looked up SNPs associated in our analysis, for association with height, bone mineral density (BMD), and fracture. In meta-analysis (combined Discovery and Replication stages), GWS associations were found at 5p15 (IRX1 and ADAMTS16); 5q35 near FGFR4; at 12p11 (in CCDC91); 11q13 (near LRP5 and PPP6R3 (rs7102273)). Several hip geometry signals overlapped with BMD, including LRP5 (chr. 11). Chr. 11 SNP rs7102273 was associated with any-type fracture (p = 7.5 × 10-5 ). We used bone transcriptome data and discovered several significant eQTLs, including rs7102273 and PPP6R3 expression (p = 0.0007), and rs6556301 (intergenic, chr.5 near FGFR4) and PDLIM7 expression (p = 0.005). In conclusion, we found associations between several genes and hip geometry measures that explained 12% to 22% of heritability at different sites. The results provide a defined set of genes related to biological pathways relevant to BMD and etiology of bone fragility. © 2019 American Society for Bone and Mineral Research.

Efficient region-based test strategy uncovers genetic risk factors for functional outcome in bipolar disorder.

Genome-wide association studies of case-control status have advanced the understanding of the genetic basis of psychiatric disorders. Further progress may be gained by increasing sample size but also by new analysis strategies that advance the exploitation of existing data, especially for clinically important quantitative phenotypes. The functionally-informed efficient region-based test strategy (FIERS) introduced herein uses prior knowledge on biological function and dependence of genotypes within a powerful statistical framework with improved sensitivity and specificity for detecting consistent genetic effects across studies. As proof of concept, FIERS was used for the first genome-wide single nucleotide polymorphism (SNP)-based investigation on bipolar disorder (BD) that focuses on an important aspect of disease course, the functional outcome. FIERS identified a significantly associated locus on chromosome 15 (hg38: chr15:48965004 - 49464789 bp) with consistent effect strength between two independent studies (GAIN/TGen: European Americans, BOMA: Germans; n = 1592 BD patients in total). Protective and risk haplotypes were found on the most strongly associated SNPs. They contain a CTCF binding site (rs586758); CTCF sites are known to regulate sets of genes within a chromatin domain. The rs586758 - rs2086256 - rs1904317 haplotype is located in the promoter flanking region of the COPS2 gene, close to microRNA4716, and the EID1, SHC4, DTWD1 genes as plausible biological candidates. While implication with BD is novel, COPS2, EID1, and SHC4 are known to be relevant for neuronal differentiation and function and DTWD1 for psychopharmacological side effects. The test strategy FIERS that enabled this discovery is equally applicable for tag SNPs and sequence data.

Erratum: Large meta-analysis of genome-wide association studies identifies five loci for lean body mass.

A correction to this article has been published and is linked from the HTML version of this article.

Large meta-analysis of genome-wide association studies identifies five loci for lean body mass.

Lean body mass, consisting mostly of skeletal muscle, is important for healthy aging. We performed a genome-wide association study for whole body (20 cohorts of European ancestry with n = 38,292) and appendicular (arms and legs) lean body mass (n = 28,330) measured using dual energy X-ray absorptiometry or bioelectrical impedance analysis, adjusted for sex, age, height, and fat mass. Twenty-one single-nucleotide polymorphisms were significantly associated with lean body mass either genome wide (p < 5 × 10-8) or suggestively genome wide (p < 2.3 × 10-6). Replication in 63,475 (47,227 of European ancestry) individuals from 33 cohorts for whole body lean body mass and in 45,090 (42,360 of European ancestry) subjects from 25 cohorts for appendicular lean body mass was successful for five single-nucleotide polymorphisms in/near HSD17B11, VCAN, ADAMTSL3, IRS1, and FTO for total lean body mass and for three single-nucleotide polymorphisms in/near VCAN, ADAMTSL3, and IRS1 for appendicular lean body mass. Our findings provide new insight into the genetics of lean body mass.Lean body mass is a highly heritable trait and is associated with various health conditions. Here, Kiel and colleagues perform a meta-analysis of genome-wide association studies for whole body lean body mass and find five novel genetic loci to be significantly associated.

Genetic Variants Associated with Circulating Parathyroid Hormone.

Parathyroid hormone (PTH) is a primary calcium regulatory hormone. Elevated serum PTH concentrations in primary and secondary hyperparathyroidism have been associated with bone disease, hypertension, and in some studies, cardiovascular mortality. Genetic causes of variation in circulating PTH concentrations are incompletely understood. We performed a genome-wide association study of serum PTH concentrations among 29,155 participants of European ancestry from 13 cohort studies (n=22,653 and n=6502 in discovery and replication analyses, respectively). We evaluated the association of single nucleotide polymorphisms (SNPs) with natural log-transformed PTH concentration adjusted for age, sex, season, study site, and principal components of ancestry. We discovered associations of SNPs from five independent regions with serum PTH concentration, including the strongest association with rs6127099 upstream of CYP24A1 (P=4.2 × 10-53), a gene that encodes the primary catabolic enzyme for 1,25-dihydroxyvitamin D and 25-dihydroxyvitamin D. Each additional copy of the minor allele at this SNP associated with 7% higher serum PTH concentration. The other SNPs associated with serum PTH concentration included rs4074995 within RGS14 (P=6.6 × 10-17), rs219779 adjacent to CLDN14 (P=3.5 × 10-16), rs4443100 near RTDR1 (P=8.7 × 10-9), and rs73186030 near CASR (P=4.8 × 10-8). Of these five SNPs, rs6127099, rs4074995, and rs219779 replicated. Thus, common genetic variants located near genes involved in vitamin D metabolism and calcium and renal phosphate transport associated with differences in circulating PTH concentrations. Future studies could identify the causal variants at these loci, and the clinical and functional relevance of these variants should be pursued.

NIPT and Informed Consent: an Assessment of Patient Understanding of a Negative NIPT Result.

Since becoming clinically available in 2011, the use of noninvasive prenatal testing (NIPT) to screen for fetal aneuploidy has continued to increase. However, it has been questioned whether the education of patients undergoing NIPT consistently meets informed consent standards. We sought to evaluate patients' basic understanding of NIPT, such as conditions assessed and accuracy. In addition, we investigated patient self-assessment of NIPT knowledge and satisfaction with the testing process. We distributed an anonymous paper survey to pregnant women during prenatal visits following a negative NIPT result. The survey assessed patient NIPT knowledge, gathered pregnancy-specific and demographic information, and allowed respondents to rank their basic understanding of NIPT and provide written feedback about the testing process. A total of 95 completed and 3 partially completed surveys were returned. Participants scored lowest on knowledge questions involving whether a negative NIPT result ensures a healthy baby or eliminates the possibility of Down syndrome. Most perceived themselves to have a good basic understanding of NIPT and two-thirds of the written feedback proposed no changes to NIPT administration. Overall, most patients appear satisfied with their understanding of NIPT and the testing process, yet they may not fully appreciate the limitations of this screening method.

Genome-wide association study of serum iron phenotypes in premenopausal women of European descent.

A genome-wide association study was performed on 1130 premenopausal women to detect common variants associated with three serum iron-related phenotypes. Total iron binding capacity was strongly associated (p=10(-14)) with variants in and near the TF gene (transferrin), the serum iron transporting protein, and with variants in HFE (p=4×10(-7)), which encodes the human hemochromatosis gene. Association was also detected between percent iron saturation (p=10(-8)) and variants in the chromosome 6 region containing both HFE and SLC17A2, which encodes a phosphate transport protein. No significant associations were detected with serum iron, but variants in HFE were suggestive (p=10(-6)). Our results corroborate prior studies in older subjects and demonstrate that the association of these genetic variants with iron phenotypes can be detected in premenopausal women.

Fine mapping of bone structure and strength QTLs in heterogeneous stock rat.

We previously demonstrated that skeletal structure and strength phenotypes vary considerably in heterogeneous stock (HS) rats. These phenotypes were found to be strongly heritable, suggesting that the HS rat model represents a unique genetic resource for dissecting the complex genetic etiology underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone structure and strength phenotypes using 1524 adult male and female HS rats between 17 to 20 weeks of age. Structure measures included femur length, neck width, head width; femur and lumbar spine (L3-5) areas obtained by DXA; and cross-sectional areas (CSA) at the midshaft, distal femur and femoral neck, and the 5th lumbar vertebra measured by CT. In addition, measures of strength of the whole femur and femoral neck were obtained. Approximately 70,000 polymorphic SNPs distributed throughout the rat genome were selected for genotyping, with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent at each locus from each of the 8 HS founder strains. The haplotypes were then tested for association with each structure and strength phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for structure phenotypes on chromosomes 3, 8, 10, 12, 17 and 20, and QTLs for strength phenotypes on chromosomes 5, 10 and 11 that met a conservative genome-wide empiric significance threshold (FDR=5%; P<3×10(-6)). Importantly, most QTLs were localized to very narrow genomic regions (as small as 0.3 Mb and up to 3 Mb), each harboring a small set of candidate genes, both novel and previously shown to have roles in skeletal development and homeostasis.

Assessment of first and second degree relatives of individuals with bipolar disorder shows increased genetic risk scores in both affected relatives and young At-Risk Individuals.

Recent studies have revealed the polygenic nature of bipolar disorder (BP), and identified common risk variants associated with illness. However, the role of common polygenic risk in multiplex families has not previously been examined. The present study examined 249 European-ancestry families from the NIMH Genetics Initiative sample, comparing subjects with narrowly defined BP (excluding bipolar II and recurrent unipolar depression; n = 601) and their adult relatives without BP (n = 695). Unrelated adult controls (n = 266) were from the NIMH TGEN control dataset. We also examined a prospective cohort of young (12-30 years) offspring and siblings of individuals with BPI and BPII disorder (at risk; n = 367) and psychiatrically screened controls (n = 229), ascertained from five sites in the US and Australia and assessed with standardized clinical protocols. Thirty-two disease-associated SNPs from the PGC-BP Working Group report (2011) were genotyped and additive polygenic risk scores (PRS) derived. We show increased PRS in adult cases compared to unrelated controls (P = 3.4 × 10(-5) , AUC = 0.60). In families with a high-polygenic load (PRS score ≥32 in two or more subjects), PRS distinguished cases with BPI/SAB from other relatives (P = 0.014, RR = 1.32). Secondly, a higher PRS was observed in at-risk youth, regardless of affected status, compared to unrelated controls (GEE-χ(2) = 5.15, P = 0.012). This report is the first to explore common polygenic risk in multiplex families, albeit using only a small number of robustly associated risk variants. We show that individuals with BP have a higher load of common disease-associated variants than unrelated controls and first-degree relatives, and illustrate the potential utility of PRS assessment in a family context.

Genetic and childhood trauma interaction effect on age of onset in bipolar disorder: An exploratory analysis.

INTRODUCTION: This study investigated whether early life trauma mediates genetic effects on the age at onset (AAO) of bipolar disorder. METHOD: Data from the BiGS Consortium case samples (N=1119) were used. Childhood traumatic events were documented using the Childhood Life Events Scale (CLES). Interaction between occurrence of childhood trauma and common genetic variants throughout the genome was tested to identify single nucleotide polymorphic gene variants (SNPs) whose effects on bipolar AAO differ between individuals clearly exposed (CLES≥2) and not exposed (CLES=0) to childhood trauma. RESULTS: The modal response to the CLES was 0 (N=480), but an additional 276 subjects had CLES=1, and 363 subjects reported 2 or more traumatic lifetime events. The distribution of age at onset showed a broad peak between ages 12 and 18, with the majority of subjects having onset during that period, and a significant decrease in age of onset with the number of traumatic events. No single SNP showed a statistically significant interaction with the presence of traumatic events to impact bipolar age at onset. However, SNPs in or near genes coding for calcium channel activity-related proteins (Gene Ontology: 0005262) were found to be more likely than other SNPs to show evidence of interaction using the INRICH method (p<0.001). LIMITATIONS: Retrospective ascertainment of trauma and AAO. CONCLUSION: Interaction effects of early life trauma with genotype may have a significant effect on the development and manifestation of bipolar disorder. These effects may be mediated in part by genes involved in calcium signaling.

Characteristics of Bipolar I patients grouped by externalizing disorders.

BACKGROUND: Bipolar disorder co-occurs with a number of disorders with externalizing features. The aim of this study is to determine whether Bipolar I (BPI) subjects with comorbid externalizing disorders and a subgroup with externalizing symptoms prior to age 15 have different clinical features than those without externalizing disorders and whether these could be attributed to specific genetic variations. METHODS: A large cohort (N=2505) of Bipolar I subjects was analyzed. Course of illness parameters were compared between an Externalizing Group, an Early-Onset Subgroup and a Non-Externalizing Group in the Discovery sample (N=1268). Findings were validated using an independent set of 1237 BPI subjects (Validation sample). Genetic analyses were carried out. RESULTS: Subjects in the Externalizing Group (and Early-Onset Subgroup) tended to have a more severe clinical course, even in areas specifically related to mood disorder such as cycling frequency and rapid mood switching. Regression analysis showed that the differences are not completely explainable by substance use. Genetic analyses identified nominally associated SNPs; calcium channel genes were not enriched in the gene variants identified. LIMITATIONS: Validation in independent samples is needed to confirm the genetic findings in the present study. CONCLUSIONS: Our findings support the presence of an externalizing disorder subphenotype within BPI with greater severity of mood disorder and possible specific genetic features.

Lessons learned from whole exome sequencing in multiplex families affected by a complex genetic disorder, intracranial aneurysm.

Genetic risk factors for intracranial aneurysm (IA) are not yet fully understood. Genomewide association studies have been successful at identifying common variants; however, the role of rare variation in IA susceptibility has not been fully explored. In this study, we report the use of whole exome sequencing (WES) in seven densely-affected families (45 individuals) recruited as part of the Familial Intracranial Aneurysm study. WES variants were prioritized by functional prediction, frequency, predicted pathogenicity, and segregation within families. Using these criteria, 68 variants in 68 genes were prioritized across the seven families. Of the genes that were expressed in IA tissue, one gene (TMEM132B) was differentially expressed in aneurysmal samples (n=44) as compared to control samples (n=16) (false discovery rate adjusted p-value=0.023). We demonstrate that sequencing of densely affected families permits exploration of the role of rare variants in a relatively common disease such as IA, although there are important study design considerations for applying sequencing to complex disorders. In this study, we explore methods of WES variant prioritization, including the incorporation of unaffected individuals, multipoint linkage analysis, biological pathway information, and transcriptome profiling. Further studies are needed to validate and characterize the set of variants and genes identified in this study.

Report of new haplotype for ABCC2 gene: rs17222723 and rs8187718 in cis.

The ATP-binding cassette, subfamily C [CFTR/MRP], member 2 (ABCC2) gene is a member of the ATP-binding cassette transporters and is involved in the transport of molecules across cellular membranes. Substrates transported by ABCC2 include antiepileptics, statins, tenofovir, cisplatin, irinotecan, and carbamazepine. Because of the pharmacogenomics implications, we developed a clinical laboratory-developed assay to test for seven variants in the ABCC2 gene: c.3563T>A (p.V1188E, rs17222723), c.1249G>A (p.V417I, rs2273697), c.3972C>T (p.I1324I, rs3740066), c.2302C>T (p.R768W, rs56199535), c.2366C>T (p.S789F, rs56220353), c.-24C>T (5'UTR, rs717620), and c.4544G>A (p.C1515Y, rs8187710). During the validation process, we noted several DNA samples, obtained from the Coriell Cell Repository, that contained both c.3563T>A, c.4544G>A, and a third variant, suggesting that c.3563T>A and c.4544G>A are in cis on the chromosome in some individuals. We obtained DNA samples from a trio (father, mother, and child), tested their ABCC2 variants, and confirmed that c.3563T>A and c.4544G>A were in cis on the same chromosome. Here, we report a new haplotype in ABCC2.

Affected twins in the familial intracranial aneurysm study.

BACKGROUND AND PURPOSE: Very few cases of intracranial aneurysms (IAs) in twins have been reported. Previous work has suggested that vulnerability to IA formation is heritable. Twin studies provide an opportunity to evaluate the impact of genetics on IA characteristics, including IA location. We therefore sought to examine IA location concordance, multiplicity, and rupture status within affected twin-pairs. METHODS: The Familial Intracranial Aneurysm study was a multicenter study whose goal was to identify genetic and other risk factors for formation and rupture of IAs. The study required at least three affected family members or an affected sibling pair for inclusion. Subjects with fusiform aneurysms, an IA associated with an AVM, or a family history of conditions known to predispose to IA formation, such as polycystic kidney disease, Ehlers-Danlos syndrome, Marfan syndrome, fibromuscular dysplasia, or moyamoya syndrome were excluded. Twin-pairs were identified by birth date and were classified as monozygotic (MZ) or dizygotic (DZ) through DNA marker genotypes. In addition to zygosity, we evaluated twin-pairs by smoking status, major arterial territory of IAs, and rupture status. Location concordance was defined as the presence of an IA in the same arterial distribution (ICA, MCA, ACA, and vertebrobasilar), irrespective of laterality, in both members of a twin-pair. The Fisher exact test was used for comparisons between MZ and DZ twin-pairs. RESULTS: A total of 16 affected twin-pairs were identified. Location concordance was observed in 8 of 11 MZ twin-pairs but in only 1 of 5 DZ twin-pairs (p = 0.08). Three MZ subjects had unknown IA locations and comprised the three instances of MZ discordance. Six of the 11 MZ twin-pairs and none of the 5 DZ twin-pairs had IAs in the ICA distribution (p = 0.03). Multiple IAs were observed in 11 of 22 MZ and 5 of 10 DZ twin-pairs. Thirteen (13) of the 32 subjects had an IA rupture, including 10 of 22 MZ twins. CONCLUSIONS: We found that arterial location concordance was greater in MZ than DZ twins, which suggests a genetic influence upon aneurysm location. The 16 twin-pairs in the present study are nearly the total of affected twin-pairs that have been reported in the literature to date. Further studies are needed to determine the impact of genetics in the formation and rupture of IAs.

Association of substance dependence phenotypes in the COGA sample.

Alcohol and drug use disorders are individually heritable (50%). Twin studies indicate that alcohol and substance use disorders share common genetic influences, and therefore may represent a more heritable form of addiction and thus be more powerful for genetic studies. This study utilized data from 2322 subjects from 118 European-American families in the Collaborative Study on the Genetics of Alcoholism sample to conduct genome-wide association analysis of a binary and a continuous index of general substance dependence liability. The binary phenotype (ANYDEP) was based on meeting lifetime criteria for any DSM-IV dependence on alcohol, cannabis, cocaine or opioids. The quantitative trait (QUANTDEP) was constructed from factor analysis based on endorsement across the seven DSM-IV criteria for each of the four substances. Heritability was estimated to be 54% for ANYDEP and 86% for QUANTDEP. One single-nucleotide polymorphism (SNP), rs2952621 in the uncharacterized gene LOC151121 on chromosome 2, was associated with ANYDEP (P = 1.8 × 10(-8) ), with support from surrounding imputed SNPs and replication in an independent sample [Study of Addiction: Genetics and Environment (SAGE); P = 0.02]. One SNP, rs2567261 in ARHGAP28 (Rho GTPase-activating protein 28), was associated with QUANTDEP (P = 3.8 × 10(-8) ), and supported by imputed SNPs in the region, but did not replicate in an independent sample (SAGE; P = 0.29). The results of this study provide evidence that there are common variants that contribute to the risk for a general liability to substance dependence.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: ß = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

SIBLING family genes and bone mineral density: association and allele-specific expression in humans.

Osteoporosis is a common complex disorder with reduced bone mineral density (BMD) and increased susceptibility to fracture. Peak BMD is one of the primary determinants of osteoporotic fracture risk, and is under substantial genetic control. Extracellular matrix, a major component of the bone, influences BMD by regulating mineral deposition and maintaining cellular activity. It contains several SIBLING family proteins, null mutations of which cause mineralization defects in humans. In this study, we tested 59 single-nucleotide polymorphisms (SNPs) located in the 5 SIBLING family genes (DSPP, DMP1, IBSP, MEPE and SPP1) for association with normal variation in peak BMD in healthy men and women. We measured femoral neck (FN) and lumbar spine (LS) areal BMD by dual energy x-ray absorptiometry (DXA) in 1692 premenopausal European-American women, 512 premenopausal African-American women and 715 European-American men. SNPs were tested for association with FN and LS-BMD in the 3 subsamples. In the European-American women, we observed association (p≤0.005) with LS-BMD for SNPs in DSPP, IBSP and MEPE, and for FN-BMD with SNPs in DMP1 and IBSP. Allele-specific regulation of gene expression (ASE) is an important mechanism in which an allele giving rise to modest influence in transcript abundance might result in a predisposition to disease. To identify whether there was ASE of SIBLING family genes at these SNPs, we examined 52 human bone samples obtained from the femoral neck during surgical hip replacement (27 female, 25 male; 44 European-American and 8 African-American). We observed unidirectional ASE for the IBSP gene, with lower expression of the G allele compared to the A allele for SNP rs17013181. Our data suggest that SNPs within the SIBLING genes may contribute to normal variation of peak BMD. Further studies are necessary to identify the functional variants and to determine the mechanisms underlying the differences in ASE and how these differences relate to the pathophysiology of osteoporosis.

Identification of pathways for bipolar disorder: a meta-analysis.

IMPORTANCE: Genome-wide investigations provide systematic information regarding the neurobiology of psychiatric disorders. OBJECTIVE: To identify biological pathways that contribute to risk for bipolar disorder (BP) using genes with consistent evidence for association in multiple genome-wide association studies (GWAS). DATA SOURCES: Four independent data sets with individual genome-wide data available in July 2011 along with all data sets contributed to the Psychiatric Genomics Consortium Bipolar Group by May 2012. A prior meta-analysis was used as a source for brain gene expression data. STUDY SELECTION: The 4 published GWAS were included in the initial sample. All independent BP data sets providing genome-wide data in the Psychiatric Genomics Consortium were included as a replication sample. DATA EXTRACTION AND SYNTHESIS: We identified 966 genes that contained 2 or more variants associated with BP at P < .05 in 3 of 4 GWAS data sets (n = 12,127 [5253 cases, 6874 controls]). Simulations using 10,000 replicates of these data sets corrected for gene size and allowed the calculation of an empirical P value for each gene; empirically significant genes were entered into a pathway analysis. Each of these pathways was then tested in the replication sample (n = 8396 [3507 cases, 4889 controls]) using gene set enrichment analysis for single-nucleotide polymorphisms. The 226 genes were also compared with results from a meta-analysis of gene expression in the dorsolateral prefrontal cortex. MAIN OUTCOMES AND MEASURES: Empirically significant genes and biological pathways. RESULTS Among 966 genes, 226 were empirically significant (P < .05). Seventeen pathways were overrepresented in analyses of the initial data set. Six of the 17 pathways were associated with BP in both the initial and replication samples: corticotropin-releasing hormone signaling, cardiac ß-adrenergic signaling, phospholipase C signaling, glutamate receptor signaling, endothelin 1 signaling, and cardiac hypertrophy signaling. Among the 226 genes, 9 differed in expression in the dorsolateral prefrontal cortex in patients with BP: CACNA1C, DTNA, FOXP1, GNG2, ITPR2, LSAMP, NPAS3, NCOA2, and NTRK3. CONCLUSIONS AND RELEVANCE: Pathways involved in the genetic predisposition to BP include hormonal regulation, calcium channels, second messenger systems, and glutamate signaling. Gene expression studies implicate neuronal development pathways as well. These results tend to reinforce specific hypotheses regarding BP neurobiology and may provide clues for new approaches to treatment and prevention.

High-resolution genome screen for bone mineral density in heterogeneous stock rat.

We previously demonstrated that skeletal mass, structure, and biomechanical properties vary considerably in heterogeneous stock (HS) rat strains. In addition, we observed strong heritability for several of these skeletal phenotypes in the HS rat model, suggesting that it represents a unique genetic resource for dissecting the complex genetics underlying bone fragility. The purpose of this study was to identify and localize genes associated with bone mineral density in HS rats. We measured bone phenotypes from 1524 adult male and female HS rats between 17 and 20 weeks of age. Phenotypes included dual-energy X-ray absorptiometry (DXA) measurements for bone mineral content and areal bone mineral density (aBMD) for femur and lumbar spine (L3-L5), and volumetric BMD measurements by CT for the midshaft and distal femur, femur neck, and fifth lumbar vertebra (L5). A total of 70,000 polymorphic single-nucleotide polymorphisms (SNPs) distributed throughout the genome were selected from genotypes obtained from the Affymetrix rat custom SNPs array for the HS rat population. These SNPs spanned the HS rat genome with a mean linkage disequilibrium coefficient between neighboring SNPs of 0.95. Haplotypes were estimated across the entire genome for each rat using a multipoint haplotype reconstruction method, which calculates the probability of descent for each genotyped locus from each of the eight founder HS strains. The haplotypes were tested for association with each bone density phenotype via a mixed model with covariate adjustment. We identified quantitative trait loci (QTLs) for BMD phenotypes on chromosomes 2, 9, 10, and 13 meeting a conservative genomewide empiric significance threshold (false discovery rate [FDR] = 5%; p < 3 × 10(-6)). Importantly, most QTLs were localized to very small genomic regions (1-3 megabases [Mb]), allowing us to identify a narrow set of potential candidate genes including both novel genes and genes previously shown to have roles in skeletal development and homeostasis.