High-density fine-mapping of a chromosome 10q26 linkage peak suggests association between endometriosis and variants close to CYP2C19

OBJECTIVE To refine a previously reported linkage peak for endometriosis on chromosome 10q26, and conduct follow-up analyses and a fine-mapping association study across the region to identify new candidate genes for endometriosis. DESIGN Case-control study. SETTING Academic research. PATIENT(S) Cases=3,223 women with surgically confirmed endometriosis; controls=1,190 women without endometriosis and 7,060 population samples. INTERVENTION(S) Analysis of 11,984 single nucleotide polymorphisms on chromosome 10. MAIN OUTCOME MEASURE(S) Allele frequency differences between cases and controls. RESULT(S) Linkage analyses on families grouped by endometriosis symptoms (primarily subfertility) provided increased evidence for linkage (logarithm of odds score=3.62) near a previously reported linkage peak. Three independent association signals were found at 96.59 Mb (rs11592737), 105.63 Mb (rs1253130), and 124.25 Mb (rs2250804). Analyses including only samples from linkage families supported the association at all three regions. However, only rs11592737 in the cytochrome P450 subfamily C (CYP2C19) gene was replicated in an independent sample of 2,079 cases and 7,060 population controls. CONCLUSION(S) The role of the CYP2C19 gene in conferring risk for endometriosis warrants further investigation.

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

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.

Meta-analysis of 28,141 individuals identifies common variants within five new loci that influence uric acid concentrations

Elevated serum uric acid levels cause gout and are a risk factor for cardiovascular disease and diabetes. To investigate the polygenetic basis of serum uric acid levels, we conducted a meta-analysis of genome-wide association scans from 14 studies totalling 28,141 participants of European descent, resulting in identification of 954 SNPs distributed across nine loci that exceeded the threshold of genome-wide significance, five of which are novel. Overall, the common variants associated with serum uric acid levels fall in the following nine regions: SLC2A9 (p = 5.2x10(-201)), ABCG2 (p = 3.1x10(-26)), SLC17A1 (p = 3.0x10(-14)), SLC22A11 (p = 6.7x10(-14)), SLC22A12 (p = 2.0x10(-9)), SLC16A9 (p = 1.1x10(-8)), GCKR (p = 1.4x10(-9)), LRRC16A (p = 8.5x10(-9)), and near PDZK1 (p = 2.7x10(-9)). Identified variants were analyzed for gender differences. We found that the minor allele for rs734553 in SLC2A9 has greater influence in lowering uric acid levels in women and the minor allele of rs2231142 in ABCG2 elevates uric acid levels more strongly in men compared to women. To further characterize the identified variants, we analyzed their association with a panel of metabolites. rs12356193 within SLC16A9 was associated with DL-carnitine (p = 4.0x10(-26)) and propionyl-L-carnitine (p = 5.0x10(-8)) concentrations, which in turn were associated with serum UA levels (p = 1.4x10(-57) and p = 8.1x10(-54), respectively), forming a triangle between SNP, metabolites, and UA levels. Taken together, these associations highlight additional pathways that are important in the regulation of serum uric acid levels and point toward novel potential targets for pharmacological intervention to prevent or treat hyperuricemia. In addition, these findings strongly support the hypothesis that transport proteins are key in regulating serum uric acid levels.

Sequence variants in three loci influence monocyte counts and erythrocyte volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p=8.9x10(-14)). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p=1.2x10(-9), 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p=5.3x10(-9), 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values=0.001, 9.9x10(-5), and 7x10(-5), respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.

Common variants in TMPRSS6 are associated with iron status and erythrocyte volume

We report a genome-wide association study to iron status. We identify an association of SNPs in TPMRSS6 to serum iron (rs855791, combined P = 1.5 x 10(-20)), transferrin saturation (combined P = 2.2 x 10(-23)) and erythrocyte mean cell volume (MCV, combined P = 1.1 x 10(-10)). We also find suggestive evidence of association with blood hemoglobin levels (combined P = 5.3 x 10(-7)). These findings demonstrate the involvement of TMPRSS6 in control of iron homeostasis and in normal erythropoiesis.

Polymorphisms in the vascular endothelial growth factor gene and the risk of familial endometriosis

