Association of ERAP1, but not IL23R, with ankylosing spondylitis in a Han Chinese population

Objective The results of a recent genome-wide association study have shown that ERAP1 and IL23R are associated with ankylosing spondylitis (AS) in Caucasian populations from North America and the UK. Based on these findings, we undertook the current study to investigate whether single-nucleotide polymorphisms (SNPs) covering the genes ERAP1 and IL23R are associated with AS in a Han Chinese population. Methods A case-control study was performed in Han Chinese patients with AS (n = 527) and controls (n = 945) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The Sequenom iPlex platform was used to genotype cases and controls for 21 tag SNPs covering IL23R and 38 tag SNPs covering ERAP1. Statistical analysis was performed using the Cochran-Armitage test for trend. Results Multiple SNPs in ERAP1 were significantly associated with AS (for rs27980, P = 0.0048; for rs7711564, P = 0.0081). However, no association was observed between IL23R and AS (for all SNPs, P > 0.1). The nonsynonymous SNP in IL23R, rs11209026, widely thought to be the primary AS-associated SNP in IL23R in Europeans, was found not to be polymorphic in Chinese. Conclusion Our results demonstrate that genetic polymorphisms in ERAP1 are associated with AS in Han Chinese, suggesting a common pathogenic mechanism for the disease in Chinese and Caucasian populations, and that IL23R is not associated with AS in Chinese, indicating a difference in the mechanism of disease pathogenesis between Chinese and Caucasian populations. This may result from the fact that rs11209026, the nonsynonymous SNP in IL23R, is not polymorphic in Chinese patients, providing further evidence that rs11209026 is the key polymorphism associated with AS (and likely inflammatory bowel disease and psoriasis) in this gene.

New England Youth Speeding Project: A pilot study of young drivers and risky driving: Effective strategies for preventing speeding amongst young Australians aged 16-25 years in the New England Region

Young Australian drivers aged 17 – 25 years are overwhelmingly represented in road fatalities where speed is a factor. In the combined LGAs of Armidale Dumaresq, Guyra, Uralla and Walcha in the 5 years 1999-2003 inclusive, 43% of speeding related casualty crashes involved a young driver aged less than 25 years. This is despite the fact that the 17-25 age group account for only 25% of the driving population in this area. Young male drivers account for the majority of these crashes and also tend to have a higher number of driving offences and accrue more penalties for road traffic offences, especially speeding. By analysing data from questionnaires by male and female participants this research project has been able to evaluate road safety advertisements to determine which ones are most effective to young drivers, what features of these advertisements are effective, how males differ from females in their receptiveness and preferences for road safety advertisements and specifically how to target young people especially young men in conveying road safety messages. Finally this research project has identified factors that are important in the production of media road safety advertisements and has made recommendations for how best to convey effective road safety messages to young Australian drivers in rural areas.

Genome-wide association study for radiographic vertebral fractures: A potential role for the 16q24 BMD locus

Vertebral fracture risk is a heritable complex trait. The aim of this study was to identify genetic susceptibility factors for osteoporotic vertebral fractures applying a genome-wide association study (GWAS) approach. The GWAS discovery was based on the Rotterdam Study, a population-based study of elderly Dutch individuals aged >55years; and comprising 329 cases and 2666 controls with radiographic scoring (McCloskey-Kanis) and genetic data. Replication of one top-associated SNP was pursued by de-novo genotyping of 15 independent studies across Europe, the United States, and Australia and one Asian study. Radiographic vertebral fracture assessment was performed using McCloskey-Kanis or Genant semi-quantitative definitions. SNPs were analyzed in relation to vertebral fracture using logistic regression models corrected for age and sex. Fixed effects inverse variance and Han-Eskin alternative random effects meta-analyses were applied. Genome-wide significance was set at p<5×10-8. In the discovery, a SNP (rs11645938) on chromosome 16q24 was associated with the risk for vertebral fractures at p=4.6×10-8. However, the association was not significant across 5720 cases and 21,791 controls from 14 studies. Fixed-effects meta-analysis summary estimate was 1.06 (95% CI: 0.98-1.14; p=0.17), displaying high degree of heterogeneity (I2=57%; Qhet p=0.0006). Under Han-Eskin alternative random effects model the summary effect was significant (p=0.0005). The SNP maps to a region previously found associated with lumbar spine bone mineral density (LS-BMD) in two large meta-analyses from the GEFOS consortium. A false positive association in the GWAS discovery cannot be excluded, yet, the low-powered setting of the discovery and replication settings (appropriate to identify risk effect size >1.25) may still be consistent with an effect size <1.10, more of the type expected in complex traits. Larger effort in studies with standardized phenotype definitions is needed to confirm or reject the involvement of this locus on the risk for vertebral fractures.

