New susceptibility locus for coronary artery disease on chromosome 3q22.3
We present a three-stage analysis of genome-wide SNP data in 1,222 German individuals with myocardial infarction and 1,298 controls, in silico replication in three additional genome-wide datasets of coronary artery disease (CAD) and subsequent replication in approximately 25,000 subjects. We identified one new CAD risk locus on 3q22.3 in MRAS (P = 7.44 x 10(-13); OR = 1.15, 95% CI = 1.11-1.19), and suggestive association with a locus on 12q24.31 near HNF1A-C12orf43 (P = 4.81 x 10(-7); OR = 1.08, 95% CI = 1.05-1.11).
Novel multiple sclerosis susceptibility loci implicated in epigenetic regulation
Genome-wide study in Germans identifies four novel multiple sclerosis risk genes and confirms already known gene loci.
Removing Batch Effects from Longitudinal Gene Expression - Quantile Normalization Plus ComBat as Best Approach for Microarray Transcriptome Data
International audience; 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 fol…
Testing for goodness rather than lack of fit of an X–chromosomal SNP to the Hardy-Weinberg model
The problem of checking the genotype distribution obtained for some diallelic marker for compatibility with the Hardy-Weinberg equilibrium (HWE) condition arises also for loci on the X chromosome. The possible genotypes depend on the sex of the individual in this case: for females, the genotype distribution is trinomial, as in the case of an autosomal locus, whereas a binomial proportion is observed for males. Like in genetic association studies with autosomal SNPs, interest is typically in establishing approximate compatibility of the observed genotype frequencies with HWE. This requires to replace traditional methods tailored for detecting lack of fit to the model with an equivalence test…
Genome-Wide Haplotype Analysis of Cis Expression Quantitative Trait Loci in Monocytes
In order to assess whether gene expression variability could be influenced by several SNPs acting in cis, either through additive or more complex haplotype effects, a systematic genome-wide search for cis haplotype expression quantitative trait loci (eQTL) was conducted in a sample of 758 individuals, part of the Cardiogenics Transcriptomic Study, for which genome-wide monocyte expression and GWAS data were available. 19,805 RNA probes were assessed for cis haplotypic regulation through investigation of ∼2,1×109 haplotypic combinations. 2,650 probes demonstrated haplotypic p-values >104-fold smaller than the best single SNP p-value. Replication of significant haplotype effects were tested f…
A trans-acting locus regulates an anti-viral expression network and type 1 diabetes risk
Combined analyses of gene networks and DNA sequence variation can provide new insights into the aetiology of common diseases that may not be apparent from genome-wide association studies alone. Recent advances in rat genomics are facilitating systems-genetics approaches. Here we report the use of integrated genome-wide approaches across seven rat tissues to identify gene networks and the loci underlying their regulation. We defined an interferon regulatory factor 7 (IRF7)-driven inflammatory network (IDIN) enriched for viral response genes, which represents a molecular biomarker for macrophages and which was regulated in multiple tissues by a locus on rat chromosome 15q25. We show that Epst…
Adaptive linear rank tests for eQTL studies
Expression quantitative trait loci (eQTL) studies are performed to identify single-nucleotide polymorphisms that modify average expression values of genes, proteins, or metabolites, depending on the genotype. As expression values are often not normally distributed, statistical methods for eQTL studies should be valid and powerful in these situations. Adaptive tests are promising alternatives to standard approaches, such as the analysis of variance or the Kruskal-Wallis test. In a two-stage procedure, skewness and tail length of the distributions are estimated and used to select one of several linear rank tests. In this study, we compare two adaptive tests that were proposed in the literatur…
Meta-analysis of genome-wide association studies from the CHARGE consortium identifies common variants associated with carotid intima media thickness and plaque
Carotid intima media thickness (cIMT) and plaque determined by ultrasonography are established measures of subclinical atherosclerosis that each predicts future cardiovascular disease events. We conducted a meta-analysis of genome-wide association data in 31,211 participants of European ancestry from nine large studies in the setting of the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. We then sought additional evidence to support our findings among 11,273 individuals using data from seven additional studies. In the combined meta-analysis, we identified three genomic regions associated with common carotid intima media thickness and two different regions a…
The Choice of the Filtering Method in Microarrays Affects the Inference Regarding Dosage Compensation of the Active X-Chromosome
BackgroundThe hypothesis of dosage compensation of genes of the X chromosome, supported by previous microarray studies, was recently challenged by RNA-sequencing data. It was suggested that microarray studies were biased toward an over-estimation of X-linked expression levels as a consequence of the filtering of genes below the detection threshold of microarrays.Methodology/principal findingsTo investigate this hypothesis, we used microarray expression data from circulating monocytes in 1,467 individuals. In total, 25,349 and 1,156 probes were unambiguously assigned to autosomes and the X chromosome, respectively. Globally, there was a clear shift of X-linked expressions toward lower levels…
Genetic variants associated with cardiac structure and function: a meta-analysis and replication of genome-wide association data.
Context Echocardiographic measures of left ventricular (LV) structure and function are heritable phenotypes of cardiovascular disease. Objective To identify common genetic variants associated with cardiac structure and function by conducting a meta-analysis of genome-wide association data in 5 population-based cohort studies (stage 1) with replication (stage 2) in 2 other community-based samples. Design, Setting, and Participants Within each of 5 community-based cohorts comprising the EchoGen consortium (stage 1; n = 12 612 individuals of European ancestry; 55% women, aged 26-95 years; examinations between 1978-2008), we estimated the association between approximately 2.5 million single-nuc…
Comprehensive exploration of the effects of miRNA SNPs on monocyte gene expression.
We aimed to assess whether pri-miRNA SNPs (miSNPs) could influence monocyte gene expression, either through marginal association or by interacting with polymorphisms located in 3'UTR regions (3utrSNPs). We then conducted a genome-wide search for marginal miSNPs effects and pairwise miSNPs × 3utrSNPs interactions in a sample of 1,467 individuals for which genome-wide monocyte expression and genotype data were available. Statistical associations that survived multiple testing correction were tested for replication in an independent sample of 758 individuals with both monocyte gene expression and genotype data. In both studies, the hsa-mir-1279 rs1463335 was found to modulate in cis the expres…
Successful Replication of GWAS Hits for Multiple Sclerosis in 10,000 Germans Using the Exome Array
Genome-wide association studies (GWAS) successfully identified various chromosomal regions to be associated with multiple sclerosis (MS). The primary aim of this study was to replicate reported associations from GWAS using an exome array in a large German study. German MS cases (n = 4,476) and German controls (n = 5,714) were genotyped using the Illumina HumanExome v1-Chip. Genotype calling was performed with the Illumina Genome Studio(TM) Genotyping Module, followed by zCall. Single-nucleotide polymorphisms (SNPs) in seven regions outside the human leukocyte antigen (HLA) region showed genome-wide significant associations with MS (P values < 5 × 10(-8) ). These associations have been repor…
Large-scale association analysis identifies new risk loci for coronary artery disease
Coronary artery disease (CAD) is the commonest cause of death. Here, we report an association analysis in 63,746 CAD cases and 130,681 controls identifying 15 loci reaching genome-wide significance, taking the number of susceptibility loci for CAD to 46, and a further 104 independent variants (r 2 < 0.2) strongly associated with CAD at a 5% false discovery rate (FDR). Together, these variants explain approximately 10.6% of CAD heritability. Of the 46 genome-wide significant lead SNPs, 12 show a significant association with a lipid trait, and 5 show a significant association with blood pressure, but none is significantly associated with diabetes. Network analysis with 233 candidate genes …
Towards a Physical Map of the HLA Complex
The human major histocompatibility (HLA) complex is located on the short arm of chromosome 6 in the 6p21.31→6p21.33 region (1,2). There are three clusters of genes, the HLA class I, II, and III regions. Whereas the class III loci are only moderately polymorphic (see (3) for review), the class I and II genes which encode cell surface glycoproteins show an extreme degree of polymorphism. There are a minimum of 17 class I loci (4) and at least 15 genes for class II alpha and beta chains (5). In addition, we (6) and others (7) have recently demonstrated that also the loci for tumor necrosis factor (TNFA) and lymphotoxin (TNFB) are part of the HLA complex [see also Ragoussis et al., this volume]…
Association Between Chromosome 9p21 Variants and the Ankle-Brachial Index Identified by a Meta-Analysis of 21 Genome-Wide Association Studies
Background— Genetic determinants of peripheral arterial disease (PAD) remain largely unknown. To identify genetic variants associated with the ankle-brachial index (ABI), a noninvasive measure of PAD, we conducted a meta-analysis of genome-wide association study data from 21 population-based cohorts. Methods and Results— Continuous ABI and PAD (ABI ≤0.9) phenotypes adjusted for age and sex were examined. Each study conducted genotyping and imputed data to the ≈2.5 million single nucleotide polymorphisms (SNPs) in HapMap. Linear and logistic regression models were used to test each SNP for association with ABI and PAD using additive genetic models. Study-specific data were combined using fi…
Influence of sex and genetic variability on expression of X-linked genes in human monocytes
Abstract In humans, the fraction of X-linked genes with higher expression in females has been estimated to be 5% from microarray studies, a proportion lower than the 25% of genes thought to escape X inactivation. We analyzed 715 X-linked transcripts in circulating monocytes from 1,467 subjects and found an excess of female-biased transcripts on the X compared to autosomes (9.4% vs 5.5%, p −5 ). Among the genes not previously known to escape inactivation, the most significant one was EFHC2 whose 20% of variability was explained by sex. We also investigated cis expression quantitative trait loci (eQTLs) by analyzing 15,703 X-linked SNPs. The frequency and magnitude of X-linked cis eQTLs were…
Transcriptome-Wide Analysis Identifies Novel Associations With Blood Pressure.
Hypertension represents a major cardiovascular risk factor. The pathophysiology of increased blood pressure (BP) is not yet completely understood. Transcriptome profiling offers possibilities to uncover genetics effects on BP. Based on 2 populations including 2549 individuals, a meta-analyses of monocytic transcriptome-wide profiles were performed to identify transcripts associated with BP. Replication was performed in 2 independent studies of whole-blood transcriptome data including 1990 individuals. For identified candidate genes, a direct link between long-term changes in BP and gene expression over time and by treatment with BP-lowering therapy was assessed. The predictive value of pro…
Treatment choices and neuropsychological symptoms of a large cohort of early MS
ObjectiveTo assess clinical characteristics, distribution of disease-modifying treatments (DMTs), and neuropsychological symptoms in a large cohort of patients with early-stage MS.MethodsThe German National MS Cohort is a multicenter prospective longitudinal cohort study that has recruited DMT-naive patients with clinically isolated syndrome (CIS) and relapsing-remitting MS (RRMS) since 2010. We evaluated their baseline characteristics and the prevalence of neuropsychological symptoms.ResultsOf 1,124 patients, with a 2.2:1 female-to-male ratio and median age at onset of 31.71 years (interquartile range [IQR]: 26.06–40.33), 44.6% and 55.3% had CIS and RRMS, respectively. The median Expanded …
Low frequency and rare coding variation contributes to multiple sclerosis risk
AbstractMultiple sclerosis is a common, complex neurological disease, where almost 20% of risk heritability can be attributed to common genetic variants, including >230 identified by genome-wide association studies (Patsopoulos et al., 2017). Multiple strands of evidence suggest that the majority of the remaining heritability is also due to the additive effects of individual variants, rather than epistatic interactions between these variants, or mutations exclusive to individual families. Here, we show in 68,379 cases and controls that as much as 5% of this heritability is explained by low-frequency variation in gene coding sequence. We identify four novel genes driving MS risk independe…
Genetics and Beyond – The Transcriptome of Human Monocytes and Disease Susceptibility
BACKGROUND: Variability of gene expression in human may link gene sequence variability and phenotypes; however, non-genetic variations, alone or in combination with genetics, may also influence expression traits and have a critical role in physiological and disease processes. METHODOLOGY/PRINCIPAL FINDINGS: To get better insight into the overall variability of gene expression, we assessed the transcriptome of circulating monocytes, a key cell involved in immunity-related diseases and atherosclerosis, in 1,490 unrelated individuals and investigated its association with >675,000 SNPs and 10 common cardiovascular risk factors. Out of 12,808 expressed genes, 2,745 expression quantitative trait …
Integrating genome-wide genetic variations and monocyte expression data reveals trans-regulated gene modules in humans.
One major expectation from the transcriptome in humans is to characterize the biological basis of associations identified by genome-wide association studies. So far, few cis expression quantitative trait loci (eQTLs) have been reliably related to disease susceptibility. Trans-regulating mechanisms may play a more prominent role in disease susceptibility. We analyzed 12,808 genes detected in at least 5% of circulating monocyte samples from a population-based sample of 1,490 European unrelated subjects. We applied a method of extraction of expression patterns—independent component analysis—to identify sets of co-regulated genes. These patterns were then related to 675,350 SNPs to identify maj…
Discovery and Fine Mapping of Serum Protein Loci through Transethnic Meta-analysis
Many disorders are associated with altered serum protein concentrations, including malnutrition, cancer, and cardiovascular, kidney, and inflammatory diseases. Although these protein concentrations are highly heritable, relatively little is known about their underlying genetic determinants. Through transethnic meta-analysis of European-ancestry and Japanese genome-wide association studies, we identified six loci at genome-wide significance (p −8 ) for serum albumin ( HPN-SCN1B , GCKR-FNDC4 , SERPINF2-WDR81 , TNFRSF11A-ZCCHC2 , FRMD5-WDR76 , and RPS11-FCGRT , in up to 53,190 European-ancestry and 9,380 Japanese individuals) and three loci for total protein ( TNFRS13B , 6q21.3, and ELL2 , in …
Plasma HDL cholesterol and risk of myocardial infarction: a mendelian randomisation study
BACKGROUND: High plasma HDL cholesterol is associated with reduced risk of myocardial infarction, but whether this association is causal is unclear. Exploiting the fact that genotypes are randomly assigned at meiosis, are independent of non-genetic confounding, and are unmodified by disease processes, mendelian randomisation can be used to test the hypothesis that the association of a plasma biomarker with disease is causal. METHODS: We performed two mendelian randomisation analyses. First, we used as an instrument a single nucleotide polymorphism (SNP) in the endothelial lipase gene (LIPG Asn396Ser) and tested this SNP in 20 studies (20,913 myocardial infarction cases, 95,407 controls). Se…
Design of the Coronary ARtery DIsease Genome-Wide Replication And Meta-Analysis (CARDIoGRAM) Study
Background— Recent genome-wide association studies (GWAS) of myocardial infarction (MI) and other forms of coronary artery disease (CAD) have led to the discovery of at least 13 genetic loci. In addition to the effect size, power to detect associations is largely driven by sample size. Therefore, to maximize the chance of finding novel susceptibility loci for CAD and MI, the Coronary ARtery DIsease Genome-wide Replication And Meta-analysis (CARDIoGRAM) consortium was formed. Methods and Results— CARDIoGRAM combines data from all published and several unpublished GWAS in individuals with European ancestry; includes >22 000 cases with CAD, MI, or both and >60 000 controls; and unifies …
Linkage disequilibrium of MTHFR genotypes 677C/T-1298A/C in the German population and association studies in probands with neural tube defects(NTD).
A number of studies have demonstrated that the common polymorphism 677CT in the gene encoding 5, 10-methylenetetrahydrofolate reductase (MTHFR) leads to a thermolabile variant with decreased enzyme activity and to mildly elevated plasma homocysteine. 677TT homozygosity was shown to be more frequent in NTD probands compared with controls in some studies. Recently, another polymorphism, 1298AC, in the MTHFR gene was described and combined heterozygosity 677CT/1298AC was suggested to be an additional risk factor for NTD. The present study examines the genotype and haplotype distribution of the two polymorphisms in the German population and evaluates the impact on NTD individuals and their rela…
Large-scale genome-wide analysis identifies genetic variants associated with cardiac structure and function
BACKGROUND: Understanding the genetic architecture of cardiac structure and function may help to prevent and treat heart disease. This investigation sought to identify common genetic variations associated with inter-individual variability in cardiac structure and function. METHODS: A GWAS meta-analysis of echocardiographic traits was performed, including 46,533 individuals from 30 studies (EchoGen consortium). The analysis included 16 traits of left ventricular (LV) structure, and systolic and diastolic function. RESULTS: The discovery analysis included 21 cohorts for structural and systolic function traits (n = 32,212) and 17 cohorts for diastolic function traits (n = 21,852). Replication …
Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function
Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, 19 associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation o…
Analyzing Illumina Gene Expression Microarray Data from Different Tissues: Methodological Aspects of Data Analysis in the MetaXpress Consortium
Microarray profiling of gene expression is widely applied in molecular biology and functional genomics. Experimental and technical variations make meta-analysis of different studies challenging. In a total of 3358 samples, all from German population-based cohorts, we investigated the effect of data preprocessing and the variability due to sample processing in whole blood cell and blood monocyte gene expression data, measured on the Illumina HumanHT-12 v3 BeadChip array. Gene expression signal intensities were similar after applying the log(2) or the variance-stabilizing transformation. In all cohorts, the first principal component (PC) explained more than 95% of the total variation. Technic…
New loci associated with kidney function and chronic kidney disease
Chronic kidney disease (CKD) is a significant public health problem, and recent genetic studies have identified common CKD susceptibility variants. The CKDGen consortium performed a meta-analysis of genome-wide association data in 67,093 individuals of European ancestry from 20 predominantly population-based studies in order to identify new susceptibility loci for reduced renal function as estimated by serum creatinine (eGFRcrea), serum cystatin c (eGFRcys) and CKD (eGFRcrea 60 ml/min/1.73 m 2; n = 5,807 individuals with CKD (cases)). Follow-up of the 23 new genome-wide-significant loci (P 5 × 10 8) in 22,982 replication samples identified 13 new loci affecting renal function and CKD (in or…