Pseudo-marginal algorithms with multiple CPUs

The emergence of pseudo-marginal algorithms has led to improved computational efficiency for dealing with complex Bayesian models with latent variables. Here an unbiased estimator of the likelihood replaces the true likelihood in order to produce a Bayesian algorithm that remains on the marginal space of the model parameter (with latent variables integrated out), with a target distribution that is still the correct posterior distribution. Very efficient proposal distributions can be developed on the marginal space relative to the joint space of model parameter and latent variables. Thus psuedo-marginal algorithms tend to have substantially better mixing properties. However, for pseudo-marginal approaches to perform well, the likelihood has to be estimated rather precisely. This can be difficult to achieve in complex applications. In this paper we propose to take advantage of multiple central processing units (CPUs), that are readily available on most standard desktop computers. Here the likelihood is estimated independently on the multiple CPUs, with the ultimate estimate of the likelihood being the average of the estimates obtained from the multiple CPUs. The estimate remains unbiased, but the variability is reduced. We compare and contrast two different technologies that allow the implementation of this idea, both of which require a negligible amount of extra programming effort. The superior performance of this idea over the standard approach is demonstrated on simulated data from a stochastic volatility model.

The institutional limits to, and possibility of, understanding young people diagnosed with 'mental disorders', within inclusive teacher-student relationships

A persistent pattern of exclusion of young people with ‘mental disorders’ from school systems, despite the best intentions of schools and teachers, has prompted a call for a more reflexive understanding of their behaviours. This thesis, by describing how institutionally recognised ways of understanding can result in otherwise avoidable moral collisions and exclusion, produces new insights into the nature and processes of understanding required to promote inclusion. These insights were produced through an intensive qualitative examination of a violent classroom episode, identifying key points in the interaction that could make the difference between misrecognition and recognition, turning exclusion into inclusion.

Ingenuity pathways analysis of urine metabonomics phenotypes toxicity of gentamicin in multiple organs

We introduce the use of Ingenuity Pathway Analysis to analyzing global metabonomics in order to characterize phenotypically biochemical perturbations and the potential mechanisms of the gentamicin-induced toxicity in multiple organs. A single dose of gentamicin was administered to Sprague Dawley rats (200 mg/kg, n = 6) and urine samples were collected at -24-0 h pre-dosage, 0-24, 24-48, 48-72 and 72-96 h post-dosage of gentamicin. The urine metabonomics analysis was performed by UPLC/MS, and the mass spectra signals of the detected metabolites were systematically deconvoluted and analyzed by pattern recognition analyses (Heatmap, PCA and PLS-DA), revealing a time-dependency of the biochemical perturbations induced by gentamicin toxicity. As result, the holistic metabolome change induced by gentamicin toxicity in the animal's organisms was characterized. Several metabolites involved in amino acid metabolism were identified in urine, and it was confirmed that gentamicin biochemical perturbations can be foreseen from these biomarkers. Notoriously, it was found that gentamicin induced toxicity in multiple organs system in the laboratory rats. The proof-of-knowledge based Ingenuity Pathway Analysis revealed gentamicin induced liver and heart toxicity, along with the previously known toxicity in kidney. The metabolites creatine, nicotinic acid, prostaglandin E2, and cholic acid were identified and validated as phenotypic biomarkers of gentamicin induced toxicity. Altogether, the significance of the use of metabonomics analyses in the assessment of drug toxicity is highlighted once more; furthermore, this work demonstrated the powerful predictive potential of the Ingenuity Pathway Analysis to study of drug toxicity and its valuable complementation for metabonomics based assessment of the drug toxicity.

Postpartum care for Aboriginal and non-Aboriginal women with Gestational Diabetes Mellitus across urban, rural and remote locations : a protocol for a cohort linkage study

