Honghuang Lin

Meta-analysis identifies six new susceptibility loci for atrial fibrillation

Atrial fibrillation is a highly prevalent arrhythmia and a major risk factor for stroke, heart failure and death. We conducted a genome-wide association study (GWAS) in individuals of European ancestry, including 6,707 with and 52,426 without atrial fibrillation. Six new atrial fibrillation susceptibility loci were identified and replicated in an additional sample of individuals of European ancestry, including 5,381 subjects with and 10,030 subjects without atrial fibrillation (P < 5 × 10−8). Four of the loci identified in Europeans were further replicated in silico in a GWAS of Japanese individuals, including 843 individuals with and 3,350 individuals without atrial fibrillation. The identified loci implicate candidate genes that encode transcription factors related to cardiopulmonary development, cardiac-expressed ion channels and cell signaling molecules.

Evaluation of MHC-II peptide binding prediction servers: applications for vaccine research.

BackgroundInitiation and regulation of immune responses in humans involves recognition of peptides presented by human leukocyte antigen class II (HLA-II) molecules. These peptides (HLA-II T-cell epitopes) are increasingly important as research targets for the development of vaccines and immunotherapies. HLA-II peptide binding studies involve multiple overlapping peptides spanning individual antigens, as well as complete viral proteomes. Antigen variation in pathogens and tumor antigens, and extensive polymorphism of HLA molecules increase the number of targets for screening studies. Experimental screening methods are expensive and time consuming and reagents are not readily available for many of the HLA class II molecules. Computational prediction methods complement experimental studies, minimize the number of validation experiments, and significantly speed up the epitope mapping process. We collected test data from four independent studies that involved 721 peptide binding assays. Full overlapping studies of four antigens identified binding affinity of 103 peptides to seven common HLA-DR molecules (DRB1*0101, 0301, 0401, 0701, 1101, 1301, and 1501). We used these data to analyze performance of 21 HLA-II binding prediction servers accessible through the WWW.ResultsBecause not all servers have predictors for all tested HLA-II molecules, we assessed a total of 113 predictors. The length of test peptides ranged from 15 to 19 amino acids. We tried three prediction strategies – the best 9-mer within the longer peptide, the average of best three 9-mer predictions, and the average of all 9-mer predictions within the longer peptide. The best strategy was the identification of a single best 9-mer within the longer peptide. Overall, measured by the receiver operating characteristic method (AROC), 17 predictors showed good (AROC > 0.8), 41 showed marginal (AROC > 0.7), and 55 showed poor performance (AROC < 0.7). Good performance predictors included HLA-DRB1*0101 (seven), 1101 (six), 0401 (three), and 0701 (one). The best individual predictor was NETMHCIIPAN, closely followed by PROPRED, IEDB (Consensus), and MULTIPRED (SVM). None of the individual predictors was shown to be suitable for prediction of promiscuous peptides. Current predictive capabilities allow prediction of only 50% of actual T-cell epitopes using practical thresholds.ConclusionThe available HLA-II servers do not match prediction capabilities of HLA-I predictors. Currently available HLA-II prediction servers offer only a limited prediction accuracy and the development of improved predictors is needed for large-scale studies, such as proteome-wide epitope mapping. The requirements for accuracy of HLA-II binding predictions are stringent because of the substantial effect of false positives.

Evaluation of MHC class I peptide binding prediction servers: Applications for vaccine research

BackgroundProtein antigens and their specific epitopes are formulation targets for epitope-based vaccines. A number of prediction servers are available for identification of peptides that bind major histocompatibility complex class I (MHC-I) molecules. The lack of standardized methodology and large number of human MHC-I molecules make the selection of appropriate prediction servers difficult. This study reports a comparative evaluation of thirty prediction servers for seven human MHC-I molecules.ResultsOf 147 individual predictors 39 have shown excellent, 47 good, 33 marginal, and 28 poor ability to classify binders from non-binders. The classifiers for HLA-A*0201, A*0301, A*1101, B*0702, B*0801, and B*1501 have excellent, and for A*2402 moderate classification accuracy. Sixteen prediction servers predict peptide binding affinity to MHC-I molecules with high accuracy; correlation coefficients ranging from r = 0.55 (B*0801) to r = 0.87 (A*0201).ConclusionNon-linear predictors outperform matrix-based predictors. Most predictors can be improved by non-linear transformations of their raw prediction scores. The best predictors of peptide binding are also best in prediction of T-cell epitopes. We propose a new standard for MHC-I binding prediction – a common scale for normalization of prediction scores, applicable to both experimental and predicted data. The results of this study provide assistance to researchers in selection of most adequate prediction tools and selection criteria that suit the needs of their projects.

Atrial Fibrillation Current Knowledge and Future Directions in Epidemiology and Genomics

