Using genome-Wwide complex trait analysis to quantify 'missing heritability' in Parkinson's disease

Genome-wide association studies (GWASs) have been successful at identifying single-nucleotide polymorphisms (SNPs) highly associated with common traits; however, a great deal of the heritable variation associated with common traits remains unaccounted for within the genome. Genome-wide complex trait analysis (GCTA) is a statistical method that applies a linear mixed model to estimate phenotypic variance of complex traits explained by genome-wide SNPs, including those not associated with the trait in a GWAS. We applied GCTA to 8 cohorts containing 7096 case and 19 455 control individuals of European ancestry in order to examine the missing heritability present in Parkinson's disease (PD). We meta-analyzed our initial results to produce robust heritability estimates for PD types across cohorts. Our results identify 27% (95% CI 17-38, P = 8.08E - 08) phenotypic variance associated with all types of PD, 15% (95% CI -0.2 to 33, P = 0.09) phenotypic variance associated with early-onset PD and 31% (95% CI 17-44, P = 1.34E - 05) phenotypic variance associated with late-onset PD. This is a substantial increase from the genetic variance identified by top GWAS hits alone (between 3 and 5%) and indicates there are substantially more risk loci to be identified. Our results suggest that although GWASs are a useful tool in identifying the most common variants associated with complex disease, a great deal of common variants of small effect remain to be discovered. © Published by Oxford University Press 2012.

Common variants near ATM are associated with glycemic response to metformin in type 2 diabetes

Metformin is the most commonly used pharmacological therapy for type 2 diabetes. We report a genome-wide association study for glycemic response to metformin in 1,024 Scottish individuals with type 2 diabetes with replication in two cohorts including 1,783 Scottish individuals and 1,113 individuals from the UK Prospective Diabetes Study. In a combined meta-analysis, we identified a SNP, rs11212617, associated with treatment success (n = 3,920, P = 2.9 P×-9, odds ratio = 1.35, 95% CI 1.22-1.49) at a locus containing ATM, the ataxia telangiectasia mutated gene. In a rat hepatoma cell line, inhibition of ATM with KU-55933 attenuated the phosphorylation and activation of AMP-activated protein kinase in response to metformin. We conclude that ATM, a gene known to be involved in DNA repair and cell cycle control, plays a role in the effect of metformin upstream of AMP-activated protein kinase, and variation in this gene alters glycemic response to metformin. © 2011 Nature America, Inc. All rights reserved.

Dissection of the genetics of Parkinson's disease identifies an additional association 5' of SNCA and multiple associated haplotypes at 17q21

We performed a genome-wide association study (GWAS) in 1705 Parkinson's disease (PD) UK patients and 5175 UK controls, the largest sample size so far for a PD GWAS. Replication was attempted in an additional cohort of 1039 French PD cases and 1984 controls for the 27 regions showing the strongest evidence of association (P < 10 4). We replicated published associations in the 4q22/SNCA and 17q21/MAPT chromosome regions (P < 10 10) and found evidence for an additional independent association in 4q22/SNCA.A detailed analysis of the haplotype structure at 17q21 showed that there are three separate risk groups within this region. We found weak but consistent evidence of association for common variants located in three previously published associated regions (4p15/BST1, 4p16/GAK and 1q32/PARK16). We found no support for the previously reported SNP association in 12q12/LRRK2. We also found an association of the two SNPs in 4q22/SNCA with the age of onset of the disease. © The Author 2010. Published by Oxford University Press.

Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis

Multiple sclerosis is a common disease of the central nervous system in which the interplay between inflammatory and neurodegenerative processes typically results in intermittent neurological disturbance followed by progressive accumulation of disability. Epidemiological studies have shown that genetic factors are primarily responsible for the substantially increased frequency of the disease seen in the relatives of affected individuals, and systematic attempts to identify linkage in multiplex families have confirmed that variation within the major histocompatibility complex (MHC) exerts the greatest individual effect on risk. Modestly powered genome-wide association studies (GWAS) have enabled more than 20 additional risk loci to be identified and have shown that multiple variants exerting modest individual effects have a key role in disease susceptibility. Most of the genetic architecture underlying susceptibility to the disease remains to be defined and is anticipated to require the analysis of sample sizes that are beyond the numbers currently available to individual research groups. In a collaborative GWAS involving 9,772 cases of European descent collected by 23 research groups working in 15 different countries, we have replicated almost all of the previously suggested associations and identified at least a further 29 novel susceptibility loci. Within the MHC we have refined the identity of the HLA-DRB1 risk alleles and confirmed that variation in the HLA-A gene underlies the independent protective effect attributable to the class I region. Immunologically relevant genes are significantly overrepresented among those mapping close to the identified loci and particularly implicate T-helper-cell differentiation in the pathogenesis of multiple sclerosis. © 2011 Macmillan Publishers Limited. All rights reserved.

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10-8 in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10-6 overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.

