Genetics and beyond--the transcriptome of human monocytes and disease susceptibility

Coronary Artery Disease-Associated LIPA Coding Variant rs1051338 Reduces Lysosomal Acid Lipase Levels and Activity in Lysosomes

Supplemental Digital Content is available in the text.

Objective—

Genome-wide association studies have linked variants at chromosome 10q23 with increased coronary artery disease risk. The disease-associated variants fall in LIPA, which encodes lysosomal acid lipase (LAL), the enzyme responsible for lysosomal cholesteryl ester hydrolysis. Loss-of-function mutations in LIPA result in accelerated atherosclerosis. Surprisingly, the coronary artery disease variants are associated with increased LIPA expression in some cell types. In this study, we address this apparent contradiction.

Approach and Results—

We investigated a coding variant rs1051338, which is in high linkage disequilibrium (r²=0.89) with the genome-wide association study lead–associated variant rs2246833 and causes a nonsynonymous threonine to proline change within the signal peptide of LAL. Transfection of allele-specific expression constructs showed that the risk allele results in reduced lysosomal LAL protein (P=0.004) and activity (P=0.005). Investigation of LAL localization and turnover showed the risk LAL protein is degraded more quickly. This mechanism was confirmed in disease-relevant macrophages from individuals homozygous for either the nonrisk or risk allele. There was no difference in LAL protein or activity in whole macrophage extracts; however, we found reduced LAL protein (P=0.02) and activity (P=0.026) with the risk genotype in lysosomal extracts, suggesting that the risk genotype affects lysosomal LAL activity. Inhibition of the proteasome resulted in equal amounts of lysosomal LAL protein in risk and nonrisk macrophages.

Conclusions—

Our findings show that the coronary artery disease–associated coding variant rs1051338 causes reduced lysosomal LAL protein and activity because of increased LAL degradation, providing a plausible causal mechanism of increased coronary artery disease risk.

Reference intervals of plasma homoarginine from the German Gutenberg Health Study

BACKGROUND

Low circulating homoarginine has been associated with adverse cardiovascular (CV) outcome and mortality in patients at risk and in the general population. The present study aimed to define plasma homoarginine reference intervals from a representative population sample to improve risk stratification between healthy individuals and individuals at risk.

METHODS

We determined age- and sex-specific reference intervals for circulating plasma homoarginine in a subgroup of 786 healthy participants (no CV disease or risk factors) of the Gutenberg Health Study. Homoarginine concentrations were measured using a validated liquid chromatography-tandem mass spectrometry method.

RESULTS

Median EDTA plasma homoarginine concentration was 1.88 [25th; 75th percentile, 1.47; 2.41] μmol/L, with lower concentrations in women (1.77 [1.38; 2.26] μmol/L) than in men (2.01 [1.61; 2.56] μmol/L; p<0.001). Sex-specific 2.5th and 97.5th percentiles of reference intervals were 0.84 and 3.89 μmol/L in women and 0.98 and 4.10 μmol/L in men, respectively. Homoarginine concentrations also depended on age and single nucleotide polymorphisms related to the L-arginine:glycine amidinotransferase gene.

CONCLUSIONS

We provide plasma homoarginine reference intervals in men and women of the general population. The determination of homoarginine levels might be favorable for individual risk stratification.

Cis-eQTL-based trans-ethnic meta-analysis reveals novel genes associated with breast cancer risk

Breast cancer is the most common solid organ malignancy and the most frequent cause of cancer death among women worldwide. Previous research has yielded insights into its genetic etiology, but there remains a gap in the understanding of genetic factors that contribute to risk, and particularly in the biological mechanisms by which genetic variation modulates risk. The National Cancer Institute's "Up for a Challenge" (U4C) competition provided an opportunity to further elucidate the genetic basis of the disease. Our group leveraged the seven datasets made available by the U4C organizers and data from the publicly available UK Biobank cohort to examine associations between imputed gene expression and breast cancer risk. In particular, we used reference datasets describing the breast tissue and whole blood transcriptomes to impute expression levels in breast cancer cases and controls. In trans-ethnic meta-analyses of U4C and UK Biobank data, we found significant associations between breast cancer risk and the expression of RCCD1 (joint p-value: 3.6x10-06) and DHODH (p-value: 7.1x10-06) in breast tissue, as well as a suggestive association for ANKLE1 (p-value: 9.3x10-05). Expression of RCCD1 in whole blood was also suggestively associated with disease risk (p-value: 1.2x10-05), as were expression of ACAP1 (p-value: 1.9x10-05) and LRRC25 (p-value: 5.2x10-05). While genome-wide association studies (GWAS) have implicated RCCD1 and ANKLE1 in breast cancer risk, they have not identified the remaining three genes. Among the genetic variants that contributed to the predicted expression of the five genes, we found 23 nominally (p-value < 0.05) associated with breast cancer risk, among which 15 are not in high linkage disequilibrium with risk variants previously identified by GWAS. In summary, we used a transcriptome-based approach to investigate the genetic underpinnings of breast carcinogenesis. This approach provided an avenue for deciphering the functional relevance of genes and genetic variants involved in breast cancer.

Age-Dependent Association of TNFSF15/TNFSF8 Variants and Leprosy Type 1 Reaction

A current major challenge in leprosy control is the prevention of permanent disabilities. Host pathological inflammatory responses termed type 1 reaction (T1R) are a leading cause of nerve damage for leprosy patients. The environmental or inherited factors that predispose leprosy cases to undergo T1R are not known. However, studies have shown an important contribution of host genetics for susceptibility to T1R. We have previously identified variants encompassing the TNFSF15/TNFSF8 genes as T1R risk factors in a Vietnamese sample and replicated this association in a Brazilian sample. However, we failed to validate in Brazilian patients the strong association of TNFSF15/TNFSF8 markers rs6478108 and rs7863183 with T1R that we had observed in Vietnamese patients. Here, we investigated if the lack of validation of these variants was due to age-dependent effects on association using four independent population samples, two from Brazil and two from Vietnam. In the combined analysis across the four samples, we observed a strong association of the TNFSF15/TNFSF8 variants rs6478108, rs7863183, and rs3181348 with T1R (p_combined = 1.5E−05, p_combined = 1.8E−05, and p_combined = 6.5E−06, respectively). However, the association of rs6478108 with T1R was more pronounced in leprosy cases under 30 years of age compared to the global sample [odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.54–2.46, p_combined = 2.5E−08 versus OR = 1.46, 95% CI = 1.23–1.73, p_combined = 1.5E−05]. A multivariable analysis indicated that the association of rs6478108 with T1R was independent of either rs7863183 or rs3181348. These three variants are known regulators of the TNFSF8 gene transcription level in multiple tissues. The age dependency of association of rs6478108 and T1R suggests that the genetic control of gene expression varies across the human life span.

Fine-mapping of antipsychotic response genome-wide association studies reveals novel regulatory mechanisms

Aim: Noncoding variation has demonstrated regulatory effects on disease treatment outcomes. This study investigated the potential functionality of previously implicated noncoding variants on schizophrenia treatment response. Materials & methods: Predicted regulatory potential of variation identified from antipsychotic response genome-wide association studies was determined. Prioritized variants were assessed for association(s) with treatment outcomes in a South African first episode schizophrenia cohort (n = 103). Results: Bioinformatic and association results implicated a relationship between regulatory variants, expression of MANBA, COL9A2 and NFKB1, and treatment response. Three SNPs were associated with poor outcomes (rs230493: p = 1.88 × 10-6; rs3774959: p = 1.75 × 10-5; and rs230504: p = 1.48 × 10-4). Conclusion: This study has thoroughly investigated previous GWAS to pinpoint variants that may play a causal role in poor schizophrenia treatment outcomes, and provides potential candidate genes for further study in the field of antipsychotic response.

Identifying Cell Type-Specific Transcription Factors by Integrating ChIP-seq and eQTL Data-Application to Monocyte Gene Regulation

We describe a novel computational approach to identify transcription factors (TFs) that are candidate regulators in a human cell type of interest. Our approach involves integrating cell type-specific expression quantitative trait locus (eQTL) data and TF data from chromatin immunoprecipitation-to-tag-sequencing (ChIP-seq) experiments in cell lines. To test the method, we used eQTL data from human monocytes in order to screen for TFs. Using a list of known monocyte-regulating TFs, we tested the hypothesis that the binding sites of cell type-specific TF regulators would be concentrated in the vicinity of monocyte eQTLs. For each of 397 ChIP-seq data sets, we obtained an enrichment ratio for the number of ChIP-seq peaks that are located within monocyte eQTLs. We ranked ChIP-seq data sets according to their statistical significances for eQTL overlap, and from this ranking, we observed that monocyte-regulating TFs are more highly ranked than would be expected by chance. We identified 27 TFs that had significant monocyte enrichment scores and mapped them into a protein interaction network. Our analysis uncovered two novel candidate monocyte-regulating TFs, BCLAF1 and SIN3A. Our approach is an efficient method to identify candidate TFs that can be used for any cell/tissue type for which eQTL data are available.

