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.

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected

Rank regression analysis of correlated water quality data from South East Queensland

With growing population and fast urbanization in Australia, it is a challenging task to maintain our water quality. It is essential to develop an appropriate statistical methodology in analyzing water quality data in order to draw valid conclusions and hence provide useful advices in water management. This paper is to develop robust rank-based procedures for analyzing nonnormally distributed data collected over time at different sites. To take account of temporal correlations of the observations within sites, we consider the optimally combined estimating functions proposed by Wang and Zhu (Biometrika, 93:459-464, 2006) which leads to more efficient parameter estimation. Furthermore, we apply the induced smoothing method to reduce the computational burden. Smoothing leads to easy calculation of the parameter estimates and their variance-covariance matrix. Analysis of water quality data from Total Iron and Total Cyanophytes shows the differences between the traditional generalized linear mixed models and rank regression models. Our analysis also demonstrates the advantages of the rank regression models for analyzing nonnormal data.

Study of Pinna nobilis growth from inner record: How biased are posterior adductor muscle scars estimates?

Previous studies have shown that the external growth records of the posterior adductor muscle scar (PAMS) of the bivalve Pinna nobilis are incomplete and do not produce accurate age estimations. We have developed a new methodology to study age and growth using the inner record of the PAMS, which avoids the necessity of costly in situ shell measurements or isotopic studies. Using the inner record we identified the positions of PAMS previously obscured by nacre and estimated the number of missing records in adult specimens with strong abrasion of the calcite layer in the anterior portion of the shell. The study of the PAMS and inner record of two shells that were 6 years old when collected showed that only 2 and 3 PAMS were observed, while 6 inner records could be counted, thus confirming our working methodology. Growth parameters of a P. nobilis population located in Moraira, Spain (western Mediterranean) were estimated with the new methodology and compared to those obtained using PAMS data and in situ measurements. For the comparisons, we applied different models considering the data alternatively as length-at-age (LA) and tag-recapture (TR). Among every method we tested to fit the Von Bertalanffy growth model, we observed that LA data from inner record fitted to the model using non-linear mixed effects and the estimation of missing records using the calcite width was the most appropriate. The equation obtained with this method, L = 573*(1 - e(-0.16(t-0.02))), is very similar to that calculated previously from in situ measurements for the same population.

Patterns in the distribution and abundance of female eastern king prawns, Melicertus plebejus (Hess, 1865), capable of spawning and reproductive potential in waters off

Linear mixed models were used to test the null hypothesis that there were no differences between seasons and locations in the reproductive potential of female eastern king prawns, Melicertus plebejus along the east coast of Australia. Three samples were collected in each season between autumn 1991 and winter 1992 (inclusive). Females capable of spawning were found at all locations but proportions were greater in lower than higher latitudes. Females capable of spawning were not found at the southern (highest latitude) most location in all seasons. There was a significant interaction in reproductive potential between seasons and locations suggesting that patterns among seasons differed between locations and vice versa. Reproductive potential was greatest amongst the northern (lower latitudes) most locations and was greatest in autumn at these locations. Seasonal patterns were less pronounced further south (higher latitudes). The length composition of females in catches differed between locations with more larger prawns found in samples from northern locations. The challenge that remains is to quantify the oceanic sources of larvae that contribute to recruitment in each nursery area and the estuarine sources of juveniles that contribute adults back to the effective spawning stock. Maintaining the effective spawning stock and important nursery areas are crucial to the sustainability of this resource.

Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer's disease, multiple sclerosis and endometriosis

Common diseases such as endometriosis (ED), Alzheimer's disease (AD) and multiple sclerosis (MS) account for a significant proportion of the health care burden in many countries. Genome-wide association studies (GWASs) for these diseases have identified a number of individual genetic variants contributing to the risk of those diseases. However, the effect size for most variants is small and collectively the known variants explain only a small proportion of the estimated heritability. We used a linear mixed model to fit all single nucleotide polymorphisms (SNPs) simultaneously, and estimated genetic variances on the liability scale using SNPs from GWASs in unrelated individuals for these three diseases. For each of the three diseases, case and control samples were not all genotyped in the same laboratory. We demonstrate that a careful analysis can obtain robust estimates, but also that insufficient quality control (QC) of SNPs can lead to spurious results and that too stringent QC is likely to remove real genetic signals. Our estimates show that common SNPs on commercially available genotyping chips capture significant variation contributing to liability for all three diseases. The estimated proportion of total variation tagged by all SNPs was 0.26 (SE 0.04) for ED, 0.24 (SE 0.03) for AD and 0.30 (SE 0.03) for MS. Further, we partitioned the genetic variance explained into five categories by a minor allele frequency (MAF), by chromosomes and gene annotation. We provide strong evidence that a substantial proportion of variation in liability is explained by common SNPs, and thereby give insights into the genetic architecture of the diseases.

