Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10−8). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10−4, Bonferroni corrected), of which six reached P < 5 × 10−8, including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.

Brain-specific epigenetic markers of schizophrenia

Epigenetics plays a crucial role in schizophrenia susceptibility. In a previous study, we identified over 4500 differentially methylated sites in prefrontal cortex (PFC) samples from schizophrenia patients. We believe this was the first genome-wide methylation study performed on human brain tissue using the Illumina Infinium HumanMethylation450 Bead Chip. To understand the biological significance of these results, we sought to identify a smaller number of differentially methylated regions (DMRs) of more functional relevance compared with individual differentially methylated sites. Since our schizophrenia whole genome methylation study was performed, another study analysing two separate data sets of post-mortem tissue in the PFC from schizophrenia patients has been published. We analysed all three data sets using the bumphunter function found in the Bioconductor package minfi to identify regions that are consistently differentially methylated across distinct cohorts. We identified seven regions that are consistently differentially methylated in schizophrenia, despite considerable heterogeneity in the methylation profiles of patients with schizophrenia. The regions were near CERS3, DPPA5, PRDM9, DDX43, REC8, LY6G5C and a region on chromosome 10. Of particular interest is PRDM9 which encodes a histone methyltransferase that is essential for meiotic recombination and is known to tag genes for epigenetic transcriptional activation. These seven DMRs are likely to be key epigenetic factors in the aetiology of schizophrenia and normal brain neurodevelopment.

A Phenomic Scan of the Norfolk Island Genetic Isolate Identifies a Major Pleiotropic Effect Locus Associated with Metabolic and Renal Disorder Markers

Multiphenotype genome-wide association studies (GWAS) may reveal pleiotropic genes, which would remain undetected using single phenotype analyses. Analysis of large pedigrees offers the added advantage of more accurately assessing trait heritability, which can help prioritise genetically influenced phenotypes for GWAS analysis. In this study we performed a principal component analysis (PCA), heritability (h2) estimation and pedigree-based GWAS of 37 cardiovascular disease -related phenotypes in 330 related individuals forming a large pedigree from the Norfolk Island genetic isolate. PCA revealed 13 components explaining >75% of the total variance. Nine components yielded statistically significant h2 values ranging from 0.22 to 0.54 (P<0.05). The most heritable component was loaded with 7 phenotypic measures reflecting metabolic and renal dysfunction. A GWAS of this composite phenotype revealed statistically significant associations for 3 adjacent SNPs on chromosome 1p22.2 (P<1x10-8). These SNPs form a 42kb haplotype block and explain 11% of the genetic variance for this renal function phenotype. Replication analysis of the tagging SNP (rs1396315) in an independent US cohort supports the association (P = 0.000011). Blood transcript analysis showed 35 genes were associated with rs1396315 (P<0.05). Gene set enrichment analysis of these genes revealed the most enriched pathway was purine metabolism (P = 0.0015). Overall, our findings provide convincing evidence for a major pleiotropic effect locus on chromosome 1p22.2 influencing risk of renal dysfunction via purine metabolism pathways in the Norfolk Island population. Further studies are now warranted to interrogate the functional relevance of this locus in terms of renal pathology and cardiovascular disease risk.

Prolonged hyperinsulinemia affects metabolic signal transduction markers in a tissue specific manner

Insulin dysregulation is common in horses although the mechanisms of metabolic dysfunction are poorly understood. We hypothesized that insulin signaling in striated (cardiac and skeletal) muscle and lamellae may be mediated through different receptors as a result of receptor content, and that transcriptional regulation of downstream signal transduction and glucose transport may also differ between tissues sites during hyperinsulinemia. Archived samples from horses treated with a prolonged insulin infusion or a balanced electrolyte solution were used. All treated horses developed marked hyperinsulinemia and clinical laminitis. Protein expression was compared across tissues for the insulin receptor and insulin-like growth factor 1 receptor (IGF-1R) by immunoblotting. Gene expression of metabolic insulin-signaling markers (insulin receptor substrate 1, Akt2, and glycogen synthase kinase 3 beta [GSK-3β]) and glucose transport (basal glucose transporter 1 and insulin-sensitive glucose transporter 4) was evaluated using real-time reverse transcription polymerase chain reaction. Lamellar tissue contained significantly more IGF-1R protein than skeletal muscle, indicating the potential significance of IGF-1R signaling for this tissue. Gene expression of the selected markers of insulin signaling and glucose transport in skeletal muscle and lamellar tissues was unaffected by prolonged hyperinsulinemia. In contrast, the significant upregulation of Akt2, GSK-3β, GLUT1, and GLUT4 gene expression in cardiac tissue suggested that the prolonged hyperinsulinemia induced an increase in insulin sensitivity and a transcriptional activation of glucose transport. Responses to insulin are tissue-specific, and extrapolation of data across tissue sites is inappropriate.