A genome-wide asociation study identifies new psoriasis susceptibility loci and an interaction betwEn HLA-C and ERAP1

To identify new susceptibility loci for psoriasis, we undertOk a genome-wide asociation study of 594,224 SNPs in 2,622 individuals with psoriasis and 5,667 controls. We identified asociations at eight previously unreported genomic loci. Seven loci harbored genes with recognized iMune functions (IL28RA, REL, IFIH1, ERAP1, TRAF3IP2, NFKBIA and TYK2). These asociations were replicated in 9,079 European samples (six loci with a combined P < 5-10 -8 and two loci with a combined P < 5-10-7). We also report compeLing evidence for an interaction betwEn the HLA-C and ERAP1 loci (combined P = 6.95-10-6). ERAP1 plays an important role in MHC claS I peptide proceSing. ERAP1 variants only influenced psoriasis susceptibility in individuals carrying the HLA-C risk aLele. Our findings implicate pathways that integrate epidermal barrier dysfunction with iNate and adaptive iMune dysregulation in psoriasis pathogenesis.

Genome-wide association study of ulcerative colitis identifies three new susceptibility loci, including the HNF4A region

Ulcerative colitis is a common form of inflammatory bowel disease with a complex etiology. As part of the Wellcome Trust Case Control Consortium 2, we performed a genome-wide association scan for ulcerative colitis in 2,361 cases and 5,417 controls. Loci showing evidence of association at P 1 × 10 5 were followed up by genotyping in an independent set of 2,321 cases and 4,818 controls. We find genome-wide significant evidence of association at three new loci, each containing at least one biologically relevant candidate gene, on chromosomes 20q13 (HNF4A; P = 3.2 × 10 17), 16q22 (CDH1 and CDH3; P = 2.8 × 10 8) and 7q31 (LAMB1; P = 3.0 × 10 8). Of note, CDH1 has recently been associated with susceptibility to colorectal cancer, an established complication of longstanding ulcerative colitis. The new associations suggest that changes in the integrity of the intestinal epithelial barrier may contribute to the pathogenesis of ulcerative colitis. © 2009 Nature America, Inc. All rights reserved.

Combined analysis of three whole genome linkage scans for Ankylosing Spondylitis

Objective. Ankylosing spondylitis (AS) is a debilitating chronic inflammatory condition with a high degree of familiality (λs=82) and heritability (>90%) that primarily affects spinal and sacroiliac joints. Whole genome scans for linkage to AS phenotypes have been conducted, although results have been inconsistent between studies and all have had modest sample sizes. One potential solution to these issues is to combine data from multiple studies in a retrospective meta-analysis. Methods: The International Genetics of Ankylosing Spondylitis Consortium combined data from three whole genome linkage scans for AS (n=3744 subjects) to determine chromosomal markers that show evidence of linkage with disease. Linkage markers typed in different centres were integrated into a consensus map to facilitate effective data pooling. We performed a weighted meta-analysis to combine the linkage results, and compared them with the three individual scans and a combined pooled scan. Results: In addition to the expected region surrounding the HLA-B27 gene on chromosome 6, we determined that several marker regions showed significant evidence of linkage with disease status. Regions on chromosome 10q and 16q achieved 'suggestive' evidence of linkage, and regions on chromosomes 1q, 3q, 5q, 6q, 9q, 17q and 19q showed at least nominal linkage in two or more scans and in the weighted meta-analysis. Regions previously associated with AS on chromosome 2q (the IL-1 gene cluster) and 22q (CYP2D6) exhibited nominal linkage in the meta-analysis, providing further statistical support for their involvement in susceptibility to AS. Conclusion: These findings provide a useful guide for future studies aiming to identify the genes involved in this highly heritable condition. . Published by on behalf of the British Society for Rheumatology.

The v-MFG test: Investigating maternal, offspring and maternal-fetal genetic incompatibility effects on disease and viability

