Exploring causal networks of bovine milk fatty acids in a multivariate mixed model context

Background

Knowledge regarding causal relationships among traits is important to understand complex biological systems. Structural equation models (SEM) can be used to quantify the causal relations between traits, which allow prediction of outcomes to interventions applied to such a network. Such models are fitted conditionally on a causal structure among traits, represented by a directed acyclic graph and an Inductive Causation (IC) algorithm can be used to search for causal structures. The aim of this study was to explore the space of causal structures involving bovine milk fatty acids and to select a network supported by data as the structure of a SEM.

Results

The IC algorithm adapted to mixed models settings was applied to study 14 correlated bovine milk fatty acids, resulting in an undirected network. The undirected pathway from C4:0 to C12:0 resembled the de novo synthesis pathway of short and medium chain saturated fatty acids. By using prior knowledge, directions were assigned to that part of the network and the resulting structure was used to fit a SEM that led to structural coefficients ranging from 0.85 to 1.05. The deviance information criterion indicated that the SEM was more plausible than the multi-trait model.

Conclusions

The IC algorithm output pointed towards causal relations between the studied traits. This changed the focus from marginal associations between traits to direct relationships, thus towards relationships that may result in changes when external interventions are applied. The causal structure can give more insight into underlying mechanisms and the SEM can predict conditional changes due to such interventions.     

Targeted resequencing of GWAS loci reveals novel genetic variants for milk production traits

Background

Genome wide association study (GWAS) has been proven to be a powerful tool for detecting genomic variants associated with complex traits. However, the specific genes and causal variants underlying these traits remain unclear.

Results

Here, we used target-enrichment strategy coupled with next generation sequencing technique to study target regions which were found to be associated with milk production traits in dairy cattle in our previous GWAS. Among the large amount of novel variants detected by targeted resequencing, we selected 200 SNPs for further association study in a population consisting of 2634 cows. Sixty six SNPs distributed in 53 genes were identified to be associated significantly with on milk production traits. Of the 53 genes, 26 were consistent with our previous GWAS results. We further chose 20 significant genes to analyze their mRNA expression in different tissues of lactating cows, of which 15 were specificly highly expressed in mammary gland.

Conclusions

Our study illustrates the potential for identifying causal mutations for milk production traits using target-enrichment resequencing and extends the results of GWAS by discovering new and potentially functional mutations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1105) contains supplementary material, which is available to authorized users.     

Presence of Hydrogen Peroxide,a Source of Hydroxyl Radicals,in Acid Electrolyzed Water

Background

Acid electrolyzed water (AEW), which is produced through the electrolysis of dilute sodium chloride (NaCl) or potassium chloride solution, is used as a disinfectant in various fields because of its potent antimicrobial activity. The hydroxyl radical, an oxygen radical species, is often suggested as a putative active ingredient for AEW antimicrobial activity.

Methodology/Principal Findings

The aim of the present study is to detect hydroxyl radicals in AEW. The hydroxyl radicals in AEW prepared under different conditions were determined using an electron spin resonance (ESR) technique. A signal from 5,5-dimethyl-1-pyrroline N-oxide (DMPO)-OH, an adduct of DMPO and the hydroxyl radical, was detected in AEW prepared by double or triple electrolyses of 1% NaCl but not of 0.1% NaCl solution. Then the presence of hydrogen peroxide as a proposed source of hydroxyl radicals was examined using a combination of ESR and a Fenton reaction. The DMPO-OH signal was clearly detected, even in AEW prepared by single electrolysis of 0.1% NaCl solution, when ferrous sulfate was added to induce a Fenton reaction, indicating the presence of hydrogen peroxide in the AEW. Since sodium formate, a hydroxyl radical scavenger, did not affect the bactericidal activity of AEW, it is concluded that the radical is unlikely to contribute to the antimicrobial activity of AEW, although a small amount of the radical is produced from hydrogen peroxide. Dimethyl sulfoxide, the other hydroxyl radical scavenger used in the present study, canceled the bactericidal activity of AEW, accompanied by complete depletion of free available chlorine, suggesting that hypochlorous acid is probably a major contributor to the antimicrobial activity.

Conclusions

