A versatile gene-based test for genome-wide association studies

We have derived a versatile gene-based test for genome-wide association studies (GWAS). Our approach, called VEGAS (versatile gene-based association study), is applicable to all GWAS designs, including family-based GWAS, meta-analyses of GWAS on the basis of summary data, and DNA-pooling-based GWAS, where existing approaches based on permutation are not possible, as well as singleton data, where they are. The test incorporates information from a full set of markers (or a defined subset) within a gene and accounts for linkage disequilibrium between markers by using simulations from the multivariate normal distribution. We show that for an association study using singletons, our approach produces results equivalent to those obtained via permutation in a fraction of the computation time. We demonstrate proof-of-principle by using the gene-based test to replicate several genes known to be associated on the basis of results from a family-based GWAS for height in 11,536 individuals and a DNA-pooling-based GWAS for melanoma in approximately 1300 cases and controls. Our method has the potential to identify novel associated genes; provide a basis for selecting SNPs for replication; and be directly used in network (pathway) approaches that require per-gene association test statistics. We have implemented the approach in both an easy-to-use web interface, which only requires the uploading of markers with their association p-values, and a separate downloadable application.

A simple and fast two-locus quality control test to detect false positives due to batch effects in genome-wide association studies

The impact of erroneous genotypes having passed standard quality control (QC) can be severe in genome-wide association studies, genotype imputation, and estimation of heritability and prediction of genetic risk based on single nucleotide polymorphisms (SNP). To detect such genotyping errors, a simple two-locus QC method, based on the difference in test statistic of association between single SNPs and pairs of SNPs, was developed and applied. The proposed approach could detect many problematic SNPs with statistical significance even when standard single SNP QC analyses fail to detect them in real data. Depending on the data set used, the number of erroneous SNPs that were not filtered out by standard single SNP QC but detected by the proposed approach varied from a few hundred to thousands. Using simulated data, it was shown that the proposed method was powerful and performed better than other tested existing methods. The power of the proposed approach to detect erroneous genotypes was approximately 80% for a 3% error rate per SNP. This novel QC approach is easy to implement and computationally efficient, and can lead to a better quality of genotypes for subsequent genotype-phenotype investigations.

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.

Interactions between a Trawl fishery and spatial closures for biodiversity conservation in the Great Barrier Reef World Heritage Area, Australia.

Background: The Queensland East Coast Otter Trawl Fishery (ECOTF) for penaeid shrimp fishes within Australia's Great Barrier Reef World Heritage Area (GBRWHA). The past decade has seen the implementation of conservation and fisheries management strategies to reduce the impact of the ECOTF on the seabed and improve biodiversity conservation. New information from electronic vessel location monitoring systems (VMS) provides an opportunity to review the interactions between the ECOTF and spatial closures for biodiversity conservation. Methodology and Results: We used fishing metrics and spatial information on the distribution of closures and modelled VMS data in a geographical information system (GIS) to assess change in effort of the trawl fishery from 2001-2009 and to quantify the exposure of 70 reef, non-reef and deep water bioregions to trawl fishing. The number of trawlers and the number of days fished almost halved between 2001 and 2009 and new spatial closures introduced in 2004 reduced the area zoned available for trawl fishing by 33%. However, we found that there was only a relatively minor change in the spatial footprint of the fishery as a result of new spatial closures. Non-reef bioregions benefited the most from new spatial closures followed by deep and reef bioregions. Conclusions/Significance: Although the catch of non target species remains an issue of concern for fisheries management, the small spatial footprint of the ECOTF relative to the size of the GBRWHA means that the impact on benthic habitats is likely to be negligible. The decline in effort as a result of fishing industry structural adjustment, increasing variable costs and business decisions of fishers is likely to continue a trend to fish only in the most productive areas. This will provide protection for most benthic habitats without any further legislative or management intervention.

Economic case study of ABCD cane management practices in the Burdekin River Irrigation Area (BRIA) region

