The application of network label propagation to rank biomarkers in genome-wide Alzheimer's data

Background: Ranking and identifying biomarkers that are associated with disease from genome-wide measurements holds significant promise for understanding the genetic basis of common diseases. The large number of single nucleotide polymorphisms (SNPs) in genome-wide studies (GWAS), however, makes this task computationally challenging when the ranking is to be done in a multivariate fashion. This paper evaluates the performance of a multivariate graph-based method called label propagation (LP) that efficiently ranks SNPs in genome-wide data.Results: The performance of LP was evaluated on a synthetic dataset and two late onset Alzheimer's disease (LOAD) genome-wide datasets, and the performance was compared to that of three control methods. The control methods included chi squared, which is a commonly used univariate method, as well as a Relief method called SWRF and a sparse logistic regression (SLR) method, which are both multivariate ranking methods. Performance was measured by evaluating the top-ranked SNPs in terms of classification performance, reproducibility between the two datasets, and prior evidence of being associated with LOAD.On the synthetic data LP performed comparably to the control methods. On GWAS data, LP performed significantly better than chi squared and SWRF in classification performance in the range from 10 to 1000 top-ranked SNPs for both datasets, and not significantly different from SLR. LP also had greater ranking reproducibility than chi squared, SWRF, and SLR. Among the 25 top-ranked SNPs that were identified by LP, there were 14 SNPs in one dataset that had evidence in the literature of being associated with LOAD, and 10 SNPs in the other, which was higher than for the other methods.Conclusion: LP performed considerably better in ranking SNPs in two high-dimensional genome-wide datasets when compared to three control methods. It had better performance in the evaluation measures we used, and is computationally efficient to be applied practically to data from genome-wide studies. These results provide support for including LP in the methods that are used to rank SNPs in genome-wide datasets. © 2014 Stokes et al.; licensee BioMed Central Ltd

Genetic contribution of SCARB1 variants to lipid traits in African Blacks: A candidate gene association study

Background: High-density lipoprotein cholesterol (HDL-C) exerts many anti-atherogenic properties including its role in reverse cholesterol transport (RCT). Scavenger receptor class B member 1 (SCARB1) plays a key role in RCT by selective uptake of HDL cholesteryl esters. We aimed to explore the genetic contribution of SCARB1 to affecting lipid levels in African Blacks from Nigeria. Methods: We resequenced 13 exons and exon-intron boundaries of SCARB1 in 95 individuals with extreme HDL-C levels using Sanger method. Then, we genotyped 147 selected variants (78 sequence variants, 69 HapMap tagSNPs, and 2 previously reported relevant variants) in the entire sample of 788 African Blacks using either the iPLEX Gold or TaqMan methods. A total of 137 successfully genotyped variants were further evaluated for association with major lipid traits. Results: The initial gene-based analysis demonstrated evidence of association with HDL-C and apolipoprotein A-I (ApoA-I). The follow-up single-site analysis revealed nominal evidence of novel associations of nine common variants with HDL-C and/or ApoA-I (P < 0.05). The strongest association was between rs11057851 and HDL-C (P = 0.0043), which remained significant after controlling for multiple testing using false discovery rate. Rare variant association testing revealed a group of 23 rare variants (frequencies ≤1 %) associated with HDL-C (P = 0.0478). Haplotype analysis identified four SCARB1 regions associated with HDL-C (global P < 0.05). Conclusions: To our knowledge, this is the first report of a comprehensive association study of SCARB1 variations with lipid traits in an African Black population. Our results showed the consistent association of SCARB1 variants with HDL-C across various association analyses, supporting the role of SCARB1 in lipoprotein-lipid regulatory mechanism

A rare duplication on chromosome 16p11.2 is identified in patients with psychosis in Alzheimer's disease