J, Chu AY, Bragg-Gresham JL, van Dongen J, Huffman JE, Ahluwalia TS, Cadby G, Eklund N, Eriksson J, Esko T, Feitosa MF, Goel A, Gorski M, Hayward C, Heard-Costa NL, Jackson AU, Jokinen E, Kanoni S, Kristiansson K, Kutalik Z, Lahti J, Luan J, Mägi R, Mahajan A, Mangino M, Medina-Gomez C, Monda KL, Nolte IM, Pérusse L, Prokopenko I, Qi L, Rose LM, Salvi E, Smith MT, Snieder H, Stančáková A, Ju Sung Y, Tachmazidou I, Teumer A, Thorleifsson G, van der Harst P, Walker RW, Wang SR, Wild SH, Willems SM, Wong A, Zhang W, Albrecht E, Couto Alves A, Bakker SJL, Barlassina C, Bartz TM, Beilby J, Bellis C, Bergman RN, Bergmann S, Blangero J, Blüher M, Boerwinkle E, Bonnycastle LL, Bornstein SR, Bruinenberg M, Campbell H, Chen YDI, Chiang CWK, Chines PS, Collins FS, Cucca F, Cupples LA, D'Avila F, de Geus EJ.C, Dedoussis G, Dimitriou M, Döring A, Eriksson JG, Farmaki AE, Farrall M, Ferreira T, Fischer K, Forouhi NG, Friedrich N, Gjesing AP, Glorioso N, Graff M, Grallert H, Grarup N, Gräßler J, Grewal J, Hamsten A, Harder MN, Hartman CA, Hassinen M, Hastie N, Hattersley AT, Havulinna AS, Heliövaara M, Hillege H, Hofman A, Holmen O, Homuth G, Hottenga JJ, Hui J, Husemoen LL, Hysi PG, Isaacs A, Ittermann T, Jalilzadeh S, James AL, Jørgensen T, Jousilahti P, Jula A, Marie Justesen J, Justice AE, Kähönen M, Karaleftheri M, Tee Khaw K, Keinanen-Kiukaanniemi SM, Kinnunen L, Knekt PB, Koistinen HA, Kolcic I, Kooner IK, Koskinen S, Kovacs P, Kyriakou T, Laitinen T, Langenberg C, Lewin AM, Lichtner P, Lindgren CM, Lindström J, Linneberg A, Lorbeer R, Lorentzon M, Luben R, Lyssenko V, Männistö S, Manunta P, Leach IM, McArdle WL, Mcknight B, Mohlke KL, Mihailov E, Milani L, Mills R, Montasser ME, Morris AP, Müller G, Musk AW, Narisu N, Ong KK, Oostra BA, Osmond C, Palotie A, Pankow JS, Paternoster L, Penninx BW, Pichler I, Pilia MG, Polašek O, Pramstaller PP, Raitakari OT, Rankinen T, Rao DC, Rayner NW, Ribel-Madsen R, Rice TK, Richards M, Ridker PM, Rivadeneira F, Ryan KA, Sanna S, Sarzynski MA, Scholtens S, Scott RA, Sebert S, Southam L, Sparsø TH, Steinthorsdottir V, Stirrups K, Stolk RP, Strauch K, Stringham HM, Swertz MA, Swift AJ, Tönjes A, Tsafantakis E, van der Most PJ, Van Vliet-Ostaptchouk JV, Vandenput L, Vartiainen E, Venturini C, Verweij N, Viikari JS, Vitart V, Vohl MC, Vonk JM, Waeber G, Widén E, Willemsen G, Wilsgaard T, Winkler TW, Wright AF, Yerges-Armstrong LM, Hua Zhao J, Carola Zillikens M, Boomsma DI, Bouchard C, Chambers JC, Chasman DI, Cusi D, Gansevoort RT, Gieger C, Hansen T, Hicks AA, Hu F, Hveem K, Jarvelin MR, Kajantie E, Kooner JS, Kuh D, Kuusisto J, Laakso M, Lakka TA, Lehtimäki T, Metspalu A, Njølstad I, Ohlsson C, Oldehinkel AJ, Palmer LJ, Pedersen O, Perola M, Peters A, Psaty BM, Puolijoki H, Rauramaa R, Rudan I, Salomaa V, Schwarz PEH, Shudiner AR, Smit JH, Sørensen TIA, Spector TD, Stefansson K, Stumvoll M, Tremblay A, Tuomilehto J, Uitterlinden AG, Uusitupa M, Völker U, Vollenweider P, Wareham NJ, Watkins H, Wilson JF, Zeggini E, Abecasis GR, Boehnke M, Borecki IB, Deloukas P, van Duijn CM, Fox C, Groop LC, Heid IM, Hunter DJ, Kaplan RC, McCarthy MI, North KE, O'Connell JR, Schlessinger D, Thorsteinsdottir U, Strachan DP, Frayling T, Hirschhorn JN, Müller-Nurasyid M, Loos RJF. A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

The precision relationships between eight GWAS-identified genetic variants and breast cancer in a Chinese population

Some of the new breast cancer susceptibility loci discovered in recent Genome-wide association studies (GWASs) have not been confirmed in Chinese populations. To determine whether eight novel Single-Nucleotide Polymorphisms (SNPs) have associations with breast cancer risk in women from southeast China, we conducted a case-control study of 1,156 breast cancer patients and 1,256 healthy controls. We first validated that the SNPs rs12922061, rs2290203, and rs2981578 were associated with overall breast cancer risk in southeast Chinese women, with the per-allele OR of 1.209 (95%CI: 1.064-1.372), 1.176 (95%CI: 1.048-1.320), and 0.852 (95%CI: 0.759-0.956), respectively. Rs12922061 and rs2290203 even passed the threshold for Bonferroni correction (P value: 0.00625). In stratified analysis, we found another three SNPs were significantly associated within different subgroups. However, after Bonferroni correction (P value: 0.000446), there were no statistically significant was observed. In gene-environment interaction analysis, we observed gene-environment interactions played a potential role of in the risk of breast cancer. These findings provide new insight into the associations between the genetic susceptibility and fine classifications of breast cancer. Based on these results, we encourage further large series studies and functional research to confirm these finding.

Targeted high-throughput sequencing of candidate genes for chronic obstructive pulmonary disease

Background

Reduced lung function in patients with chronic obstructive pulmonary disease (COPD) is likely due to both environmental and genetic factors. We report here a targeted high-throughput DNA sequencing approach to identify new and previously known genetic variants in a set of candidate genes for COPD.

Methods

Exons in 22 genes implicated in lung development as well as 61 genes and 10 genomic regions previously associated with COPD were sequenced using individual DNA samples from 68 cases with moderate or severe COPD and 66 controls matched for age, gender and smoking. Cases and controls were selected from the Obstructive Lung Disease in Northern Sweden (OLIN) studies.

Results

In total, 37 genetic variants showed association with COPD (p < 0.05, uncorrected). Several variants previously discovered to be associated with COPD from genetic genome-wide analysis studies were replicated using our sample. Two high-risk variants were followed-up for functional characterization in a large eQTL mapping study of 1,111 human lung specimens. The C allele of a synonymous variant, rs8040868, predicting a p.(S45=) in the gene for cholinergic receptor nicotinic alpha 3 (CHRNA3) was associated with COPD (p = 8.8 x 10⁻³). This association remained (p = 0.003 and OR = 1.4, 95 % CI 1.1-1.7) when analysing all available cases and controls in OLIN (n = 1,534). The rs8040868 variant is in linkage disequilibrium with rs16969968 previously associated with COPD and altered expression of the CHRNA5 gene. A follow-up analysis for detection of expression quantitative trait loci revealed that rs8040868-C was found to be significantly associated with a decreased expression of the nearby gene cholinergic receptor, nicotinic, alpha 5 (CHRNA5) in lung tissue.

Conclusion

Our data replicate previous result suggesting CHRNA5 as a candidate gene for COPD and rs8040868 as a risk variant for the development of COPD in the Swedish population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12890-016-0309-y) contains supplementary material, which is available to authorized users.

Human cells involved in atherosclerosis have a sex

The influence of sex has been largely described in cardiovascular diseases. Atherosclerosis is a complex process that involves many cell types such as vessel cells, immune cells and endothelial progenitor cells; however, many, if not all, studies do not report the sex of the cells. This review focuses on sex differences in human cells involved in the atherosclerotic process, emphasizing the role of sex hormones. Furthermore, we report sex differences and issues related to the processes that determine the fate of the cells such as apoptotic and autophagic mechanisms. The analysis of the data reveals that there are still many gaps in our knowledge regarding sex influences in atherosclerosis, largely for the cell types that have not been well studied, stressing the urgent need for a clear definition of experimental conditions and the inclusion of both sexes in preclinical studies.

Genome-Wide Interaction Analyses between Genetic Variants and Alcohol Consumption and Smoking for Risk of Colorectal Cancer