The MFG test is a family-based association test that detects genetic effects contributing to disease in offspring, including offspring allelic effects, maternal allelic effects and MFG incompatibility effects. Like many other family-based association tests, it assumes that the offspring survival and the offspring-parent genotypes are conditionally independent provided the offspring is affected. However, when the putative disease-increasing locus can affect another competing phenotype, for example, offspring viability, the conditional independence assumption fails and these tests could lead to incorrect conclusions regarding the role of the gene in disease. We propose the v-MFG test to adjust for the genetic effects on one phenotype, e.g., viability, when testing the effects of that locus on another phenotype, e.g., disease. Using genotype data from nuclear families containing parents and at least one affected offspring, the v-MFG test models the distribution of family genotypes conditional on offspring phenotypes. It simultaneously estimates genetic effects on two phenotypes, viability and disease. Simulations show that the v-MFG test produces accurate genetic effect estimates on disease as well as on viability under several different scenarios. It generates accurate type-I error rates and provides adequate power with moderate sample sizes to detect genetic effects on disease risk when viability is reduced. We demonstrate the v-MFG test with HLA-DRB1 data from study participants with rheumatoid arthritis (RA) and their parents, we show that the v-MFG test successfully detects an MFG incompatibility effect on RA while simultaneously adjusting for a possible viability loss.

miR-514a regulates the tumour suppressor NF1 and modulates BRAFi sensitivity in melanoma

To identify ‘melanoma-specific’ microRNAs (miRNAs) we used an unbiased microRNA profiling approach to comprehensively study cutaneous melanoma in relation to other solid malignancies, which revealed 233 differentially expressed (≥ 2 fold, p < 0.05) miRNAs. Among the top 20 most significantly different miRNAs was hsa-miR-514a-3p. miR-514a is a member of a cluster of miRNAs (miR-506-514) involved in initiating melanocyte transformation and promotion of melanoma growth. We found miR-514a was expressed in 38/55 (69%) melanoma cell lines but in only 1/34 (3%) other solid cancers. To identify miR-514a regulated targets we conducted a miR-514a-mRNA ‘pull-down’ experiment, which revealed hundreds of genes, including: CTNNB1, CDK2, MC1R, and NF1, previously associated with melanoma. NF1 was selected for functional validation because of its recent implication inacquired resistance to BRAFV600E-targeted therapy. Luciferase-reporter assays confirmed NF1 as a direct target of miR-514a and over-expression of miR-514a in melanoma cell lines inhibited NF1 expression, which correlated with increased survival of BRAFV600E cells treated with PLX4032. These data provide another mechanism for the dysregulation of the MAPK pathway which may contribute to the profound resistance associated with current RAF-targeted therapies.

Validation of a Competitive Enzyme-Linked Immunosorbent Assay for Detection of Babesia bigemina Antibodies in Cattle

A competitive enzyme-linked immunosorbent assay (cELISA) based on a broadly conserved, species-specific, B-cell epitope within the C terminus of Babesia bigemina rhoptry-associated protein 1a was validated for international use. Receiver operating characteristic analysis revealed 16% inhibition as the threshold for a negative result, with an associated specificity of 98.3% and sensitivity of 94.7%. Increasing the threshold to 21% increased the specificity to 100% but modestly decreased the sensitivity to 87.2%. By using 21% inhibition, the positive predictive values ranged from 90.7% (10% prevalence) to 100% (95% prevalence) and the negative predictive values ranged from 97.0% (10% prevalence) to 48.2% (95% prevalence). The assay was able to detect serum antibody as early as 7 days after intravenous inoculation. The cELISA was distributed to five different laboratories along with a reference set of 100 defined bovine serum samples, including known positive, known negative, and field samples. The pairwise concordance among the five laboratories ranged from 100% to 97%, and all kappa values were above 0.8, indicating a high degree of reliability. Overall, the cELISA appears to have the attributes necessary for international application.

Developments in controlled green-water larval culture technologies for estuarine fishes in Queensland, Australia and elsewhere