Fishing power and standardised catch rates: Implications of missing vessel-characteristic data from the Australian eastern king prawn fishery

Standardised time series of fishery catch rates require collations of fishing power data on vessel characteristics. Linear mixed models were used to quantify fishing power trends and study the effect of missing data encountered when relying on commercial logbooks. For this, Australian eastern king prawn (Melicertus plebejus) harvests were analysed with historical (from vessel surveys) and current (from commercial logbooks) vessel data. Between 1989 and 2010, fishing power increased up to 76%. To date, both forward-filling and, alternatively, omitting records with missing vessel information from commercial logbooks produce broadly similar fishing power increases and standardised catch rates, due to the strong influence of years with complete vessel data (16 out of 23 years of data). However, if gaps in vessel information had not originated randomly and skippers from the most efficient vessels were the most diligent at filling in logbooks, considerable errors would be introduced. Also, the buffering effect of complete years would be short lived as years with missing data accumulate. Given ongoing changes in fleet profile with high-catching vessels fishing proportionately more of the fleet’s effort, compliance with logbook completion, or alternatively ongoing vessel gear surveys, is required for generating accurate estimates of fishing power and standardised catch rates.

Soil carbon stock in the tropical rangelands of Australia: Effects of soil type and grazing pressure, and determination of sampling requirement.

On-going, high-profile public debate about climate change has focussed attention on how to monitor the soil organic carbon stock (C(s)) of rangelands (savannas). Unfortunately, optimal sampling of the rangelands for baseline C(s) - the critical first step towards efficient monitoring - has received relatively little attention to date. Moreover, in the rangelands of tropical Australia relatively little is known about how C(s) is influenced by the practice of cattle grazing. To address these issues we used linear mixed models to: (i) unravel how grazing pressure (over a 12-year period) and soil type have affected C(s) and the stable carbon isotope ratio of soil organic carbon (delta(13)C) (a measure of the relative contributions of C(3) and C(4) vegetation to C(s)); (ii) examine the spatial covariation of C(s) and delta(13)C; and, (iii) explore the amount of soil sampling required to adequately determine baseline C(s). Modelling was done in the context of the material coordinate system for the soil profile, therefore the depths reported, while conventional, are only nominal. Linear mixed models revealed that soil type and grazing pressure interacted to influence C(s) to a depth of 0.3 m in the profile. At a depth of 0.5 m there was no effect of grazing on C(s), but the soil type effect on C(s) was significant. Soil type influenced delta(13)C to a soil depth of 0.5 m but there was no effect of grazing at any depth examined. The linear mixed model also revealed the strong negative correlation of C(s) with delta(13)C, particularly to a depth of 0.1 m in the soil profile. This suggested that increased C(s) at the study site was associated with increased input of C from C(3) trees and shrubs relative to the C(4) perennial grasses; as the latter form the bulk of the cattle diet, we contend that C sequestration may be negatively correlated with forage production. Our baseline C(s) sampling recommendation for cattle-grazing properties of the tropical rangelands of Australia is to: (i) divide the property into units of apparently uniform soil type and grazing management; (ii) use stratified simple random sampling to spread at least 25 soil sampling locations about each unit, with at least two samples collected per stratum. This will be adequate to accurately estimate baseline mean C(s) to within 20% of the true mean, to a nominal depth of 0.3 m in the profile.

Fishing power increases continue in Queensland's east coast trawl fishery, Australia

The Queensland east coast trawl fishery is by far the largest prawn and scallop otter trawl fleet in Australia in terms of number of vessels, with 504 vessels licensed to fish for species including tiger prawns, endeavour prawns, red spot king prawns, eastern king prawns and saucer scallops by the end of 2004. The vessel fleet has gradually upgraded characteristics such as engine power and use of propeller nozzles, quad nets, global positioning systems (GPS) and computer mapping software. These changes, together with the ever-changing profile of the fleet, were analysed by linear mixed models to quantify annual efficiency increases of an average vessel at catching prawns or scallops. The analyses included vessel characteristics (treated as fixed effects) and vessel identifier codes (treated as random effects). For the period from 1989 to 2004 the models estimated overall fishing power increases of 6% in the northern tiger, 6% in the northern endeavour, 12% in the southern tiger, 18% in the red spot king, 46% in the eastern king prawn and 15% in the saucer scallop sector. The results illustrate the importance of ongoing monitoring of vessel and fleet characteristics and the need to use this information to standardise catch rate indices used in stock assessment and management.