Genome-wide association studies (GWAS) have identified many genetic susceptibility loci for colorectal cancer (CRC). However, variants in these loci explain only a small proportion of familial aggregation, and there are likely additional variants that are associated with CRC susceptibility. Genome-wide studies of gene-environment interactions may identify variants that are not detected in GWAS of marginal gene effects. To study this, we conducted a genome-wide analysis for interaction between genetic variants and alcohol consumption and cigarette smoking using data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Interactions were tested using logistic regression. We identified interaction between CRC risk and alcohol consumption and variants in the 9q22.32/HIATL1 (P_interaction = 1.76×10⁻⁸; permuted p-value 3.51x10^-8) region. Compared to non-/occasional drinking light to moderate alcohol consumption was associated with a lower risk of colorectal cancer among individuals with rs9409565 CT genotype (OR, 0.82 [95% CI, 0.74–0.91]; P = 2.1×10⁻⁴) and TT genotypes (OR,0.62 [95% CI, 0.51–0.75]; P = 1.3×10⁻⁶) but not associated among those with the CC genotype (p = 0.059). No genome-wide statistically significant interactions were observed for smoking. If replicated our suggestive finding of a genome-wide significant interaction between genetic variants and alcohol consumption might contribute to understanding colorectal cancer etiology and identifying subpopulations with differential susceptibility to the effect of alcohol on CRC risk.

Alcohol consumption and smoking are associated with CRC risk. We performed a genome-wide analysis for interaction between genetic variants and alcohol consumption and cigarette smoking to identify potential new genetic regions associated with CRC. About 8,000 CRC cases and 8,800 controls were included in alcohol-related analysis and over 11,000 cases and 11,000 controls were involved in smoking-related analysis. We identified interaction between variants at 9q22.32/HIATL1 and alcohol consumption in relation to CRC risk (P_interaction = 1.76×10⁻⁸). If replicated our suggested finding of the interaction between genetic variants and alcohol consumption might contribute to understanding colorectal cancer etiology and identifying subpopulations with differential susceptible to the effect of alcohol on CRC risk.

Meta-analysis identifies common and rare variants influencing blood pressure and overlapping with metabolic trait loci

Meta-analyses of association results for blood pressure using exome-centric single-variant and gene-based tests identified 31 new loci in a discovery stage among 146,562 individuals, with follow-up and meta-analysis in 180,726 additional individuals (total n = 327,288). These blood pressure-associated loci are enriched for known variants for cardiometabolic traits. Associations were also observed for the aggregation of rare and low-frequency missense variants in three genes, NPR1, DBH, and PTPMT1. In addition, blood pressure associations at 39 previously reported loci were confirmed. The identified variants implicate biological pathways related to cardiometabolic traits, vascular function, and development. Several new variants are inferred to have roles in transcription or as hubs in protein-protein interaction networks. Genetic risk scores constructed from the identified variants were strongly associated with coronary disease and myocardial infarction. This large collection of blood pressure-associated loci suggests new therapeutic strategies for hypertension, emphasizing a link with cardiometabolic risk.

FAM26F: An Enigmatic Protein Having a Complex Role in the Immune System

Mammalian immune system is a complex amalgam of diverse cellular and noncellular components such as cytokines, receptors and co-receptors. FAM26F (family with sequence similarity 26, member F) is a recently identified tetraspanin-like membrane glycoprotein which is predicted to make homophilic interactions and potential synapses between several immune cells including CD4⁺, CD8⁺, NK, dendritic cells and macrophages. Various whole transcriptome analyses have demonstrated the differential expression of FAM26F in several bacterial, viral and parasitic infections, in certain pathophysiological conditions such as liver and heart transplantation, and in various cancers. The complete understanding of transcriptional regulation of FAM26F is in its infancy however it is up regulated by various stimulants such as polyI:C, LPS, INF gamma and TNF alpha, and via various proposed pathways including TLR3, TLR4 IFN-β and Dectin-1. These pathways can merge in STAT1 activation. The synergistic expression of FAM26F on both NK-cells and myeloid dendritic cells is required to activate NK-cells against tumors via its cytoplasmic tail, thus emphasizing therapeutic potential of FAM26F for NK sensitive tumors. Current review provides a comprehensive basis to propose that FAM26F expression level is at least a hallmark for IFN-γ-lead immune responses and thus can proficiently be regarded as an early diagnostic marker. Future investigation dissecting the role of FAM26F in activation of various immune cell populations in local amplification by cell-cell contact is crucial to provide the missing link imperative for elucidating the relevance of this protein in immune responses.

Recommendations for genetic testing to reduce the incidence of anthracycline-induced cardiotoxicity

AIMS

Anthracycline-induced cardiotoxicity (ACT) occurs in 57% of treated patients and remains an important limitation of anthracycline-based chemotherapy. In various genetic association studies, potential genetic risk markers for ACT have been identified. Therefore, we developed evidence-based clinical practice recommendations for pharmacogenomic testing to further individualize therapy based on ACT risk.

METHODS

We followed a standard guideline development process, including a systematic literature search, evidence synthesis and critical appraisal, and the development of clinical practice recommendations with an international expert group.

RESULTS

RARG rs2229774, SLC28A3 rs7853758 and UGT1A6 rs17863783 variants currently have the strongest and the most consistent evidence for association with ACT. Genetic variants in ABCC1, ABCC2, ABCC5, ABCB1, ABCB4, CBR3, RAC2, NCF4, CYBA, GSTP1, CAT, SULT2B1, POR, HAS3, SLC22A7, SCL22A17, HFE and NOS3 have also been associated with ACT, but require additional validation. We recommend pharmacogenomic testing for the RARG rs2229774 (S427L), SLC28A3 rs7853758 (L461L) and UGT1A6*4 rs17863783 (V209V) variants in childhood cancer patients with an indication for doxorubicin or daunorubicin therapy (Level B - moderate). Based on an overall risk stratification, taking into account genetic and clinical risk factors, we recommend a number of management options including increased frequency of echocardiogram monitoring, follow-up, as well as therapeutic options within the current standard of clinical practice.

CONCLUSIONS

Existing evidence demonstrates that genetic factors have the potential to improve the discrimination between individuals at higher and lower risk of ACT. Genetic testing may therefore support both patient care decisions and evidence development for an improved prevention of ACT.

2014 Jeffrey M. Hoeg Award Lecture: Transcriptional Control of Monocyte Development

Monocytes and macrophages are key immune cells involved in the early progression of atherosclerosis. Transcription factors that control their development in the bone marrow are important therapeutic targets to control the numbers and functions of these cells in disease. This review highlights what is currently known about the transcription factors that are critical for monocyte development.

Screening and identification of microRNA involved in unstable angina using gene-chip analysis

Increasing evidence has suggested that microRNA (miRNA) may play a role in the pathogenesis of cardiovascular disease, which has led to a greater understanding of the complex pathophysiological processes underlying unstable angina (UA). The present study aimed to investigate changes in the miRNA expression profiles of patients with UA using gene-chip analysis, in order to further elucidate the pathogenesis of UA. Total RNA was extracted and purified from plasma samples collected from patients with UA and healthy controls. The samples underwent microarray analysis using an Exiqon miRCURY LNA™ microRNA Array. Differentially expressed miRNAs were identified by volcano plot filtering, and were validated using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). In addition, functional annotation of the differentially expressed miRNAs involved gene ontology analyses. Among the 212 miRNAs differentially expressed between the two groups, 82 were upregulated and 130 were downregulated. Notably, the results of the RT-qPCR were consistent with the gene-chip results. The miRNAs identified in the present study may be potential novel biomarkers for the prevention and early diagnosis of UA. Furthermore, the results of the present study suggested that UA occurs as a result of complex and dynamic processes regulated by numerous factors, including multiple miRNAs.

Gene-based analysis of regulatory variants identifies 4 putative novel asthma risk genes related to nucleotide synthesis and signaling

BACKGROUND

Hundreds of genetic variants are thought to contribute to variation in asthma risk by modulating gene expression. Methods that increase the power of genome-wide association studies (GWASs) to identify risk-associated variants are needed.

OBJECTIVE

We sought to develop a method that aggregates the evidence for association with disease risk across expression quantitative trait loci (eQTLs) of a gene and use this approach to identify asthma risk genes.

METHODS

We developed a gene-based test and software package called EUGENE that (1) is applicable to GWAS summary statistics; (2) considers both cis- and trans-eQTLs; (3) incorporates eQTLs identified in different tissues; and (4) uses simulations to account for multiple testing. We applied this approach to 2 published asthma GWASs (combined n = 46,044) and used mouse studies to provide initial functional insights into 2 genes with novel genetic associations.

RESULTS

We tested the association between asthma and 17,190 genes that were found to have cis- and/or trans-eQTLs across 16 published eQTL studies. At an empirical FDR of 5%, 48 genes were associated with asthma risk. Of these, for 37, the association was driven by eQTLs located in established risk loci for allergic disease, including 6 genes not previously implicated in disease cause (eg, LIMS1, TINF2, and SAFB). The remaining 11 significant genes represent potential novel genetic associations with asthma. The association with 4 of these replicated in an independent GWAS: B4GALT3, USMG5, P2RY13, and P2RY14, which are genes involved in nucleotide synthesis or nucleotide-dependent cell activation. In mouse studies, P2ry13 and P2ry14-purinergic receptors activated by adenosine 5-diphosphate and UDP-sugars, respectively-were upregulated after allergen challenge, notably in airway epithelial cells, eosinophils, and neutrophils. Intranasal exposure with receptor agonists induced the release of IL-33 and subsequent eosinophil infiltration into the lungs.

