Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

Comparison of pre-processing methods for multiplex bead-based immunoassays.

BACKGROUND: High throughput protein expression studies can be performed using bead-based protein immunoassays, such as the Luminex® xMAP® technology. Technical variability is inherent to these experiments and may lead to systematic bias and reduced power. To reduce technical variability, data pre-processing is performed. However, no recommendations exist for the pre-processing of Luminex® xMAP® data. RESULTS: We compared 37 different data pre-processing combinations of transformation and normalization methods in 42 samples on 384 analytes obtained from a multiplex immunoassay based on the Luminex® xMAP® technology. We evaluated the performance of each pre-processing approach with 6 different performance criteria. Three performance criteria were plots. All plots were evaluated by 15 independent and blinded readers. Four different combinations of transformation and normalization methods performed well as pre-processing procedure for this bead-based protein immunoassay. CONCLUSIONS: The following combinations of transformation and normalization were suitable for pre-processing Luminex® xMAP® data in this study: weighted Box-Cox followed by quantile or robust spline normalization (rsn), asinh transformation followed by loess normalization and Box-Cox followed by rsn.

Joint analysis of multiple phenotypes: summary of results and discussions from the Genetic Analysis Workshop 19.

For Genetic Analysis Workshop 19, 2 extensive data sets were provided, including whole genome and whole exome sequence data, gene expression data, and longitudinal blood pressure outcomes, together with nongenetic covariates. These data sets gave researchers the chance to investigate different aspects of more complex relationships within the data, and the contributions in our working group focused on statistical methods for the joint analysis of multiple phenotypes, which is part of the research field of data integration. The analysis of data from different sources poses challenges to researchers but provides the opportunity to model the real-life situation more realistically.Our 4 contributions all used the provided real data to identify genetic predictors for blood pressure. In the contributions, novel multivariate rare variant tests, copula models, structural equation models and a sparse matrix representation variable selection approach were applied. Each of these statistical models can be used to investigate specific hypothesized relationships, which are described together with their biological assumptions.The results showed that all methods are ready for application on a genome-wide scale and can be used or extended to include multiple omics data sets. The results provide potentially interesting genetic targets for future investigation and replication. Furthermore, all contributions demonstrated that the analysis of complex data sets could benefit from modeling correlated phenotypes jointly as well as by adding further bioinformatics information.

GUESS-ing polygenic associations with multiple phenotypes using a GPU-based evolutionary stochastic search algorithm.

Genome-wide association studies (GWAS) yielded significant advances in defining the genetic architecture of complex traits and disease. Still, a major hurdle of GWAS is narrowing down multiple genetic associations to a few causal variants for functional studies. This becomes critical in multi-phenotype GWAS where detection and interpretability of complex SNP(s)-trait(s) associations are complicated by complex Linkage Disequilibrium patterns between SNPs and correlation between traits. Here we propose a computationally efficient algorithm (GUESS) to explore complex genetic-association models and maximize genetic variant detection. We integrated our algorithm with a new Bayesian strategy for multi-phenotype analysis to identify the specific contribution of each SNP to different trait combinations and study genetic regulation of lipid metabolism in the Gutenberg Health Study (GHS). Despite the relatively small size of GHS (n  =  3,175), when compared with the largest published meta-GWAS (n > 100,000), GUESS recovered most of the major associations and was better at refining multi-trait associations than alternative methods. Amongst the new findings provided by GUESS, we revealed a strong association of SORT1 with TG-APOB and LIPC with TG-HDL phenotypic groups, which were overlooked in the larger meta-GWAS and not revealed by competing approaches, associations that we replicated in two independent cohorts. Moreover, we demonstrated the increased power of GUESS over alternative multi-phenotype approaches, both Bayesian and non-Bayesian, in a simulation study that mimics real-case scenarios. We showed that our parallel implementation based on Graphics Processing Units outperforms alternative multi-phenotype methods. Beyond multivariate modelling of multi-phenotypes, our Bayesian model employs a flexible hierarchical prior structure for genetic effects that adapts to any correlation structure of the predictors and increases the power to identify associated variants. This provides a powerful tool for the analysis of diverse genomic features, for instance including gene expression and exome sequencing data, where complex dependencies are present in the predictor space.

