What Can Causal Networks Tell Us about Metabolic Pathways?

Graphical models describe the linear correlation structure of data and have been used to establish causal relationships among phenotypes in genetic mapping populations. Data are typically collected at a single point in time. Biological processes on the other hand are often non-linear and display time varying dynamics. The extent to which graphical models can recapitulate the architecture of an underlying biological processes is not well understood. We consider metabolic networks with known stoichiometry to address the fundamental question: “What can causal networks tell us about metabolic pathways?”. Using data from an Arabidopsis BaySha population and simulated data from dynamic models of pathway motifs, we assess our ability to reconstruct metabolic pathways using graphical models. Our results highlight the necessity of non-genetic residual biological variation for reliable inference. Recovery of the ordering within a pathway is possible, but should not be expected. Causal inference is sensitive to subtle patterns in the correlation structure that may be driven by a variety of factors, which may not emphasize the substrate-product relationship. We illustrate the effects of metabolic pathway architecture, epistasis and stochastic variation on correlation structure and graphical model-derived networks. We conclude that graphical models should be interpreted cautiously, especially if the implied causal relationships are to be used in the design of intervention strategies

An integrated workflow for robust alignment and simplified quantitative analysis of NMR spectrometry data

<p>Abstract</p> <p>Background</p> <p>Nuclear magnetic resonance spectroscopy (NMR) is a powerful technique to reveal and compare quantitative metabolic profiles of biological tissues. However, chemical and physical sample variations make the analysis of the data challenging, and typically require the application of a number of preprocessing steps prior to data interpretation. For example, noise reduction, normalization, baseline correction, peak picking, spectrum alignment and statistical analysis are indispensable components in any NMR analysis pipeline.</p> <p>Results</p> <p>We introduce a novel suite of informatics tools for the quantitative analysis of NMR metabolomic profile data. The core of the processing cascade is a novel peak alignment algorithm, called hierarchical Cluster-based Peak Alignment (CluPA). The algorithm aligns a target spectrum to the reference spectrum in a top-down fashion by building a hierarchical cluster tree from peak lists of reference and target spectra and then dividing the spectra into smaller segments based on the most distant clusters of the tree. To reduce the computational time to estimate the spectral misalignment, the method makes use of Fast Fourier Transformation (FFT) cross-correlation. Since the method returns a high-quality alignment, we can propose a simple methodology to study the variability of the NMR spectra. For each aligned NMR data point the ratio of the between-group and within-group sum of squares (BW-ratio) is calculated to quantify the difference in variability between and within predefined groups of NMR spectra. This differential analysis is related to the calculation of the F-statistic or a one-way ANOVA, but without distributional assumptions. Statistical inference based on the BW-ratio is achieved by bootstrapping the null distribution from the experimental data.</p> <p>Conclusions</p> <p>The workflow performance was evaluated using a previously published dataset. Correlation maps, spectral and grey scale plots show clear improvements in comparison to other methods, and the down-to-earth quantitative analysis works well for the CluPA-aligned spectra. The whole workflow is embedded into a modular and statistically sound framework that is implemented as an R package called "speaq" ("spectrum alignment and quantitation"), which is freely available from <url>http://code.google.com/p/speaq/</url>.</p

C2 and CFB Genes in Age-Related Maculopathy and Joint Action with CFH and LOC387715 Genes

Background: Age-related maculopathy (ARM) is a common cause of visual impairment in the elderly populations of industrialized countries and significantly affects the quality of life of those suffering from the disease. Variants within two genes, the complement factor H (CFH) and the poorly characterized LOC387715 (ARMS2), are widely recognized as ARM risk factors. CFH is important in regulation of the alternative complement pathway suggesting this pathway is involved in ARM pathogenesis. Two other complement pathway genes, the closely linked complement component receptor (C2) and complement factor B (CFB), were recently shown to harbor variants associated with ARM. Methods/Principal Findings: We investigated two SNPs in C2 and two in CFB in independent case-control and family cohorts of white subjects and found rs547154, an intronic SNP in C2, to be significantly associated with ARM in both our case-control (P-value 0.00007) and family data (P-value 0.00001). Logistic regression analysis suggested that accounting for the effect at this locus significantly (P-value 0.002) improves the fit of a genetic risk model of CFH and LOC387715 effects only. Modeling with the generalized multifactor dimensionality reduction method showed that adding C2 to the two-factor model of CFH and LOC387715 increases the sensitivity (from 63% to 73%). However, the balanced accuracy increases only from 71% to 72%, and the specificity decreases from 80% to 72%. Conclusions/Significance: C2/CFB significantly influences AMD susceptibility and although accounting for effects at this locus does not dramatically increase the overall accuracy of the genetic risk model, the improvement over the CFH-LOC387715 model is statistically significant. © 2008 Jakobsdottir et al