Suggestive linkage of 2p22-25 and 11q12-13 with low bone mineral density at the lumbar spine in the Irish population

Osteoporosis is a disease characterized by low bone mineral density (BMD) and poor bone quality. Peak bone density is achieved by the third decade of life, after which bone is maintained by a balanced cycle of bone resorption and synthesis. Age-related bone loss occurs as the bone resorption phase outweighs the bone synthesis phase of bone metabolism. Heritability accounts for up to 90% of the variability in BMD. Chromosomal loci including 1p36, 2p22-25, 11q12-13, parathyroid hormone receptor type 1 (PTHR1), interleukin-6 (IL-6), interleukin 1 alpha (IL-1α) and type II collagen A1/vitamin D receptor (COL11A1/VDR) have been linked or shown suggestive linkage with BMD in other populations. To determine whether these loci predispose to low BMD in the Irish population, we investigated 24 microsatellite markers at 7 chromosomal loci by linkage studies in 175 Irish families of probands with primary low BMD (T-score ≤ -1.5). Nonparametric analysis was performed using the maximum likelihood variance estimation and traditional Haseman-Elston tests on the Mapmaker/Sibs program. Suggestive evidence of linkage was observed with lumbar spine BMD at 2p22-25 (maximum LOD score 2.76) and 11q12-13 (MLS 2.55). One region, 1p36, approached suggestive linkage with femoral neck BMD (MLS 2.17). In addition, seven markers achieved LOD scores > 1.0, D2S149, D11S1313, D11S987, D11S1314 including those encompassing the PTHR1 (D3S3559, D3S1289) for lumbar spine BMD and D2S149 for femoral neck BMD. Our data suggest that genes within a these chromosomal regions are contributing to a predisposition to low BMD in the Irish population.

Association of IL23R and ERAP1 genes with ankylosing spondylitis in a Portuguese population

Objective: Association between ankylosing spondylitis (AS) and two genes, ERAP1 and IL23R, has recently been reported in North American and British populations. The population attributable risk fraction for ERAP1 in this study was 25%, and for IL23R, 9%. Confirmation of these findings to ERAP1 in other ethnic groups has not yet been demonstrated. We sought to test the association between single nucleotide polymorphisms (SNPs) in these genes and susceptibility to AS among a Portuguese population. We also investigated the role of these genes in clinical manifestations of AS, including age of symptom onset, the Bath Ankylosing Spondylitis Disease Activity, Metrology and Functional Indices, and the modified Stoke Ankylosing Spondylitis Spinal Score. Methods: The study was conducted on 358 AS cases and 285 ethnically matched Portuguese healthy controls. AS was defined according to the modified New York Criteria. Genotyping of IL23R and ERAP1 allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests of genotypes as implemented in PLINK for dichotomous and quantitative variables respectively. A meta-analysis for Portuguese and previously published Spanish IL23R data was performed using the StatsDirect® Statistical tools, by fixed and random effects models. Results: A total of 14 nsSNPs markers (8 for IL23R, 5 for ERAPl, 1 for LN-PEP) were analysed. Three markers (2 for IL23R and 1 for ERAP1) showed significant single-locus disease associations, confirming that the association of these genes with AS in the Portuguese population. The strongest associated SNP in IL23R was rs1004819 (OR=1.4, p=0.0049), and in ERAP1 was rs30187 (OR=1.26, p=0.035). The population attributable risk fractions in the Portuguese population for these SNPs are 11% and 9.7% respectively. No association was seen with any SNP in LN-PEP, which flanks ERAP1 and was associated with AS in the British population. No association was seen with clinical manifestations of AS. Conclusions: These results show that IL23R and ERAP1 genes are also associated with susceptibility to AS in the Portuguese population, and that they contribute a significant proportion of the population risk for this disease.

Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation

Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K(+), inward rectifying K(+), L-type Ca(2+), and Na(+)/K(+) pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.

The return to work experiences of middle-aged Australian workers diagnosed with colorectal cancer: a matched cohort study