Since 1989, researchers with the Department of Primary Industries and Fisheries (DPI&F) in Queensland, Australia, have successfully used controlled low-water exchange green-water cultures to rear the larvae of estuarine fishes and crustaceans through to metamorphosis. High survivals and excellent fry condition have been achieved for several commercially important endemic species produced for various projects. They include barramundi or sea bass, Lates calcarifer, Australian bass, Macquaria novemaculeata, dusky flathead, Platycephalus fuscus, sand whiting, Sillago ciliata, red sea bream or snapper, Pagrus auratus, banana prawn, Fenneropenaeus merguiensis, and others. The consistent success of our standardised and relatively simple approach at different localities has led to it being incorporated into general fingerling production practices at several establishments in Australia. Although post-metamorphosis rearing methods have differed for each species investigated, due to various biological and behavioural traits and project requirements, these larval rearing methods have been successful with few species-specific modifications. Initially modelled on the Taiwanese approach to rearing Penaeids in aerated low-water exchange cultures, the approach similarly appears to rely on a beneficial assemblage of micro-organisms. Conceptually, these micro-organisms may include a mixture of the air-borne primary invaders of pure phytoplankton cultures when exposed to outdoor conditions. Whilst this would vary with different sites, our experiences with these methods have consistently been favourable. Mass microalgal cultures with eco-physiological youth are used to regularly augment larval fish cultures so that rearing conditions simulate an exponential growth-phase microalgal bloom. Moderate to heavy aeration prevents settlement of particulate matter and encourages aerobic bacterial decomposition of wastes. The green-water larval rearing approach described herein has demonstrated high practical utility in research and commercial applications, and has greatly simplified marine finfish hatchery operations whilst generally lifting production capacities for metamorphosed fry in Australia. Its potential uses in areas of aquaculture other than larviculture are also discussed.

The introduction and release of Chiasmia inconspicua and C. assimilis (Lepidoptera: Geometridae) for the biological control of Acacia nilotica in Australia

Two geometrid moths Chiasmia inconspicua and Chiasmia assimilis, identified as potential biological control agents for prickly acacia Acacia nilotica subsp. indica, were collected in Kenya and imported into quarantine facilities in Australia where laboratory cultures were established. Aspects of the biologies of both insects were studied and CLIMEX® models indicating the climatically favourable areas of Australia were developed. Host range tests were conducted using an approved test list of 74 plant species and no-choice tests of neonate larvae placed on both cut foliage and potted plants. C. inconspicua developed through to adult on prickly acacia and, in small numbers, Acacia pulchella. C. assimilis developed through to adult on prickly acacia and also in very small numbers on A. pulchella, A. deanei, A. decurrens, and A. mearnsii. In all experiments, the response on prickly acacia could be clearly differentiated from the responses on the non-target species. Both insects were approved for release in Australia. Over a three-year period releases were made at multiple sites in north Queensland, almost all in inland areas. There was no evidence of either insect's establishment and both colonies were terminated. A new colony of C. assimilis was subsequently established from insects collected in South Africa and releases of C. assimilis from this new colony were made into coastal and inland infestations of prickly acacia. Establishment was rapid at one coastal site and the insect quickly spread to other infestations. Establishment at one inland area was also confirmed in early 2006. The establishment in coastal areas supported a CLIMEX model that indicated that the climate of coastal areas was more suitable than inland areas.

Tree-ring based hydroclimate reconstruction over northern Vietnam from Fokienia hodginsii: eighteenth century mega-drought and tropical Pacific influence

We present here the first statistically calibrated and verified tree-ring reconstruction of climate from continental Southeast Asia.The reconstructed variable is March-May (MAM) Palmer Drought Severity Index (PDSI) based on ring widths from 22 trees (42 radial cores) of rare and long-lived conifer, Fokienia hodginsii (Po Mu as locally called) from northern Vietnam. This is the first published tree ring chronology from Vietnam as well as the first for this species. Spanning 535 years, this is the longest cross-dated tree-ring series yet produced from continental Southeast Asia. Response analysis revealed that the annual growth of Fokienia at this site was mostly governed by soil moisture in the pre-monsoon season. The reconstruction passed the calibration-verification tests commonly used in dendroclimatology, and revealed two prominent periods of drought in the mid-eighteenth and late-nineteenth enturies. The former lasted nearly 30 years and was concurrent with a similar drought over northwestern Thailand inferred from teak rings, suggesting a ``mega-drought'' extending across Indochina in the eighteenth century. Both of our reconstructed droughts are consistent with the periods of warm sea surface temperature (SST)anomalies in the tropical Pacific. Spatial correlation analyses with global SST indicated that ENSO-like anomalies might play a role in modulating droughts over the region, with El Nio (warm) phases resulting in reduced rainfall. However, significant correlation was also seen with SST over the Indian Ocean and the north Pacific,suggesting that ENSO is not the only factor affecting the climate of the area. Spectral analyses revealed significant peaks in the range of 53.9-78.8 years as well as in the ENSO-variability range of 2.0 to 3.2 years.