Genome-wide association study for refractive astigmatism reveals genetic co-determination with spherical equivalent refractive error: the CREAM consortium.

To identify genetic variants associated with refractive astigmatism in the general population, meta-analyses of genome-wide association studies were performed for: White Europeans aged at least 25 years (20 cohorts, N = 31,968); Asian subjects aged at least 25 years (7 cohorts, N = 9,295); White Europeans aged <25 years (4 cohorts, N = 5,640); and all independent individuals from the above three samples combined with a sample of Chinese subjects aged <25 years (N = 45,931). Participants were classified as cases with refractive astigmatism if the average cylinder power in their two eyes was at least 1.00 diopter and as controls otherwise. Genome-wide association analysis was carried out for each cohort separately using logistic regression. Meta-analysis was conducted using a fixed effects model. In the older European group the most strongly associated marker was downstream of the neurexin-1 (NRXN1) gene (rs1401327, P = 3.92E-8). No other region reached genome-wide significance, and association signals were lower for the younger European group and Asian group. In the meta-analysis of all cohorts, no marker reached genome-wide significance: The most strongly associated regions were, NRXN1 (rs1401327, P = 2.93E-07), TOX (rs7823467, P = 3.47E-07) and LINC00340 (rs12212674, P = 1.49E-06). For 34 markers identified in prior GWAS for spherical equivalent refractive error, the beta coefficients for genotype versus spherical equivalent, and genotype versus refractive astigmatism, were highly correlated (r = -0.59, P = 2.10E-04). This work revealed no consistent or strong genetic signals for refractive astigmatism; however, the TOX gene region previously identified in GWAS for spherical equivalent refractive error was the second most strongly associated region. Analysis of additional markers provided evidence supporting widespread genetic co-susceptibility for spherical and astigmatic refractive errors.

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 < 5 × 10(-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 up to 25,539 European-ancestry and 10,168 Japanese individuals). We observed little evidence of heterogeneity in allelic effects at these loci between groups of European and Japanese ancestry but obtained substantial improvements in the resolution of fine mapping of potential causal variants by leveraging transethnic differences in the distribution of linkage disequilibrium. We demonstrated a functional role for the most strongly associated serum albumin locus, HPN, for which Hpn knockout mice manifest low plasma albumin concentrations. Other loci associated with serum albumin harbor genes related to ribosome function, protein translation, and proteasomal degradation, whereas those associated with serum total protein include genes related to immune function. Our results highlight the advantages of transethnic meta-analysis for the discovery and fine mapping of complex trait loci and have provided initial insights into the underlying genetic architecture of serum protein concentrations and their association with human disease.

Nerve conduction velocity is regulated by the inositol polyphosphate-4-phosphatase II gene.

Impairment of nerve conduction is common in neurodegenerative and neuroinflammatory diseases such as multiple sclerosis (MS), and measurement of evoked potentials (visual, motor, or sensory) has been widely used for diagnosis and recently also as a prognostic marker for MS. We used a classical genetic approach to identify novel genes controlling nerve conduction. First, we used quantitative trait mapping in F2 progeny of B10/SJL mice to identify EAE31, a locus controlling latency of motor evoked potentials (MEPs) and clinical onset of experimental autoimmune encephalomyelitis. Then, by combining congenic mapping, in silico haplotype analyses, and comparative genomics we identified inositol polyphosphate-4-phosphatase, type II (Inpp4b) as the quantitative trait gene for EAE31. Sequence variants of Inpp4b (C/A, exon 13; A/C, exon 14) were identified as differing among multiple mouse strains and correlated with individual cortical MEP latency differences. To evaluate the functional relevance of the amino acid exchanges at positions S474R and H548P, we generated transgenic mice carrying the longer-latency allele (Inpp4b(474R/548P)) in the C57BL/6J background. Inpp4b(474R/548P) mice exhibited significantly longer cortical MEP latencies (4.5 ± 0.22 ms versus 3.7 ± 0.13 ms; P = 1.04 × 10(-9)), indicating that INPP4B regulates nerve conduction velocity. An association of an INPP4B polymorphism (rs13102150) with MS was observed in German and Spanish MS cohorts (3676 controls and 911 cases) (P = 8.8 × 10(-3)).