CONCLUSION

We identified novel associations between asthma and eQTLs for 4 genes related to nucleotide synthesis/signaling and demonstrated the power of gene-based analyses of GWASs.

Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies

BACKGROUND & AIMS

A genome-wide association study (GWAS) of 280 cases identified the hepatic cholesterol transporter ABCG8 as a locus associated with risk for gallstone disease, but findings have not been reported from any other GWAS of this phenotype. We performed a large-scale, meta-analysis of GWASs of individuals of European ancestry with available prior genotype data, to identify additional genetic risk factors for gallstone disease.

METHODS

We obtained per-allele odds ratio (OR) and standard error estimates using age- and sex-adjusted logistic regression models within each of the 10 discovery studies (8720 cases and 55,152 controls). We performed an inverse variance weighted, fixed-effects meta-analysis of study-specific estimates to identify single-nucleotide polymorphisms that were associated independently with gallstone disease. Associations were replicated in 6489 cases and 62,797 controls.

RESULTS

We observed independent associations for 2 single-nucleotide polymorphisms at the ABCG8 locus: rs11887534 (OR, 1.69; 95% confidence interval [CI], 1.54-1.86; P = 2.44 × 10(-60)) and rs4245791 (OR, 1.27; P = 1.90 × 10(-34)). We also identified and/or replicated associations for rs9843304 in TM4SF4 (OR, 1.12; 95% CI, 1.08-1.16; P = 6.09 × 10(-11)), rs2547231 in SULT2A1 (encodes a sulfoconjugation enzyme that acts on hydroxysteroids and cholesterol-derived sterol bile acids) (OR, 1.17; 95% CI, 1.12-1.21; P = 2.24 × 10(-10)), rs1260326 in glucokinase regulatory protein (OR, 1.12; 95% CI, 1.07-1.17; P = 2.55 × 10(-10)), and rs6471717 near CYP7A1 (encodes an enzyme that catalyzes conversion of cholesterol to primary bile acids) (OR, 1.11; 95% CI, 1.08-1.15; P = 8.84 × 10(-9)). Among individuals of African American and Hispanic American ancestry, rs11887534 and rs4245791 were associated positively with gallstone disease risk, whereas the association for the rs1260326 variant was inverse.

CONCLUSIONS

In this large-scale GWAS of gallstone disease, we identified 4 loci in genes that have putative functions in cholesterol metabolism and transport, and sulfonylation of bile acids or hydroxysteroids.

Targeted genomic analysis reveals widespread autoimmune disease association with regulatory variants in the TNF superfamily cytokine signalling network

BACKGROUND

Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored.

METHODS

We took a hypothesis-driven approach using available genome-wide datasets to identify genetic variants regulating gene expression in the TNF superfamily cytokine signalling network and the association of these variants with autoimmune and autoinflammatory disease. Using paired gene expression and genetic data, we identified genetic variants associated with gene expression, expression quantitative trait loci (eQTLs), in four peripheral blood cell subsets. We then examined whether eQTLs were dependent on gene expression level or the presence of active enhancer chromatin marks. Using these eQTLs as genetic markers of the TNF superfamily signalling network, we performed targeted gene set association analysis in eight autoimmune and autoinflammatory disease genome-wide association studies.

RESULTS

Comparison of TNF superfamily network gene expression and regulatory variants across four leucocyte subsets revealed patterns that differed between cell types. eQTLs for genes in this network were not dependent on absolute gene expression levels and were not enriched for chromatin marks of active enhancers. By examining autoimmune disease risk variants among our eQTLs, we found that risk alleles can be associated with either increased or decreased expression of co-stimulatory TNF superfamily cytokines, receptors or downstream signalling molecules. Gene set disease association analysis revealed that eQTLs for genes in the TNF superfamily pathway were associated with six of the eight autoimmune and autoinflammatory diseases examined, demonstrating associations beyond single genome-wide significant hits.

CONCLUSIONS

This systematic analysis of the influence of regulatory genetic variants in the TNF superfamily network reveals widespread and diverse roles for these cytokines in susceptibility to a number of immune-mediated diseases.

RNA-Seq of Huntington's disease patient myeloid cells reveals innate transcriptional dysregulation associated with proinflammatory pathway activation

Innate immune activation beyond the central nervous system is emerging as a vital component of the pathogenesis of neurodegeneration. Huntington's disease (HD) is a fatal neurodegenerative disorder caused by a CAG repeat expansion in the huntingtin gene. The systemic innate immune system is thought to act as a modifier of disease progression; however, the molecular mechanisms remain only partially understood. Here we use RNA-sequencing to perform whole transcriptome analysis of primary monocytes from thirty manifest HD patients and thirty-three control subjects, cultured with and without a proinflammatory stimulus. In contrast with previous studies that have required stimulation to elicit phenotypic abnormalities, we demonstrate significant transcriptional differences in HD monocytes in their basal, unstimulated state. This includes previously undetected increased resting expression of genes encoding numerous proinflammatory cytokines, such as IL6 Further pathway analysis revealed widespread resting enrichment of proinflammatory functional gene sets, while upstream regulator analysis coupled with Western blotting suggests that abnormal basal activation of the NFĸB pathway plays a key role in mediating these transcriptional changes. That HD myeloid cells have a proinflammatory phenotype in the absence of stimulation is consistent with a priming effect of mutant huntingtin, whereby basal dysfunction leads to an exaggerated inflammatory response once a stimulus is encountered. These data advance our understanding of mutant huntingtin pathogenesis, establish resting myeloid cells as a key source of HD immune dysfunction, and further demonstrate the importance of systemic immunity in the potential treatment of HD and the wider study of neurodegeneration.

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.

Identification of STOML2 as a putative novel asthma risk gene associated with IL6R

BACKGROUND

Functional variants in the interleukin-6 receptor gene (IL6R) are associated with asthma risk. We hypothesized that genes co-expressed with IL6R might also be regulated by genetic polymorphisms that are associated with asthma risk. The aim of this study was to identify such genes.

METHODS

To identify genes whose expression was correlated with that of IL6R, we analyzed gene expression levels generated for 373 human lymphoblastoid cell lines by the Geuvadis consortium and for 38 hematopoietic cell types by the Differentiation Map Portal (DMAP) project. Genes correlated with IL6R were then screened for nearby single nucleotide polymorphisms (SNPs) that were significantly associated with both variation in gene expression levels (eSNPs) and asthma risk.

RESULTS

We identified 90 genes with expression levels correlated with those of IL6R and that also had a nearby eSNP associated with disease risk in a published asthma GWAS (N = 20 776). For 16 (18%) genes, the association between the eSNP and asthma risk replicated with the same direction of effect in a further independent published asthma GWAS (N = 27 378). Among the top replicated associations (FDR < 0.05) were eSNPs for four known (IL18R1, IL18RAP, BCL6, and STAT6) and one putative novel asthma risk gene, stomatin-like protein 2 (STOML2). The expression of STOML2 was negatively correlated with IL6R, while eSNPs that increased the expression of STOML2 were associated with an increased asthma risk.

CONCLUSION

The expression of STOML2, a gene that plays a key role in mitochondrial function and T-cell activation, is associated with both IL-6 signaling and asthma risk.

From Loci to Biology: Functional Genomics of Genome-Wide Association for Coronary Disease

Genome-wide association studies have provided a rich collection of ≈ 58 coronary artery disease (CAD) loci that suggest the existence of previously unsuspected new biology relevant to atherosclerosis. However, these studies only identify genomic loci associated with CAD, and many questions remain even after a genomic locus is definitively implicated, including the nature of the causal variant(s) and the causal gene(s), as well as the directionality of effect. There are several tools that can be used for investigation of the functional genomics of these loci, and progress has been made on a limited number of novel CAD loci. New biology regarding atherosclerosis and CAD will be learned through the functional genomics of these loci, and the hope is that at least some of these new pathways relevant to CAD pathogenesis will yield new therapeutic targets for the prevention and treatment of CAD.

Genome-wide association study in East Asians identifies two novel breast cancer susceptibility loci

Breast cancer is one of the most common malignancies among women worldwide. Genetic factors have been shown to play an important role in breast cancer aetiology. We conducted a two-stage genome-wide association study (GWAS) including 14 224 cases and 14 829 controls of East Asian women to search for novel genetic susceptibility loci for breast cancer. Single nucleotide polymorphisms (SNPs) in two loci were found to be associated with breast cancer risk at the genome-wide significance level. The first locus, represented by rs12118297 at 1p22.3 (near the LMO4 gene), was associated with breast cancer risk with odds ratio (OR) and (95% confidence interval (CI)) of 0.91 (0.88-0.94) and a P-value of 4.48 × 10^{-  8} This association was replicated in another study, DRIVE GAME-ON Consortium, including 16 003 cases and 41 335 controls of European ancestry (OR = 0.95, 95% CI = 0.91-0.99, P-value = 0.019). The second locus, rs16992204 at 21q22.12 (near the LINC00160 gene), was associated with breast cancer risk with OR (95% CI) of 1.13 (1.07-1.18) and a P-value of 4.63 × 10 ^{-  8} The risk allele frequency for this SNP is zero in European-ancestry populations in 1000 Genomes Project and thus its association with breast cancer risk cannot be assessed in DRIVE GAME-ON Consortium. Functional annotation using the ENCODE data indicates that rs12118297 might be located in a repressed element and locus 21q22.12 may affect breast cancer risk through regulating LINC00160 expressions and interaction with oestrogen receptor signalling. Our findings provide additional insights into the genetics of breast cancer.