A case study was undertaken to determine the economic impact of a change in management class as detailed in the A, B, C and D management class framework. This document focuses on the implications of changing from D to C, C to B and B to A class management in the Burdekin River irrigation area (BRIA) and if the change is worthwhile from an economic perspective. This report provides a guide to the economic impact that may be expected when undertaking a particular change in farming practices and will ultimately lead to more informed decisions being made by key industry stakeholders. It is recognised that these management classes have certain limitations and in many cases the grouping of practices may not be reflective of the real situation. The economic case study is based on the A, B, C and D management class framework for water quality improvement developed in 2007/2008 for the Burdekin natural resource management region. The framework for the Burdekin is currently being updated to clarify some issues and incorporate new knowledge since the earlier version of the framework. However, this updated version is not yet complete and so the Paddock to Reef project has used the most current available version of the framework for the modelling and economics. As part of the project specification, sugarcane crop production data for the BRIA was provided by the APSIM model. The information obtained from the APSIM crop modelling programme included sugarcane yields and legume grain yield (legume grain yield only applies to A class management practice). Because of the complexity involved in the economic calculations, a combination of the FEAT, PiRisk and a custom made spreadsheet was used for the economic analysis. Figures calculated in the FEAT program were transferred to the custom made spreadsheet to develop a discounted cash flow analysis. The marginal cash flow differences for each farming system were simulated over a 5-year and 10-year planning horizon to determine the net present value of changing across different management practices. PiRisk was used to test uncertain parameters in the economic analysis and the potential risk associated with a change in value.

Paddock to Reef Monitoring & Evaluation - Economic analysis of ‘ABCD’ cane management for the Burdekin River Irrigation Area.

The economic analysis is based on the A, B, C and D management practice framework for water quality improvement developed in 2007/2008 by the respective natural resource management region. This document focuses on the economic implications of these management practices in the Burdekin River Irrigation Area (BRIA). A review of the management practices is currently being undertaken to clarify some issues and incorporate new knowledge since the earlier version of the framework. However, this updated version is not yet complete and so the Paddock to Reef project has used the most current available version of the framework for the modelling and economics.

Model of Strangriede Cemetery during World War II when it was a holding area for Jews to be deported; Hannover Cemeteries

Digital Image

Estimating winter crop area across seasons and regions using time-sequential MODIS imagery.

The wheat grain industry is Australia's second largest agricultural export commodity. There is an increasing demand for accurate, objective and near real-time crop production information by industry. The advent of the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite platform has augmented the capability of satellite-based applications to capture reflectance over large areas at acceptable pixel scale, cost and accuracy. The use of multi-temporal MODIS-enhanced vegetation index (EVI) imagery to determine crop area was investigated in this article. Here the rigour of the harmonic analysis of time-series (HANTS) and early-season metric approaches was assessed when extrapolating over the entire Queensland (QLD) cropping region for the 2005 and 2006 seasons. Early-season crop area estimates, at least 4 months before harvest, produced high accuracy at pixel and regional scales with percent errors of -8.6% and -26% for the 2005 and 2006 seasons, respectively. In discriminating among crops at pixel and regional scale, the HANTS approach showed high accuracy. The errors for specific area estimates for wheat, barley and chickpea were 9.9%, -5.2% and 10.9% (for 2005) and -2.8%, -78% and 64% (for 2006), respectively. Area estimates of total winter crop, wheat, barley and chickpea resulted in coefficient of determination (R(2)) values of 0.92, 0.89, 0.82 and 0.52, when contrasted against the actual shire-scale data. A significantly high coefficient of determination (0.87) was achieved for total winter crop area estimates in August across all shires for the 2006 season. Furthermore, the HANTS approach showed high accuracy in discriminating cropping area from non-cropping area and highlighted the need for accurate and up-to-date land use maps. The extrapolability of these approaches to determine total and specific winter crop area estimates, well before flowering, showed good utility across larger areas and seasons. Hence, it is envisaged that this technology might be transferable to different regions across Australia.

Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.

Genome-wide association studies of quantitative traits with related individuals: Little (power) lost but much to be gained