No evidence for genetic association with the let-7 microRNA-binding site or other common KRAS variants in risk of endometriosis

STUDY QUESTION Is there a contribution of the minor allele at the KRAS single nucleotide polymorphism (SNP) rs61764370 in the let-7 microRNA-binding site to endometriosis risk? SUMMARY ANSWER We found no evidence for association between endometriosis risk and rs61764370 or any other SNPs in KRAS. WHAT IS KNOWN ALREADY The rs61764370 SNP in the 3' untranslated region of the KRAS gene is predicted to disrupt a complementary binding site (LCS6) for the let-7 microRNA, and was recently reported to be at a high frequency (31%) in 132 women of varying ancestry with endometriosis compared with frequencies in a database of population controls (up to 7.6% depending on ancestry), suggesting a strong effect of this KRAS SNP in the aetiology of endometriosis. STUDY DESIGN, SIZE AND DURATION This was a case-control study with a total of 11 206 subjects. The study was performed between February 2012 and July 2012. PARTICIPANTS/MATERIALS, SETTINGAND METHODS We first investigated a possible association between common markers in KRAS and endometriosis risk from our genome-wide association (GWA) data in 3194 surgically confirmed endometriosis cases and 7060 controls of European ancestry. Although rs61764370 was not genotyped on the GWA arrays, five SNPs typed in the study were highly correlated with this variant. The rs61764370 and two SNPs highly correlated with rs61764370 were then genotyped in 933 endometriosis cases and 952 controls using the Sequenom MassARRAY platform. MAIN RESULTS AND THE ROLE OF CHANCE There was no evidence for an association between rs61764370 and endometriosis risk P = 0.411 and odds ratio = 1.10 (95% confidence intervals: 0.88-1.36). We also found no evidence for an association between the highly correlated SNP rs17387019 and endometriosis. Their minor allele frequencies in cases and controls were of 0.087-0.091 similar to the population frequency reported previously for this variant in controls. Analyses of endometriosis cases with revised American Fertility Society stage III/IV disease also showed no evidence for an association between these SNPs and endometriosis risk. LIMITATIONS AND REASONS FOR CAUTION The GWA and genotyped data sets were not independent since individuals and cases from some families overlap. Controls in our GWA study were not screened for endometriosis. WIDER IMPLICATIONS OF THE FINDINGS The key SNP, rs61764370, was genotyped in a subset of samples. Our results do not support the suggestion that carrying the minor allele at rs61764370 contributes to a significant number of endometriosis cases and rs61764370 is, therefore, unlikely to be a useful marker of endometriosis risk. STUDY FUNDING/COMPETING INTEREST(S) The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Common variants in the trichohyalin gene are associated with straight hair in Europeans

Hair morphology is highly differentiated between populations and among people of European ancestry. Whereas hair morphology in East Asian populations has been studied extensively, relatively little is known about the genetics of this trait in Europeans. We performed a genome-wide association scan for hair morphology (straight, wavy, curly) in three Australian samples of European descent. All three samples showed evidence of association implicating the Trichohyalin gene (TCHH), which is expressed in the developing inner root sheath of the hair follicle, and explaining approximately 6% of variance (p=1.5x10(-31)). These variants are at their highest frequency in Northern Europeans, paralleling the distribution of the straight-hair EDAR variant in Asian populations.

A reanalysis of 409 European-Ancestry and African American schizophrenia pedigrees reveals significant linkage to 8p23.3 with evidence of locus heterogeneity

The detection and replication of schizophrenia risk loci can require substantial sample sizes, which has prompted various collaborative efforts for combining multiple samples. However, pooled samples may comprise sub-samples with substantial population genetic differences, including allele frequency differences. We investigated the impact of population differences via linkage reanalysis of Molecular Genetics of Schizophrenia 1 (MGS1) affected sibling-pair data, comprising two samples of distinct ancestral origin: European (EA: 263 pedigrees) and African-American (AA: 146 pedigrees). To exploit the linkage information contained within these distinct continental samples, we performed separate analyses of the individual samples, allowing for within-sample locus heterogeneity, and the pooled sample, allowing for both within-sample and between-sample heterogeneity. Significance levels, corrected for the multiple tests, were determined empirically. For all suggestive peaks, stronger linkage evidence was obtained in either the EA or AA sample than the combined sample, regardless of how heterogeneity was modeled for the latter. Notably, we report genomewide significant linkage of schizophrenia to 8p23.3 and evidence for a second, independent susceptibility locus, reaching suggestive linkage, 29 cM away on 8p21.3. We also detected suggestive linkage on chromosomes 5p13.3 and 7q36.2. Many regions showed pronounced differences in the extent of linkage between the EA and AA samples. This reanalysis highlights the potential impact of population differences upon linkage evidence in pooled data and demonstrates a useful approach for the analysis of samples drawn from distinct continental groups.