Genome-wide association study identifies new susceptibility loci for cutaneous lupus erythematosus.

Cutaneous lupus erythematosus (CLE) is a chronic autoimmune disease of the skin with typical clinical manifestations. Here, we genotyped 906 600 single nucleotide polymorphisms (SNPs) in 183 CLE cases and 1288 controls of Central European ancestry. Replication was performed for 13 SNPs in 219 case subjects and 262 controls from Finland. Association was particularly pronounced at 4 loci, all with genomewide significance (P < 5 × 10(-8) ): rs2187668 (PGWAS  = 1.4 × 10(-12) ), rs9267531 (PGWAS  = 4.7 × 10(-10) ), rs4410767 (PGWAS  = 1.0 × 10(-9) ) and rs3094084 (PGWAS  = 1.1 × 10(-9) ). All mentioned SNPs are located within the major histocompatibility complex (MHC) region of chromosome 6 and near genes of known immune functions or associations with other autoimmune diseases such as HLA-DQ alpha chain 1 (HLA-DQA1), MICA, MICB, MSH5, TRIM39 and RPP21. For example, TRIM39/RPP21 read through transcript is a known mediator of the interferon response, a central pathway involved in the pathogenesis of CLE and systemic lupus erythematosus (SLE). Taken together, this genomewide analysis of disease association of CLE identified candidate genes and genomic regions that may contribute to pathogenic mechanisms in CLE via dysregulated antigen presentation (HLA-DQA1), apoptosis regulation, RNA processing and interferon response (MICA, MICB, MSH5, TRIM39 and RPP21).

MicroRNA profiling of low-grade and transformed nodal marginal zone lymphoma reveals a similar signature pattern distinct from diffuse large B cell lymphoma.

BACKGROUND/AIMS: As critical post-transcriptional regulators of gene expression, microRNAs are involved in several cellular processes of vital impact including cell growth and apoptosis. Many hematologic malignancies exhibit distinct microRNA signatures. MicroRNA implication in the pathogenesis of nodal marginal zone lymphoma (NMZL), however, remains widely elusive. METHODS: Comprehensive morphologic, immunophenotypic and cytogenetic studies were carried out on a cohort of NMZL (n = 30) incorporating indolent as well as transformed MZL. In addition, microRNA signatures were generated, employing a quantitative real-time polymerase chain reaction approach. These were then compared to signatures from cases of diffuse large B cell lymphoma (DLBCL) alongside reactive lymph node controls. RESULTS: While microRNA signatures of low-grade and transformed NMZL did not differ significantly, several microRNAs were differentially expressed between transformed NMZL and DLBCL, hinting at molecularly distinct mechanisms of lymphomagenesis and indicating the biological disparity of transformed NMZL from DLBCL. CONCLUSION: In the light of the unresolved issue regarding the classification of marginal zone-derived transformed B-cell neoplasms, microRNAs may be a valuable aid in discriminating NMZL from DLBCL.

Genome-wide association study in musician's dystonia: a risk variant at the arylsulfatase G locus?