Background Few studies have been undertaken to understand the employment impact in patients with colorectal cancer and none in middle-aged individuals with cancer. This study described transitions in, and key factors influencing, work participation during the 12 months following a diagnosis of colorectal cancer. Methods We enrolled 239 adults during 2010 and 2011who were employed at the time of their colorectal cancer diagnosis and were prospectively followed over 12 months. They were compared to an age- and gender-matched general population group of 717 adults from the Household, Income and Labour Dynamics in Australia (HILDA) Survey. Data were collected using telephone and postal surveys. Primary outcomes included work participation at 12 months, changes in hours worked and time to work re-entry. Multivariable logistic and Cox proportional hazards models were undertaken. Results A significantly higher proportion of participants with colorectal cancer (27%) had stopped working at 12 months than participants from the comparison group (8%) (p < 0.001). Participants with cancer who returned to work took a median of 91 days off work (25–75 percentiles: 14–183 days). For participants with cancer, predictors of not working at 12 months included: being older, lower BMI and lower physical well-being. Factors related to delayed work re-entry included not being university-educated, working for an employer with more than 20 employees in a non-professional or managerial role, longer hospital stay, poorer perceived financial status and having or had chemotherapy. Conclusions In middle-adulthood, those working and diagnosed with colorectal cancer can expect to take around three months off work. Individuals treated with chemotherapy, without a university degree and from large employers could be targeted for specific assistance for a more timely work entry.

Dynamic pattern of suicide in Australia, 1986-2005: A descriptive-analytic study

Objective This study explores the spatiotemporal variations of suicide across Australia from 1986 to 2005, discusses the reasons for dynamic changes, and considers future suicide research and prevention strategies. Design Suicide (1986–2005) and population data were obtained from the Australian Bureau of Statistics. A series of analyses were conducted to examine the suicide pattern by sex, method and age group over time and geography. Results Differences in suicide rates across sex, age groups and suicide methods were found across geographical areas. Male suicides were mainly completed by hanging, firearms, gases and self-poisoning. Female suicides were primarily completed by hanging and self-poisoning. Suicide rates were higher in rural areas than in urban areas (capital cities and regional centres). Suicide rates by firearms were higher in rural areas than in urban areas, while the pattern for self-poisoning showed the reverse trend. Suicide rates had relatively stable trend for the total population and those aged between 15 and 54, while suicide decreased among 55 years and over during the study period. There was a decrease in suicides by firearms during the study period especially after 1996 when a new firearm control law was implemented, while suicide by hanging continued to increase. Areas with a high proportion of indigenous population (eg, northwest of Queensland and top north of the Northern Territory) had shown a substantial increase in suicide incidence after 1995. Conclusions Suicide rates varied over time and space and across sexes, age groups and suicide methods. This study provides detailed patterns of suicide to inform suicide control and prevention strategies for specific subgroups and areas of high and increased risk.

Comparing multilevel and Bayesian spatial random effects survival models to assess geographical inequalities in colorectal cancer survival: a case study

Background Multilevel and spatial models are being increasingly used to obtain substantive information on area-level inequalities in cancer survival. Multilevel models assume independent geographical areas, whereas spatial models explicitly incorporate geographical correlation, often via a conditional autoregressive prior. However the relative merits of these methods for large population-based studies have not been explored. Using a case-study approach, we report on the implications of using multilevel and spatial survival models to study geographical inequalities in all-cause survival. Methods Multilevel discrete-time and Bayesian spatial survival models were used to study geographical inequalities in all-cause survival for a population-based colorectal cancer cohort of 22,727 cases aged 20–84 years diagnosed during 1997–2007 from Queensland, Australia. Results Both approaches were viable on this large dataset, and produced similar estimates of the fixed effects. After adding area-level covariates, the between-area variability in survival using multilevel discrete-time models was no longer significant. Spatial inequalities in survival were also markedly reduced after adjusting for aggregated area-level covariates. Only the multilevel approach however, provided an estimation of the contribution of geographical variation to the total variation in survival between individual patients. Conclusions With little difference observed between the two approaches in the estimation of fixed effects, multilevel models should be favored if there is a clear hierarchical data structure and measuring the independent impact of individual- and area-level effects on survival differences is of primary interest. Bayesian spatial analyses may be preferred if spatial correlation between areas is important and if the priority is to assess small-area variations in survival and map spatial patterns. Both approaches can be readily fitted to geographically enabled survival data from international settings

