Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

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.

Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 x 10(-)(9), odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 x 10(-)(1)(1) (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 x 10(-)(5), permuted threshold for genome-wide significance 7.7 x 10(-)(5). To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Common genetic influences underlie comorbidity of migraine and endometriosis

We examined the co-occurrence of migraine and endometriosis within the largest known collection of families containing multiple women with surgically confirmed endometriosis and in an independent sample of 815 monozygotic and 457 dizygotic female twin pairs. Within the endometriosis families, a significantly increased risk of migrainous headache was observed in women with endometriosis compared to women without endometriosis (odds ratio [OR] 1.57, 95% confidence interval [CI]: 1.12-2.21, P=0.009). Bivariate heritability analyses indicated no evidence for common environmental factors influencing either migraine or endometriosis but significant genetic components for both traits, with heritability estimates of 69 and 49%, respectively. Importantly, a significant additive genetic correlation (r(G) = 0.27, 95% CI: 0.06-0.47) and bivariate heritability (h(2)=0.17, 95% CI: 0.08-0.27) was observed between migraine and endometriosis. Controlling for the personality trait neuroticism made little impact on this association. These results confirm the previously reported comorbidity between migraine and endometriosis and indicate common genetic influences completely explain their co-occurrence within individuals. Given pharmacological treatments for endometriosis typically target hormonal pathways and a number of findings provide support for a relationship between hormonal variations and migraine, hormone-related genes and pathways are highly plausible candidates for both migraine and endometriosis. Therefore, taking into account the status of both migraine and endometriosis may provide a novel opportunity to identify the genes underlying them. Finally, we propose that the analysis of such genetically correlated comorbid traits can increase power to detect genetic risk loci through the use of more specific, homogenous and heritable phenotypes.

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.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Wild dogma II: the role and implications of wild dogma for wild dog management in Australia

The studies of Allen (2011) and Allen et al. (2011) recently examined the methodology underpinning claims that dingoes provide net benefits to biodiversity by suppressing foxes and cats. They found most studies to have design flaws and/or observational methods that preclude valid interpretations from the data, describing most of the current literature as ‘wild dogma’. In this short supplement, we briefly highlight the roles and implications of wild dogma for wild dog management in Australia. We discuss nomenclature, and the influence that unreliable science can have on policy and practice changes related to apex predator management.

Fumarate Hydratase and Succinate Dehydrogenase in Neoplasia

Germline mutations in fumarate hydratase (FH) cause hereditary leiomyomatosis and renal cell cancer (HLRCC). FH is a nuclear encoded enzyme which functions in the Krebs tricarboxylic acid cycle, and homozygous mutation in FH lead to severe developmental defects. Both uterine and cutaneous leiomyomas are components of the HLRCC phenotype. Most of these tumours show loss of the wild-type allele and, also, the mutations reduce FH enzyme activity, which indicate that FH is a tumour suppressor gene. The renal cell cancers associated with HLRCC are of rare papillary type 2 histology. Other genes involved in the Krebs cycle, which are also implicated in neoplasia are 3 of the 4 subunits encoding succinate dehydrogenase (SDH); mutations in SHDB, SDHC, and SDHD predispose to paraganglioma and phaeochromocytoma. Although uterine leiomyomas (or fibroids) are very common, the estimations of affected women ranging from 25% to 77%, not much is known about their genetic background. Cytogenetic studies have revealed that rearrangements involving chromosomes 6, 7, 12 and 14 are most commonly seen in fibroids. Deletions on the long arm of chromosome 7 have been reported to be involved in about 17 to 34 % of leiomyomas and the small commonly deleted region on 7q22 suggests that there might be an underlying tumour suppressor gene in that region. The purpose of this study was to investigate the genetic mechanisms behind the development of tumours associated with HLRCC, both renal cell cancer and uterine fibroids. Firstly, a database search at the Finnish cancer registry was conducted in order to identify new families with early-onset RCC and to test if the family history was compatible with HLRCC. Secondly, sporadic uterine fibroids were tested for deletions on 7q in order to define the minimal deleted 7q-region, followed by mutation analysis of the candidate genes. Thirdly, oligonucleotide chips were utilised to study the global gene expression profiles of uterine fibroids in order to test whether 7q-deletions and FH mutations significantly affected fibroid biology. In the screen for early-onset RCC, 214 families were identified. Subsequently, the pedigrees were constructed and clinical data obtained. One of the index cases (RCC at the age of 28) had a mother who had been diagnosed with a heart tumour, which in further investigation turned out to be a paraganglioma. This lead to an alternative hypothesis that SDH, instead of FH, could be involved. SDHA, SDHB, SDHC and SDHD were sequenced from these individuals; a germline SDHB R27X mutation was detected with loss of the wild-type allele in both tumours. These results suggest that germline mutations in the SDHB gene predispose to early-onset RCC establishing a novel form of hereditary RCC. This has immediate clinical implications in the surveillance of patients suffering from early-onset RCC and phaeochromocytoma/paraganglioma. For the studies on sporadic uterine fibroids, a set of 166 fibroids from 51 individuals were collected. The 7q LOH mapping defined a commonly deleted region of about 3.2 mega bases in 11 of the 166 tumours. The deletion was consistent with previously reported allelotyping studies of leiomyomas and it therefore suggested the presence of a tumour suppressor gene in the deleted region. Furthermore, the high-resolution aCGH-chip analysis refined the deleted region to only 2.79Mb. When combined with previous data, the commonly deleted region was only 2.3Mb. The mutation screening of the known genes within the commonly deleted region did not reveal pathogenic mutations, however. The expression microarray analysis revealed that FH-deficient fibroids, both sporadic and familial, had their distinct gene expression profile as they formed their own group in the unsupervised clustering. On the other hand, the presence or absence of 7q-deletions did not significantly alter the global gene expression pattern of fibroids, suggesting that these two groups do not have different biological backgrounds. Multiple differentially expressed genes were identified between FH wild-type and FH-mutant fibroids, and the most significant increase was seen in the expression of carbohydrate metabolism-related and hypoxia inducible factor (HIF) target genes.

Hybridisation of Australian chickpea cultivars with wild Cicer spp. increases resistance to root-lesion nematodes (Pratylenchus thornei and P. neglectus).

Australian and international chickpea (Cicer arietinum) cultivars and germplasm accessions, and wild annual Cicer spp. in the primary and secondary gene pools, were assessed in glasshouse experiments for levels of resistance to the root-lesion nematodes Pratylenchus thornei and P. neglectus. Lines were grown in replicated experiments in pasteurised soil inoculated with a pure culture of either P. thornei or P. neglectus and the population density of the nematodes in the soil plus roots after 16 weeks growth was used as a measure of resistance. Combined statistical analyses of experiments (nine for P. thornei and four for P. neglectus) were conducted and genotypes were assessed using best linear unbiased predictions. Australian and international chickpea cultivars possessed a similar range of susceptibilities through to partial resistance. Wild relatives from both the primary (C. reticulatum and C. echinospermum) and secondary (C. bijugum) gene pools of chickpea were generally more resistant than commercial chickpea cultivars to either P. thornei or P. neglectus or both. Wild relatives of chickpea have probably evolved to have resistance to endemic root-lesion nematodes whereas modern chickpea cultivars constitute a narrower gene pool with respect to nematode resistance. Resistant accessions of C. reticulatum and C. echinospermum were crossed and topcrossed with desi chickpea cultivars and resistant F(4) lines were obtained. Development of commercial cultivars with the high levels of resistance to P. thornei and P. neglectus in these hybrids will be most valuable for areas of the Australian grain region and other parts of the world where alternating chickpea and wheat crops are the preferred rotation.

Regulation of tumor suppressor protein p53 in cellular stress and tumorigenesis

Functional loss of tumor suppressor protein p53 is a common feature in diverse human cancers. The ability of this protein to sense cellular damage and halt the progression of the cell cycle or direct the cells to apoptosis is essential in preventing tumorigenesis. Tumors having wild-type p53 also respond better to current chemotherapies. The loss of p53 function may arise from TP53 mutations or dysregulation of factors controlling its levels and activity. Probably the most significant inhibitor of p53 function is Mdm2, a protein mediating its degradation and inactivation. Clearly, the maintenance of a strictly controlled p53-Mdm2 route is of great importance in preventing neoplastic transformation. Moreover, impairing Mdm2 function could be a nongenotoxic way to increase p53 levels and activity. Understanding the precise molecular mechanisms behind p53-Mdm2 relationship is thus essential from a therapeutic point of view. The aim of this thesis study was to discover factors affecting the negative regulation of p53 by Mdm2, causing activation of p53 in stressed cells. As a model of cellular damage, we used UVC radiation, inducing a complex cellular stress pathway. Exposure to UVC, as well as to several chemotherapeutic drugs, causes robust transcriptional stress in the cells and leads to activation of p53. By using this model of cellular stress, our goal was to understand how and by which proteins p53 is regulated. Furthermore, we wanted to address whether these pathways affecting p53 function could be altered in human cancers. In the study, two different p53 pathway proteins, nucleophosmin (NPM) and promyelocytic leukemia protein (PML), were found to participate in the p53 stress response following UV stress. Subcellular translocations of these proteins were discovered rapidly after exposure to UV. The alterations in the cellular localizations were connected to transient interactions with p53 and Mdm2, implicating their significance in the regulation of p53 stress response. NPM was shown to control Mdm2-p53 interface and mediate p53 stabilization by blocking the ability of Mdm2 to promote p53 degradation. Furthermore, NPM mediated p53 stabilization upon viral insult. We further detected a connection between cellular pathways of NPM and PML, as PML was found to associate with NPM in UV-radiated cells. The observed temporal UV-induced interactions strongly imply existence of a multiprotein complex participating in the p53 response. In addition, PML controlled the UV response of NPM, its localization and complex formation with chromatin associated factors. The relevance of the UV-promoted interactions was demonstrated in studies in a human leukemia cell line, being under abnormal transcriptional repression due to expression of oncogenic PML-RARa fusion protein. Reversing the leukemic phenotype with a therapeutically significant drug was associated with similar complex formation between p53 and its partners as following UV. In conclusion, this thesis study identifies novel p53 pathway interactions associated with the recovery from UV-promoted as well as oncogenic transcriptional repression.

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.

Common variation in the CYP17A1 and IFIT1 genes on chromosome 10 does not contribute to the risk of endometriosis

Endometriosis is a complex disease involving multiple susceptibility genes and environmental factors. Our previous studies on endometriosis identified a region of significant linkage on chromosome 10q. Two biological candidate genes (CYP17A1 and IFIT1) located on chromosome 10q, have previously been implicated in endometriosis and/or uterine function. We hypothesized that variation in CYP17A1 and/or IFIT1 could contribute to the risk of endometriosis and may account for some of the linkage signal on chromosome 10q. We genotyped 17 single nucleotide polymorphisms (SNPs) in the CYP17A1 and IFIT1 genes including SNP rs743572 previously associated with endometriosis in 768 endometriosis cases and 768 unrelated controls. We found no evidence for association between endometriosis and individual SNPs or SNP haplotypes in CYP17A1 and IFIT1. Common variation in these genes does not appear to be a major contributor to endometriosis susceptibility in our Australian sample.

Corymbia leaf oils, latitude, hybrids and herbivory: a test using common-garden field trials

Foliar oils, particularly monoterpenes, can influence the susceptibility of plants to herbivory. In plants, including eucalypts, monoterpenes are often associated with plant defence. A recent analysis revealed an increase in foliar oil content with increasing latitudinal endemism, and we tested this pattern using three eucalypt taxa comprising a latitudinal replacement cline. We also examined the relative concentrations of two monoterpenes (alpha-pinene and 1,8-cineole), for which meta-analyses also showed latitudinal variation, using hybrids of these three taxa with Corymbia torelliana. These, and pure C. torelliana, were then assessed in common-garden field plots for the abundance and distribution of herbivory by four distinct herbivore taxa. Differing feeding strategies among these herbivores allowed us to test hypotheses regarding heritability of susceptibility and relationships to alpha-pinene and 1,8-cineole. We found no support for an increase in foliar oil content with increasing latitude, nor did our analysis support predictions for consistent variation in alpha-pinene and 1,8-cineole contents with latitude. However, herbivore species showed differential responses to different taxa and monoterpene contents. For example, eriophyid mites, the most monophagous of our censused herbivores, avoided the pure species, but fed on hybrid taxa, supporting hypotheses on hybrid susceptibility. The most polyphagous herbivore (leaf blister sawfly Phylacteophaga froggatti) showed no evidence of response to plant secondary metabolites, while the distribution and abundance patterns of Paropsis atomaria showed some relationship to monoterpene yields.