Musician's dystonia (MD) affects 1% to 2% of professional musicians and frequently terminates performance careers. It is characterized by loss of voluntary motor control when playing the instrument. Little is known about genetic risk factors, although MD or writer's dystonia (WD) occurs in relatives of 20% of MD patients. We conducted a 2-stage genome-wide association study in whites. Genotypes at 557,620 single-nucleotide polymorphisms (SNPs) passed stringent quality control for 127 patients and 984 controls. Ten SNPs revealed P < 10(-5) and entered the replication phase including 116 MD patients and 125 healthy musicians. A genome-wide significant SNP (P < 5 × 10(-8) ) was also genotyped in 208 German or Dutch WD patients, 1,969 Caucasian, Spanish, and Japanese patients with other forms of focal or segmental dystonia as well as in 2,233 ethnically matched controls. Genome-wide significance with MD was observed for an intronic variant in the arylsulfatase G (ARSG) gene (rs11655081; P = 3.95 × 10(-9) ; odds ratio [OR], 4.33; 95% confidence interval [CI], 2.66-7.05). rs11655081 was also associated with WD (P = 2.78 × 10(-2) ) but not with any other focal or segmental dystonia. The allele frequency of rs11655081 varies substantially between different populations. The population stratification in our sample was modest (λ = 1.07), but the effect size may be overestimated. Using a small but homogenous patient sample, we provide data for a possible association of ARSG with MD. The variant may also contribute to the risk of WD, a form of dystonia that is often found in relatives of MD patients.

BiomarCaRE: rationale and design of the European BiomarCaRE project including 300,000 participants from 13 European countries.