Epidemiological and genetic studies suggest that schizophrenia and autism may share genetic links. Besides common single nucleotide polymorphisms, recent data suggest that some rare copy number variants (CNVs) are risk factors for both disorders. Because we have previously found that schizophrenia and psychosis in Alzheimer's disease (AD+P) share some genetic risk, we investigated whether CNVs reported in schizophrenia and autism are also linked to AD+P. We searched for CNVs associated with AD+P in 7 recurrent CNV regions that have been previously identified across autism and schizophrenia, using the Illumina HumanOmni1-Quad BeadChip. A chromosome 16p11.2 duplication CNV (chr16: 29,554,843-30,105,652) was identified in 2 of 440 AD+P subjects, but not in 136 AD subjects without psychosis, or in 593 AD subjects with intermediate psychosis status, or in 855 non-AD individuals. The frequency of this duplication CNV in AD+P (0.46%) was similar to that reported previously in schizophrenia (0.46%). This duplication CNV was further validated using the NanoString nCounter CNV Custom CodeSets. The 16p11.2 duplication has been associated with developmental delay, intellectual disability, behavioral problems, autism, schizophrenia (SCZ), and bipolar disorder. These two AD+P patients had no personal of, nor any identified family history of, SCZ, bipolar disorder and autism. To the best of our knowledge, our case report is the first suggestion that 16p11.2 duplication is also linked to AD+P. Although rare, this CNV may have an important role in the development of psychosis

Identifying genetic interactions associated with late-onset Alzheimer's disease

Background: Identifying genetic interactions in data obtained from genome-wide association studies (GWASs) can help in understanding the genetic basis of complex diseases. The large number of single nucleotide polymorphisms (SNPs) in GWASs however makes the identification of genetic interactions computationally challenging. We developed the Bayesian Combinatorial Method (BCM) that can identify pairs of SNPs that in combination have high statistical association with disease. Results: We applied BCM to two late-onset Alzheimer's disease (LOAD) GWAS datasets to identify SNPs that interact with known Alzheimer associated SNPs. We also compared BCM with logistic regression that is implemented in PLINK. Gene Ontology analysis of genes from the top 200 dataset SNPs for both GWAS datasets showed overrepresentation of LOAD-related terms. Four genes were common to both datasets: APOE and APOC1, which have well established associations with LOAD, and CAMK1D and FBXL13, not previously linked to LOAD but having evidence of involvement in LOAD. Supporting evidence was also found for additional genes from the top 30 dataset SNPs. Conclusion: BCM performed well in identifying several SNPs having evidence of involvement in the pathogenesis of LOAD that would not have been identified by univariate analysis due to small main effect. These results provide support for applying BCM to identify potential genetic variants such as SNPs from high dimensional GWAS datasets

Comprehensive evaluation of the association of APOE genetic variation with plasma lipoprotein traits in U.S. Whites and African Blacks

Although common APOE genetic variation has a major influence on plasma LDL-cholesterol, its role in affecting HDL-cholesterol and triglycerides is not well established. Recent genome-wide association studies suggest that APOE also affects plasma variation in HDL-cholesterol and triglycerides. It is thus important to resequence the APOE gene to identify both common and uncommon variants that affect plasma lipid profile. Here, we have sequenced the APOE gene in 190 subjects with extreme HDL-cholesterol levels selected from two well-defined epidemiological samples of U.S. non-Hispanic Whites (NHWs) and African Blacks followed by genotyping of identified variants in the entire datasets (623 NHWs, 788 African Blacks) and association analyses with major lipid traits. We identified a total of 40 sequence variants, of which 10 are novel. A total of 32 variants, including common tagSNPs (≥5% frequency) and all uncommon variants (<5% frequency) were successfully genotyped and considered for genotype-phenotype associations. Other than the established associations of APOE∗2 and APOE∗4 with LDL-cholesterol, we have identified additional independent associations with LDL-cholesterol. We have also identified multiple associations of uncommon and common APOE variants with HDL-cholesterol and triglycerides. Our comprehensive sequencing and genotype-phenotype analyses indicate that APOE genetic variation impacts HDL-cholesterol and triglycerides in addition to affecting LDL-cholesterol

Needs assessment for research use of high- throughput sequencing at a large academic medical center