It is strongly suggested that although hydrogen peroxide is present in AEW as a source of hydroxyl radicals, the antimicrobial activity of AEW does not depend on these radicals.     

Partitioning of genomic variance reveals biological pathways associated with udder health and milk production traits in dairy cattle

Background

We have used a linear mixed model (LMM) approach to examine the joint contribution of genetic markers associated with a biological pathway. However, with these markers being scattered throughout the genome, we are faced with the challenge of modelling the contribution from several, sometimes even all, chromosomes at once. Due to linkage disequilibrium (LD), all markers may be assumed to account for some genomic variance; but the question is whether random sets of markers account for the same genomic variance as markers associated with a biological pathway?

Results

We applied the LMM approach to identify biological pathways associated with udder health and milk production traits in dairy cattle. A random gene sampling procedure was applied to assess the biological pathways in a dataset that has an inherently complex genetic correlation pattern due to the population structure of dairy cattle, and to linkage disequilibrium within the bovine genome and within the genes associated to the biological pathway.

Conclusions

Several biological pathways that were significantly associated with health and production traits were identified in dairy cattle; i.e. the markers linked to these pathways explained more of the genomic variance and provided a better model fit than 95 % of the randomly sampled gene groups. Our results show that immune related pathways are associated with production traits, and that pathways that include a causal marker for production traits are identified with our procedure.We are confident that the LMM approach provides a general framework to exploit and integrate prior biological information and could potentially lead to improved understanding of the genetic architecture of complex traits and diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0132-6) contains supplementary material, which is available to authorized users.     

Inhibition of allergic airway responses by heparin derived oligosaccharides: identification of a tetrasaccharide sequence

Background

Previous studies showed that heparin''s anti-allergic activity is molecular weight dependent and resides in oligosaccharide fractions of <2500 daltons.

Objective

To investigate the structural sequence of heparin''s anti-allergic domain, we used nitrous acid depolymerization of porcine heparin to prepare an oligosaccharide, and then fractionated it into disaccharide, tetrasaccharide, hexasaccharide, and octasaccharide fractions. The anti-allergic activity of each oligosaccharide fraction was tested in allergic sheep.

Methods

Allergic sheep without (acute responder) and with late airway responses (LAR; dual responder) were challenged with Ascaris suum antigen with and without inhaled oligosaccharide pretreatment and the effects on specific lung resistance and airway hyperresponsiveness (AHR) to carbachol determined. Additional inflammatory cell recruitment studies were performed in immunized ovalbumin-challenged BALB/C mice with and without treatment.

Results

The inhaled tetrasaccharide fraction was the minimal effective chain length to show anti-allergic activity. This fraction showed activity in both groups of sheep; it was also effective in inhibiting LAR and AHR, when administered after the antigen challenge. Tetrasaccharide failed to modify the bronchoconstrictor responses to airway smooth muscle agonists (histamine, carbachol and LTD₄), and had no effect on antigen-induced histamine release in bronchoalveolar lavage fluid in sheep. In mice, inhaled tetrasaccharide also attenuated the ovalbumin-induced peribronchial inflammatory response and eosinophil influx in the bronchoalveolar lavage fluid. Chemical analysis identified the active structure to be a pentasulfated tetrasaccharide ([IdoU2S (1→4)GlcNS6S (1→4) IdoU2S (1→4) AMan-6S]) which lacked anti-coagulant activity.

Conclusions

These results demonstrate that heparin tetrasaccharide possesses potent anti-allergic and anti-inflammatory properties, and that the domains responsible for anti-allergic and anti-coagulant activity are distinctly different.     

Causal relationships between milk quality and coagulation properties in Italian Holstein-Friesian dairy cattle

Background

Recently, selection for milk technological traits was initiated in the Italian dairy cattle industry based on direct measures of milk coagulation properties (MCP) such as rennet coagulation time (RCT) and curd firmness 30 min after rennet addition (a₃₀) and on some traditional milk quality traits that are used as predictors, such as somatic cell score (SCS) and casein percentage (CAS). The aim of this study was to shed light on the causal relationships between traditional milk quality traits and MCP. Different structural equation models that included causal effects of SCS and CAS on RCT and a₃₀ and of RCT on a₃₀ were implemented in a Bayesian framework.

