Discovery of candidate serum proteomic and metabolomic biomarkers in ankylosing spondylitis

Ankylosing Spondylitis (AS) is a common inflammatory rheumatic disease with a predilection for the axial skeleton, affecting 0.2% of the population. Current diagnostic criteria rely on a composite of clinical and radiological changes, with a mean time to diagnosis of 5 to 10 years. In this study we employed nano liquid-chromatography mass spectrometry analysis to detect and quantify proteins and small compounds including endogenous peptides and metabolites in serum from 18 AS patients and nine healthy individuals. We identified a total of 316 proteins in serum, of which 22 showed significant up- or down-regulation (p < 0.05) in AS patients. Receiver operating characteristic analysis of combined levels of serum amyloid P component and inter-α-trypsin inhibitor heavy chain 1 revealed high diagnostic value for Ankylosing Spondylitis (area under the curve = 0.98). We also depleted individual sera of proteins to analyze endogenous peptides and metabolic compounds. We detected more than 7000 molecular features in patients and healthy individuals. Quantitative MS analysis revealed compound profiles that correlate with the clinical assessment of disease activity. One molecular feature identified as a Vitamin D3 metabolite-(23S,25R)-25-hydroxyvitamin D3 26,23-peroxylactone-was down-regulated in AS. The ratio of this vitamin D metabolite versus vitamin D binding protein serum levels was also altered in AS as compared with controls. These changes may contribute to pathological skeletal changes in AS. Our study is the first example of an integration of proteomic and metabolomic techniques to find new biomarker candidates for the diagnosis of Ankylosing Spondylitis

Convergent evolution of marine mammals is associated with distinct substitutions in common genes

Phenotypic convergence is thought to be driven by parallel substitutions coupled with natural selection at the sequence level. Multiple independent evolutionary transitions of mammals to an aquatic environment offer an opportunity to test this thesis. Here, whole genome alignment of coding sequences identified widespread parallel amino acid substitutions in marine mammals; however, the majority of these changes were not unique to these animals. Conversely, we report that candidate aquatic adaptation genes, identified by signatures of likelihood convergence and/or elevated ratio of nonsynonymous to synonymous nucleotide substitution rate, are characterized by very few parallel substitutions and exhibit distinct sequence changes in each group. Moreover, no significant positive correlation was found between likelihood convergence and positive selection in all three marine lineages. These results suggest that convergence in protein coding genes associated with aquatic lifestyle is mainly characterized by independent substitutions and relaxed negative selection.

Modeling strategies in longitudinal data analysis: Covariate, variance function and correlation structure selection

A modeling paradigm is proposed for covariate, variance and working correlation structure selection for longitudinal data analysis. Appropriate selection of covariates is pertinent to correct variance modeling and selecting the appropriate covariates and variance function is vital to correlation structure selection. This leads to a stepwise model selection procedure that deploys a combination of different model selection criteria. Although these criteria find a common theoretical root based on approximating the Kullback-Leibler distance, they are designed to address different aspects of model selection and have different merits and limitations. For example, the extended quasi-likelihood information criterion (EQIC) with a covariance penalty performs well for covariate selection even when the working variance function is misspecified, but EQIC contains little information on correlation structures. The proposed model selection strategies are outlined and a Monte Carlo assessment of their finite sample properties is reported. Two longitudinal studies are used for illustration.

Working correlation structure misspecification, estimation and covariate design: Implications for generalised estimating equations performance

The method of generalised estimating equations for regression modelling of clustered outcomes allows for specification of a working matrix that is intended to approximate the true correlation matrix of the observations. We investigate the asymptotic relative efficiency of the generalised estimating equation for the mean parameters when the correlation parameters are estimated by various methods. The asymptotic relative efficiency depends on three-features of the analysis, namely (i) the discrepancy between the working correlation structure and the unobservable true correlation structure, (ii) the method by which the correlation parameters are estimated and (iii) the 'design', by which we refer to both the structures of the predictor matrices within clusters and distribution of cluster sizes. Analytical and numerical studies of realistic data-analysis scenarios show that choice of working covariance model has a substantial impact on regression estimator efficiency. Protection against avoidable loss of efficiency associated with covariance misspecification is obtained when a 'Gaussian estimation' pseudolikelihood procedure is used with an AR(1) structure.

Spectral analysis of pair-correlation bandwidth: Application to cell biology images

