Zhenwei Shang

Meta-Analysis of 125 Rheumatoid Arthritis-Related Single Nucleotide Polymorphisms Studied in the Past Two Decades

Objective Candidate gene association studies and genome-wide association studies (GWAs) have identified a large number of single nucleotide polymorphisms (SNPs) loci affecting susceptibility to rheumatoid arthritis (RA). However, for the same locus, some studies have yielded inconsistent results. To assess all the available evidence for association, we performed a meta-analysis on previously published case-control studies investigating the association between SNPs and RA. Methods Two hundred and sixteen studies, involving 125 SNPs, were reviewed. For each SNP, three genetic models were considered: the allele, dominant and recessive effects models. For each model, the effect summary odds ratio (OR) and 95% CIs were calculated. Cochran’s Q-statistics were used to assess heterogeneity. If the heterogeneity was high, a random effects model was used for meta-analysis, otherwise a fixed effects model was used. Results The meta-analysis results showed that: (1) 30, 28 and 26 SNPs were significantly associated with RA (P<0.01) for the allele, dominant, and recessive models, respectively. (2) rs2476601 (PTPN22) showed the strongest association for all the three models: OR = 1.605, 95% CI: 1.540–1.672, P<1.00E−15 for the T-allele; OR = 1.638, 95% CI: 1.565–1.714, P<1.00E−15 for the T/T+T/C genotype and OR = 2.544, 95% CI: 2.173–2.978, P<1.00E−15 for the T/T genotype. (3) Only 23 (18.4%), 13 (10.4%) and 15 (12.0%) SNPs had high heterogeneity (P<0.01) for the three models, respectively. (4) For some of the SNPs, there was no publication bias according to Funnel plots and Egger’s regression tests (P<0.01). For the other SNPs, the associations were tested in only a few studies, and may have been subject to publication bias. More studies on these loci are required. Conclusion Our meta-analysis provides a comprehensive evaluation of the RA association studies from the past two decades. The detailed meta-analysis results are available at: http://210.46.85.180/DRAP/index.php/Metaanalysis/index.

Predict potential drug targets from the ion channel proteins based on SVM.

The identification of molecular targets is a critical step in the drug discovery and development process. Ion channel proteins represent highly attractive drug targets implicated in a diverse range of disorders, in particular in the cardiovascular and central nervous systems. Due to the limits of experimental technique and low-throughput nature of patch-clamp electrophysiology, they remain a target class waiting to be exploited. In our study, we combined three types of protein features, primary sequence, secondary structure and subcellular localization to predict potential drug targets from ion channel proteins applying classical support vector machine (SVM) method. In addition, our prediction comprised two stages. In stage 1, we predicted ion channel target proteins based on whole-genome target protein characteristics. Firstly, we performed feature selection by Mann-Whitney U test, then made predictions to identify potential ion channel targets by SVM and designed a new evaluating indicator Q to prioritize results. In stage 2, we made a prediction based on known ion channel target protein characteristics. Genetic algorithm was used to select features and SVM was used to predict ion channel targets. Then, we integrated results of two stages, and found that five ion channel proteins appeared in both prediction results including CGMP-gated cation channel beta subunit and Gamma-aminobutyric acid receptor subunit alpha-5, etc., and four of which were relative to some nerve diseases. It suggests that these five proteins are potential targets for drug discovery and our prediction strategies are effective.

Prioritization of Candidate Genes for Periodontitis Using Multiple Computational Tools

BACKGROUND Both genetic and environmental factors contribute to the development of periodontitis. Genetic studies identified a variety of candidate genes for periodontitis. The aim of the present study is to identify the most promising candidate genes for periodontitis using an integrative gene ranking method. METHODS Seed genes that were confirmed to be associated with periodontitis were identified using text mining. Three types of candidate genes were then extracted from different resources (expression profiles, genome-wide association studies). Combining the seed genes, four freely available bioinformatics tools (ToppGene, DIR, Endeavour, and GPEC) were integrated for prioritization of candidate genes. Candidate genes that identified with at least three programs and ranked in the top 20 by each program were considered the most promising. RESULTS Prioritization analysis resulted in 21 promising genes involved or potentially involved in periodontitis. Among them, IL18 (interleukin 18), CD44 (CD44 molecule), CXCL1 (chemokine [CXC motif] ligand 1), IL6ST (interleukin 6 signal transducer), MMP3 (matrix metallopeptidase 3), MMP7, CCR1 (chemokine [C-C motif] receptor 1), MMP13, and TLR9 (Toll-like receptor 9) had been associated with periodontitis. However, the roles of other genes, such as CSF3 (colony stimulating factor 3 receptor), CD40, TNFSF14 (tumor necrosis factor receptor superfamily, member 14), IFNB1 (interferon-β1), TIRAP (toll-interleukin 1 receptor domain containing adaptor protein), IL2RA (interleukin 2 receptor α), ETS1 (v-ets avian erythroblastosis virus E26 oncogene homolog 1), GADD45B (growth arrest and DNA-damage-inducible 45 β), BIRC3 (baculoviral IAP repeat containing 3), VAV1 (vav 1 guanine nucleotide exchange factor), COL5A1 (collagen, type V, α1), and C3 (complement component 3), have not been investigated thoroughly in the process of periodontitis. These genes are mainly involved in bacterial infection, immune response, and inflammatory reaction, suggesting that further characterizing their roles in periodontitis will be important. CONCLUSIONS A combination of computational tools will be useful in mining candidate genes for periodontitis. These theoretical results provide new clues for experimental biologists to plan targeted experiments.