Results

Our results indicate a non-zero magnitude of the causal relationships between the traits studied. Causal effects of SCS and CAS on RCT and a₃₀ were observed, which suggests that the relationship between milk coagulation ability and traditional milk quality traits depends more on phenotypic causal pathways than directly on common genetic influence. While RCT does not seem to be largely controlled by SCS and CAS, some of the variation in a₃₀ depends on the phenotypes of these traits. However, a₃₀ depends heavily on coagulation time. Our results also indicate that, when direct effects of SCS, CAS and RCT are considered simultaneously, most of the overall genetic variability of a₃₀ is mediated by other traits.

Conclusions

This study suggests that selection for RCT and a₃₀ should not be performed on correlated traits such as SCS or CAS but on direct measures because the ability of milk to coagulate is improved through the causal effect that the former play on the latter, rather than from a common source of genetic variation. Breaking the causal link (e.g. standardizing SCS or CAS before the milk is processed into cheese) would reduce the impact of the improvement due to selective breeding. Since a₃₀ depends heavily on RCT, the relative emphasis that is put on this trait should be reconsidered and weighted for the fact that the pure measure of a₃₀ almost double-counts RCT.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0123-7) contains supplementary material, which is available to authorized users.     

Association between floral traits and rewards in Erysimum mediohispanicum (Brassicaceae)

Background and Aims

Floral rewards may be associated with certain morphological floral traits and thus act as underlying factors promoting selection on these traits. This study investigates whether some traits that are under pollinator-mediated selection (flower number, stalk height, corolla diameter, corolla tube length and corolla tube width) in the Mediterranean herb E. mediohispanicum (Brassicaceae) are associated with rewards (pollen and nectar).

Methods

During 2005 the phenotypic traits and the visitation rate of the main pollinator functional groups were quantified in 720 plants belonging to eight populations in south-east Spain, and during 2006 the same phenotypic traits and the reward production were quantified in 400 additional plants from the same populations.

Key Results

A significant correlation was found between nectar production rate and corolla tube length, and between pollen production and corolla diameter. Visitation rates of large bees and butterflies were significantly higher in plants exhibiting larger flowers with longer corolla tubes.

Conclusions

The association between reward production and floral traits may be a factor underlying the pattern of visitation rate displayed by some pollinators.Key words: Erysimum, floral traits, nectar, pollen, pollinator visitation rate, reward     

DIC Score in Pregnant Women – A Population Based Modification of the International Society on Thrombosis and Hemostasis Score

Objectives

The objectives of this study were: 1) To determine the component needed to generate a validated DIC score during pregnancy. 2) To validate such scoring system in the identification of patients with clinical diagnosis of DIC.

Material and Methods

This is a population based retrospective study, including all women who gave birth at the ‘Soroka University Medical Center’ during the study period, and have had blood coagulation tests including complete blood cell count, prothrombin time (PT)(seconds), partial thromboplastin time (aPTT), fibrinogen, and D-dimers. Nomograms for pregnancy were established, and DIC score was constructed based on ROC curve analyses.

Results

1) maternal plasma fibrinogen concentrations increased during pregnancy; 2) maternal platelet count decreased gradually during gestation; 3) the PT and PTT values did not change with advancing gestation; 4) PT difference had an area under the curve (AUC) of 0.96 (p<0.001), and a PT difference ≥1.55 had an 87% sensitivity and 90% specificity for the diagnosis of DIC; 5) the platelet count had an AUC of 0.87 (p<0.001), an 86% sensitivity and 71% specificity for the diagnosis of DIC; 6) fibrinogen concentrations had an AUC of 0.95 (p<0.001) and a cutoff point ≤3.9 g/L had a sensitivity of 87% and a specificity of 92% for the development of DIC; and 7) The pregnancy adjusted DIC score had an AUC of 0.975 (p<0.001) and at a cutoff point of ≥26 had a sensitivity of 88%, a specificity of 96%, a LR(+) of 22 and a LR(−) of 0.125 for the diagnosis of DIC.

Conclusion

We could establish a sensitive and specific pregnancy adjusted DIC score. The positive likelihood ratio of this score suggests that a patient with a score of ≥26 has a high probability to have DIC.     