Images from cell biology experiments often indicate the presence of cell clustering, which can provide insight into the mechanisms driving the collective cell behaviour. Pair-correlation functions provide quantitative information about the presence, or absence, of clustering in a spatial distribution of cells. This is because the pair-correlation function describes the ratio of the abundance of pairs of cells, separated by a particular distance, relative to a randomly distributed reference population. Pair-correlation functions are often presented as a kernel density estimate where the frequency of pairs of objects are grouped using a particular bandwidth (or bin width), Δ>0. The choice of bandwidth has a dramatic impact: choosing Δ too large produces a pair-correlation function that contains insufficient information, whereas choosing Δ too small produces a pair-correlation signal dominated by fluctuations. Presently, there is little guidance available regarding how to make an objective choice of Δ. We present a new technique to choose Δ by analysing the power spectrum of the discrete Fourier transform of the pair-correlation function. Using synthetic simulation data, we confirm that our approach allows us to objectively choose Δ such that the appropriately binned pair-correlation function captures known features in uniform and clustered synthetic images. We also apply our technique to images from two different cell biology assays. The first assay corresponds to an approximately uniform distribution of cells, while the second assay involves a time series of images of a cell population which forms aggregates over time. The appropriately binned pair-correlation function allows us to make quantitative inferences about the average aggregate size, as well as quantifying how the average aggregate size changes with time.

The designed environment and how it affects brain morphology and mental health

Background The environment is inextricably related to mental health. Recent research replicates findings of a significant, linear correlation between a childhood exposure to the urban environment and psychosis. Related studies also correlate the urban environment and aberrant brain morphologies. These findings challenge common beliefs that the mind and brain remain neutral in the face of worldly experience. Aim There is a signature within these neurological findings that suggests that specific features of design cause and trigger mental illness. The objective in this article is to work backward from the molecular dynamics to identify features of the designed environment that may either trigger mental illness or protect against it. Method This review analyzes the discrete functions putatively assigned to the affected brain areas and a neurotransmitter called dopamine, which is the primary target of most antipsychotic medications. The intention is to establish what the correlations mean in functional terms, and more specifically, how this relates to the phenomenology of urban experience. In doing so, environmental mental illness risk factors are identified. Conclusions Having established these relationships, the review makes practical recommendations for those in public health who wish to use the environment itself as a tool to improve the mental health of a community through design.

Shared genetic basis for migraine and ischemic stroke: A genome-wide analysis of common variants

OBJECTIVE To quantify genetic overlap between migraine and ischemic stroke (IS) with respect to common genetic variation. METHODS We applied 4 different approaches to large-scale meta-analyses of genome-wide data on migraine (23,285 cases and 95,425 controls) and IS (12,389 cases and 62,004 controls). First, we queried known genome-wide significant loci for both disorders, looking for potential overlap of signals. We then analyzed the overall shared genetic load using polygenic scores and estimated the genetic correlation between disease subtypes using data derived from these models. We further interrogated genomic regions of shared risk using analysis of covariance patterns between the 2 phenotypes using cross-phenotype spatial mapping. RESULTS We found substantial genetic overlap between migraine and IS using all 4 approaches. Migraine without aura (MO) showed much stronger overlap with IS and its subtypes than migraine with aura (MA). The strongest overlap existed between MO and large artery stroke (LAS; p = 6.4 x 10(-28) for the LAS polygenic score in MO) and between MO and cardioembolic stroke (CE; p = 2.7 x 10(-20) for the CE score in MO). CONCLUSIONS Our findings indicate shared genetic susceptibility to migraine and IS, with a particularly strong overlap between MO and both LAS and CE pointing towards shared mechanisms. Our observations on MA are consistent with a limited role of common genetic variants in this subtype.

Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 x 10(-)(9), odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 x 10(-)(1)(1) (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 x 10(-)(5), permuted threshold for genome-wide significance 7.7 x 10(-)(5). To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Common genetic influences underlie comorbidity of migraine and endometriosis

We examined the co-occurrence of migraine and endometriosis within the largest known collection of families containing multiple women with surgically confirmed endometriosis and in an independent sample of 815 monozygotic and 457 dizygotic female twin pairs. Within the endometriosis families, a significantly increased risk of migrainous headache was observed in women with endometriosis compared to women without endometriosis (odds ratio [OR] 1.57, 95% confidence interval [CI]: 1.12-2.21, P=0.009). Bivariate heritability analyses indicated no evidence for common environmental factors influencing either migraine or endometriosis but significant genetic components for both traits, with heritability estimates of 69 and 49%, respectively. Importantly, a significant additive genetic correlation (r(G) = 0.27, 95% CI: 0.06-0.47) and bivariate heritability (h(2)=0.17, 95% CI: 0.08-0.27) was observed between migraine and endometriosis. Controlling for the personality trait neuroticism made little impact on this association. These results confirm the previously reported comorbidity between migraine and endometriosis and indicate common genetic influences completely explain their co-occurrence within individuals. Given pharmacological treatments for endometriosis typically target hormonal pathways and a number of findings provide support for a relationship between hormonal variations and migraine, hormone-related genes and pathways are highly plausible candidates for both migraine and endometriosis. Therefore, taking into account the status of both migraine and endometriosis may provide a novel opportunity to identify the genes underlying them. Finally, we propose that the analysis of such genetically correlated comorbid traits can increase power to detect genetic risk loci through the use of more specific, homogenous and heritable phenotypes.