For complex disease genetics research in human populations, remarkable progress has been made in recent times with the publication of a number of genome-wide association scans (GWAS) and subsequent statistical replications. These studies have identified new genes and pathways implicated in disease, many of which were not known before. Given these early successes, more GWAS are being conducted and planned, both for disease and quantitative phenotypes. Many researchers and clinicians have DNA samples available on collections of families, including both cases and controls. Twin registries around the world have facilitated the collection of large numbers of families, with DNA and multiple quantitative phenotypes collected on twin pairs and their relatives. In the design of a new GWAS with a fixed budget for the number of chips, the question arises whether to include or exclude related individuals. It is commonly believed to be preferable to use unrelated individuals in the first stage of a GWAS because relatives are 'over-matched' for genotypes. In this study, we quantify that for GWAS of a quantitative phenotype, relative to a sample of unrelated individuals surprisingly little power is lost when using relatives. The advantages of using relatives are manifold, including the ability to perform more quality control, the choice to perform within-family tests of association that are robust to population stratification, and the ability to perform joint linkage and association analysis. Therefore, the advantages of using relatives in GWAS for quantitative traits may well outweigh the small disadvantage in terms of statistical power.

A genome-wide linkage scan provides evidence for both new and previously reported loci influencing common migraine

Latent class analysis was performed on migraine symptom data collected in a Dutch population sample (N = 12,210, 59% female) in order to obtain empirical groupings of individuals suffering from symptoms of migraine headache. Based on these heritable groupings (h(2) = 0.49, 95% CI: 0.41-0.57) individuals were classified as affected (migrainous headache) or unaffected. Genome-wide linkage analysis was performed using genotype data from 105 families with at least 2 affected siblings. In addition to this primary phenotype, linkage analyses were performed for the individual migraine symptoms. Significance levels, corrected for the analysis of multiple traits, were determined empirically via a novel simulation approach. Suggestive linkage for migrainous headache was found on chromosomes 1 (LOD = 1.63; pointwise P = 0.0031), 13 (LOD = 1.63; P = 0.0031), and 20 (LOD = 1.85; P = 0.0018). Interestingly, the chromosome 1 peak was located close to the ATP1A2 gene, associated with familial hemiplegic migraine type 2 (FHM2). Individual symptom analysis produced a LOD score of 1.97 (P = 0.0013) on chromosome 5 (photo/phonophobia), a LOD score of 2.13 (P = 0.0009) on chromosome 10 (moderate/severe pain intensity) and a near significant LOD score of 3.31 (P = 0.00005) on chromosome 13 (pulsating headache). These peaks were all located near regions previously reported in migraine linkage studies. Our results provide important replication and support for the presence of migraine susceptibility genes within these regions, and further support the utility of an LCA-based phenotyping approach and analysis of individual symptoms in migraine genetic research. Additionally, our novel "2-step" analysis and simulation approach provides a powerful means to investigate linkage to individual trait components.

Effects of living crown reduction on needle element status of Scots pine

A wide range of biotic and abiotic factors, operating over different time perspectives and intensities, cause defoliation and a rapid decrease in the crown size of trees. Scleroderris canker disease [Gremmeniella abietina (Lagerb.) Morelet] has caused widespread crown reduction and tree mortality in Scots pine (Pinus sylvestris L) in forests in Scandinavia during the last three decades. In the 1980's, attempts were made to show, on the basis of the higher foliar N and S concentrations of affected pines in the diseased area, that sulphur and nitrogen deposition predispose trees to G. abietina. Unfortunately, in many studies on defoliated trees, exceptionally high or low needle mineral nutrient concentrations are still often interpreted as one of the causes of tree injury and not, conversely, as the result. In this thesis, three different field experiments, with foliar analysis as the main study method, were conducted in order to asses the possible long-term effects of living crown reduction on the needle nutrient concentrations of Scots pine trees in southern Finland. The crown ratio and length of the living crown were used to estimate the amount of defoliation in the reduced canopies. The material for the partial studies was collected and a total of 968 foliar samples were analysed individually (15-17 elements/sample) on a total of 488 sample trees (140 diseased, 116 pruned and 232 control trees) during the years 1987-1996 in 13 Scots pine stands. All the three experiments of this thesis provided significant evidence that severe, disease-induced defoliation or artificial pruning of the living branches can induce long-lasting nutritional changes in the foliage of the recovering trees under the typical growing conditions for Scots pine. The foliar concentrations of all the 17 mineral nutrients/elements analysed were affected, to a varying degree, by artificial pruning during the following three years. Although Scots pine, as an evergreen conifer, is considered to have low induced chemical responses to defoliation, this study proved experimentally under natural forest conditions that severe artificial pruning or disease-induced defoliation of Scots pine trees may induce biologically significant changes in the concentrations of most of the important macro- and micronutrients, as well as of carbon, in refoliated needles. Concerning the studies in this thesis, I find the results significant in providing new information about the long-term effects of rapid living crown reduction on the foliar nutrient and element status of Scots pine trees. Key words: Foliar analysis, defoliation, needle loss, pruning, nutrients, Pinus sylvestris, Gremmeniella abietina