Xeno-Free Defined Conditions for Culture of Human Embryonic Stem Cells,Neural Stem Cells and Dopaminergic Neurons Derived from Them

Background

Human embryonic stem cells (hESCs) may provide an invaluable resource for regenerative medicine. To move hESCs towards the clinic it is important that cells with therapeutic potential be reproducibly generated under completely defined conditions.

Methodology/Principal Findings

Here we report a four-step scalable process that is readily transferable to a Good Manufacture Practice (GMP) facility for the production of functional dopaminergic neurons from hESCs for potential clinical uses. We show that each of the steps (propagation of ESC→generation of neural stem cells (NSC)→induction of dopaminergic precursors→maturation of dopaminergic neurons) could utilize xeno-free defined media and substrate, and that cells could be stored at intermediate stages in the process without losing their functional ability. Neurons generated by this process expressed midbrain and A9 dopaminergic markers and could be transplanted at an appropriate time point in development to survive after transplant.

Conclusions/Significance

hESCs and NSCs can be maintained in xeno-free defined media for a prolonged period of time while retaining their ability to differentiate into authentic dopaminergic neurons. Our defined medium system provides a path to a scalable GMP-applicable process of generation of dopaminergic neurons from hESCs for therapeutic applications, and a ready source of large numbers of neurons for potential screening applications.     

Inherited variability in multiple traits determines fitness in populations of an annual legume from contrasting latitudinal origins

Background and Aims

Variation in fitness depends on corresponding variation in multiple traits which have both genetically controlled and plastic components. These traits are subjected to varying degrees of local adaptation in specific populations and, consequently, are genetically controlled to different extents. In this study it is hypothesized that modulation of different traits would have contrasting relevance for the fitness of populations of diverse origins. Specifically, assuming that environmental pressures vary across a latitudinal gradient, it is suggested that inherited variation in traits differentially determines fitness in annual Lupinus angustifolius populations from contrasting latitudinal origins in western Spain.

Methods

Seeds of L. angustifolius from three contrasting origins were grown in a common garden. Traits related to more plastic vegetative growth and more genetically conserved phenology were measured, together with estimates of reproductive success. Fitness was estimated by the number of viable seeds per plant. Structural Equation Models were used to infer causal relationships among multiple traits and fitness, separating the direct and indirect effects of morphological, phenological and reproductive traits.

Key Results

Phenological, vegetative and reproductive traits accounted for most of the fitness variation. Fitness was highest in plants of southernmost origin, mainly due to earlier flowering. Fitness within each seed origin was controlled by variation in different traits. Southern origin plants that grew to a larger size achieved higher fitness. However, plant size in plants of northernmost origin was irrelevant, but early flowering promoted higher fitness. Variation in fruit and seed set had a greater effect on the fitness of plants of central origin than phenological and size variation.

Conclusions

It is concluded that modulation of a functional trait can be relevant to fitness in a given population (i.e. affecting intensity and direction), but irrelevant in other populations. This points to the need to consider integrated phenotypes when trying to unravel local adaptation effects over single traits.Key words: Lupinus, Structural Equation Models, fitness, phenology, functional traits, reproductive success, SLA, seed size     

Tyrosine sulfation of native mouse Psgl-1 is required for optimal leukocyte rolling on P-selectin in vivo

Background

We recently demonstrated that tyrosine sulfation is an important contributor to monocyte recruitment and retention in a mouse model of atherosclerosis. P-selectin glycoprotein ligand-1 (Psgl-1) is tyrosine-sulfated in mouse monocyte/macrophages and its interaction with P-selectin is important in monocyte recruitment in atherosclerosis. However, whether tyrosine sulfation is required for the P-selectin binding function of mouse Psgl-1 is unknown. Here we test the function of native Psgl-1 expressed in leukocytes lacking endogenous tyrosylprotein sulfotransferase (TPST) activity.

Methodology/Principal Findings