Background: Gestational diabetes mellitus (GDM) is increasing, along with obesity and type 2 diabetes (T2DM), with Aboriginal and Torres Strait Islander people* in Australia particularly affected. GDM causes serious complications in pregnancy, birth, and the longer term, for both women and their infants. Women diagnosed with GDM have an eightfold risk of developing T2DM after pregnancy, compared with women who have not had GDM. Indigenous women have an even higher risk, at a younger age, and progress more quickly from GDM to T2DM, compared to non-Indigenous women. If left undetected and untreated, T2DM can lead to heart disease, stroke, renal disease, kidney failure, amputations and blindness. A GDM diagnosis offers a ‘window of opportunity’ for diabetes health interventions and it is vital that acceptable and effective prevention, treatment, and post-pregnancy care are provided. Low rates of post-pregnancy screening for T2DM are reported among non-Aboriginal women in Australia and among Indigenous women in other countries, however data for Aboriginal women are scarce. Breastfeeding, a healthy diet, and exercise can also help to prevent T2DM, and together with T2DM screening are recommended elements of ‘post-pregnancy care’ for women with GDM, This paper describes methods for a data linkage study to investigate rates of post-pregnancy care among women with GDM. Methods/Design: This retrospective cohort includes all women who gave birth at Cairns Base Hospital in Far North Queensland, Australia, from 2004 to 2010, coded as having GDM in the Cairns Base Hospital Clinical Coding system. Data linkage is being conducted with the Queensland Perinatal Data Collection, and three laboratories. Hospital medical records are being reviewed to validate the accuracy of GDM case ascertainment, and gather information on breastfeeding and provision of dietary advice. Multiple logistic regression is being used to compare post-pregnancy care between Aboriginal and non-Aboriginal women, while adjusting for other factors may impact on post-pregnancy care. Survival analysis is being used to estimate the rates of progression from GDM to T2DM. Discussion: There are challenges to collecting post-pregnancy data for women with GDM. However, research is urgently needed to ensure adequate post-pregnancy care is provided for women with GDM in Australia.

Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10−11, OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

The evolution and biogeography of the austral horse fly tribe Scionini (Diptera: Tabanidae: Pangoniinae) inferred from multiple mitochondrial and nuclear genes

Phylogenetic relationships within the Tabanidae are largely unknown, despite their considerable medical and ecological importance. The first robust phylogenetic hypothesis for the horse fly tribe Scionini is provided, completing the systematic placement of all tribes in the subfamily Pangoniinae. The Scionini consists of seven mostly southern hemisphere genera distributed in Australia, New Guinea, New Zealand and South America. A 5757. bp alignment of 6 genes, including mitochondrial (COI and COII), ribosomal (28S) and nuclear (AATS and CAD regions 1, 3 and 4) genes, was analysed for 176 taxa using both Bayesian and maximum likelihood approaches. Results indicate the Scionini are strongly monophyletic, with the exclusion of the only northern hemisphere genus Goniops. The South American genera Fidena, Pityocera and Scione were strongly monophyletic, corresponding to current morphology-based classification schemes. The most widespread genus Scaptia was paraphyletic and formed nine strongly supported monophyletic clades, each corresponding to either the current subgenera or several previously synonymised genera that should be formally resurrected. Molecular results also reveal a newly recognised genus endemic to New Zealand, formerly placed within Scaptia. Divergence time estimation was employed to assess the global biogeographical patterns in the Pangoniinae. These analyses demonstrated that the Scionini are a typical Gondwanan group whose diversification was influenced by the fragmentation of that ancient land mass. Furthermore, results indicate that the Scionini most likely originated in Australia and subsequently radiated to New Zealand and South American by both long distance dispersal and vicariance. The phylogenetic framework of the Scionini provided herein will be valuable for taxonomic revisions of the Tabanidae.

Circulating microRNAs involved in multiple sclerosis

Multiple sclerosis (MS) is an immune-mediated, demyelinating and neurodegenerative disease of the central nervous system. After traumatic brain injury, it is the leading cause of neurology disability in young adults. Considerable advances have been made in identifying genes involved in MS but the genetic and phenotypic complexity associated with this disease significantly hinders any progress. A novel class of small RNA molecules, microRNAs (miRNAs) has acquired much attention because they regulate the expression of up to 30% of protein-coding genes and may play a pivotal role in the development of many, if not all, complex diseases. Seven published studies investigated miRNAs from peripheral blood mononuclear cells, CD4+, CD8+ T cell, B lymphocytes, peripheral blood leukocytes, whole blood and brain astrocytes with MS risk. The absence of MS studies investigating plasma miRNA prompted the current investigation of identifying a circulating miRNA signature in MS. We conducted a microarray analysis of over 900 known miRNA transcripts from plasma samples collected from four MS individuals and four sex-aged and ethnicity matched healthy controls. We identified six plasma miRNA (miR-614, miR-572, miR-648, miR-1826, miR-422a and miR-22) that were significantly up-regulated and one plasma miRNA (miR-1979) that was significantly down-regulated in MS individuals. Both miR-422a and miR-22 have previously been implicated in MS. The present study is the first to show a circulating miRNA signature involved in MS that could serve as a potential prognostic and diagnostic biomarker for MS.