Vascular endothelial growth factor (VEGF) is an endothelial cell-specific angiogenic protein suspected to be involved in the pathogenesis of endometriosis by establishing a new blood supply to the human exfoliated endometrium. Several transcription factor-binding sites are found in the VEGF 5'-untranslated region and variation within the region increases the transcriptional activity. Six previous studies which tested between one and three single nucleotide polymorphisms (SNPs) in samples comprising 105-215 cases and 100-219 controls have produced conflicting evidence for association between the SNPs in the VEGF region and endometriosis. To further investigate the reported association between VEGF variants and endometriosis, we tested the four VEGF polymorphisms (-2578 A/C, rs699947; -460 T/C, rs833061; +405 G/C, rs2010963 and +936 C/T, rs3025039) in a large Australian sample of 958 familial endometriosis cases and 959 controls. We also conducted a literature-based review of all relevant association studies of these VEGF SNPs in endometriosis and performed a meta-analysis. There was no evidence for association between endometriosis and the VEGF polymorphisms genotyped in our study. Combined association results from a meta-analysis did not provide any evidence for either genotypic or allelic association with endometriosis. Our detailed review and meta-analysis of the VEGF polymorphisms suggests that genotyping assay problems may underlie the previously reported associations between VEGF variants and endometriosis.

Susceptibility variants for male-pattern baldness on chromosome 20p11

We carried out a genome-wide association study in 296 individuals with male-pattern baldness (androgenetic alopecia) and 347 controls. We then investigated the 30 best SNPs in an independent replication sample and found highly significant association for five SNPs on chromosome 20p11 (rs2180439 combined P = 2.7 x 10(-15)). No interaction was detected with the X-chromosomal androgen receptor locus, suggesting that the 20p11 locus has a role in a yet-to-be-identified androgen-independent pathway.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.

Multi-species sequence comparison reveals conservation of ghrelin gene-derived splice variants encoding a truncated ghrelin peptide

The peptide hormone ghrelin is a potent orexigen produced predominantly in the stomach. It has a number of other biological actions, including roles in appetite stimulation, energy balance, the stimulation of growth hormone release and the regulation of cell proliferation. Recently, several ghrelin gene splice variants have been described. Here, we attempted to identify conserved alternative splicing of the ghrelin gene by cross-species sequence comparisons. We identified a novel human exon 2-deleted variant and provide preliminary evidence that this splice variant and in1-ghrelin encode a C-terminally truncated form of the ghrelin peptide, termed minighrelin. These variants are expressed in humans and mice, demonstrating conservation of alternative splicing spanning 90 million years. Minighrelin appears to have similar actions to full-length ghrelin, as treatment with exogenous minighrelin peptide stimulates appetite and feeding in mice. Forced expression of the exon 2-deleted preproghrelin variant mirrors the effect of the canonical preproghrelin, stimulating cell proliferation and migration in the PC3 prostate cancer cell line. This is the first study to characterise an exon 2-deleted preproghrelin variant and to demonstrate sequence conservation of ghrelin gene-derived splice variants that encode a truncated ghrelin peptide. This adds further impetus for studies into the alternative splicing of the ghrelin gene and the function of novel ghrelin peptides in vertebrates.

Molecular dynamics simulations of CXCL-8 and its interactions with a receptor peptide, heparin fragments, and sulfated linked cyclitols

CXCL-8 (Interleukin 8) is a CXC chemokine with a central role in the human immune response. We have undertaken extensive in silico analyses to elucidate the interactions of CXCL-8 with its various binding partners, which are crucial for its biological function. Sequence and structure analyses showed that residues in the thirdq β-sheet and basic residues in the heparin binding site are highly variable, while residues in the second β-sheet are highly conserved. Molecular dynamics simulations in aqueous solution of dimeric CXCL-8 have been performed with starting geometries from both X-ray and NMR structures showed shearing movements between the two antiparallel C-terminal helices. Dynamic conservation analyses of these simulations agreed with experimental data indicating that structural differences between the two structures at quaternary level arise from changes in the secondary structure of the N-terminal loop, the 310-helix, the 30s, 40s, and 50s loops and the third β-sheet, resulting in a different interhelical separation. Nevertheless, the observation of these different states indicates that CXCL-8 has the potential to undergo conformational changes, and it seems likely that this feature is relevant to the mode of binding of glycosaminoglycan (GAG) mimetics such as cyclitols. Simulations of the receptor peptide fragment−CXCL-8 complex identified several specific interactions of the receptor peptide with CXCL-8 that could be exploited in the structure-based design of competitive peptides and nonpeptidic molecules targeting CXCL-8 for combating inflammatory diseases. Simulations of the CXCL-8 dimer complexed with a 24-mer heparin fragment and of the CXCL-8−receptor peptide complex revealed that Arg60, Lys64, and Arg68 in the dimer bind to cyclitols in a horseshoe pattern, defining a region which is spatially distinct from the receptor binding site. There appears to be an optimum number of sulfates and an optimum length of alkyl spacers required for the interaction of cyclitol inhibitors with the dimeric form of CXCL-8. Calculation of the binding affinities of cyclitol inhibitors reflected satisfactorily the ranking of experimentally determined inhibitory potencies. The findings of these molecular modeling studies will help in the search for inhibitors which can modulate various CXCL-8 biological activities and serve as an excellent model system to study CXC-inhibitor interactions.