Comprehensive promoter level expression quantitative trait loci analysis of the human frontal lobe

BACKGROUND

Expression quantitative trait loci (eQTL) analysis is a powerful method to detect correlations between gene expression and genomic variants and is widely used to interpret the biological mechanism underlying identified genome wide association studies (GWAS) risk loci. Numerous eQTL studies have been performed on different cell types and tissues of which the majority has been based on microarray technology.

METHODS

We present here an eQTL analysis based on cap analysis gene expression sequencing (CAGEseq) data created from human postmortem frontal lobe tissue combined with genotypes obtained through genotyping arrays, exome sequencing, and CAGEseq. Using CAGEseq as an expression profiling technique combined with these different genotyping techniques allows measurement of the molecular effect of variants on individual transcription start sites and increases the resolution of eQTL analysis by also including the non-annotated parts of the genome.

RESULTS

We identified 2410 eQTLs and show that non-coding transcripts are more likely to contain an eQTL than coding transcripts, in particular antisense transcripts. We provide evidence for how previously identified GWAS loci for schizophrenia (NRGN), Parkinson's disease, and Alzheimer's disease (PARK16 and MAPT loci) could increase the risk for disease at a molecular level. Furthermore, we demonstrate that CAGEseq improves eQTL analysis because variants obtained from CAGEseq are highly enriched for having a functional effect and thus are an efficient method towards the identification of causal variants.

CONCLUSION

Our data contain both coding and non-coding transcripts and has the added value that we have identified eQTLs for variants directly adjacent to TSS. Future eQTL studies would benefit from combining CAGEseq with RNA sequencing for a more complete interpretation of the transcriptome and increased understanding of eQTL signals.

Endophenotype Network Models: Common Core of Complex Diseases

Historically, human diseases have been differentiated and categorized based on the organ system in which they primarily manifest. Recently, an alternative view is emerging that emphasizes that different diseases often have common underlying mechanisms and shared intermediate pathophenotypes, or endo(pheno)types. Within this framework, a specific disease’s expression is a consequence of the interplay between the relevant endophenotypes and their local, organ-based environment. Important examples of such endophenotypes are inflammation, fibrosis, and thrombosis and their essential roles in many developing diseases. In this study, we construct endophenotype network models and explore their relation to different diseases in general and to cardiovascular diseases in particular. We identify the local neighborhoods (module) within the interconnected map of molecular components, i.e., the subnetworks of the human interactome that represent the inflammasome, thrombosome, and fibrosome. We find that these neighborhoods are highly overlapping and significantly enriched with disease-associated genes. In particular they are also enriched with differentially expressed genes linked to cardiovascular disease (risk). Finally, using proteomic data, we explore how macrophage activation contributes to our understanding of inflammatory processes and responses. The results of our analysis show that inflammatory responses initiate from within the cross-talk of the three identified endophenotypic modules.

Removing Batch Effects from Longitudinal Gene Expression - Quantile Normalization Plus ComBat as Best Approach for Microarray Transcriptome Data

Technical variation plays an important role in microarray-based gene expression studies, and batch effects explain a large proportion of this noise. It is therefore mandatory to eliminate technical variation while maintaining biological variability. Several strategies have been proposed for the removal of batch effects, although they have not been evaluated in large-scale longitudinal gene expression data. In this study, we aimed at identifying a suitable method for batch effect removal in a large study of microarray-based longitudinal gene expression. Monocytic gene expression was measured in 1092 participants of the Gutenberg Health Study at baseline and 5-year follow up. Replicates of selected samples were measured at both time points to identify technical variability. Deming regression, Passing-Bablok regression, linear mixed models, non-linear models as well as ReplicateRUV and ComBat were applied to eliminate batch effects between replicates. In a second step, quantile normalization prior to batch effect correction was performed for each method. Technical variation between batches was evaluated by principal component analysis. Associations between body mass index and transcriptomes were calculated before and after batch removal. Results from association analyses were compared to evaluate maintenance of biological variability. Quantile normalization, separately performed in each batch, combined with ComBat successfully reduced batch effects and maintained biological variability. ReplicateRUV performed perfectly in the replicate data subset of the study, but failed when applied to all samples. All other methods did not substantially reduce batch effects in the replicate data subset. Quantile normalization plus ComBat appears to be a valuable approach for batch correction in longitudinal gene expression data.

Association of MHC region SNPs with irritant susceptibility in healthcare workers

Irritant contact dermatitis is the most common work-related skin disease, especially affecting workers in "wet-work" occupations. This study was conducted to investigate the association between single nucleotide polymorphisms (SNPs) within the major histocompatibility complex (MHC) and skin irritant response in a group of healthcare workers. 585 volunteer healthcare workers were genotyped for MHC SNPs and patch tested with three different irritants: sodium lauryl sulfate (SLS), sodium hydroxide (NaOH) and benzalkonium chloride (BKC). Genotyping was performed using Illumina Goldengate MHC panels. A number of SNPs within the MHC Class I (OR2B3, TRIM31, TRIM10, TRIM40 and IER3), Class II (HLA-DPA1, HLA-DPB1) and Class III (C2) genes were associated (p < 0.001) with skin response to tested irritants in different genetic models. Linkage disequilibrium patterns and functional annotations identified two SNPs in the TRIM40 (rs1573298) and HLA-DPB1 (rs9277554) genes, with a potential impact on gene regulation. In addition, SNPs in PSMB9 (rs10046277 and ITPR3 (rs499384) were associated with hand dermatitis. The results are of interest as they demonstrate that genetic variations in inflammation-related genes within the MHC can influence chemical-induced skin irritation and may explain the connection between inflamed skin and propensity to subsequent allergic contact sensitization.

Genetic Basis of Irritant Susceptibility in Health Care Workers

OBJECTIVE

The aim of this study was to investigate the association of single nucleotide polymorphisms (SNPs) within genes involved in inflammation, skin barrier integrity, signaling/pattern recognition, and antioxidant defense with irritant susceptibility in a group of health care workers.

METHODS

The 536 volunteer subjects were genotyped for selected SNPs and patch tested with three model irritants: sodium lauryl sulfate (SLS), sodium hydroxide (NaOH), and benzalkonium chloride (BKC). Genotyping was performed on genomic DNA using Illumina Goldengate custom panels.

RESULTS

The ACACB (rs2268387, rs16934132, rs2284685), NTRK2 (rs10868231), NTRK3 (rs1347424), IL22 (rs1179251), PLAU (rs2227564), EGFR (rs6593202), and FGF2 (rs308439) SNPs showed an association with skin response to tested irritants in different genetic models (all at P < 0.001). Functional annotations identified two SNPs in PLAU (rs2227564) and ACACB (rs2284685) genes with a potential impact on gene regulation. In addition, EGF (rs10029654), EGFR (rs12718939), CXCL12 (rs197452), and VCAM1 (rs3917018) genes showed an association with hand dermatitis (P < 0.005).

CONCLUSIONS

The results demonstrate that genetic variations in genes related to inflammation and skin homeostasis can influence responses to irritants and may explain inter-individual variation in the development of subsequent contact dermatitis.

SUGP1 is a novel regulator of cholesterol metabolism

A large haplotype on chromosome 19p13.11 tagged by rs10401969 in intron 8 of SURP and G patch domain containing 1 (SUGP1) is associated with coronary artery disease (CAD), plasma LDL cholesterol levels, and other energy metabolism phenotypes. Recent studies have suggested that TM6SF2 is the causal gene within the locus, but we postulated that this locus could harbor additional CAD risk genes, including the putative splicing factor SUGP1. Indeed, we found that rs10401969 regulates SUGP1 exon 8 skipping, causing non-sense-mediated mRNA decay. Hepatic Sugp1 overexpression in CD1 male mice increased plasma cholesterol levels 20–50%. In human hepatoma cell lines, SUGP1 knockdown stimulated 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) alternative splicing and decreased HMGCR transcript stability, thus reducing cholesterol synthesis and increasing LDL uptake. Our results strongly support a role for SUGP1 as a novel regulator of cholesterol metabolism and suggest that it contributes to the relationship between rs10401969 and plasma cholesterol.

Multivariate eQTL mapping uncovers functional variation on the X-chromosome associated with complex disease traits

