Weighted Interaction SNP Hub (WISH) network method for building genetic networks for complex diseases and traits using whole genome genotype data

Background

High-throughput genotype (HTG) data has been used primarily in genome-wide association (GWA) studies; however, GWA results explain only a limited part of the complete genetic variation of traits. In systems genetics, network approaches have been shown to be able to identify pathways and their underlying causal genes to unravel the biological and genetic background of complex diseases and traits, e.g., the Weighted Gene Co-expression Network Analysis (WGCNA) method based on microarray gene expression data. The main objective of this study was to develop a scale-free weighted genetic interaction network method using whole genome HTG data in order to detect biologically relevant pathways and potential genetic biomarkers for complex diseases and traits.

Results

We developed the Weighted Interaction SNP Hub (WISH) network method that uses HTG data to detect genome-wide interactions between single nucleotide polymorphism (SNPs) and its relationship with complex traits. Data dimensionality reduction was achieved by selecting SNPs based on its: 1) degree of genome-wide significance and 2) degree of genetic variation in a population. Network construction was based on pairwise Pearson's correlation between SNP genotypes or the epistatic interaction effect between SNP pairs. To identify modules the Topological Overlap Measure (TOM) was calculated, reflecting the degree of overlap in shared neighbours between SNP pairs. Modules, clusters of highly interconnected SNPs, were defined using a tree-cutting algorithm on the SNP dendrogram created from the dissimilarity TOM (1-TOM). Modules were selected for functional annotation based on their association with the trait of interest, defined by the Genome-wide Module Association Test (GMAT). We successfully tested the established WISH network method using simulated and real SNP interaction data and GWA study results for carcass weight in a pig resource population; this resulted in detecting modules and key functional and biological pathways related to carcass weight.

Conclusions

We developed the WISH network method which is a novel 'systems genetics' approach to study genetic networks underlying complex trait variation. The WISH network method reduces data dimensionality and statistical complexity in associating genotypes with phenotypes in GWA studies and enables researchers to identify biologically relevant pathways and potential genetic biomarkers for any complex trait of interest.

     

Intradiscal application of rhBMP-7 does not induce regeneration in a canine model of spontaneous intervertebral disc degeneration

IntroductionStrategies for biological repair and regeneration of the intervertebral disc (IVD) by cell and tissue engineering are promising, but few have made it into a clinical setting. Recombinant human bone morphogenetic protein 7 (rhBMP-7) has been shown to stimulate matrix production by IVD cells in vitro and in vivo in animal models of induced IVD degeneration. The aim of this study was to determine the most effective dose of an intradiscal injection of rhBMP-7 in a spontaneous canine IVD degeneration model for translation into clinical application for patients with low back pain.MethodsCanine nucleus pulposus cells (NPCs) were cultured with rhBMP-7 to assess the anabolic effect of rhBMP-7 in vitro, and samples were evaluated for glycosaminoglycan (GAG) and DNA content, histology, and matrix-related gene expression. Three different dosages of rhBMP-7 (2.5 μg, 25 μg, and 250 μg) were injected in vivo into early degenerated IVDs of canines, which were followed up for six months by magnetic resonance imaging (T2-weighted images, T1rho and T2 maps). Post-mortem, the effects of rhBMP-7 were determined by radiography, computed tomography, and macroscopy, and by histological, biochemical (GAG, DNA, and collagen), and biomolecular analyses of IVD tissue.ResultsIn vitro, rhBMP-7 stimulated matrix production of canine NPCs as GAG deposition was enhanced, DNA content was maintained, and gene expression levels of ACAN and COL2A1 were significantly upregulated. Despite the wide dose range of rhBMP-7 (2.5 to 250 μg) administered in vivo, no regenerative effects were observed at the IVD level. Instead, extensive extradiscal bone formation was noticed after intradiscal injection of 25 μg and 250 μg of rhBMP-7.ConclusionsAn intradiscal bolus injection of 2.5 μg, 25 μg, and 250 μg rhBMP-7 showed no regenerative effects in a spontaneous canine IVD degeneration model. In contrast, intradiscal injection of 250 μg rhBMP-7, and to a lesser extent 25 μg rhBMP-7, resulted in extensive extradiscal bone formation, indicating that a bolus injection of rhBMP-7 alone cannot be used for treatment of IVD degeneration in human or canine patients.