Genome-wide haplotype association study identify TNFRSF1A, CASP7, LRP1B, CDH1 and TG genes associated with Alzheimer’s disease in Caribbean Hispanic individuals

Alzheimers disease (AD) is an acquired disorder of cognitive and behavioral impairment. It is considered to be caused by variety of factors, such as age, environment and genetic factors. In order to identify the genetic affect factors of AD, we carried out a bioinformatic approach which combined genome-wide haplotype-based association study with gene prioritization. The raw SNP genotypes data was downloaded from GEO database (GSE33528). It contains 615 AD patients and 560 controls of Caribbean Hispanic individuals. Firstly, we identified the linkage disequilibrium (LD) haplotype blocks and performed genome-wide haplotype association study to screen significant haplotypes that were associated with AD. Then we mapped these significant haplotypes to genes and obtained candidate genes set for AD. At last, we prioritized AD candidate genes based on their similarity with 36 known AD genes, so as to identify AD related genes. The results showed that 141 haplotypes on 134 LD blocks were significantly associated with AD (P<1E-4), and these significant haplotypes were mapped to 132 AD candidate genes. After prioritizing these candidate genes, we found seven AD related genes: APOE, APOC1, TNFRSF1A, LRP1B, CDH1, TG and CASP7. Among these genes, APOE and APOC1 are known AD risk genes. For the other five genes TNFRSF1A, CDH1, CASP7, LRP1B and TG, this is the first genetic association study which showed the significant association between these five genes and AD susceptibility in Caribbean Hispanic individuals. We believe that our findings can provide a new perspective to understand the genetic affect factors of AD.

The drug target genes show higher evolutionary conservation than non-target genes

Although evidence indicates that drug target genes share some common evolutionary features, there have been few studies analyzing evolutionary features of drug targets from an overall level. Therefore, we conducted an analysis which aimed to investigate the evolutionary characteristics of drug target genes. We compared the evolutionary conservation between human drug target genes and non-target genes by combining both the evolutionary features and network topological properties in human protein-protein interaction network. The evolution rate, conservation score and the percentage of orthologous genes of 21 species were included in our study. Meanwhile, four topological features including the average shortest path length, betweenness centrality, clustering coefficient and degree were considered for comparison analysis. Then we got four results as following: compared with non-drug target genes, 1) drug target genes had lower evolutionary rates; 2) drug target genes had higher conservation scores; 3) drug target genes had higher percentages of orthologous genes and 4) drug target genes had a tighter network structure including higher degrees, betweenness centrality, clustering coefficients and lower average shortest path lengths. These results demonstrate that drug target genes are more evolutionarily conserved than non-drug target genes. We hope that our study will provide valuable information for other researchers who are interested in evolutionary conservation of drug targets.

Genome-wide pathway-based association analysis identifies risk pathways associated with Parkinson’s disease

Parkinsons disease (PD) is the second most common neurodegenerative disease. It is generally believed that it is influenced by both genetic and environmental factors, but the precise pathogenesis of PD is unknown to date. In this study, we performed a pathway analysis based on genome-wide association study (GWAS) to detect risk pathways of PD in three GWAS datasets. We first mapped all SNP markers to autosomal genes in each GWAS dataset. Then, we evaluated gene risk values using the minimum P-value of the tagSNPs. We took a pathway as a unit to identify the risk pathways based on the cumulative risks of the genes in the pathway. Finally, we combine the analysis results of the three datasets to detect the high risk pathways associated with PD. We found there were five same pathways in the three datasets. Besides, we also found there were five pathways which were shared in two datasets. Most of these pathways are associated with nervoussystem. Five pathways had been reported to be PD-related pathways in the previous literature. Our findings also implied that there was a close association between immune response and PD. Continued investigation of these pathways will further help us explain the pathogenesis of PD.

Genes with stable DNA methylation levels show higher evolutionary conservation than genes with fluctuant DNA methylation levels