Very few studies have investigated the associations between genetic polymorphisms and gene expression on the X-chromosome. This is a major bottleneck when conducting functional follow-up studies of trait-associated variants, as those identified in genome-wide association studies (GWAS). We used a multivariate approach to test the association between individual single nucleotide polymorphisms (SNPs) and exon expression levels measured in 356 Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) from the Geuvadis RNA sequencing project to identify SNPs associated with variation in gene expression on the X-chromosome, which we refer to as eSNPs. At an FDR of 5 %, we discovered 548 independent [linkage disequilibrium (LD) r (2) < 0.1] eSNPs on the X-chromosome. Of these, 35 were in LD (r (2) > 0.8) with previously published disease- or trait-associated variants identified through GWAS. One of the strongest eSNPs identified was rs35975601, which was associated with F8A1 expression (p value = 3 × 10(-20)) and was in LD with a type 1 diabetes risk variant. Additionally, we identified a number of genes for which eSNPs were in LD with multiple diseases or traits, including DNASE1L1 which was mapped to bilirubin levels, type 1 diabetes and schizophrenia. Our results also indicate that multivariate exon-level analysis provides a more powerful approach than univariate gene-level analysis, particularly when SNPs influence the expression of different exons with different magnitude and/or direction of effect. The associations identified in our study may provide new insights into the molecular process by which gene expression may contribute to trait variation or disease risk in humans.

Fine mapping under linkage peaks for symptomatic or asymptomatic outcomes of Leishmania infantum infection in Brazil

Infection with the protozoan Leishmania infantum can lead to asymptomatic infection and protective immunity, or to the progressive and potentially fatal disease visceral leishmaniasis (VL). Published studies show host genetic background determines in part whether infected individuals will develop a symptomatic or asymptomatic outcome. The purpose of the current study was to fine map chromosome regions previously linked with risk for symptomatic (chromosome 9) or asymptomatic (chromosomes 15 and 19) manifestations of L. infantum infection. We conducted a family-based genetic study of VL and asymptomatic infection (detected by a DTH skin test) with a final post quality control sample of 961 individuals with full genotype and phenotype information from highly endemic neighborhoods of northeast Brazil. A total of 5485 SNPs under the linkage peaks on chromosomes 9, 15 and 19 were genotyped. No strong SNP associations were observed for the DTH phenotype. The most significant associations with the VL phenotype were with SNP rs1470217 (p=5.9e-05; pcorrected=0.057) on chromosome 9, and with SNP rs8107014 (p=1.4e-05; pcorrected=0.013) on chromosome 19. SNP rs1470217 is situated in a 180kb intergenic region between TMEM215 (Transmembrane protein 215) and APTX (Aprataxin). SNP rs8107014 lies in the intron between exons 26 and 27 of a 34 exon transcript (ENST00000204005) of LTBP4, (Latent transforming growth factor-beta-binding protein 4a). The latter supports growing evidence that the transforming growth factor-beta pathway is important in the immunopathogenesis of VL.

Identification of Susceptibility Genes of Adult Asthma in French Canadian Women

Susceptibility genes of asthma may be more successfully identified by studying subgroups of phenotypically similar asthma patients. This study aims to identify single nucleotide polymorphisms (SNPs) associated with asthma in French Canadian adult women. A pooling-based genome-wide association study was performed in 240 allergic asthmatic and 120 allergic nonasthmatic women. The top associated SNPs were selected for individual genotyping in an extended cohort of 349 asthmatic and 261 nonasthmatic women. The functional impact of asthma-associated SNPs was investigated in a lung expression quantitative trait loci (eQTL) mapping study (n = 1035). Twenty-one of the 38 SNPs tested by individual genotyping showed P values lower than 0.05 for association with asthma. Cis-eQTL analyses supported the functional contribution of rs17801353 associated with C3AR1 (P = 7.90E - 10). The asthma risk allele for rs17801353 is associated with higher mRNA expression levels of C3AR1 in lung tissue. In silico functional characterization of the asthma-associated SNPs also supported the contribution of C3AR1 and additional genes including SYNE1, LINGO2, and IFNG-AS1. This pooling-based GWAS in French Canadian adult women followed by lung eQTL mapping suggested C3AR1 as a functional locus associated with asthma. Additional susceptibility genes were suggested in this homogenous subgroup of asthma patients.

Genetic and environmental risk factors for atherosclerosis regulate transcription of phosphatase and actin regulating gene PHACTR1

Background and aims

Coronary artery disease (CAD) risk is associated with non-coding genetic variants at the phosphatase and actin regulating protein 1(PHACTR1) gene locus. The PHACTR1 gene encodes an actin-binding protein with phosphatase regulating activity. The mechanism whereby PHACTR1 influences CAD risk is unknown. We hypothesized that PHACTR1 would be expressed in human cell types relevant to CAD and regulated by atherogenic or genetic factors.

Methods and results

Using immunohistochemistry, we demonstrate that PHACTR1 protein is expressed strongly in human atherosclerotic plaque macrophages, lipid-laden foam cells, adventitial lymphocytes and endothelial cells. Using a combination of genomic analysis and molecular techniques, we demonstrate that PHACTR1 is expressed as multiple previously uncharacterized transcripts in macrophages, foam cells, lymphocytes and endothelial cells. Immunoblotting confirmed a total absence of PHACTR1 in vascular smooth muscle cells. Real-time quantitative PCR showed that PHACTR1 is regulated by atherogenic and inflammatory stimuli. In aortic endothelial cells, oxLDL and TNF-alpha both upregulated an intermediate length transcript. A short transcript expressed only in immune cells was upregulated in macrophages by oxidized low-density lipoprotein, and oxidized phospholipids but suppressed by lipopolysaccharide or TNF-alpha. In primary human macrophages, we identified a novel expression quantitative trait locus (eQTL) specific for this short transcript, whereby the risk allele at CAD risk SNP rs9349379 is associated with reduced PHACTR1 expression, similar to the effect of an inflammatory stimulus.

Conclusions

Our data demonstrate that PHACTR1 is a key atherosclerosis candidate gene since it is regulated by atherogenic stimuli in macrophages and endothelial cells and we identify an effect of the genetic risk variant on PHACTR1 expression in macrophages that is similar to that of an inflammatory stimulus.

•

PHACTR1 is expressed as two transcripts in both immune and endothelial cells in human atherosclerotic plaque.

•

Oxidized-LDL upregulates a short PHACTR1 transcript, but suppresses an intermediate length transcript in macrophages.

•

Lipopolysaccharide and TNF-alpha cause the opposite effect with strong suppression of the short transcript in macrophages.

•

The coronary artery disease risk SNP, rs9349379, is associated with expression of the short transcript in macrophages.

•

The effect of the coronary artery disease risk allele on PHACTR1 mirrors that of inflammatory stimuli.

Characterization of Expression Quantitative Trait Loci in Pedigrees from Colombia and Costa Rica Ascertained for Bipolar Disorder

The observation that variants regulating gene expression (expression quantitative trait loci, eQTL) are at a high frequency among SNPs associated with complex traits has made the genome-wide characterization of gene expression an important tool in genetic mapping studies of such traits. As part of a study to identify genetic loci contributing to bipolar disorder and other quantitative traits in members of 26 pedigrees from Costa Rica and Colombia, we measured gene expression in lymphoblastoid cell lines derived from 786 pedigree members. The study design enabled us to comprehensively reconstruct the genetic regulatory network in these families, provide estimates of heritability, identify eQTL, evaluate missing heritability for the eQTL, and quantify the number of different alleles contributing to any given locus. In the eQTL analysis, we utilize a recently proposed hierarchical multiple testing strategy which controls error rates regarding the discovery of functional variants. Our results elucidate the heritability and regulation of gene expression in this unique Latin American study population and identify a set of regulatory SNPs which may be relevant in future investigations of complex disease in this population. Since our subjects belong to extended families, we are able to compare traditional kinship-based estimates with those from more recent methods that depend only on genotype information.

Antiphospholipid antibodies in a large population-based cohort: genome-wide associations and effects on monocyte gene expression

The antiphospholipid syndrome (APS) is characterised by venous and/or arterial thrombosis and pregnancy morbidity in women combined with the persistent presence of antiphospholipid antibodies (aPL). We aimed to identify genetic factors associated with the presence of aPL in a population based cohort. Furthermore, we wanted to clarify if the presence of aPL affects gene expression in circulating monocytes. Titres of IgG and IgM against cardiolipin, β2glycoprotein 1 (anti-β2GPI), and IgG against domain 1 of β2GPI (anti-domain 1) were determined in approx. 5,000 individuals from the Gutenberg Health Study (GHS) a population based cohort of German descent. Genotyping was conducted on Affymetrix Genome-Wide Human SNP 6.0 arrays. Monocyte gene expression was determined in a subgroup of 1,279 individuals by using the Illumina HT-12 v3 BeadChip. Gene expression data were confirmed in vitro and ex vivo by qRT-PCR. Genome wide analysis revealed significant associations of anti-β2GPI IgG and APOH on chromosome 17, which had been previously identified by candidate gene approaches, and of anti-domain1 and MACROD2 on chromosome 20 which has been listed in a previous GWAS as a suggestive locus associated with the occurrence of anti-β2GPI antibodies. Expression analysis confirmed increased expression of TNFα in monocytes and identified and confirmed neuron navigator 3 (NAV3) as a novel gene induced by aPL. In conclusion, MACROD2 represents a novel genetic locus associated with aPL. Furthermore, we show that aPL induce the expression of NAV3 in monocytes and endothelial cells. This will stimulate further research into the role of these genes in the APS.