Psgl-1 function was assessed by examining P-selectin dependent leukocyte rolling in post-capillary venules of C57BL6 mice transplanted with hematopoietic progenitors from wild type (WT→B6) or Tpst1;Tpst2 double knockout mice (Tpst DKO→B6) which lack TPST activity. We observed that rolling flux fractions were lower and leukocyte rolling velocities were higher in Tpst DKO→B6 venules compared to WT→B6 venules. Similar results were observed on immobilized P-selectin in vitro. Finally, Tpst DKO leukocytes bound less P-selectin than wild type leukocytes despite equivalent surface expression of Psgl-1.

Conclusions/Significance

These findings provide direct and convincing evidence that tyrosine sulfation is required for optimal function of mouse Psgl-1 in vivo and suggests that tyrosine sulfation of Psgl-1 contributes to the development of atherosclerosis.     

Quantitative trait loci mapping for canine hip dysplasia and its related traits in UK Labrador Retrievers

Background

Canine hip dysplasia (CHD) is characterised by a malformation of the hip joint, leading to osteoarthritis and lameness. Current breeding schemes against CHD have resulted in measurable but moderate responses. The application of marker-assisted selection, incorporating specific markers associated with the disease, or genomic selection, incorporating genome-wide markers, has the potential to dramatically improve results of breeding schemes. Our aims were to identify regions associated with hip dysplasia or its related traits using genome and chromosome-wide analysis, study the linkage disequilibrium (LD) in these regions and provide plausible gene candidates. This study is focused on the UK Labrador Retriever population, which has a high prevalence of the disease and participates in a recording program led by the British Veterinary Association (BVA) and The Kennel Club (KC).

Results

Two genome-wide and several chromosome-wide QTLs affecting CHD and its related traits were identified, indicating regions related to hip dysplasia.

Conclusion

Consistent with previous studies, the genetic architecture of CHD appears to be based on many genes with small or moderate effect, suggesting that genomic selection rather than marker-assisted selection may be an appropriate strategy for reducing this disease.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-833) contains supplementary material, which is available to authorized users.     

Associations of Children’s Appetitive Traits with Weight and Dietary Behaviours in the Context of General Parenting

Background

Individual variations in child weight can be explained by genetic and behavioural susceptibility to obesity. Behavioural susceptibility can be expressed in appetite-related traits, e.g. food responsiveness. Research into such behavioural factors is important, as it can provide starting points for (preventive) interventions.

Objectives

To examine associations of children’s appetitive traits with weight and with fruit, snack and sugar-sweetened beverage intake, and to examine whether parenting style interacts with appetite in determining child weight/intake.

Methods

Data were used from 1275 children participating in the INPACT study in 2009–2010, with a mean age of 9 years in 2009. Their height and weight were measured to calculate body mass index (BMI). Parents completed a questionnaire to measure children’s appetitive traits, children’s dietary intake and parenting style. Child BMI z-scores, fruit, snack and sugar-sweetened beverage intake were regressed on appetitive traits. Moderation by parenting style was tested by adding interaction terms to the regression analyses.

Results

Food-approaching appetitive traits were positively, and food-avoidant appetitive traits were negatively related to child BMI z-scores and to child fruit intake. There were no or less consistent associations for snack and sugar-sweetened beverage intake. Authoritative parenting voided the negative association between food fussiness and fruit intake, while neglecting parenting strengthened the positive association between food-approaching appetitive traits and weight.

Conclusions

Early assessment of appetitive traits could be used to identify children at risk for overweight. As parenting style can moderate the associations between appetitive traits and weight/intake in a favourable way, parents are a promising target group for preventive interventions aimed at influencing the effect of appetitive traits on children.     

Immune Sculpting of Norepinephrine on MHC-I,B7-1, IDO and B7-H1 Expression and Regulation of Proliferation and Invasion in Pancreatic Carcinoma Cells

Background

The sympathetic neurotransmitter Norepinephrine (NE) contributes to tumorigenesis and cancer progression. This study aims to investigate the role of NE in modulating the immune phenotype and allowing pancreatic carcinoma (PC) cells to escape the immune response.

Methods

Varied concentrations of NE and interferon-gamma (IFN-γ) were administrated to MIA PaCa-2 and BxPC-3 cell lines for 48 hours. Proliferation and invasion were then investigated using an MTT assay and a membrane invasion culture system respectively. MHC-I, B7-1, IDO and B7-H1 expression were measured using real-time quantitative RT-PCR, western blotting and immunocytochemistry. The synergistic and time-dependent effects of NE/IFN-γ were also investigated. Adrenergic antagonists were used to identify the relevant target receptor of NE.