Different human genes often exhibit different degrees of stability in their DNA methylation levels between tissues, samples or cell types. This may be related to the evolution of human genome. Thus, we compared the evolutionary conservation between two types of genes: genes with stable DNA methylation levels (SM genes) and genes with fluctuant DNA methylation levels (FM genes). For long-term evolutionary characteristics between species, we compared the percentage of the orthologous genes, evolutionary rate dn/ds and protein sequence identity. We found that the SM genes had greater percentages of the orthologous genes, lower dn/ds, and higher protein sequence identities in all the 21 species. These results indicated that the SM genes were more evolutionarily conserved than the FM genes. For short-term evolutionary characteristics among human populations, we compared the single nucleotide polymorphism (SNP) density, and the linkage disequilibrium (LD) degree in HapMap populations and 1000 genomes project populations. We observed that the SM genes had lower SNP densities, and higher degrees of LD in all the 11 HapMap populations and 13 1000 genomes project populations. These results mean that the SM genes had more stable chromosome genetic structures, and were more conserved than the FM genes.

DBGSA: a novel method of distance-based gene set analysis

When compared with single gene functional analysis, gene set analysis (GSA) can extract more information from gene expression profiles. Currently, several gene set methods have been proposed, but most of the methods cannot detect gene sets with a large number of minor-effect genes. Here, we propose a novel distance-based gene set analysis method. The distance between two groups of genes with different phenotypes based on gene expression should be larger if a certain gene set is significantly associated with the given phenotype. We calculated the distance between two groups with different phenotypes, estimated the significant P-values using two permutation methods and performed multiple hypothesis testing adjustments. This method was performed on one simulated data set and three real data sets. After a comparison and literature verification, we determined that the gene resampling-based permutation method is more suitable for GSA, and the centroid statistical and average linkage statistical distance methods are efficient, especially in detecting gene sets containing more minor-effect genes. We believe that this distance-based method will assist us in finding functional gene sets that are significantly related to a complex trait. Additionally, we have prepared a simple and publically available Perl and R package (http://bioinfo.hrbmu.edu.cn/dbgsa or http://cran.r-project.org/web/packages/DBGSA/).

The shared and specific mechanism of four autoimmune diseases

Interaction between genetic and epigenetic mechanisms may lead to autoimmune diseases. The features of these diseases show familial aggregation. The generality and specificity are keys to studying pathogenesis and etiology of them. This research integrated data of genetics and epigenetics, to find disease-related genes based on the levels of expression and regulation, and explored then to the shared and specific mechanism of them by analyzing shared and specific pathways of common four autoimmune diseases, including Type 1 Diabetes Mellitus (T1D), Multiple Sclerosis (MS), Rheumatoid Arthritis (RA) and Systemic Lupus Erythematosus (SLE). The results showed that Lysosome and Fc gamma R-mediated phagocytosis are shared pathways of the four diseases. It means that the occurrence and development of them may associate with lysosomes and phagocytosis. And there were 2 pathways are shared pathways of three diseases, ribosome pathway associated with susceptibility to MS, RA and SLE, and Pathogenic Escherichia coli infection associated with susceptibility to T1D, MS and RA; 9 pathways are shared pathways of two diseases. The corporate underlying causes of these diseases may be these shared pathways activated. Furthermore, we found that T1D-related specific pathways (Insulin signaling,etc.) were 9, MS (Proteasome,etc.) is also 9, RA and SLE is 10 and 6 respectively. These pathways could help us to reveal shared and specific mechanisms of the four diseases.

Genome-wide haplotype association study identify the FGFR2 gene as a risk gene for acute myeloid leukemia.

Acute myeloid leukemia (AML) is a cancer of the myeloid line of blood cells, and generally considered to be caused by environment and genetic factors. In this study, we combined a genome-wide haplotype association study (GWHAS) and gene prioritization strategy to mine AML-related genetic affect factors and understand its pathogenesis. A total of 175 AML patients were downloaded from the public GEO database (GSE32462) and 218 matched Caucasian controls were from the HapMap Project. We first identified the linkage disequilibrium (LD) blocks and performed a GWHAS to scan AML-related haplotypes. Then we mapped these haplotypes to the corresponding genes as candidate. And finally, we prioritized all the AML candidate genes based on the similarity with 38 known AML susceptibility genes. The results showed that 1754 haplotypes were significant associated with AML (P<1E-5) and mapped to 591 candidate genes. After prioritizing all 591 AML candidate genes, we obtained four genes ranking at the front as AML risk genes: RUNX1, JAK1, PDGFRA, and FGFR2. Among them, RUNX1, JAK1 and PDGFRA had been confirmed as AML risk genes. In particular, we found that the gene FGFR2 was a novel AML susceptibility gene with a haplotype TT (rs7090018 and rs2912759) showed significant association with AML (P-value = 7.07E-06). In a word, our findings might provide a new perspective to understand the pathogenesis of AML.