Familial recurrence risks for multiple sclerosis in Australia

BACKGROUND: Genetic susceptibility to multiple sclerosis (MS) has been recognised for many years. Considerable data exist from the northern hemisphere regarding the familial recurrence risks for MS, but there are few data for the southern hemisphere and regions at lower latitude such as Australia. To investigate the interaction between environmental and genetic causative factors in MS, the authors undertook a familial recurrence risk study in three latitudinally distinct regions of Australia. METHODS: Immediate and extended family pedigrees have been collected for three cohorts of people with MS in Queensland, Victoria and Tasmania spanning 15° of latitude. Age of onset data from Queensland were utilised to estimate age-adjusted recurrence rates. RESULTS: Recurrence risks in Australia were significantly lower than in studies from northern hemisphere populations. The age-adjusted risk for siblings across Australia was 2.13% compared with 3.5% for the northern hemisphere. A similar pattern was seen for other relatives. The risks to relatives were proportional to the population risks for each site, and hence the sibling recurrence-risk ratio (λ(s)) was similar across all sites. DISCUSSION: The familial recurrence risk of MS in Australia is lower than in previously reported studies. This is directly related to the lower population prevalence of MS. The overall genetic susceptibility in Australia as measured by the λ(s) is similar to the northern hemisphere, suggesting that the difference in population risk is explained largely by environmental factors rather than by genetic admixture.

Genome-wide meta-analysis identifies novel multiple sclerosis susceptibility loci

Objective: To perform a 1-stage meta-analysis of genome-wide association studies (GWAS) of multiple sclerosis (MS) susceptibility and to explore functional consequences of new susceptibility loci. Methods: We synthesized 7 MS GWAS. Each data set was imputed using HapMap phase II, and a per single nucleotide polymorphism (SNP) meta-analysis was performed across the 7 data sets. We explored RNA expression data using a quantitative trait analysis in peripheral blood mononuclear cells (PBMCs) of 228 subjects with demyelinating disease. Results: We meta-analyzed 2,529,394 unique SNPs in 5,545 cases and 12,153 controls. We identified 3 novel susceptibility alleles: rs170934T at 3p24.1 (odds ratio [OR], 1.17; p ¼ 1.6 � 10�8) near EOMES, rs2150702G in the second intron of MLANA on chromosome 9p24.1 (OR, 1.16; p ¼ 3.3 � 10�8), and rs6718520A in an intergenic region on chromosome 2p21, with THADA as the nearest flanking gene (OR, 1.17; p ¼ 3.4 � 10�8). The 3 new loci do not have a strong cis effect on RNA expression in PBMCs. Ten other susceptibility loci had a suggestive p < 1 � 10�6, some of these loci have evidence of association in other inflammatory diseases (ie, IL12B, TAGAP, PLEK, and ZMIZ1). Interpretation: We have performed a meta-analysis of GWAS in MS that more than doubles the size of previous gene discovery efforts and highlights 3 novel MS susceptibility loci. These and additional loci with suggestive evidence of association are excellent candidates for further investigations to refine and validate their role in the genetic architecture of MS.

Modeling the cumulative genetic risk for multiple sclerosis from genome-wide association data