Results

The results showed that NE had dose-dependent and time-dependent effects on cell biological processes as well as on the expression of MHC-I, B7-1, IDO and B7-H1. These effects occurred mainly via the β₂-adrenergic receptor. Long-term NE treatment was able to antagonize some of the effects of IFN-γ (after 2 weeks of treatment), but NE and IFN-γ had significant synergistic stimulatory effects on IDO and B7-H1 expression. The residual effects on biological activities lasted for 2 weeks, while the immunophenotypic changes decreased at early time points after treatment.

Conclusions

NE plays important roles in modulating PC cell biological activities and affecting MHC-I, B7-1, IDO and B7-H1 expression in vitro, mainly via the β2-adrenergic receptor (β2-AR) in a time- and dose-dependent fashion. Only at extended treatment durations could NE affect PC cell progression and immune evasion.     

Immune Modulation by Different Types of β2→1-Fructans Is Toll-Like Receptor Dependent

Introduction

β2→1-fructans are dietary fibers. Main objectives of this study were 1) to demonstrate direct signalling of β2→1-fructans on immune cells, 2) to study whether this is mediated by the pattern recognition receptors Toll-like receptors (TLRs) and nucleotide-binding oligomerisation domain-containing proteins (NODs), and 3) to relate the observed effects to the chain length differences in β2→1-fructans.

Methods

Four different β2→1-fructan formulations were characterised for their chain length profile. Human peripheral blood mononuclear cells (PBMCs) were stimulated in vitro with β2→1-fructans, and production of IL-1Ra, IL-1β, IL-6, IL-10, IL-12p70, and TNF-α was analysed. Reporter cells for TLRs and NODs were incubated with β2→1-fructans and analysed for NF-κB/AP-1 activation.

Results

Cytokine production in human PBMCs was dose- and chain length-dependent. Strikingly, short chain enriched β2→1-fructans induced a regulatory cytokine balance compared to long chain enriched β2→1-fructans as measured by IL-10/IL-12 ratios. Activation of reporter cells showed that signalling was highly dependent on TLRs and their adapter, myeloid differentiation primary response protein 88 (MyD88). In human embryonic kidney reporter cells, TLR2 was prominently activated, while TLR4, 5, 7, 8, and NOD2 were mildly activated.

Conclusions

β2→1-fructans possess direct signalling capacity on human immune cells. By activating primarily TLR2, and to a lesser extent TLR4, 5, 7, 8, and NOD2, β2→1-fructan stimulation results in NF-κB/AP-1 activation. Chain length of β2→1-fructans is important for the induced activation pattern and IL-10/IL-12 ratios.     

The adaptive value of functional and life-history traits across fertility treatments in an annual plant

Background and Aims

Plant functional traits are assumed to be adaptive. As selection acts on individuals and not on traits, interpreting the adaptive value of a trait not may be straightforward. For example, productive leaves are associated with fertile environments. However, it is not clear if productive leaves confer an advantage in these habitats, or if they are an advantage as part of a suite of coordinated traits.

Methods

Genotypes of Arabidopsis thaliana were grown in high and low nutrient treatments and low, neutral and high pH treatments. Nutrient availability is reduced in acidic or basic soils relative to neutral pH soils. pH treatments were used to alter the availability of resources rather than the amount of resources.

Key Results

Leaf function (specific leaf area, SLA) and life history (size at reproduction, age at reproduction) were variable across genotypes and were plastic. High nutrient availability induced higher SLA and larger size at reproduction. Genotypes that reproduced at large size in high nutrient conditions at neutral pH had the greatest fruit production. SLA was only indirectly related to fruit production through a causal relationship with rosette size; in high nutrient conditions, plants with high SLA were large at reproduction and had higher fruit production. In high nutrient and high pH treatments, plants were large at reproduction, but large size at reproduction was associated with low fecundity. This suggests that large size is adaptive under high nutrient availability.

Conclusions