Next Generation Sequencing (NGS) methods are driving profound changes in biomedical research, with a growing impact on patient care. Many academic medical centers are evaluating potential models to prepare for the rapid increase in NGS information needs. This study sought to investigate (1) how and where sequencing data is generated and analyzed, (2) research objectives and goals for NGS, (3) workforce capacity and unmet needs, (4) storage capacity and unmet needs, (5) available and anticipated funding resources, and (6) future challenges. As a precursor to informed decision making at our institution, we undertook a systematic needs assessment of investigators using survey methods. We recruited 331 investigators from over 60 departments and divisions at the University of Pittsburgh Schools of Health Sciences and had 140 respondents, or a 42% response rate. Results suggest that both sequencing and analysis bottlenecks currently exist. Significant educational needs were identified, including both investigator-focused needs, such as selection of NGS methods suitable for specific research objectives, and program-focused needs, such as support for training an analytic workforce. The absence of centralized infrastructure was identified as an important institutional gap. Key principles for organizations managing this change were formulated based on the survey responses. This needs assessment provides an in-depth case study which may be useful to other academic medical centers as they identify and plan for future needs

Mechanisms of CFTR Functional Variants That Impair Regulated Bicarbonate Permeation and Increase Risk for Pancreatitis but Not for Cystic Fibrosis

CFTR is a dynamically regulated anion channel. Intracellular WNK1-SPAK activation causes CFTR to change permeability and conductance characteristics from a chloride-preferring to bicarbonate-preferring channel through unknown mechanisms. Two severe CFTR mutations (CFTRsev) cause complete loss of CFTR function and result in cystic fibrosis (CF), a severe genetic disorder affecting sweat glands, nasal sinuses, lungs, pancreas, liver, intestines, and male reproductive system. We hypothesize that those CFTR mutations that disrupt the WNK1-SPAK activation mechanisms cause a selective, bicarbonate defect in channel function (CFTRBD) affecting organs that utilize CFTR for bicarbonate secretion (e.g. the pancreas, nasal sinus, vas deferens) but do not cause typical CF. To understand the structural and functional requirements of the CFTR bicarbonate-preferring channel, we (a) screened 984 well-phenotyped pancreatitis cases for candidate CFTRBD mutations from among 81 previously described CFTR variants; (b) conducted electrophysiology studies on clones of variants found in pancreatitis but not CF; (c) computationally constructed a new, complete structural model of CFTR for molecular dynamics simulation of wild-type and mutant variants; and (d) tested the newly defined CFTRBD variants for disease in non-pancreas organs utilizing CFTR for bicarbonate secretion. Nine variants (CFTR R74Q, R75Q, R117H, R170H, L967S, L997F, D1152H, S1235R, and D1270N) not associated with typical CF were associated with pancreatitis (OR 1.5, p = 0.002). Clones expressed in HEK 293T cells had normal chloride but not bicarbonate permeability and conductance with WNK1-SPAK activation. Molecular dynamics simulations suggest physical restriction of the CFTR channel and altered dynamic channel regulation. Comparing pancreatitis patients and controls, CFTRBD increased risk for rhinosinusitis (OR 2.3, p<0.005) and male infertility (OR 395, p≪0.0001). WNK1-SPAK pathway-activated increases in CFTR bicarbonate permeability are altered by CFTRBD variants through multiple mechanisms. CFTRBD variants are associated with clinically significant disorders of the pancreas, sinuses, and male reproductive system. © 2014 Whitcomb et al

Involvement of Endoplasmic Reticulum Stress in Inflammatory Bowel Disease: A Different Implication for Colonic and Ileal Disease?