A missense variant in CST3 exerts a recessive effect on susceptibility to age-related macular degeneration resembling its association with Alzheimer’s disease

Age-related macular degeneration (AMD) and Alzheimer’s disease (AD) are degenerative, multifactorial diseases involving age-related accumulation of extracellular deposits linked to dysregulation of protein homeostasis. Here, we strengthen the evidence that an nsSNP (p.Ala25Thr) in the cysteine proteinase inhibitor cystatin C gene CST3, previously confirmed by meta-analysis to be associated with AD, is associated with exudative AMD. To our knowledge, this is the first report highlighting a genetic variant that increases the risk of developing both AD and AMD. Furthermore, we demonstrate that the risk associated with the mutant allele follows a recessive model for both diseases. We perform an AMD-CST3 case–control study genotyping 350 exudative AMD Caucasian individuals. Bringing together our data with the previously reported AMD-CST3 association study, the evidence of a recessive effect on AMD risk is strengthened (OR = 1.89, P = 0.005). This effect closely resembles the AD-CST3 recessive effect (OR = 1.73, P = 0.005) previously established by meta-analysis. This resemblance is substantiated by the high correlation between CST3 genotype and effect size across the two diseases (R2 = 0.978). A recessive effect is in line with the known function of cystatin C, a potent enzyme inhibitor. Its potency means that, in heterozygous individuals, a single functional allele is sufficient to maintain its inhibitory function; only homozygous individuals will lack this form of proteolytic regulation. Our findings support the hypothesis that recessively acting variants account for some of the missing heritability of multifactorial diseases. Replacement therapy represents a translational opportunity for individuals homozygous for the mutant allele

A Non Membrane-Targeted Human Soluble CD59 Attenuates Choroidal Neovascularization in a Model of Age Related Macular Degeneration

Age related macular degeneration (AMD) is the most common cause of blindness amongst the elderly. Approximately 10% of AMD patients suffer from an advanced form of AMD characterized by choroidal neovascularization (CNV). Recent evidence implicates a significant role for complement in the pathogenesis of AMD. Activation of complement terminates in the incorporation of the membrane attack complex (MAC) in biological membranes and subsequent cell lysis. Elevated levels of MAC have been documented on choroidal blood vessels and retinal pigment epithelium (RPE) of AMD patients. CD59 is a naturally occurring membrane bound inhibitor of MAC formation. Previously we have shown that membrane bound human CD59 delivered to the RPE cells of mice via an adenovirus vector can protect those cells from human complement mediated lysis ex vivo. However, application of those observations to choroidal blood vessels are limited because protection from MAC- mediated lysis was restricted only to the cells originally transduced by the vector. Here we demonstrate that subretinal delivery of an adenovirus vector expressing a transgene for a soluble non-membrane binding form of human CD59 can attenuate the formation of laser-induced choroidal neovascularization and murine MAC formation in mice even when the region of vector delivery is distal to the site of laser induced CNV. Furthermore, this same recombinant transgene delivered to the intravitreal space of mice by an adeno-associated virus vector (AAV) can also attenuate laser-induced CNV. To our knowledge, this is the first demonstration of a non-membrane targeting CD59 having biological potency in any animal model of disease in vivo. We propose that the above approaches warrant further exploration as potential approaches for alleviating complement mediated damage to ocular tissues in AMD

Analysis of copy number variation at DMBT1 and age-related macular degeneration

BACKGROUND: DMBT1 is a gene that shows extensive copy number variation (CNV) that alters the number of bacteria-binding domains in the protein and has been shown to activate the complement pathway. It lies next to the ARMS2/HTRA1 genes in a region of chromosome 10q26, where single nucleotide variants have been strongly associated with age-related macular degeneration (AMD), the commonest cause of blindness in Western populations. Complement activation is thought to be a key factor in the pathogenesis of this condition. We sought to investigate whether DMBT1 CNV plays any role in the susceptibility to AMD. METHODS: We analysed long-range linkage disequilibrium of DMBT1 CNV1 and CNV2 with flanking single nucleotide polymorphisms (SNPs) using our previously published CNV and HapMap Phase 3 SNP data in the CEPH Europeans from Utah (CEU). We then typed a large cohort of 860 AMD patients and 419 examined age-matched controls for copy number at DMBT1 CNV1 and CNV2 and combined these data with copy numbers from a further 480 unexamined controls. RESULTS: We found weak linkage disequilibrium between DMBT1 CNV1 and CNV2 with the SNPs rs1474526 and rs714816 in the HTRA1/ARMS2 region. By directly analysing copy number variation, we found no evidence of association of CNV1 or CNV2 with AMD. CONCLUSIONS: We have shown that copy number variation at DMBT1 does not affect risk of developing age-related macular degeneration and can therefore be ruled out from future studies investigating the association of structural variation at 10q26 with AMD