Atrial fibrillation (AF) is of public health importance and profoundly increases morbidity, mortality and health-related expenditures. Morbidities include the increased risks of cardiovascular outcomes such as heart failure and stroke, and the deleterious effects on quality of life, functional status and cognition. The clinical epidemiology of AF, its risk factors and outcomes, have been extensively investigated. Genetic advances over the last decade have facilitated the identification of mutations and common polymorphisms associated with AF. Metabolomics, proteomics and other “omics” technologies have only recently been applied to the study of AF, and have not yet been systematically investigated. Systems biology approaches, while still in their infancy, offer the promise of providing novel insights into pathways influencing AF risk. In the present review, we address the current state of the epidemiology and genomics of AF. We seek to emphasize how epidemiology and “omic” advances will contribute towards a systems biology approach that will help to unravel the pathogenesis, risk stratification, and novel targets for AF therapies. Our purpose is to articulate questions and challenges that hinge on integrating novel scientific advances in the epidemiology and genomics of AF. As a reference we have provided a glossary in the inset box.

Integrating Genetic, Transcriptional, and Functional Analyses to Identify 5 Novel Genes for Atrial Fibrillation

Background— Atrial fibrillation (AF) affects >30 million individuals worldwide and is associated with an increased risk of stroke, heart failure, and death. AF is highly heritable, yet the genetic basis for the arrhythmia remains incompletely understood. Methods and Results— To identify new AF-related genes, we used a multifaceted approach, combining large-scale genotyping in 2 ethnically distinct populations, cis-eQTL (expression quantitative trait loci) mapping, and functional validation. Four novel loci were identified in individuals of European descent near the genes NEURL (rs12415501; relative risk [RR]=1.18; 95% confidence interval [CI], 1.13–1.23; P=6.5×10−16), GJA1 (rs13216675; RR=1.10; 95% CI, 1.06–1.14; P=2.2×10−8), TBX5 (rs10507248; RR=1.12; 95% CI, 1.08–1.16; P=5.7×10−11), and CAND2 (rs4642101; RR=1.10; 95% CI, 1.06–1.14; P=9.8×10−9). In Japanese, novel loci were identified near NEURL (rs6584555; RR=1.32; 95% CI, 1.26–1.39; P=2.0×10−25) and CUX2 (rs6490029; RR=1.12; 95% CI, 1.08–1.16; P=3.9×10−9). The top single-nucleotide polymorphisms or their proxies were identified as cis-eQTLs for the genes CAND2 (P=2.6×10−19), GJA1 (P=2.66×10−6), and TBX5 (P=1.36×10−5). Knockdown of the zebrafish orthologs of NEURL and CAND2 resulted in prolongation of the atrial action potential duration (17% and 45%, respectively). Conclusions— We have identified 5 novel loci for AF. Our results expand the diversity of genetic pathways implicated in AF and provide novel molecular targets for future biological and pharmacological investigation.

Novel genetic markers associate with atrial fibrillation risk in Europeans and Japanese

OBJECTIVES This study sought to identify nonredundant atrial fibrillation (AF) genetic susceptibility signals and examine their cumulative relations with AF risk. BACKGROUND AF-associated loci span broad genomic regions that may contain multiple susceptibility signals. Whether multiple signals exist at AF loci has not been systematically explored. METHODS We performed association testing conditioned on the most significant, independently associated genetic markers at 9 established AF loci using 2 complementary techniques in 64,683 individuals of European ancestry (3,869 incident and 3,302 prevalent AF cases). Genetic risk scores were created and tested for association with AF in Europeans and an independent sample of 11,309 individuals of Japanese ancestry (7,916 prevalent AF cases). RESULTS We observed at least 4 distinct AF susceptibility signals on chromosome 4q25 upstream of PITX2, but not at the remaining 8 AF loci. A multilocus score comprised 12 genetic markers demonstrated an estimated 5-fold gradient in AF risk. We observed a similar spectrum of risk associated with these markers in Japanese. Regions containing AF signals on chromosome 4q25 displayed a greater degree of evolutionary conservation than the remainder of the locus, suggesting that they may tag regulatory elements. CONCLUSIONS The chromosome 4q25 AF locus is architecturally complex and harbors at least 4 AF susceptibility signals in individuals of European ancestry. Similar polygenic AF susceptibility exists between Europeans and Japanese. Future work is necessary to identify causal variants, determine mechanisms by which associated loci predispose to AF, and explore whether AF susceptibility signals classify individuals at risk for AF and related morbidity.

A support vector machines approach for virtual screening of active compounds of single and multiple mechanisms from large libraries at an improved hit-rate and enrichment factor.

Support vector machines (SVM) and other machine-learning (ML) methods have been explored as ligand-based virtual screening (VS) tools for facilitating lead discovery. While exhibiting good hit selection performance, in screening large compound libraries, these methods tend to produce lower hit-rate than those of the best performing VS tools, partly because their training-sets contain limited spectrum of inactive compounds. We tested whether the performance of SVM can be improved by using training-sets of diverse inactive compounds. In retrospective database screening of active compounds of single mechanism (HIV protease inhibitors, DHFR inhibitors, dopamine antagonists) and multiple mechanisms (CNS active agents) from large libraries of 2.986 million compounds, the yields, hit-rates, and enrichment factors of our SVM models are 52.4-78.0%, 4.7-73.8%, and 214-10,543, respectively, compared to those of 62-95%, 0.65-35%, and 20-1200 by structure-based VS and 55-81%, 0.2-0.7%, and 110-795 by other ligand-based VS tools in screening libraries of >or=1 million compounds. The hit-rates are comparable and the enrichment factors are substantially better than the best results of other VS tools. 24.3-87.6% of the predicted hits are outside the known hit families. SVM appears to be potentially useful for facilitating lead discovery in VS of large compound libraries.