Background: Multiple sclerosis (MS) is the most common cause of chronic neurologic disability beginning in early to middle adult life. Results from recent genome-wide association studies (GWAS) have substantially lengthened the list of disease loci and provide convincing evidence supporting a multifactorial and polygenic model of inheritance. Nevertheless, the knowledge of MS genetics remains incomplete, with many risk alleles still to be revealed. Methods: We used a discovery GWAS dataset (8,844 samples, 2,124 cases and 6,720 controls) and a multi-step logistic regression protocol to identify novel genetic associations. The emerging genetic profile included 350 independent markers and was used to calculate and estimate the cumulative genetic risk in an independent validation dataset (3,606 samples). Analysis of covariance (ANCOVA) was implemented to compare clinical characteristics of individuals with various degrees of genetic risk. Gene ontology and pathway enrichment analysis was done using the DAVID functional annotation tool, the GO Tree Machine, and the Pathway-Express profiling tool. Results: In the discovery dataset, the median cumulative genetic risk (P-Hat) was 0.903 and 0.007 in the case and control groups, respectively, together with 79.9% classification sensitivity and 95.8% specificity. The identified profile shows a significant enrichment of genes involved in the immune response, cell adhesion, cell communication/ signaling, nervous system development, and neuronal signaling, including ionotropic glutamate receptors, which have been implicated in the pathological mechanism driving neurodegeneration. In the validation dataset, the median cumulative genetic risk was 0.59 and 0.32 in the case and control groups, respectively, with classification sensitivity 62.3% and specificity 75.9%. No differences in disease progression or T2-lesion volumes were observed among four levels of predicted genetic risk groups (high, medium, low, misclassified). On the other hand, a significant difference (F = 2.75, P = 0.04) was detected for age of disease onset between the affected misclassified as controls (mean = 36 years) and the other three groups (high, 33.5 years; medium, 33.4 years; low, 33.1 years). Conclusions: The results are consistent with the polygenic model of inheritance. The cumulative genetic risk established using currently available genome-wide association data provides important insights into disease heterogeneity and completeness of current knowledge in MS genetics.

A polymorphism in the HLA-DPB1 gene is associated with susceptibility to multiple sclerosis

We conducted an association study across the human leukocyte antigen (HLA) complex to identify loci associated with multiple sclerosis (MS). Comparing 1927 SNPs in 1618 MS cases and 3413 controls of European ancestry, we identified seven SNPs that were independently associated with MS conditional on the others (each ). All associations were significant in an independent replication cohort of 2212 cases and 2251 controls () and were highly significant in the combined dataset (). The associated SNPs included proxies for HLA-DRB1*15:01 and HLA-DRB1*03:01, and SNPs in moderate linkage disequilibrium (LD) with HLA-A*02:01, HLA-DRB1*04:01 and HLA-DRB1*13:03. We also found a strong association with rs9277535 in the class II gene HLA-DPB1 (discovery set , replication set , combined ). HLA-DPB1 is located centromeric of the more commonly typed class II genes HLA-DRB1, -DQA1 and -DQB1. It is separated from these genes by a recombination hotspot, and the association is not affected by conditioning on genotypes at DRB1, DQA1 and DQB1. Hence rs9277535 represents an independent MS-susceptibility locus of genome-wide significance. It is correlated with the HLA-DPB1*03:01 allele, which has been implicated previously in MS in smaller studies. Further genotyping in large datasets is required to confirm and resolve this association.

Investigation of the [−/A]8 and C1236T genetic variations within the human toll-like receptor 3 gene for association with multiple sclerosis

Multiple sclerosis (MS) is a serious cause of neurological disability among young adults. The clinical course remains difficult to predict, and the pathogenesis of the disease is still modestly understood. Autoimmunity is thought to be a key aspect of the disease, with autoreactive T cells thought to mediate central nervous system (CNS) inflammation to some extent. Toll-like receptors are known to mediate cellular recognition of pathogens by way of patterns of molecular presentation. Toll-like receptor 3 is coded by the gene TLR3 and is recognized as an important factor in virus recognition and is known to be involved in the expression of neuroprotective mediators. We set out to investigate two variations within the TLR3 gene, an 8 bp insertion-deletion \[-/A](8) and a single base-pair variation C1236T, in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We used capillary gel electrophoresis and TaqMan genotyping assay techniques to resolve genotypes for each marker, respectively. Our work found no significant difference between frequencies for TLR3 \[-/A](8) by genotype (chi(2)=1.03, p=0.60) or allele (chi(2)=1.09, p=0.30). Similarly, we found no evidence for the association of TLR3 C1236T by genotype (chi(2)=0.35, p=0.84) or allele frequency (chi(2)=0.31, p=0.58). This work reveals no evidence to suggest that these markers are associated with MS in the tested population. Although the role of TLR3 and the wider toll-like receptor family remain significant in neurological and CNS inflammatory disorders, our current work does not support a role for the two tested variants in this gene with regard to MS susceptibility.