A Chaperone Trap Contributes to the Onset of Cystic Fibrosis

Protein folding is the primary role of proteostasis network (PN) where chaperone interactions with client proteins determine the success or failure of the folding reaction in the cell. We now address how the Phe508 deletion in the NBD1 domain of the cystic fibrosis (CF) transmembrane conductance regulator (CFTR) protein responsible for cystic fibrosis (CF) impacts the binding of CFTR with cellular chaperones. We applied single ion reaction monitoring mass spectrometry (SRM-MS) to quantitatively characterize the stoichiometry of the heat shock proteins (Hsps) in CFTR folding intermediates in vivo and mapped the sites of interaction of the NBD1 domain of CFTR with Hsp90 in vitro. Unlike folding of WT-CFTR, we now demonstrate the presence of ΔF508-CFTR in a stalled folding intermediate in stoichiometric association with the core Hsps 40, 70 and 90, referred to as a ‘chaperone trap’. Culturing cells at 30 C resulted in correction of ΔF508-CFTR trafficking and function, restoring the sub-stoichiometric association of core Hsps observed for WT-CFTR. These results support the interpretation that ΔF508-CFTR is restricted to a chaperone-bound folding intermediate, a state that may contribute to its loss of trafficking and increased targeting for degradation. We propose that stalled folding intermediates could define a critical proteostasis pathway branch-point(s) responsible for the loss of function in misfolding diseases as observed in CF

A 32 kb Critical Region Excluding Y402H in CFH Mediates Risk for Age-Related Macular Degeneration

Complement factor H shows very strong association with Age-related Macular Degeneration (AMD), and recent data suggest that multiple causal variants are associated with disease. To refine the location of the disease associated variants, we characterized in detail the structural variation at CFH and its paralogs, including two copy number polymorphisms (CNP), CNP147 and CNP148, and several rare deletions and duplications. Examination of 34 AMD-enriched extended families (N = 293) and AMD cases (White N = 4210 Indian = 134; Malay = 140) and controls (White N = 3229; Indian = 117; Malay = 2390) demonstrated that deletion CNP148 was protective against AMD, independent of SNPs at CFH. Regression analysis of seven common haplotypes showed three haplotypes, H1, H6 and H7, as conferring risk for AMD development. Being the most common haplotype H1 confers the greatest risk by increasing the odds of AMD by 2.75-fold (95% CI = [2.51, 3.01]; p = 8.31×10−109); Caucasian (H6) and Indian-specific (H7) recombinant haplotypes increase the odds of AMD by 1.85-fold (p = 3.52×10−9) and by 15.57-fold (P = 0.007), respectively. We identified a 32-kb region downstream of Y402H (rs1061170), shared by all three risk haplotypes, suggesting that this region may be critical for AMD development. Further analysis showed that two SNPs within the 32 kb block, rs1329428 and rs203687, optimally explain disease association. rs1329428 resides in 20 kb unique sequence block, but rs203687 resides in a 12 kb block that is 89% similar to a noncoding region contained in ΔCNP148. We conclude that causal variation in this region potentially encompasses both regulatory effects at single markers and copy number

Modelling the Genetic Risk in Age-Related Macular Degeneration

Late-stage age-related macular degeneration (AMD) is a common sight-threatening disease of the central retina affecting approximately 1 in 30 Caucasians. Besides age and smoking, genetic variants from several gene loci have reproducibly been associated with this condition and likely explain a large proportion of disease. Here, we developed a genetic risk score (GRS) for AMD based on 13 risk variants from eight gene loci. The model exhibited good discriminative accuracy, area-under-curve (AUC) of the receiver-operating characteristic of 0.820, which was confirmed in a cross-validation approach. Noteworthy, younger AMD patients aged below 75 had a significantly higher mean GRS (1.87, 95% CI: 1.69–2.05) than patients aged 75 and above (1.45, 95% CI: 1.36–1.54). Based on five equally sized GRS intervals, we present a risk classification with a relative AMD risk of 64.0 (95% CI: 14.11–1131.96) for individuals in the highest category (GRS 3.44–5.18, 0.5% of the general population) compared to subjects with the most common genetic background (GRS −0.05–1.70, 40.2% of general population). The highest GRS category identifies AMD patients with a sensitivity of 7.9% and a specificity of 99.9% when compared to the four lower categories. Modeling a general population around 85 years of age, 87.4% of individuals in the highest GRS category would be expected to develop AMD by that age. In contrast, only 2.2% of individuals in the two lowest GRS categories which represent almost 50% of the general population are expected to manifest AMD. Our findings underscore the large proportion of AMD cases explained by genetics particularly for younger AMD patients. The five-category risk classification could be useful for therapeutic stratification or for diagnostic testing purposes once preventive treatment is available