Disease-associated polymorphisms in ERAP1 do not alter endoplasmic reticulum stress in patients with ankylosing spondylitis

The mechanism by which human leukocyte antigen B27 (HLA-B27) contributes to ankylosing spondylitis (AS) remains unclear. Genetic studies demonstrate that association with and interaction between polymorphisms of endoplasmic reticulum aminopeptidase 1 (ERAP1) and HLA-B27 influence the risk of AS. It has been hypothesised that ERAP1-mediated HLA-B27 misfolding increases endoplasmic reticulum (ER) stress, driving an interleukin (IL) 23-dependent, pro-inflammatory immune response. We tested the hypothesis that AS-risk ERAP1 variants increase ER-stress and concomitant pro-inflammatory cytokine production in HLA-B27 + but not HLA-B27-AS patients or controls. Forty-nine AS cases and 22 healthy controls were grouped according to HLA-B27 status and AS-associated ERAP1 rs30187 genotypes: HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective, HLA-B27-ERAP1 risk and HLA-B27-ERAP1 protective. Expression levels of ER-stress markers GRP78 (8 kDa glucose-regulated protein), CHOP (C/EBP-homologous protein) and inflammatory cytokines were determined in peripheral blood mononuclear cell and ileal biopsies. We found no differences in ER-stress gene expression between HLA-B27 + and HLA-B27-cases or healthy controls, or between cases or controls stratified by carriage of ERAP1 risk or protective alleles in the presence or absence of HLA-B27. No differences were observed between expression of IL17A or TNF (tumour necrosis factor) in HLA-B27 + ERAP1 risk, HLA-B27 + ERAP1 protective and HLA-B27-ERAP1 protective cases. These data demonstrate that aberrant ERAP1 activity and HLA-B27 carriage does not alter ER-stress levels in AS, suggesting that ERAP1 and HLA-B27 may influence disease susceptibility through other mechanisms. © 2015 Macmillan Publishers Limited.

Association of FOXE1 Polyalanine Repeat Region with Papillary Thyroid Cancer

CONTEXT: Polyalanine tract variations in transcription factors have been identified for a wide spectrum of developmental disorders. The thyroid transcription factor forkhead factor E1 (FOXE1) contains a polymorphic polyalanine tract with 12-22 alanines. Single-nucleotide polymorphisms (SNP) close to this locus are associated with papillary thyroid cancer (PTC), and a strong linkage disequilibrium block extends across this region. OBJECTIVE: The objective of the study was to assess whether the FOXE1 polyalanine repeat region was associated with PTC and to assess the effect of polyalanine repeat region variants on protein expression, DNA binding, and transcriptional function on FOXE1-responsive promoters. DESIGN: This was a case-control study. SETTING: The study was conducted at a tertiary referral hospital. PATIENTS AND METHODS: The FOXE1 polyalanine repeat region and tag SNP were genotyped in 70 PTC, with a replication in a further 92 PTC, and compared with genotypes in 5767 healthy controls (including 5667 samples from the Wellcome Trust Case Control Consortium). In vitro studies were performed to examine the protein expression, DNA binding, and transcriptional function for FOXE1 variants of different polyalanine tract lengths. RESULTS: All the genotyped SNP were in tight linkage disequilibrium, including the FOXE1 polyalanine repeat region. We confirmed the strong association of rs1867277 with PTC (overall P = 1 × 10(-7), odds ratio 1.84, confidence interval 1.31-2.57). rs1867277 was in tight linkage disequilibrium with the FOXE1 polyalanine repeat region (r(2) = 0.95). FOXE1(16Ala) was associated with PTC with an odds ratio of 2.23 (confidence interval 1.42-3.50; P = 0.0005). Functional studies in vitro showed that FOXE1(16Ala) was transcriptionally impaired compared with FOXE1(14Ala), which was not due to differences in protein expression or DNA binding. CONCLUSIONS: We have confirmed the previous association of FOXE1 with PTC. Our data suggest that the coding polyalanine expansion in FOXE1 may be responsible for the observed association between FOXE1 and PTC.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10-8) level: 14q24.2 (rs227425, P-value 3.98 × 10-13, SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10-9, CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n 5 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis. © The Author 2013. Published by Oxford University Press.