Multiple sclerosis susceptibility-associated SNPs do not influence disease severity measures in a cohort of Australian MS patients

Recent association studies in multiple sclerosis (MS) have identified and replicated several single nucleotide polymorphism (SNP) susceptibility loci including CLEC16A, IL2RA, IL7R, RPL5, CD58, CD40 and chromosome 12q13–14 in addition to the well established allele HLA-DR15. There is potential that these genetic susceptibility factors could also modulate MS disease severity, as demonstrated previously for the MS risk allele HLA-DR15. We investigated this hypothesis in a cohort of 1006 well characterised MS patients from South-Eastern Australia. We tested the MS-associated SNPs for association with five measures of disease severity incorporating disability, age of onset, cognition and brain atrophy. We observed trends towards association between the RPL5 risk SNP and time between first demyelinating event and relapse, and between the CD40 risk SNP and symbol digit test score. No associations were significant after correction for multiple testing. We found no evidence for the hypothesis that these new MS disease risk-associated SNPs influence disease severity.

An investigation of the C77G and C772T variations within the human protein tyrosine phosphatase receptor type C gene for association with multiple sclerosis in an Australian population

Multiple sclerosis (MS) is a common cause of neurological disability in young adults. The disease generally manifests in early to middle adulthood and causes various neurological deficits. Autoreactive T lymphocytes and their associated antigens have long been presumed important features of MS pathogenesis. The Protein tyrosine phosphatase receptor type C gene (PTPRC) encodes the T-cell receptor CD45. Variations within PTPRC have been previously associated with diseases of autoimmune origin such as type 1 diabetes mellitus and Graves' disease. We set out to investigate two variants within the PTPRC gene, C77G and C772T in subjects with MS and matched healthy controls to determine whether significant differences exist in these markers in an Australian population. We employed high resolution melt analysis (HRM) and restriction length polymorphism (RFLP) techniques to determine genotypic and allelic frequencies. Our study found no significant difference between frequencies for PTPRC C77G by either genotype (Χ2 = 0.65, P = 0.72) or allele (Χ2 = 0.48, P = 0.49). Similarly, we did not find evidence to suggest an association between PTPRC C772T by genotype (Χ2 = 1.06, P = 0.59) or allele (Χ2 = 0.20, P = 0.66). Linkage disequilibrium (LD) analysis showed strong linkage disequilibrium between the two tested markers (D' = 0.9970, SD = 0.0385). This study reveals no evidence to suggest that these markers are associated with MS in the tested Australian Caucasian population. Although the PTPRC gene has a significant role in regulating CD4+ and CD8+ autoreactive T-cells, interferon-beta responsiveness, and potentially other important processes, our study does not support a role for the two tested variants of this gene in MS susceptibility in the Australian population.

Genome-wide association study identifies new multiple sclerosis susceptibility loci on chromosomes 12 and 20

To identify multiple sclerosis (MS) susceptibility loci, we conducted a genome-wide association study (GWAS) in 1,618 cases and used shared data for 3,413 controls. We performed replication in an independent set of 2,256 cases and 2,310 controls, for a total of 3,874 cases and 5,723 controls. We identified risk-associated SNPs on chromosome 12q13-14 (rs703842, P = 5.4 x 10(-11); rs10876994, P = 2.7 x 10(-10); rs12368653, P = 1.0 x 10(-7)) and upstream of CD40 on chromosome 20q13 (rs6074022, P = 1.3 x 10(-7); rs1569723, P = 2.9 x 10(-7)). Both loci are also associated with other autoimmune diseases. We also replicated several known MS associations (HLA-DR15, P = 7.0 x 10(-184); CD58, P = 9.6 x 10(-8); EVI5-RPL5, P = 2.5 x 10(-6); IL2RA, P = 7.4 x 10(-6); CLEC16A, P = 1.1 x 10(-4); IL7R, P = 1.3 x 10(-3); TYK2, P = 3.5 x 10(-3)) and observed a statistical interaction between SNPs in EVI5-RPL5 and HLA-DR15 (P = 0.001).