Risk factors associated with corneal nerve alteration in type 1 diabetes in the absence of neuropathy: A longitudinal in vivo corneal confocal microscopy study

Purpose The aim of this study was to determine alterations to the corneal subbasal nerve plexus (SNP) over four years using in vivo corneal confocal microscopy (IVCM) in participants with type 1 diabetes and to identify significant risk factors associated with these alterations. Methods A cohort of 108 individuals with type 1 diabetes and no evidence of peripheral neuropathy at enrollment underwent laser-scanning IVCM, ocular screening, and health and metabolic assessment at baseline and the examinations continued for four subsequent annual visits. At each annual visit, eight central corneal images of the SNP were selected and analyzed to quantify corneal nerve fiber density (CNFD), branch density (CNBD) and fiber length (CNFL). Linear mixed model approaches were fitted to examine the relationship between risk factors and corneal nerve parameters. Results A total of 96 participants completed the final visit and 91 participants completed all visits. No significant relationships were found between corneal nerve parameters and time, sex, duration of diabetes, smoking, alcohol consumption, blood pressure or BMI. However, CNFD was negatively associated with HbA1c (β=-0.76, P<0.01) and age (β=-0.13, P<0.01) and positively related to high density lipids (HDL) (β=2.01, P=0.03). Higher HbA1c (β=-1.58, P=0.04) and age (β=-0.23, P<0.01) also negatively impacted CNBD. CNFL was only affected by higher age (β=-0.06, P<0.01). Conclusions Glycemic control, HDL and age have significant effects on SNP structure. These findings highlight the importance of diabetic management to prevent corneal nerve damage as well as the capability of IVCM for monitoring subclinical alterations in the corneal SNP in diabetes.

Generalised growth models for aquatic species with an application to blacklip abalone (Haliotis rubra)

This paper presents a maximum likelihood method for estimating growth parameters for an aquatic species that incorporates growth covariates, and takes into consideration multiple tag-recapture data. Individual variability in asymptotic length, age-at-tagging, and measurement error are also considered in the model structure. Using distribution theory, the log-likelihood function is derived under a generalised framework for the von Bertalanffy and Gompertz growth models. Due to the generality of the derivation, covariate effects can be included for both models with seasonality and tagging effects investigated. Method robustness is established via comparison with the Fabens, improved Fabens, James and a non-linear mixed-effects growth models, with the maximum likelihood method performing the best. The method is illustrated further with an application to blacklip abalone (Haliotis rubra) for which a strong growth-retarding tagging effect that persisted for several months was detected. (C) 2013 Elsevier B.V. All rights reserved.

Resistance to Fusarium oxysporum f. sp fragariae and Predicted Breeding Values in Strawberry

Fusarium wilt of strawberry, incited by Fusarium oxysporum f. sp. fragariae (Fof), is a major disease of the cultivated strawberry (Fragaria xananassa) worldwide. An increase in disease outbreaks of the pathogen in Western Australia and Queensland plus the search for alternative disease management strategies place emphasis on the development of resistant cultivars. In response, a partial incomplete diallel cross involving four parents was performed for use in glasshouse resistance screenings. The resulting progeny were evaluated for their susceptibility to Fof. Best-performing progeny and suitability of progenies as parents were determined using data from disease severity ratings and analyzed using a linear mixed model incorporating a pedigree to produce best linear unbiased predictions of breeding values. Variation in disease response, ranging from highly susceptible to resistant, indicates a quantitative effect. The estimate of the narrow-sense heritability was 0.49 +/- 0.04 (SE), suggesting the population should be responsive to phenotypic recurrent selection. Several progeny genotypes have predicted breeding values higher than any of the parents. Knowledge of Fof resistance derived from this study can help select best parents for future crosses for the development of new strawberry cultivars with Fof resistance.

Effect of Broadcast Baiting on Abundance Patterns of Red Imported Fire Ants (Hymenoptera: Formicidae) and Key Local Ant Genera at Long-Term Monitoring Sites in Brisbane, Australia