Electronic supplementary material

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

Hypercholesterolemia boosts joint destruction in chronic arthritis. An experimental model aggravated by foam macrophage infiltration

Objective

The aim of this study was to determine whether hypercholesterolemia increases articular damage in a rabbit model of chronic arthritis.

Methods

Hypercholesterolemia was induced in 18 rabbits by administrating a high-fat diet (HFD). Fifteen rabbits were fed normal chow as controls. Chronic antigen-induced arthritis (AIA) was induced in half of the HFD and control rabbits, previously immunized, by intra-articular injections of ovalbumin. After sacrifice, lipid and systemic inflammation markers were analyzed in blood serum. Synovium was analyzed by Krenn score, multinucleated cell counting, immunohistochemistry of RAM11 and CD31, and TNF-α and macrophage chemoattractant protein-1 (MCP-1) gene expression. Active bone resorption was assessed by protein expression of receptor activator of nuclear factor kappa-B ligand (RANKL), osteoprotegerin (OPG) and quantification of cathepsin K, contact surface and the invasive area of pannus into bone.

Results

Rabbits receiving the HFD showed higher total serum cholesterol, HDL, triglycerides and CRP levels than rabbits fed a normal diet. Synovitis score was increased in HFD, and particularly in AIA and AIA + HFD groups. AIA + HFD synovium was characterized by a massive infiltration of RAM11+ cells, higher presence of multinucleated foam cells and bigger vascularization than AIA. Cathepsin K+ osteoclasts and the contact surface of bone resorbing pannus were also increased in rabbits with AIA + HFD compared with AIA alone. Synovial TNF-α and MCP-1 gene expression was increased in AIA and HFD rabbits compared with healthy animals. RANKL protein expression in AIA and AIA + HFD groups was higher compared with either HFD or normal groups.

Conclusions

This experimental model demonstrates that hypercholesterolemia increments joint tissue damage in chronic arthritis, with foam macrophages being key players in this process.     

(Cost)-effectiveness of a multi-component intervention for adults with epilepsy: study protocol of a Dutch randomized controlled trial (ZMILE study)

Background

In patients with epilepsy, poor adherence to anti-epileptic drugs has been shown to be the most important cause of poorly controlled epilepsy. Furthermore, it has been noted that the quality of life among patients with epilepsy can be improved by counseling and treatments aimed at increasing their self-efficacy and concordance, thus stimulating self-management skills. However, there is a need for evidence on the effectiveness of such programs, especially within epilepsy care. Therefore, we have developed a multi-component intervention (MCI) which combines a self-management/education program with e-Health interventions. Accordingly, the overall objective of this study is to assess the (cost)-effectiveness and feasibility of the MCI, aiming to improve self-efficacy and concordance in patients with epilepsy.

Methods

A RCT in two parallel groups will be conducted to compare the MCI with a control condition in epilepsy patients. One hundred eligible epilepsy patients will be recruited and allocated to either the intervention or control group. The intervention group will receive the MCI consisting of a self-management/education program of six meetings, including e-Health interventions, and will be followed for 12 months. The control group will receive care as usual and will be followed for 6 months, after which patients will be offered the possibility of participating in the MCI. The study will consist of three parts: 1) a clinical effectiveness study, 2) a cost-effectiveness study, and 3) process evaluation. The primary outcome will be self-efficacy. Secondary outcomes include adherence, side effects, change in seizure severity & frequency, improved quality of life, proactive coping, and societal costs. Outcome assessments will be done using questionnaires at baseline and after 3, 6, 9, and 12 months (last two applicable only for intervention group).

Discussion

In times of budget constraints, MCI could be a valuable addition to the current healthcare provision for epilepsy, as it is expected that higher concordance and self-efficacy will result in reduced use of healthcare resources and an increased QOL. Accordingly, this study is aimed helping patients to be their own provider of health care, shifting epilepsy management from professionals to self-care by patients equipped with appropriate skills and tools.

Trial registration number

NTR4484.

     

Selective use of abciximab in coronary stenting: overall outcomes can still be equivalent to those in the EPISTENT treatment group