TreeQTL: hierarchical error control for eQTL findings

: Commonly used multiplicity adjustments fail to control the error rate for reported findings in many expression quantitative trait loci (eQTL) studies. TreeQTL implements a hierarchical multiple testing procedure which allows control of appropriate error rates defined relative to a grouping of the eQTL hypotheses.

AVAILABILITY AND IMPLEMENTATION

The R package TreeQTL is available for download at http://bioinformatics.org/treeqtl

CONTACT

sabatti@stanford.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

C9orf72 is differentially expressed in the central nervous system and myeloid cells and consistently reduced in C9orf72, MAPT and GRN mutation carriers

A non-coding hexanucleotide repeat expansion (HRE) in C9orf72 is a common cause of amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) acting through a loss of function mechanism due to haploinsufficiency of C9orf72 or a gain of function mediated by aggregates of bidirectionally transcribed HRE-RNAs translated into di-peptide repeat (DPR) proteins. To fully understand regulation of C9orf72 expression we surveyed the C9orf72 locus using Cap Analysis of Gene Expression sequence data (CAGEseq). We observed C9orf72 was generally lowly expressed with the exception of a subset of myeloid cells, particularly CD14+ monocytes that showed up to seven fold higher expression as compared to central nervous system (CNS) and other tissues. The expression profile at the C9orf72 locus showed a complex architecture with differential expression of the transcription start sites (TSSs) for the annotated C9orf72 transcripts between myeloid and CNS tissues suggesting cell and/or tissue specific functions. We further detected novel TSSs in both the sense and antisense strand at the C9orf72 locus and confirmed their existence in brain tissues and CD14+ monocytes. Interestingly, our experiments showed a consistent decrease of C9orf72 coding transcripts not only in brain tissue and monocytes from C9orf72-HRE patients, but also in brains from MAPT and GRN mutation carriers together with an increase in antisense transcripts suggesting these could play a role in regulation of C9orf72. We found that the non-HRE related expression changes cannot be explained by promoter methylation but by the presence of the C9orf72-HRE risk haplotype and unknown functional interactions between C9orf72, MAPT and GRN.

Genome-wide association and Mendelian randomization study of NT-proBNP in patients with acute coronary syndrome

N-terminal pro-B-type natriuretic peptide (NT-proBNP) is a strong predictor of mortality in coronary artery disease and is widely employed as a prognostic biomarker. However, a causal relationship between NT-proBNP and clinical endpoints has not been established. We have performed a genome-wide association and Mendelian randomization study of NT-proBNP. We used a discovery set of 3740 patients from the PLATelet inhibition and patient Outcomes (PLATO) trial, which enrolled 18 624 patients with acute coronary syndrome (ACS). A further set of 5492 patients, from the same trial, was used for replication. Genetic variants at two novel loci (SLC39A8 and POC1B/GALNT4) were associated with NT-proBNP levels and replicated together with the previously known NPPB locus. The most significant SNP (rs198389, pooled P = 1.07 × 10(-15)) in NPPB interrupts an E-box consensus motif in the gene promoter. The association in SLC39A8 is driven by a deleterious variant (rs13107325, pooled P = 5.99 × 10(-10)), whereas the most significant SNP in POC1B/GALNT4 (rs11105306, pooled P = 1.02 × 10(-16)) is intronic. The SLC39A8 SNP was associated with higher risk of cardiovascular (CV) death (HR = 1.39, 95% CI: 1.08-1.79, P = 0.0095), but the other loci were not associated with clinical endpoints. We have identified two novel loci to be associated with NT-proBNP in patients with ACS. Only the SLC39A8 variant, but not the NPPB variant, was associated with a clinical endpoint. Due to pleotropic effects of SLC39A8, these results do not suggest that NT-proBNP levels have a direct effect on mortality in ACS patients. PLATO Clinical Trial Registration: www.clinicaltrials.gov; NCT00391872.

The importance of p53 pathway genetics in inherited and somatic cancer genomes

Decades of research have shown that mutations in the p53 stress response pathway affect the incidence of diverse cancers more than mutations in other pathways. However, most evidence is limited to somatic mutations and rare inherited mutations. Using newly abundant genomic data, we demonstrate that commonly inherited genetic variants in the p53 pathway also affect the incidence of a broad range of cancers more than variants in other pathways. The cancer-associated single nucleotide polymorphisms (SNPs) of the p53 pathway have strikingly similar genetic characteristics to well-studied p53 pathway cancer-causing somatic mutations. Our results enable insights into p53-mediated tumour suppression in humans and into p53 pathway-based cancer surveillance and treatment strategies.

Signatures of Evolutionary Adaptation in Quantitative Trait Loci Influencing Trace Element Homeostasis in Liver

Essential trace elements possess vital functions at molecular, cellular, and physiological levels in health and disease, and they are tightly regulated in the human body. In order to assess variability and potential adaptive evolution of trace element homeostasis, we quantified 18 trace elements in 150 liver samples, together with the expression levels of 90 genes and abundances of 40 proteins involved in their homeostasis. Additionally, we genotyped 169 single nucleotide polymorphism (SNPs) in the same sample set. We detected significant associations for 8 protein quantitative trait loci (pQTL), 10 expression quantitative trait loci (eQTLs), and 15 micronutrient quantitative trait loci (nutriQTL). Six of these exceeded the false discovery rate cutoff and were related to essential trace elements: 1) one pQTL for GPX2 (rs10133290); 2) two previously described eQTLs for HFE (rs12346) and SELO (rs4838862) expression; and 3) three nutriQTLs: The pathogenic C282Y mutation at HFE affecting iron (rs1800562), and two SNPs within several clustered metallothionein genes determining selenium concentration (rs1811322 and rs904773). Within the complete set of significant QTLs (which involved 30 SNPs and 20 gene regions), we identified 12 SNPs with extreme patterns of population differentiation (FST values in the top 5% percentile in at least one HapMap population pair) and significant evidence for selective sweeps involving QTLs at GPX1, SELENBP1, GPX3, SLC30A9, and SLC39A8. Overall, this detailed study of various molecular phenotypes illustrates the role of regulatory variants in explaining differences in trace element homeostasis among populations and in the human adaptive response to environmental pressures related to micronutrients.

Associations Between Genetic Variants of the Natriuretic Peptide System and Blood Pressure Response to Dietary Sodium Intervention: The GenSalt Study

BACKGROUND

The aim of this study was to comprehensively test the association of genetic variants in the natriuretic peptide (NP) system with blood pressure (BP) response to dietary sodium intervention in a Chinese population.

METHODS

We conducted a 7-day low-sodium intervention followed by a 7-day high-sodium intervention among 1,906 participants in rural China. BP measurements were obtained at baseline and each dietary intervention using a random-zero sphygmomanometer. Linear mixed-effect models were used to assess the associations of 48 single-nucleotide polymorphisms (SNPs) in 6 genes of NP system with BP response to dietary sodium intervention.

RESULTS

SNP rs5063 in the NPPA gene and SNP rs2077386 in the NPPC gene exhibited significant associations with BP response to low-sodium dietary intervention under recessive genetic model. For rs5063, absolute mean arterial pressure responses (95% confidence interval) to the low-sodium intervention were 1.31 (-1.08, 3.70) mm Hg for TT genotype and -3.74 (-4.01, -3.46) mm Hg for CC or TC genotype, respectively (P = 4.1 × 10(-5)). Individuals with at least one copy of the C allele of rs2077386 had significantly reduction in systolic BP during the low-sodium intervention compared to those with genotype GG with responses of -5.48 (-5.83, -5.14) vs. -2.76 (-3.52, -2.00) mm Hg, respectively (P = 1.9 × 10(-13)).

CONCLUSIONS

These novel findings suggested that genetic variants of NP system may contribute to the variation of BP response to sodium intervention in Chinese population. Certainly, replication of these results in other populations and further functional studies are warranted to clarify their role in the regulation of BP and hypertension.

Alzheimer's Disease Risk Polymorphisms Regulate Gene Expression in the ZCWPW1 and the CELF1 Loci

Late onset Alzheimer’s disease (LOAD) is a genetically complex and clinically heterogeneous disease. Recent large-scale genome wide association studies (GWAS) have identified more than twenty loci that modify risk for AD. Despite the identification of these loci, little progress has been made in identifying the functional variants that explain the association with AD risk. Thus, we sought to determine whether the novel LOAD GWAS single nucleotide polymorphisms (SNPs) alter expression of LOAD GWAS genes and whether expression of these genes is altered in AD brains. The majority of LOAD GWAS SNPs occur in gene dense regions under large linkage disequilibrium (LD) blocks, making it unclear which gene(s) are modified by the SNP. Thus, we tested for brain expression quantitative trait loci (eQTLs) between LOAD GWAS SNPs and SNPs in high LD with the LOAD GWAS SNPs in all of the genes within the GWAS loci. We found a significant eQTL between rs1476679 and PILRB and GATS, which occurs within the ZCWPW1 locus. PILRB and GATS expression levels, within the ZCWPW1 locus, were also associated with AD status. Rs7120548 was associated with MTCH2 expression, which occurs within the CELF1 locus. Additionally, expression of several genes within the CELF1 locus, including MTCH2, were highly correlated with one another and were associated with AD status. We further demonstrate that PILRB, as well as other genes within the GWAS loci, are most highly expressed in microglia. These findings together with the function of PILRB as a DAP12 receptor supports the critical role of microglia and neuroinflammation in AD risk.

