Association of variants in MMEL1 and CTLA4 with rheumatoid arthritis in the Han Chinese population

Background: The genome-wide association study era has made great progress in identifying susceptibility genes and genetic loci for rheumatoid arthritis (RA) in populations of White European ancestry. However, few studies have tried to dissect disease aetiopathogenesis in other ethnic populations. Objective: To investigate these associations in the Han Chinese population. Methods: Haplotypes from the HapMap database Chinese population were used to select tag-single-nucleotide polymorphisms (SNPs) (r2 =0.8) across 19 distinct RA genomic regions. A two phase case-control association study was performed, with 169 SNPs genotyped in phase I (n=571 cases, n=880 controls), and 64 SNPs achieving p<0.2 in the first phase being genotyped in phase II (n=464 cases, n=822 controls). Association statistics were calculated using permutation tests both unadjusted and adjusted for the number of markers studied. Results: Robust association was detected for MMEL1 and CTLA4 , and modest association was identified for another six loci: PADI4 , STAT4 , PRDM1 , CDK6 , TRAF1-C5 and KIF5A-PIP4K2C. All three markers genotyped in MMEL1 demonstrated association, with peak signal for rs3890745 (p=2.6×10 -5unadjusted, p=0.003 adjusted, OR=0.79). For CTLA4 , significance was detected for three of five variants showing association, with peak association for marker rs12992492 (p=4.3×10-5 unadjusted, p=0.0021 adjusted, OR=0.77). Lack of association of common variants in PTPN22 with RA in Han Chinese was confirmed. Conclusion: This study identifies MMEL1 and CTLA4 as RA susceptibility genes, provides suggestive evidence of association for a further six loci in the Han Chinese population and confirms lack of PTPN22 association in Asian populations. It also confirms the value of multiethnic population studies to help dissect disease aetiopathogenesis.

High-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population

Objective Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. Methods A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. Results We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly149Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly149Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). Conclusion This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

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.