BACKGROUND: The EPISTENT trial reported improved early outcomes with routine use of abciximab after coronary stenting. Increasing use of stents means that routine abciximab adds significantly to costs of percutaneous coronary intervention (PCI). This paper reports the results of a protocol encouraging restriction of abciximab therapy to high-risk patients only. METHODS: Data were collected prospectively over a 34-month period for patients undergoing PCI with stenting. In addition to those who fulfilled criteria for inclusion in the EPISTENT trial, patients treated in the setting of acute myocardial infarction (AMI) were studied. Demographic data, procedural details and early clinical outcomes were recorded. RESULTS: Of 808 patients studied, 601 fulfilled EPISTENT inclusion criteria and comprised 367 patients (45%) treated for stable angina and 234 (30%) treated for unstable or post-infarct angina. The additional 207 patients (25%) were treated during AMI. The 808 patients received a total of 981 stents. Abciximab was given in only 88 cases (10.9%). Major adverse clinical events occurred in 39 patients (4.8%). CONCLUSION: Selective use of abciximab for patients undergoing coronary stenting can be associated with outcomes equivalent to those reported for routine use, but with significant cost savings.     

A model of oxidative phosphorylation in mammalian skeletal muscle

A dynamic computer model of oxidative phosphorylation in oxidative mammalian skeletal muscle was developed. The previously published model of oxidative phosphorylation in isolated skeletal muscle mitochondria was extended by incorporation of the creatine kinase system (creatine kinase plus phosphocreatine/creatine pair), cytosolic proton production/consumption system (proton production/consumption by the creatine kinase-catalysed reaction, efflux/influx of protons), physiological size of the adenine nucleotide pool and some additional minor changes. Theoretical studies performed by means of the extended model demonstrated that the CK system, which allows for large changes in P(i) in relation to isolated mitochondria system, has no significant influence on the kinetic properties of oxidative phosphorylation, as inorganic phosphate only slightly modifies the relationship between the respiration rate and [ADP]. Computer simulations also suggested that the second-order dependence of oxidative phosphorylation on [ADP] proposed in the literature refers only to the ATP synthesis flux, but not to the oxygen consumption flux (the difference between these two fluxes being due to the proton leak). Next, time courses of changes in fluxes and metabolite concentrations during transition between different steady-states were simulated. The model suggests, in accordance with previous theoretical predictions, that activation of oxidative phosphorylation by an increase in [ADP] can (roughly) explain the behaviour of the system only at low work intensities, while at higher work intensities parallel activation of different steps of oxidative phosphorylation is involved.     

Effect of different cycling frequencies during incremental exercise on the venous plasma potassium concentration in humans

The effect of different muscle shortening velocity was studied during cycling at a pedalling rate of 60 and 120 rev.min(-1) on the [K+]v in humans. Twenty-one healthy young men aged 22.5+/-2.2 years, body mass 72.7+/-6.4 kg, VO2 max 3.720+/-0.426 l. min(-1), performed an incremental exercise test until exhaustion. The power output increased by 30 W every 3 min, using an electrically controlled ergometer Ergoline 800 S (see Zoladz et al. J. Physiol. 488: 211-217, 1995). The test was performed twice: once at a cycling frequency of 60 rev.min(-1) (test A) and a few days later at a frequency of 120 rev. min(-1) (test B). At rest and at the end of each step (i.e. the last 15 s) antecubital venous blood samples for [K+]p were taken. Gas exchange variables were measured continuously (breath-by-breath) using Oxycon Champion Jaeger. The pre-exercise [K+]v in both tests was not significantly different amounting to 4.24+/-0.36 mmol.l(-1) in test A, and 4.37+/-0.45 mmol.l(-1) in test B. However, the [K+]p during cycling at 120 rev. min(-1) was significantly higher (p<0.001, ANOVA for repeated measurements) at each power output when compared to cycling at 60 rev.min(-1). The maximal power output reached 293+/-31 W in test A which was significantly higher (p<0.001) than in test B, which amounted to 223+/-40 W. The VO2max values in both tests reached 3.720+/-0.426 l. min(-1) vs 3.777+/-0.514 l. min(-1). These values were not significantly different. When the [K+]v was measured during incremental cycling exercise, a linear increase in [K+]v was observed in both tests. However, a significant (p<0.05) upward shift in the [K+]v and a % VO2max relationship was detected during cycling at 120 rev.min(-1). The [K+]v measured at the VO2max level in tests A and B amounted to 6.00+/-0.47 mmol.l-1 vs 6.04+/-0.41 mmol.l-1, respectively. This difference was not significant. It may thus be concluded that: a) generation of the same external mechanical power output during cycling at a pedalling rate of 120 rev.min(-1) causes significantly higher [K+]v changes than when cycling at 60 rev.min(-1), b) the increase of venous plasma potassium concentration during dynamic incremental exercise is linearly related to the metabolic cost of work expressed by the percentage of VO2max (increase as reported previously by Vollestad et al. J. Physiol. 475: 359-368, 1994), c) there is a tendency towards upward up shift in the [K+]v and % VO2max relation during cycling at 120 rev.min(-1) when compared to cycling at 60 rev.min(-1).     