Life-course determinants of bone mass in young adults from a transitional rural community in India: The Andhra Pradesh Children and Parents Study (APCAPS)

Background: Undernutrition and physical inactivity are both associated with lower bone mass. Objective: This study aimed to investigate the combined effects of early-life undernutrition and urbanized lifestyles in later life on bone mass accrual in young adults from a rural community in India that is undergoing rapid socioeconomic development. Design: This was a prospective cohort study of participants of the Hyderabad Nutrition Trial (1987–1990), which offered balanced protein-calorie supplementation to pregnant women and preschool children younger than 6 y in the intervention villages. The 2009–2010 follow-up study collected data on current anthropometric measures, bone mineral density (BMD) measured by dual-energy X-ray absorptiometry, blood samples, diet, physical activity, and living standards of the trial participants (n = 1446, aged 18–23 y). Results: Participants were generally lean and had low BMD [mean hip BMD: 0.83 (women), 0.95 (men) g/cm2; lumbar spine: 0.86 (women), 0.93 (men) g/cm2]. In models adjusted for current risk factors, no strong evidence of a positive association was found between BMD and early-life supplementation. On the other hand, current lean mass and weight-bearing physical activity were positively associated with BMD. No strong evidence of an association was found between BMD and current serum 25-hydroxyvitamin D or dietary intake of calcium, protein, or calories. Conclusions: Current lean mass and weight-bearing physical activity were more important determinants of bone mass than was early-life undernutrition in this population. In transitional rural communities from low-income countries, promotion of physical activity may help to mitigate any potential adverse effects of early nutritional disadvantage.

5-year Chlamydia vaccination programme could reverse disease-related koala population decline: Predictions from a mathematical model using field data

BACKGROUND Many koala populations around Australia are in serious decline, with a substantial component of this decline in some Southeast Queensland populations attributed to the impact of Chlamydia. A Chlamydia vaccine for koalas is in development and has shown promise in early trials. This study contributes to implementation preparedness by simulating vaccination strategies designed to reverse population decline and by identifying which age and sex category it would be most effective to target. METHODS We used field data to inform the development and parameterisation of an individual-based stochastic simulation model of a koala population endemic with Chlamydia. The model took into account transmission, morbidity and mortality caused by Chlamydia infections. We calibrated the model to characteristics of typical Southeast Queensland koala populations. As there is uncertainty about the effectiveness of the vaccine in real-world settings, a variety of potential vaccine efficacies, half-lives and dosing schedules were simulated. RESULTS Assuming other threats remain constant, it is expected that current population declines could be reversed in around 5-6 years if female koalas aged 1-2 years are targeted, average vaccine protective efficacy is 75%, and vaccine coverage is around 10% per year. At lower vaccine efficacies the immunological effects of boosting become important: at 45% vaccine efficacy population decline is predicted to reverse in 6 years under optimistic boosting assumptions but in 9 years under pessimistic boosting assumptions. Terminating a successful vaccination programme at 5 years would lead to a rise in Chlamydia prevalence towards pre-vaccination levels. CONCLUSION For a range of vaccine efficacy levels it is projected that population decline due to endemic Chlamydia can be reversed under realistic dosing schedules, potentially in just 5 years. However, a vaccination programme might need to continue indefinitely in order to maintain Chlamydia prevalence at a sufficiently low level for population growth to continue.

The prevalence, predictors and associations of hypertension in Sri Lanka: A cross-sectional population based national survey

We studied the community prevalence, patterns and predictors of hypertension in a large sub-population of South Asian adults with a view of identifying differential risk factors. Data were collected between years 2005-2006 and 5000 adults were invited for the study. The sample size was 4485, and about 39.5% were males. Mean systolic and diastolic blood pressures were 127.1 ± 19.8 mmHg and 75.4 ± 11.3 mmHg, respectively. Age-adjusted prevalence in all adults, males and females was 23.7%, 23.4% and 23.8%, respectively. Urban adults had a significantly higher prevalence of hypertension than rural adults. In the binary logistic-regression analysis, male gender (OR: 1.2), increasing age, Sri Lankan Moor ethnicity (OR: 1.6), physical inactivity (OR: 1.7), presence of diabetes (OR: 2.2) and central obesity (OR: 2.3) all were significantly associated with hypertension. In conclusion, nearly one-third of the Sri Lankan adult population is hypertensive. Hence, public health initiatives should encourage healthier lifestyles with emphasis on preventing obesity and increasing physical activity.