Further evidence for linkage of Gilles de la Tourette syndrome (GTS) susceptibility loci on chromosomes 2p11, 8q22 and 11q23-24 in South African Afrikaners

Utilizing DNA samples from 91 Afrikaner nuclear families with one or more affected children, five genomic regions on chromosomes 2p, 8q, 11q, 20q, and 21q that gave evidence for association with GTS in previous case-control association studies were investigated for linkage and association with GTS. Highly polymorphic markers with mean heterozygosity of 0.77 were typed and resulting genotypes evaluated using single marker transmission disequilibrium (TDT), single marker haplotype relative risk (HRR), and multi-marker "extended" TDT and HRR methods. Single marker TDT analysis showed evidence for linkage or association, with p-values near 0.05, for markers D2S139, GATA28F12, and D11S1377 on chromosomes 2p11, 8q22 and 11q23-24, respectively. Extended, two-locus TDT and HRR analysis provided further evidence for linkage or association on chromosome 2 with p-values of 0.007 and 0.025, and chromosome 8 with p-values of 0.059 and 0.013, respectively. These results provide important additional evidence for the location of GTS susceptibility loci.

Exclusion of angiotensinogen gene in molecular basis of human hypertension: Sibpair linkage and association analyses in Australian anglo-caucasians

Linkage with essential hypertension has been claimed for a microsatellite marker near the angiotensinogen gene (AGT; chromosome 1q42), as has association for the AGT variants M235T, G(-6)A and A(-20)C. To more rigorously evaluate AGT as a candidate gene for hypertension we performed sibpair analysis with multiple microsatellite markers surrounding this locus and using more sophisticated analysis programs. We also performed an association study of the AGT variants in unrelated subjects with a strong family history (two affected parents). For the linkage study, single and multiplex polymerase chain reaction (PCRs) and automated genescan analysis were conducted on DNA from 175 Australian Anglo-Celtic Caucasian hypertensives for the following markers: D1S2880-(2.1 cM)-D1S213-(2.8 cM)-D1S251-(6.5 cM)-AGT-(2.0 cM) -D1S235. Statistical evaluation of genotype data by nonparametric methods resulted in the following scores: Single-point analysis - SPLINK, P > 0.18; APM method, P > 0.25; ASPEX, MLOD < 0.28; SIB-PAIR, P > 0. 24; Multipoint analysis - MAPMAKER/SIBS, MLOD < 0.24; GENEHUNTER, P > 0.35. Exclusion scores of Lod -4.1 to -5.1 were obtained for these markers using MAPMAKER/SIBS for a lambda(s) of 1.6. The association study of G(-6)A, A(-20)C and M235T variants in 111 hypertensives with strong family history and 190 normotensives with no family history showed significant linkage disequilibrium between particular haplotypes, but we could find no association with hypertension. The present study therefore excludes AGT in the etiology of hypertension, at least in the population of Australian Anglo-Celtic Caucasians studied.

Research update: LANDMark study: Longitudinal assessment of neuropathy in diabetes using novel ophthalmic MARKers

The eye is a simple, non-invasive location for screening, diagnosing and follow up of diabetic peripheral neuropathy.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: - 1) priming site polymorphism in the template leads to inferior or erratic amplification; - 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions, and; - 3) at least occasionally, a PCR primer straddles an exon–intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3′ end of coding sequence and 3′ UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.

The study of genetic diversity and population structure of the Korean fleshy shrimp, Fenneropenaeus chinensis, using newly developed microsatellite markers

The fleshy shrimp, Fenneropenaeus chinensis, is the family of Penaeidae and one of the most economically important marine culture species in Korea. However, its genetic characteristics have never been studied. In this study, a total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n = 60), Beopseongpo (BSP, n = 60), Chaesukpo (CSP, n = 60), and Cheonsuman (CSM, n = 60). Genetic variability and the relationships among four wild F. chinensis populations were analyzed using 13 newly developed microsatellite loci. Relatively high levels of genetic variability (mean allelic richness = 16.87; mean heterozygosity = 0.845) were found among localities. Among the 52 population loci, 13 showed significant deviation from the Hardy–Weinberg equilibrium. Neighbor-joining, principal coordinate, and molecular variance analyses revealed the presence of three subpopulations (NRD, CSM, BSP and CSP), which was consistent with clustering based on genetic distance. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These genetic variability and differentiation results of the four wild populations can be applied for future genetic improvement using selective breeding and to design suitable management guidelines for Korean F. chinensis culture.