In 2001, the red imported fire ant (Solenopsis invicta Buren) was identified in Brisbane, Australia. An eradication program involving broadcast bait treatment with two insect growth regulators and a metabolic inhibitor began in September of that year and is currently ongoing. To gauge the impacts of these treatments on local ant populations, we examined long-term monitoring data and quantified abundance patterns of S. invicta and common local ant genera using a linear mixed-effects model. For S. invicta, presence in pitfalls reduced over time to zero on every site. Significantly higher numbers of S. invicta workers were collected on high-density polygyne sites, which took longer to disinfest compared with monogyne and low-density polygyne sites. For local ants, nine genus groups of the 10 most common genera analyzed either increased in abundance or showed no significant trend. Five of these genus groups were significantly less abundant at the start of monitoring on high-density polygyne sites compared with monogyne and low-density polygyne sites. The genus Pheidole significantly reduced in abundance over time, suggesting that it was affected by treatment efforts. These results demonstrate that the treatment regime used at the time successfully removed S. invicta from these sites in Brisbane, and that most local ant genera were not seriously impacted by the treatment. These results have important implications for current and future prophylactic treatment efforts, and suggest that native ants remain in treated areas to provide some biological resistance to S. invicta.

The effect of pasture utilization rate on stocks of soil organic carbon and total nitrogen in a semi-arid tropical grassland

The influence of grazing management on total soil organic carbon (SOC) and soil total nitrogen (TN) in tropical grasslands is an issue of considerable ecological and economic interest. Here we have used linear mixed models to investigate the effect of grazing management on stocks of SOC and TN in the top 0.5 m of the soil profile. The study site was a long-term pasture utilization experiment, 26 years after the experiment was established for sheep grazing on native Mitchell grass (Astrebla spp.) pasture in northern Australia. The pasture utilization rates were between 0% (exclosure) and 80%, assessed visually. We found that a significant amount of TN had been lost from the top 0.1 m of the soil profile as a result of grazing, with 80% pasture utilization resulting in a loss of 84 kg ha−1 over the 26-year period. There was no significant effect of pasture utilization rate on TN when greater soil depths were considered. There was no significant effect of pasture utilization rate on stocks of SOC and soil particulate organic carbon (POC), or the C:N ratio at any depth; however, visual trends in the data suggested some agreement with the literature, whereby increased grazing pressure appeared to: (i) decrease SOC and POC stocks; and, (ii) increase the C:N ratio. Overall, the statistical power of the study was limited, and future research would benefit from a more comprehensive sampling scheme. Previous studies at the site have found that a pasture utilization rate of 30% is sustainable for grazing production on Mitchell grass; however, given our results, we conclude that N inputs (possibly through management of native N2-fixing pasture legumes) should be made for long-term maintenance of soil health, and pasture productivity, within this ecosystem.

Flying-Fox Roost Disturbance and Hendra Virus Spillover Risk

Bats of the genus Pteropus (flying-foxes) are the natural host of Hendra virus (HeV) which periodically causes fatal disease in horses and humans in Australia. The increased urban presence of flying-foxes often provokes negative community sentiments because of reduced social amenity and concerns of HeV exposure risk, and has resulted in calls for the dispersal of urban flying-fox roosts. However, it has been hypothesised that disturbance of urban roosts may result in a stress-mediated increase in HeV infection in flying-foxes, and an increased spillover risk. We sought to examine the impact of roost modification and dispersal on HeV infection dynamics and cortisol concentration dynamics in flying-foxes. The data were analysed in generalised linear mixed models using restricted maximum likelihood (REML). The difference in mean HeV prevalence in samples collected before (4.9%), during (4.7%) and after (3.4%) roost disturbance was small and non-significant (P = 0.440). Similarly, the difference in mean urine specific gravity-corrected urinary cortisol concentrations was small and non-significant (before = 22.71 ng/mL, during = 27.17, after = 18.39) (P= 0.550). We did find an underlying association between cortisol concentration and season, and cortisol concentration and region, suggesting that other (plausibly biological or environmental) variables play a role in cortisol concentration dynamics. The effect of roost disturbance on cortisol concentration approached statistical significance for region, suggesting that the relationship is not fixed, and plausibly reflecting the nature and timing of disturbance. We also found a small positive statistical association between HeV excretion status and urinary cortisol concentration. Finally, we found that the level of flying-fox distress associated with roost disturbance reflected the nature and timing of the activity, highlighting the need for a ‘best practice’ approach to dispersal or roost modification activities. The findings usefully inform public discussion and policy development in relation to Hendra virus and flying-fox management.