Polymorphisms in a desaturase 2 ortholog associate with cuticular hydrocarbon and male mating success variation in a natural population of Drosophila serrata

Elucidating the nature of genetic variation underlying both sexually selected traits and the fitness components of sexual selection is essential to understanding the broader consequences of sexual selection as an evolutionary process. To date, there have been relatively few attempts to connect the genetic variance in sexually selected traits with segregating DNA sequence polymorphisms. We set out to address this in a well-characterized sexual selection system - the cuticular hydrocarbons (CHCs) of Drosophila serrata - using an indirect association study design that allowed simultaneous estimation of the genetic variance in CHCs, sexual fitness and single nucleotide polymorphism (SNP) effects in an outbred population. We cloned and sequenced an ortholog of the D. melanogaster desaturase 2 gene, previously shown to affect CHC biosynthesis in D. melanogaster, and associated 36 SNPs with minor allele frequencies > 0.02 with variance in CHCs and sexual fitness. Three SNPs had significant multivariate associations with CHC phenotype (q-value < 0.05). At these loci, minor alleles had multivariate effects on CHCs that were weakly associated with the multivariate direction of sexual selection operating on these traits. Two of these SNPs had pleiotropic associations with male mating success, suggesting these variants may underlie responses to sexual selection due to this locus. There were 15 significant male mating success associations (q-value < 0.1), and interestingly, we detected a nonrandom pattern in the relationship between allele frequency and direction of effect on male mating success. The minor-frequency allele usually reduced male mating success, suggesting a positive association between male mating success and total fitness at this locus.

Systematic and Day-to-Day Effects of Chemical-Derived Population Estimates on Wastewater-Based Drug Epidemiology

Population size is crucial when estimating population-normalized drug consumption (PNDC) from wastewater-based drug epidemiology (WBDE). Three conceptually different population estimates can be used: de jure (common census, residence), de facto (all persons within a sewer catchment), and chemical loads (contributors to the sampled wastewater). De facto and chemical loads will be the same where all households contribute to a central sewer system without wastewater loss. This study explored the feasibility of determining a de facto population and its effect on estimating PNDC in an urban community over an extended period. Drugs and other chemicals were analyzed in 311 daily composite wastewater samples. The daily estimated de facto population (using chemical loads) was on average 32% higher than the de jure population. Consequently, using the latter would systemically overestimate PNDC by 22%. However, the relative day-to-day pattern of drug consumption was similar regardless of the type of normalization as daily illicit drug loads appeared to vary substantially more than the population. Using chemical loads population, we objectively quantified the total methodological uncertainty of PNDC and reduced it by a factor of 2. Our study illustrated the potential benefits of using chemical loads population for obtaining more robust PNDC data in WBDE.

A Model to Estimate the Population Contributing to the Wastewater Using Samples Collected on Census Day

An important uncertainty when estimating per capita consumption of, for example, illicit drugs by means of wastewater analysis (sometimes referred to as “sewage epidemiology”) relates to the size and variability of the de facto population in the catchment of interest. In the absence of a day-specific direct population count any indirect surrogate model to estimate population size lacks a standard to assess associated uncertainties. Therefore, the objective of this study was to collect wastewater samples at a unique opportunity, that is, on a census day, as a basis for a model to estimate the number of people contributing to a given wastewater sample. Mass loads for a wide range of pharmaceuticals and personal care products were quantified in influents of ten sewage treatment plants (STP) serving populations ranging from approximately 3500 to 500 000 people. Separate linear models for population size were estimated with the mass loads of the different chemical as the explanatory variable: 14 chemicals showed good, linear relationships, with highest correlations for acesulfame and gabapentin. De facto population was then estimated through Bayesian inference, by updating the population size provided by STP staff (prior knowledge) with measured chemical mass loads. Cross validation showed that large populations can be estimated fairly accurately with a few chemical mass loads quantified from 24-h composite samples. In contrast, the prior knowledge for small population sizes cannot be improved substantially despite the information of multiple chemical mass loads. In the future, observations other than chemical mass loads may improve this deficit, since Bayesian inference allows including any kind of information relating to population size.