Background: Endoplasmic reticulum (ER) stress has been suggested to play a role in inflammatory bowel disease (IBD). The three branches (ATF6, IRE1 and PERK) of the unfolded protein response (UPR) have different roles and are not necessarily activated simultaneously. Methodology/Principal Findings: Expression of UPR-related genes was investigated in colonic and ileal biopsies of 23 controls, 15 ulcerative colitis (UC) and 54 Crohn's disease (CD) patients. This expression was confirmed at protein level in colonic and ileal samples of five controls, UC and CD patients. HSPA5, PDIA4 and XBP1s were significantly increased in colonic IBD at mRNA and/or protein levels, indicating activation of the ATF6 and IRE1 branch. Colonic IBD was associated with increased phosphorylation of EIF2A suggesting the activation of the PERK branch, but subsequent induction of GADD34 was not observed. In ileal CD, no differential expression of the UPR-related genes was observed, but our data suggested a higher basal activation of the UPR in the ileal mucosa of controls. This was confirmed by the increased expression of 16 UPR-related genes as 12 of them were significantly more expressed in ileal controls compared to colonic controls. Tunicamycin stimulation of colonic and ileal samples of healthy individuals revealed that although the ileal mucosa is exhibiting this higher basal UPR activation, it is still responsive to ER stress, even more than colonic mucosa. Conclusions/Significance: Activation of the three UPR-related arms is seen in colonic IBD-associated inflammation. However, despite EIF2A activation, inflamed colonic tissue did not increase GADD34 expression, which is usually involved in re-establishment of ER homeostasis. This study also implies the presence of a constitutive UPR activation in healthy ileal mucosa, with no further activation during inflammation. Therefore, engagement of the UPR differs between colon and ileum and this could be a factor in the development of ileal or colonic disease

TDP-43 Identified from a Genome Wide RNAi Screen for SOD1 Regulators

Amyotrophic Lateral Sclerosis (ALS) is a late-onset, progressive neurodegenerative disease affecting motor neurons in the brain stem and spinal cord leading to loss of voluntary muscular function and ultimately, death due to respiratory failure. A subset of ALS cases are familial and associated with mutations in superoxide dismutase 1 (SOD1) that destabilize the protein and predispose it to aggregation. In spite of the fact that sporadic and familial forms of ALS share many common patho-physiological features, the mechanistic relationship between SOD1-associated and sporadic forms of the disease if any, is not well understood. To better understand any molecular connections, a cell-based protein folding assay was employed to screen a whole genome RNAi library for genes that regulate levels of soluble SOD1. Statistically significant hits that modulate SOD1 levels, when analyzed by pathway analysis revealed a highly ranked network containing TAR DNA binging protein (TDP-43), a major component of aggregates characteristic of sporadic ALS. Biochemical experiments confirmed the action of TDP-43 on SOD1. These results highlight an unexpected relationship between TDP-43 and SOD1 which may have implications in disease pathogenesis

Genetic Evidence Supporting the Association of Protease and Protease Inhibitor Genes with Inflammatory Bowel Disease: A Systematic Review

As part of the European research consortium IBDase, we addressed the role of proteases and protease inhibitors (P/PIs) in inflammatory bowel disease (IBD), characterized by chronic mucosal inflammation of the gastrointestinal tract, which affects 2.2 million people in Europe and 1.4 million people in North America. We systematically reviewed all published genetic studies on populations of European ancestry (67 studies on Crohn's disease [CD] and 37 studies on ulcerative colitis [UC]) to identify critical genomic regions associated with IBD. We developed a computer algorithm to map the 807 P/PI genes with exact genomic locations listed in the MEROPS database of peptidases onto these critical regions and to rank P/PI genes according to the accumulated evidence for their association with CD and UC. 82 P/PI genes (75 coding for proteases and 7 coding for protease inhibitors) were retained for CD based on the accumulated evidence. The cylindromatosis/turban tumor syndrome gene (CYLD) on chromosome 16 ranked highest, followed by acylaminoacyl-peptidase (APEH), dystroglycan (DAG1), macrophage-stimulating protein (MST1) and ubiquitin-specific peptidase 4 (USP4), all located on chromosome 3. For UC, 18 P/PI genes were retained (14 proteases and 4protease inhibitors), with a considerably lower amount of accumulated evidence. The ranking of P/PI genes as established in this systematic review is currently used to guide validation studies of candidate P/PI genes, and their functional characterization in interdisciplinary mechanistic studies in vitro and in vivo as part of IBDase. The approach used here overcomes some of the problems encountered when subjectively selecting genes for further evaluation and could be applied to any complex disease and gene family