Interpreting the adaptive value of functional traits will sometimes only be possible when these traits are considered as a suite of correlated and coordinated traits. Leaf functional traits may be important in defining adaptive strategies in A. thaliana but only through how they affect plant life history. Finally, manipulating soil pH can be a valuable tool in assessing adaptive plasticity on nutrient gradients.     

Performance of Polygenic Scores for Predicting Phobic Anxiety

Context

Anxiety disorders are common, with a lifetime prevalence of 20% in the U.S., and are responsible for substantial burdens of disability, missed work days and health care utilization. To date, no causal genetic variants have been identified for anxiety, anxiety disorders, or related traits.

Objective

To investigate whether a phobic anxiety symptom score was associated with 3 alternative polygenic risk scores, derived from external genome-wide association studies of anxiety, an internally estimated agnostic polygenic score, or previously identified candidate genes.

Design

Longitudinal follow-up study. Using linear and logistic regression we investigated whether phobic anxiety was associated with polygenic risk scores derived from internal, leave-one out genome-wide association studies, from 31 candidate genes, and from out-of-sample genome-wide association weights previously shown to predict depression and anxiety in another cohort.

Setting and Participants

Study participants (n = 11,127) were individuals from the Nurses'' Health Study and Health Professionals Follow-up Study.

Main Outcome Measure

Anxiety symptoms were assessed via the 8-item phobic anxiety scale of the Crown Crisp Index at two time points, from which a continuous phenotype score was derived.

Results

We found no genome-wide significant associations with phobic anxiety. Phobic anxiety was also not associated with a polygenic risk score derived from the genome-wide association study beta weights using liberal p-value thresholds; with a previously published genome-wide polygenic score; or with a candidate gene risk score based on 31 genes previously hypothesized to predict anxiety.

Conclusion

There is a substantial gap between twin-study heritability estimates of anxiety disorders ranging between 20–40% and heritability explained by genome-wide association results. New approaches such as improved genome imputations, application of gene expression and biological pathways information, and incorporating social or environmental modifiers of genetic risks may be necessary to identify significant genetic predictors of anxiety.     

Revisiting spatial distribution and biochemical composition of calcium-containing crystals in human osteoarthritic articular cartilage

Introduction

Calcium-containing (CaC) crystals, including basic calcium phosphate (BCP) and calcium pyrophosphate dihydrate (CPP), are associated with destructive forms of osteoarthritis (OA). We assessed their distribution and biochemical and morphologic features in human knee OA cartilage.

Methods

We prospectively included 20 patients who underwent total knee replacement (TKR) for primary OA. CaC crystal characterization and identification involved Fourier-transform infra-red spectrometry and scanning electron microscopy of 8 to 10 cartilage zones of each knee, including medial and lateral femoral condyles and tibial plateaux and the intercondyle zone. Differential expression of genes involved in the mineralization process between cartilage with and without calcification was assessed in samples from 8 different patients by RT-PCR. Immunohistochemistry and histology studies were performed in 6 different patients.

Results

Mean (SEM) age and body mass index of patients at the time of TKR was 74.6 (1.7) years and 28.1 (1.6) kg/m², respectively. Preoperative X-rays showed joint calcifications (chondrocalcinosis) in 4 cases only. The medial femoro-tibial compartment was the most severely affected in all cases, and mean (SEM) Kellgren-Lawrence score was 3.8 (0.1). All 20 OA cartilages showed CaC crystals. The mineral content represented 7.7% (8.1%) of the cartilage weight. All patients showed BCP crystals, which were associated with CPP crystals for 8 joints. CaC crystals were present in all knee joint compartments and in a mean of 4.6 (1.7) of the 8 studied areas. Crystal content was similar between superficial and deep layers and between medial and femoral compartments. BCP samples showed spherical structures, typical of biological apatite, and CPP samples showed rod-shaped or cubic structures. The expression of several genes involved in mineralization, including human homolog of progressive ankylosis, plasma-cell-membrane glycoprotein 1 and tissue-nonspecific alkaline phosphatase, was upregulated in OA chondrocytes isolated from CaC crystal-containing cartilages.

Conclusions

CaC crystal deposition is a widespread phenomenon in human OA articular cartilage involving the entire knee cartilage including macroscopically normal and less weight-bearing zones. Cartilage calcification is associated with altered expression of genes involved in the mineralisation process.