A variant at 9p21.3 functionally implicates CDKN2B in paediatric B-cell precursor acute lymphoblastic leukaemia aetiology

Paediatric B-cell precursor acute lymphoblastic leukaemia (BCP-ALL) is the most common cancer of childhood, yet little is known about BCP-ALL predisposition. In this study, in 2,187 cases of European ancestry and 5,543 controls, we discover and replicate a locus indexed by rs77728904 at 9p21.3 associated with BCP-ALL susceptibility (Pcombined=3.32 × 10(-15), OR=1.72) and independent from rs3731217, the previously reported ALL-associated variant in this region. Of correlated SNPs tagged by this locus, only rs662463 is significant in African Americans, suggesting it is a plausible causative variant. Functional analysis shows that rs662463 is a cis-eQTL for CDKN2B, with the risk allele associated with lower expression, and suggests that rs662463 influences BCP-ALL risk by regulating CDKN2B expression through CEBPB signalling. Functional analysis of rs3731217 suggests it is associated with BCP-ALL by acting within a splicing regulatory element determining CDKN2A exon 3 usage (P=0.01). These findings provide new insights into the critical role of the CDKN2 locus in BCP-ALL aetiology.

Interaction between variants in CLU and MS4A4E modulates Alzheimer's disease risk

INTRODUCTION

Ebbert et al. reported gene-gene interactions between rs11136000-rs670139 (CLU-MS4A4E) and rs3865444-rs670139 (CD33-MS4A4E). We evaluate these interactions in the largest data set for an epistasis study.

METHODS

We tested interactions using 3837 cases and 4145 controls from Alzheimer's Disease Genetics Consortium using meta-analyses and permutation analyses. We repeated meta-analyses stratified by apolipoprotein E (APOE) ε4 status, estimated combined odds ratio (OR) and population attributable fraction (cPAF), and explored causal variants.

RESULTS

Results support the CLU-MS4A4E interaction and a dominant effect. An association between CLU-MS4A4E and APOE ε4 negative status exists. The estimated synergy factor, OR, and cPAF for rs11136000-rs670139 are 2.23, 2.45, and 8.0, respectively. We identified potential causal variants.

DISCUSSION

We replicated the CLU-MS4A4E interaction in a large case-control series and observed APOE ε4 and possible dominant effect. The CLU-MS4A4E OR is higher than any Alzheimer's disease locus except APOE ε4, APP, and TREM2. We estimated an 8% decrease in Alzheimer's disease incidence without CLU-MS4A4E risk alleles and identified potential causal variants.

Protective effect of vitamin E against diabetes-induced oxidized LDL and aorta cell wall proliferation in rat

Hyperlipidemia and oxidized-low-density lipoproteins (Ox-LDL) are important independent cardiovascular risk factors that have been shown to stimulate vascular smooth muscle cell (VSMC) proliferation. The purpose of the present study was to investigate the effect of vitamin E on Ox-LDL, lipid profile, C-reactive protein (CRP), and VSMC proliferation of rat aorta.

METHODS

Male Wistar rats (n = 32) were divided into four groups namely: sham (SH), control (C), non-treated diabetic, and vitamin E-treated diabetic (VETD) groups. Ox-LDL, lipid profile, CRP and VSMC proliferation of aorta were measured after 42 days.

RESULTS

The results revealed that along with a significant increase in VSMC proliferation, the amount of CRP, Ox-LDL, and lipid profiles in diabetic rats. VSMC proliferation was significantly ameliorated, and elevated CRP, Ox-LDL, and lipid profiles were also restored to those of shams in VETD.

CONCLUSIONS

These findings strongly support the idea that diabetes induces Ox-LDL-mediated oxidative stress and VSMC proliferation in aorta of rat and imply that vitamin E has a strong protective effect as an antioxidant.

Non-coding genomic regions possessing enhancer and silencer potential are associated with healthy aging and exceptional survival

We have completed a genome-wide linkage scan for healthy aging using data collected from a family study, followed by fine-mapping by association in a separate population, the first such attempt reported. The family cohort consisted of parents of age 90 or above and their children ranging in age from 50 to 80. As a quantitative measure of healthy aging, we used a frailty index, called FI34, based on 34 health and function variables. The linkage scan found a single significant linkage peak on chromosome 12. Using an independent cohort of unrelated nonagenarians, we carried out a fine-scale association mapping of the region suggestive of linkage and identified three sites associated with healthy aging. These healthy-aging sites (HASs) are located in intergenic regions at 12q13-14. HAS-1 has been previously associated with multiple diseases, and an enhancer was recently mapped and experimentally validated within the site. HAS-2 is a previously uncharacterized site possessing genomic features suggestive of enhancer activity. HAS-3 contains features associated with Polycomb repression. The HASs also contain variants associated with exceptional longevity, based on a separate analysis. Our results provide insight into functional genomic networks involving non-coding regulatory elements that are involved in healthy aging and longevity.

Gene Silencing and Haploinsufficiency of Csk Increase Blood Pressure

Objective

Recent genome-wide association studies have identified 33 human genetic loci that influence blood pressure. The 15q24 locus is one such locus that has been confirmed in Asians and Europeans. There are 21 genes in the locus within a 1-Mb boundary, but a functional link of these genes to blood pressure has not been reported. We aimed to identify a causative gene for blood pressure change in the 15q24 locus.

Methods and Results

CSK and ULK3 were selected as candidate genes based on eQTL analysis studies that showed the association between gene transcript levels and the lead SNP (rs1378942). Injection of siRNAs for mouse homologs Csk, Ulk3, and Cyp1a2 (negative control) showed reduced target gene mRNA levels in vivo. However, Csk siRNA only increased blood pressure while Ulk3 and Cyp1a2 siRNA did not change it. Further, blood pressure in Csk^+/- heterozygotes was higher than in wild-type, consistent with what we observed in Csk siRNA-injected mice. We confirmed that haploinsufficiency of Csk increased the active form of Src in Csk^+/- mice aorta. We also showed that inhibition of Src by PP2, a Src inhibitor decreased high blood pressure in Csk^+/- mice and the active Src in Csk^+/- mice aorta and in Csk knock-down vascular smooth muscle cells, suggesting blood pressure regulation by Csk through Src.

Conclusions

Our study demonstrates that Csk is a causative gene in the 15q24 locus and regulates blood pressure through Src, and these findings provide a novel therapeutic target for the treatment of hypertension.

RBP-Var: a database of functional variants involved in regulation mediated by RNA-binding proteins

Transcription factors bind to the genome by forming specific contacts with the primary DNA sequence; however, RNA-binding proteins (RBPs) have greater scope to achieve binding specificity through the RNA secondary structure. It has been revealed that single nucleotide variants (SNVs) that alter RNA structure, also known as RiboSNitches, exhibit 3-fold greater local structure changes than replicates of the same DNA sequence, demonstrated by the fact that depletion of RiboSNitches could result in the alteration of specific RNA shapes at thousands of sites, including 3' UTRs, binding sites of microRNAs and RBPs. However, the network between SNVs and post-transcriptional regulation remains unclear. Here, we developed RBP-Var, a database freely available at http://www.rbp-var.biols.ac.cn/, which provides annotation of functional variants involved in post-transcriptional interaction and regulation. RBP-Var provides an easy-to-use web interface that allows users to rapidly find whether SNVs of interest can transform the secondary structure of RNA and identify RBPs whose binding may be subsequently disrupted. RBP-Var integrates DNA and RNA biology to understand how various genetic variants and post-transcriptional mechanisms cooperate to orchestrate gene expression. In summary, RBP-Var is useful in selecting candidate SNVs for further functional studies and exploring causal SNVs underlying human diseases.

R2GUESS: A Graphics Processing Unit-Based R Package for Bayesian Variable Selection Regression of Multivariate Responses

Technological advances in molecular biology over the past decade have given rise to high dimensional and complex datasets offering the possibility to investigate biological associations between a range of genomic features and complex phenotypes. The analysis of this novel type of data generated unprecedented computational challenges which ultimately led to the definition and implementation of computationally efficient statistical models that were able to scale to genome-wide data, including Bayesian variable selection approaches. While extensive methodological work has been carried out in this area, only few methods capable of handling hundreds of thousands of predictors were implemented and distributed. Among these we recently proposed GUESS, a computationally optimised algorithm making use of graphics processing unit capabilities, which can accommodate multiple outcomes. In this paper we propose R2GUESS, an R package wrapping the original C++ source code. In addition to providing a user-friendly interface of the original code automating its parametrisation, and data handling, R2GUESS also incorporates many features to explore the data, to extend statistical inferences from the native algorithm (e.g., effect size estimation, significance assessment), and to visualize outputs from the algorithm. We first detail the model and its parametrisation, and describe in details its optimised implementation. Based on two examples we finally illustrate its statistical performances and flexibility.