Cell surface heparan sulfate proteoglycans as novel markers of human neural stem cell fate determination

Multipotent neural stem cells (NSCs) provide a model to investigate neurogenesis and develop mechanisms of cell transplantation. In order to define improved markers of stemness and lineage specificity, we examined self-renewal and multi-lineage markers during long-term expansion and under neuronal and astrocyte differentiating conditions in human ESC-derived NSCs (hNSC H9 cells). In addition, with proteoglycans ubiquitous to the neural niche, we also examined heparan sulfate proteoglycans (HSPGs) and their regulatory enzymes. Our results demonstrate that hNSC H9 cells maintain self-renewal and multipotent capacity during extended culture and express HS biosynthesis enzymes and several HSPG core protein syndecans (SDCs) and glypicans (GPCs) at a high level. In addition, hNSC H9 cells exhibit high neuronal and a restricted glial differentiative potential with lineage differentiation significantly increasing expression of many HS biosynthesis enzymes. Furthermore, neuronal differentiation of the cells upregulated SDC4, GPC1, GPC2, GPC3 and GPC6 expression with increased GPC4 expression observed under astrocyte culture conditions. Finally, downregulation of selected HSPG core proteins altered hNSC H9 cell lineage potential. These findings demonstrate an involvement for HSPGs in mediating hNSC maintenance and lineage commitment and their potential use as novel markers of hNSC and neural cell lineage specification.

Ibuprofen ingestion does not affect markers of post-exercise muscle inflammation

Purpose: We investigated if oral ingestion of ibuprofen influenced leucocyte recruitment and infiltration following an acute bout of traditional resistance exercise Methods: Sixteen male subjects were divided into two groups that received the maximum over-the-counter dose of ibuprofen (1200mg d−1) or a similarly administered placebo following lower body resistance exercise. Muscle biopsies were taken from m.vastus lateralis and blood serum samples were obtained before and immediately after exercise, and at 3 and 24 h after exercise. Muscle cross-sections were stained with antibodies against neutrophils (CD66b and MPO) and macrophages (CD68). Muscle damage was assessed via creatine kinase and myoglobin in blood serum samples, and muscle soreness was rated on a ten-point pain scale. Results: The resistance exercise protocol stimulated a significant increase in the number of CD66b+ and MPO+ cells when measured 3 h post exercise. Serum creatine kinase, myoglobin and subjective muscle soreness all increased post-exercise. Muscle leucocyte infiltration, creatine kinase, myoglobin and subjective muscle soreness were unaffected by ibuprofen treatment when compared to placebo. There was also no association between increases in inflammatory leucocytes and any other marker of cellular muscle damage. Conclusion: Ibuprofen administration had no effect on the accumulation of neutrophils, markers of muscle damage or muscle soreness during the first 24 h of post-exercise muscle recovery.

Prediction of individual genetic risk to prostate cancer using a polygenic score

BACKGROUND Polygenic risk scores comprising established susceptibility variants have shown to be informative classifiers for several complex diseases including prostate cancer. For prostate cancer it is unknown if inclusion of genetic markers that have so far not been associated with prostate cancer risk at a genome-wide significant level will improve disease prediction. METHODS We built polygenic risk scores in a large training set comprising over 25,000 individuals. Initially 65 established prostate cancer susceptibility variants were selected. After LD pruning additional variants were prioritized based on their association with prostate cancer. Six-fold cross validation was performed to assess genetic risk scores and optimize the number of additional variants to be included. The final model was evaluated in an independent study population including 1,370 cases and 1,239 controls. RESULTS The polygenic risk score with 65 established susceptibility variants provided an area under the curve (AUC) of 0.67. Adding an additional 68 novel variants significantly increased the AUC to 0.68 (P = 0.0012) and the net reclassification index with 0.21 (P = 8.5E-08). All novel variants were located in genomic regions established as associated with prostate cancer risk. CONCLUSIONS Inclusion of additional genetic variants from established prostate cancer susceptibility regions improves disease prediction. Prostate 75:1467–1474, 2015. © 2015 Wiley Periodicals, Inc.