PEARLS: program for energetic analysis of receptor-ligand system.

Analysis of the energetics of small molecule ligand-protein, ligand-nucleic acid, and protein-nucleic acid interactions facilitates the quantitative understanding of molecular interactions that regulate the function and conformation of proteins. It has also been extensively used for ranking potential new ligands in virtual drug screening. We developed a Web-based software, PEARLS (Program for Energetic Analysis of Ligand-Receptor Systems), for computing interaction energies of ligand-protein, ligand-nucleic acid, protein-nucleic acid, and ligand-protein-nucleic acid complexes from their 3D structures. AMBER molecular force field, Morse potential, and empirical energy functions are used to compute the van der Waals, electrostatic, hydrogen bond, metal-ligand bonding, and water-mediated hydrogen bond energies between the binding molecules. The change in the solvation free energy of molecular binding is estimated by using an empirical solvation free energy model. Contribution from ligand conformational entropy change is also estimated by a simple model. The computed free energy for a number of PDB ligand-receptor complexes were studied and compared to experimental binding affinity. A substantial degree of correlation between the computed free energy and experimental binding affinity was found, which suggests that PEARLS may be useful in facilitating energetic analysis of ligand-protein, ligand-nucleic acid, and protein-nucleic acid interactions. PEARLS can be accessed at http://ang.cz3.nus.edu.sg/cgi-bin/prog/rune.pl.

Common genetic variation at the IL1RL1 locus regulates IL-33/ST2 signaling

The suppression of tumorigenicity 2/IL-33 (ST2/IL-33) pathway has been implicated in several immune and inflammatory diseases. ST2 is produced as 2 isoforms. The membrane-bound isoform (ST2L) induces an immune response when bound to its ligand, IL-33. The other isoform is a soluble protein (sST2) that is thought to be a decoy receptor for IL-33 signaling. Elevated sST2 levels in serum are associated with an increased risk for cardiovascular disease. We investigated the determinants of sST2 plasma concentrations in 2,991 Framingham Offspring Cohort participants. While clinical and environmental factors explained some variation in sST2 levels, much of the variation in sST2 production was driven by genetic factors. In a genome-wide association study (GWAS), multiple SNPs within IL1RL1 (the gene encoding ST2) demonstrated associations with sST2 concentrations. Five missense variants of IL1RL1 correlated with higher sST2 levels in the GWAS and mapped to the intracellular domain of ST2, which is absent in sST2. In a cell culture model, IL1RL1 missense variants increased sST2 expression by inducing IL-33 expression and enhancing IL-33 responsiveness (via ST2L). Our data suggest that genetic variation in IL1RL1 can result in increased levels of sST2 and alter immune and inflammatory signaling through the ST2/IL-33 pathway.

Pleiotropic genes for metabolic syndrome and inflammation

Metabolic syndrome (MetS) has become a health and financial burden worldwide. The MetS definition captures clustering of risk factors that predict higher risk for diabetes mellitus and cardiovascular disease. Our study hypothesis is that additional to genes influencing individual MetS risk factors, genetic variants exist that influence MetS and inflammatory markers forming a predisposing MetS genetic network. To test this hypothesis a staged approach was undertaken. (a) We analyzed 17 metabolic and inflammatory traits in more than 85,500 participants from 14 large epidemiological studies within the Cross Consortia Pleiotropy Group. Individuals classified with MetS (NCEP definition), versus those without, showed on average significantly different levels for most inflammatory markers studied. (b) Paired average correlations between 8 metabolic traits and 9 inflammatory markers from the same studies as above, estimated with two methods, and factor analyses on large simulated data, helped in identifying 8 combinations of traits for follow-up in meta-analyses, out of 130,305 possible combinations between metabolic traits and inflammatory markers studied. (c) We performed correlated meta-analyses for 8 metabolic traits and 6 inflammatory markers by using existing GWAS published genetic summary results, with about 2.5 million SNPs from twelve predominantly largest GWAS consortia. These analyses yielded 130 unique SNPs/genes with pleiotropic associations (a SNP/gene associating at least one metabolic trait and one inflammatory marker). Of them twenty-five variants (seven loci newly reported) are proposed as MetS candidates. They map to genes MACF1, KIAA0754, GCKR, GRB14, COBLL1, LOC646736-IRS1, SLC39A8, NELFE, SKIV2L, STK19, TFAP2B, BAZ1B, BCL7B, TBL2, MLXIPL, LPL, TRIB1, ATXN2, HECTD4, PTPN11, ZNF664, PDXDC1, FTO, MC4R and TOMM40. Based on large data evidence, we conclude that inflammation is a feature of MetS and several gene variants show pleiotropic genetic associations across phenotypes and might explain a part of MetS correlated genetic architecture. These findings warrant further functional investigation.