Fuel conditions and fire behavior characteristics of managed Picea abies and Pinus sylvestris forests in Finland

The objectives of this study were to analyze the impact of structural stand characteristics on ignition potential, surface fuel moisture, and fire behavior in Pinus sylvestris L. and Picea abies (L.) Karst stands in Finland and to explain stand-specific fire danger using the Canadian Fire Weather Index System and the Finnish Fire Risk Index. Additionally, the study analyzes the relationship between observed fire activity and fire weather indices at different stages of growing season. Field experiments were carried out in Pinus sylvestris or Picea abies dominated stands during fire seasons 2001 and 2002. Observations on ignition potential, fuel moisture, and fire behavior were analyzed in relation to stand structure and the outputs of the Finnish and Canadian fire weather indices. Seasonal patterns of fire activity were examined based on national fire statistics 1996 2003, effective temperature sum, and the fire weather indices. Point fire ignition potential was highest in Pinus clear-cuts and lowest in closed Picea stands. Moss-dominated surface fuels were driest in clear-cut and sapling stage stands and presented the highest moisture content under closed Picea canopy. Pinus sylvestris stands carried fire under a wide range of fire weather conditions under which Picea abies stands failed to sustain fire. In the national fire records, the daily number of reported ignitions presented its highest value during late fire season whereas the daily area burned peaked most substantially during early season. The fire weather indices correlated significantly with ignition potential and fuel moisture but were unable to explain fire behavior in the experimental fires. During the initial and final stages of the growing season, fire activity was disconnected from weather-based fire danger ratings. Information on stand structure and season stage would benefit the assessment of fire danger in Finnish forest landscape for fire suppression and controlled burning purposes.

What was behind the bark? : An assessment of decay among urban Tilia, Betula and Acer trees felled as hazardous in the Helsinki City area

Old trees growing in urban environments are often felled due to symptoms of mechanical defects that could be hazardous to people and property. The decisions concerning these removals are justified by risk assessments carried out by tree care professionals. The major motivation for this study was to determine the most common profiles of potential hazard characteristics for the three most common urban tree genera in Helsinki City: Tilia, Betula and Acer, and in this way improve management practices and protection of old amenity trees. For this research, material from approximately 250 urban trees was collected in cooperation with the City of Helsinki Public Works Department during 2001 - 2004. From the total number of trees sampled, approximately 70% were defined as hazardous. The tree species had characteristic features as potential hazard profiles. For Tilia trees, hollowed heartwood with low fungal activity and advanced decay caused by Ganoderma lipsiense were the two most common profiles. In Betula spp., the primary reason for tree removal was usually lowered amenity value in terms of decline of the crown. Internal cracks, most often due to weak fork formation, were common causes of potential failure in Acer spp. Decay caused by Rigidoporus populinus often increased the risk of stem breakage in these Acer trees. Of the decay fungi observed, G. lipsiense was most often the reason for the increased risk of stem collapse. Other fungi that also caused extensive decay were R. populinus, Inonotus obliquus, Kretzschmaria deusta and Phellinus igniarius. The most common decay fungi in terms of incidence were Pholiota spp., but decay caused by these species did not have a high potential for causing stem breakage, because it rarely extended to the cambium. The various evaluations used in the study suggested contradictions in felling decisions based on trees displaying different stages of decay. For protection of old urban trees, it is crucial to develop monitoring methods so that tree care professionals could better analyse the rate of decay progression towards the sapwood and separate those trees with decreasing amounts of sound wood from those with decay that is restricted to the heartwood area.