Biomarkers are considered as tools to enhance cardiovascular risk estimation. However, the value of biomarkers on risk estimation beyond European risk scores, their comparative impact among different European regions and their role towards personalised medicine remains uncertain. Biomarker for Cardiovascular Risk Assessment in Europe (BiomarCaRE) is an European collaborative research project with the primary objective to assess the value of established and emerging biomarkers for cardiovascular risk prediction. BiomarCaRE integrates clinical and epidemiological biomarker research and commercial enterprises throughout Europe to combine innovation in biomarker discovery for cardiovascular disease prediction with consecutive validation of biomarker effectiveness in large, well-defined primary and secondary prevention cohorts including over 300,000 participants from 13 European countries. Results from this study will contribute to improved cardiovascular risk prediction across different European populations. The present publication describes the rationale and design of the BiomarCaRE project.

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 major trans-acting regulators. We detected three genomic regions significantly associated with co-regulated gene modules. Association of these loci with multiple expression traits was replicated in Cardiogenics, an independent study in which expression profiles of monocytes were available in 758 subjects. The locus 12q13 (lead SNP rs11171739), previously identified as a type 1 diabetes locus, was associated with a pattern including two cis eQTLs, RPS26 and SUOX, and 5 trans eQTLs, one of which (MADCAM1) is a potential candidate for mediating T1D susceptibility. The locus 12q24 (lead SNP rs653178), which has demonstrated extensive disease pleiotropy, including type 1 diabetes, hypertension, and celiac disease, was associated to a pattern strongly correlating to blood pressure level. The strongest trans eQTL in this pattern was CRIP1, a known marker of cellular proliferation in cancer. The locus 12q15 (lead SNP rs11177644) was associated with a pattern driven by two cis eQTLs, LYZ and YEATS4, and including 34 trans eQTLs, several of them tumor-related genes. This study shows that a method exploiting the structure of co-expressions among genes can help identify genomic regions involved in trans regulation of sets of genes and can provide clues for understanding the mechanisms linking genome-wide association loci to disease.

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). Second, we used as an instrument a genetic score consisting of 14 common SNPs that exclusively associate with HDL cholesterol and tested this score in up to 12,482 cases of myocardial infarction and 41,331 controls. As a positive control, we also tested a genetic score of 13 common SNPs exclusively associated with LDL cholesterol. FINDINGS: Carriers of the LIPG 396Ser allele (2·6% frequency) had higher HDL cholesterol (0·14 mmol/L higher, p=8×10(-13)) but similar levels of other lipid and non-lipid risk factors for myocardial infarction compared with non-carriers. This difference in HDL cholesterol is expected to decrease risk of myocardial infarction by 13% (odds ratio [OR] 0·87, 95% CI 0·84-0·91). However, we noted that the 396Ser allele was not associated with risk of myocardial infarction (OR 0·99, 95% CI 0·88-1·11, p=0·85). From observational epidemiology, an increase of 1 SD in HDL cholesterol was associated with reduced risk of myocardial infarction (OR 0·62, 95% CI 0·58-0·66). However, a 1 SD increase in HDL cholesterol due to genetic score was not associated with risk of myocardial infarction (OR 0·93, 95% CI 0·68-1·26, p=0·63). For LDL cholesterol, the estimate from observational epidemiology (a 1 SD increase in LDL cholesterol associated with OR 1·54, 95% CI 1·45-1·63) was concordant with that from genetic score (OR 2·13, 95% CI 1·69-2·69, p=2×10(-10)). INTERPRETATION: Some genetic mechanisms that raise plasma HDL cholesterol do not seem to lower risk of myocardial infarction. These data challenge the concept that raising of plasma HDL cholesterol will uniformly translate into reductions in risk of myocardial infarction. FUNDING: US National Institutes of Health, The Wellcome Trust, European Union, British Heart Foundation, and the German Federal Ministry of Education and Research.

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06 × 10⁻6, OR  = 0.67 [95% CI 0.57-0.79] for the minor allele T). Three more SNPs showed p < 2.21 × 10⁻5. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n =564, n = 981 controls, p = 2.07 × 10⁻³, OR = 0.79 [95% CI 0.67-0.92]), France 1 (n = 433 cases, n = 395 controls, p =3.73 × 10⁻³, OR  = 0.74 [95% CI 0.60-0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26 × 10⁻4, OR  = 0.63 [95% CI 0.50-0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28 × 10⁻¹³, OR= 0.72 [95% CI 0.65-0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.

Genome-wide association study for coronary artery calcification with follow-up in myocardial infarction.

BACKGROUND: Coronary artery calcification (CAC) detected by computed tomography is a noninvasive measure of coronary atherosclerosis, which underlies most cases of myocardial infarction (MI). We sought to identify common genetic variants associated with CAC and further investigate their associations with MI. METHODS AND RESULTS: Computed tomography was used to assess quantity of CAC. A meta-analysis of genome-wide association studies for CAC was performed in 9961 men and women from 5 independent community-based cohorts, with replication in 3 additional independent cohorts (n=6032). We examined the top single-nucleotide polymorphisms (SNPs) associated with CAC quantity for association with MI in multiple large genome-wide association studies of MI. Genome-wide significant associations with CAC for SNPs on chromosome 9p21 near CDKN2A and CDKN2B (top SNP: rs1333049; P=7.58×10(-19)) and 6p24 (top SNP: rs9349379, within the PHACTR1 gene; P=2.65×10(-11)) replicated for CAC and for MI. Additionally, there is evidence for concordance of SNP associations with both CAC and MI at a number of other loci, including 3q22 (MRAS gene), 13q34 (COL4A1/COL4A2 genes), and 1p13 (SORT1 gene). CONCLUSIONS: SNPs in the 9p21 and PHACTR1 gene loci were strongly associated with CAC and MI, and there are suggestive associations with both CAC and MI of SNPs in additional loci. Multiple genetic loci are associated with development of both underlying coronary atherosclerosis and clinical events.

Population-based meta-analysis in Caucasians confirms association with COL5A1 and ZNF469 but not COL8A2 with central corneal thickness.

Central corneal thickness (CCT) has become an endophenotype of major interest for the genetically complex disorder glaucoma. CCT has a high heritability, and thin CCT is an independent risk factor for the diagnosis and progression of open-angle glaucoma. Genome-wide association studies thus provide genetic loci associated with CCT and potentially related to open-angle glaucoma. The distribution of CCT and prevalence of glaucoma in population-based studies have demonstrated ethnic differences suggesting ethnic-dependent variations in the genetic determinants of CCT. We conducted a genome-wide association study in Caucasians (n = 3,931) from the Gutenberg Health Study (Germany) followed by replication of 30 genome-wide significant SNPs or SNPs of interest (P < 10(-5)) in the Rotterdam Study (The Netherlands, n = 1,418). In a combined analysis, we confirmed quantitative trait loci on chromosomes 9q34 and 16q24 for association with CCT. On chromosome 16q24, the locus is located in an intergenic region near the ZNF469 gene (top SNP: rs9938149, P = 1.45 × 10(-12)). ZNF469 missense mutation is involved in a syndrome with very thin cornea (brittle cornea syndrome). The second locus on chromosome 9q34 represents the intergenic region between the RXRA and COL5A1 gene (top SNP: rs3132306, P = 2.71 × 10(-10)). Collagen type 5 determines the diameter of the corneal collagen fibrils. In our Caucasian population-based GWA study, we reinforce the involvement of collagen-related genes influencing CCT in Caucasians. We could not confirm the collagen type 8 locus on chromosome 1 as reported in Asian studies.

Influence of sex and genetic variability on expression of X-linked genes in human monocytes.

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<2×10(-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 quite similar in males and females. Few genes exhibited a stronger genetic effect in females than in males (ARSD, DCX, POLA1 and ITM2A). These genes would deserve further investigation since they may contribute to sex pathophysiological differences.

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-nucleotide polymorphisms (SNPs; imputed to the HapMap CEU panel) and echocardiographic traits. In stage 2, SNPs significantly associated with traits in stage 1 were tested for association in 2 other cohorts (n = 4094 people of European ancestry). Using a prespecified P value threshold of 5 x 10(-7) to indicate genome-wide significance, we performed an inverse variance-weighted fixed-effects meta-analysis of genome-wide association data from each cohort. MAIN OUTCOME MEASURES: Echocardiographic traits: LV mass, internal dimensions, wall thickness, systolic dysfunction, aortic root, and left atrial size. RESULTS: In stage 1, 16 genetic loci were associated with 5 echocardiographic traits: 1 each with LV internal dimensions and systolic dysfunction, 3 each with LV mass and wall thickness, and 8 with aortic root size. In stage 2, 5 loci replicated (6q22 locus associated with LV diastolic dimensions, explaining <1% of trait variance; 5q23, 12p12, 12q14, and 17p13 associated with aortic root size, explaining 1%-3% of trait variance). CONCLUSIONS: We identified 5 genetic loci harboring common variants that were associated with variation in LV diastolic dimensions and aortic root size, but such findings explained a very small proportion of variance. Further studies are required to replicate these findings, identify the causal variants at or near these loci, characterize their functional significance, and determine whether they are related to overt cardiovascular disease.

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 protein levels encoded by candidate genes for subsequent cardiovascular disease was investigated. Eight transcripts (CRIP1, MYADM, TIPARP, TSC22D3, CEBPA, F12, LMNA, and TPPP3) were identified jointly accounting for up to 13% (95% confidence interval, 8.7-16.2) of BP variability. Changes in CRIP1, MYADM, TIPARP, LMNA, TSC22D3, CEBPA, and TPPP3 expression associated with BP changes-among these, CRIP1 gene expression was additionally correlated to measures of cardiac hypertrophy. Assessment of circulating CRIP1 (cystein-rich protein 1) levels as biomarkers showed a strong association with increased risk for incident stroke (hazard ratio, 1.06; 95% confidence interval, 1.03-1.09; P=5.0×10-5). Our comprehensive analysis of global gene expression highlights 8 novel transcripts significantly associated with BP, providing a link between gene expression and BP. Translational approaches further established evidence for the potential use of CRIP1 as emerging disease-related biomarker.