Effect of different muscle shortening velocities during prolonged incremental cycling exercise on the plasma growth hormone, insulin, glucose, glucagon, cortisol, leptin and lactate concentrations.

BACKGROUND: Strenuous exercise was reported to involve the alteration in the release of some "stress" hormones such as growth hormone (GH), cortisol, catecholamines and appropriate adjustment of energy metabolism but the relative contribution of these hormones to metabolic response, to cycling exercise performed at different muscle shortening velocities, has not been clarified. AIMS: The purpose of this experiment was to assess the effect of applying different pedalling rates during a prolonged incremental cycling exercise test on the changes in the plasma levels of growth hormone, cortisol, insulin, glucagon and leptin in humans. Material and METHODS: Fifteen healthy non-smoking men (means +/- SD: age 22.9 +/- 2.4 years; body mass 71.9 +/- 8.2 kg; height 178 +/- 6 cm; with VO2max of 3.896 +/- 0.544 1 x min(-1), assessed in laboratory tests, were subjects in this study. The subjects performed in two different days a prolonged incremental exercise tests at two different pedalling rates, one of them at 60 and another at 120 rev x min(-1). During this tests the power output has increased by 30 W every 6 minutes. The tests were stopped when the subject reached about 70 % of the VO2max. RESULTS AND CONCLUSIONS: We have found that choosing slow or fast pedalling rates (60 or 120 rev min(-1)), while generating the same external mechanical power output, had no effect on the pattern of changes in plasma cortisol, insulin, glucagon, glucose and leptin concentrations. But, generation of the same external mechanical power output at 120 rev min(-1) causes more stepper increase (p < 0.01) in the plasma growth hormone concentration [GH]pl and plasma lactate concentrations [La]pl when compared to that observed during cycling at 60 rev x min(-1). We have also found that the onset of a significant increase in [GH]pl during cycling at 60 rev x min(-1) was not accompanied by significant increase in [La]pl. While during cycling at 120 rev x min(-1) the onset of a significant increase in [La]pl occurred without increase in [GH]pl, but with continuation of exercise when plasma [La]pl increased, there was also a parallel rise in plasma [GH]pl, as reported before. This results indicates that the increase in [GH]pl during exercise is not closely related to the increase in [La]pl.     

A rate-independent technique for analysis of nucleic acid sequences: evolutionary parsimony

The method of evolutionary parsimony--or operator invariants--is a technique of nucleic acid sequence analysis related to parsimony analysis and explicitly designed for determining evolutionary relationships among four distantly related taxa. The method is independent of substitution rates because it is derived from consideration of the group properties of substitution operators rather than from an analysis of the probabilities of substitution in branches of a tree. In both parsimony and evolutionary parsimony, three patterns of nucleotide substitution are associated one-to-one with the three topologically linked trees for four taxa. In evolutionary parsimony, the three quantities are operator invariants. These invariants are the remnants of substitutions that have occurred in the interior branch of the tree and are analogous to the substitutions assigned to the central branch by parsimony. The two invariants associated with the incorrect trees must equal zero (statistically), whereas only the correct tree can have a nonzero invariant. The chi 2-test is used to ascertain the nonzero invariant and the statistically favored tree. Examples, obtained using data calculated with evolutionary rates and branchings designed to camouflage the true tree, show that the method accurately predicts the tree, even when substitution rates differ greatly in neighboring peripheral branches (conditions under which parsimony will consistently fail). As the number of substitutions in peripheral branches becomes fewer, the parsimony and the evolutionary-parsimony solutions converge. The method is robust and easy to use.