Meta-analysis of correlated traits via summary statistics from GWASs with an application in hypertension

Genome-wide association studies (GWASs) have identified many genetic variants underlying complex traits. Many detected genetic loci harbor variants that associate with multiple-even distinct-traits. Most current analysis approaches focus on single traits, even though the final results from multiple traits are evaluated together. Such approaches miss the opportunity to systemically integrate the phenome-wide data available for genetic association analysis. In this study, we propose a general approach that can integrate association evidence from summary statistics of multiple traits, either correlated, independent, continuous, or binary traits, which might come from the same or different studies. We allow for trait heterogeneity effects. Population structure and cryptic relatedness can also be controlled. Our simulations suggest that the proposed method has improved statistical power over single-trait analysis in most of the cases we studied. We applied our method to the Continental Origins and Genetic Epidemiology Network (COGENT) African ancestry samples for three blood pressure traits and identified four loci (CHIC2, HOXA-EVX1, IGFBP1/IGFBP3, and CDH17; p < 5.0 × 10(-8)) associated with hypertension-related traits that were missed by a single-trait analysis in the original report. Six additional loci with suggestive association evidence (p < 5.0 × 10(-7)) were also observed, including CACNA1D and WNT3. Our study strongly suggests that analyzing multiple phenotypes can improve statistical power and that such analysis can be executed with the summary statistics from GWASs. Our method also provides a way to study a cross phenotype (CP) association by using summary statistics from GWASs of multiple phenotypes.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci.

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 × 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 × 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 × 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.

Impact of a cis-associated gene expression SNP on chromosome 20q11.22 on bipolar disorder susceptibility, hippocampal structure and cognitive performance.

BackgroundBipolar disorder is a highly heritable polygenic disorder. Recent enrichment analyses suggest that there may be true risk variants for bipolar disorder in the expression quantitative trait loci (eQTL) in the brain.AimsWe sought to assess the impact of eQTL variants on bipolar disorder risk by combining data from both bipolar disorder genome-wide association studies (GWAS) and brain eQTL.MethodTo detect single nucleotide polymorphisms (SNPs) that influence expression levels of genes associated with bipolar disorder, we jointly analysed data from a bipolar disorder GWAS (7481 cases and 9250 controls) and a genome-wide brain (cortical) eQTL (193 healthy controls) using a Bayesian statistical method, with independent follow-up replications. The identified risk SNP was then further tested for association with hippocampal volume (n = 5775) and cognitive performance (n = 342) among healthy individuals.ResultsIntegrative analysis revealed a significant association between a brain eQTL rs6088662 on chromosome 20q11.22 and bipolar disorder (log Bayes factor = 5.48; bipolar disorder P = 5.85×10(-5)). Follow-up studies across multiple independent samples confirmed the association of the risk SNP (rs6088662) with gene expression and bipolar disorder susceptibility (P = 3.54×10(-8)). Further exploratory analysis revealed that rs6088662 is also associated with hippocampal volume and cognitive performance in healthy individuals.ConclusionsOur findings suggest that 20q11.22 is likely a risk region for bipolar disorder; they also highlight the informative value of integrating functional annotation of genetic variants for gene expression in advancing our understanding of the biological basis underlying complex disorders, such as bipolar disorder.

Tissue-specific regulatory circuits reveal variable modular perturbations across complex diseases.

Mapping perturbed molecular circuits that underlie complex diseases remains a great challenge. We developed a comprehensive resource of 394 cell type- and tissue-specific gene regulatory networks for human, each specifying the genome-wide connectivity among transcription factors, enhancers, promoters and genes. Integration with 37 genome-wide association studies (GWASs) showed that disease-associated genetic variants-including variants that do not reach genome-wide significance-often perturb regulatory modules that are highly specific to disease-relevant cell types or tissues. Our resource opens the door to systematic analysis of regulatory programs across hundreds of human cell types and tissues (http://regulatorycircuits.org).

Evaluation of somatic cell variants deficient in glycosylphosphatidyl-inositol anchoring as candidates for genetic correction.

Many cell surface glycoproteins are anchored in the lipid bilayer by a glycosylphosphatidyl-inositol (GPI) structure. Recently, a number of cell lines which are deficient in the biosynthesis and/or addition of this anchor have been described. In this report, we summarize the current knowledge on these lines and discuss their potential use to isolate the genes involved in the GPI anchor biosynthetic pathway with a specific emphasis on L cell fibroblasts.

Impact of genetic SLC28 transporter and ITPA variants on ribavirin serum level, hemoglobin drop and therapeutic response in patients with HCV infection.

BACKGROUND & AIMS: In the last decade, pegylated interferon-α (PegIFN-α) plus ribavirin (RBV) was the standard treatment of chronic hepatitis C for genotype 1, and it remains the standard for genotypes 2 and 3. Recent studies reported associations between RBV-induced anemia and genetic polymorphisms of concentrative nucleoside transporters such as CNT3 (encoded by SLC28A3) and inosine triphosphatase (encoded by ITPA). We aimed at studying genetic determinants of RBV kinetics, efficacy and treatment-associated anemia. METHODS: We included 216 patients from two Swiss study cohorts (61% HCV genotype 1, 39% genotypes 2 or 3). Patients were analyzed for SLC28A2 single nucleotide polymorphism (SNP) rs11854484, SLC28A3 rs56350726, and SLC28A3 rs10868138 as well as ITPA SNPs rs1127354 and rs7270101, and followed for treatment-associated hemoglobin changes and sustained virological response (SVR). In 67 patients, RBV serum levels were additionally measured during treatment. RESULTS: Patients with SLC28A2 rs11854484 genotype TT had higher dosage- and body weight-adjusted RBV levels than those with genotypes TC or CC (p=0.02 and p=0.06 at weeks 4 and 8, respectively). ITPA SNP rs1127354 was associated with hemoglobin drop ≥3 g/dl during treatment, in genotype (relative risk (RR)=2.1, 95% CI 1.3-3.5) as well as allelic analyses (RR=2.0, 95%CI 1.2-3.4). SLC28A3 rs56350726 was associated with SVR in genotype (RR=2.2; 95% CI 1.1-4.3) as well as allelic analyses (RR=2.0, 95% CI 1.1-3.4). CONCLUSIONS: The newly identified association between RBV serum levels and SLC28A2 rs11854484 genotype, as well as the replicated association of ITPA and SLC28A3 genetic polymorphisms with RBV-induced anemia and treatment response, may support individualized treatment of chronic hepatitis C and warrant further investigation in larger studies.

Combined effect of 25-OH vitamin D plasma levels and genetic vitamin D receptor (NR 1I1) variants on fibrosis progression rate in HCV patients.

BACKGROUND: Decreased vitamin D levels have been described in various forms of chronic liver disease and associated with advanced fibrosis. Whether this association is a cause or consequence of advanced fibrosis remains unclear to date. AIMS: To analyse combined effects of 25-OH vitamin D plasma levels and vitamin D receptor gene (VDR; NR1I1) polymorphisms on fibrosis progression rate in HCV patients. METHODS: 251 HCV patients underwent VDR genotyping (bat-haplotype: BsmI rs1544410 C, ApaI rs7975232 A and TaqI rs731236 A). Plasma 25-OH vitamin D levels were quantified in a subgroup of 97 patients without advanced fibrosis. The VDR haplotype and genotypes as well as plasma 25-OH vitamin D levels were associated with fibrosis progression. RESULTS: The bAt[CCA]-haplotype was significantly associated with fibrosis progression >0.101 U/year (P = 0.007; OR = 2.02) and with cirrhosis (P = 0.022; OR = 1.84). Forty-five percent of bAt[CCA]-haplotype patients were rapid fibrosers, 21.1% were cirrhotic. Likewise, ApaI rs7975232 CC genotype was significantly associated with fibrosis progression and cirrhosis. Lower plasma 25-OH vitamin D levels were significantly associated with fibrosis progression >0.101 U/year in F0-2 patients (P = 0.013). Combined analysis of both variables revealed a highly significant additive effect on fibrosis progression with 45.5% rapid fibrosers for bAt[CCA]-haplotype and 25-OH vitamin D < 20 μg/L compared with only 9.1% for the most favourable combination (P = 0.006). In multivariate analysis, the bAt-haplotype was an independent risk factor for fibrosis progression (P = 0.001; OR = 2.83). CONCLUSION: Low 25-OH vitamin D plasma levels and the unfavourable VDR bAt[CCA]-haplotype are associated with rapid fibrosis progression in chronic HCV patients. In combination, both variables exert significant additive effects on fibrosis progression.

Clear Cell Papillary Renal Cell Carcinoma and Renal Angiomyoadenomatous Tumor: Two Variants of a Morphologic, Immunohistochemical, and Genetic Distinct Entity of Renal Cell Carcinoma.

Clear cell papillary renal cell carcinoma (ccpRCC) and renal angiomyoadenomatous tumor (RAT) share morphologic similarities with clear cell (ccRCC) and papillary RCC (pRCC). It is a matter of controversy whether their morphologic, immunophenotypic, and molecular features allow the definition of a separate renal carcinoma entity. The aim of our project was to investigate specific renal immunohistochemical biomarkers involved in the hypoxia-inducible factor pathway and mutations in the VHL gene to clarify the relationship between ccpRCC and RAT. We investigated 28 ccpRCC and 9 RAT samples by immunohistochemistry using 25 markers. VHL gene mutations and allele losses were investigated by Sanger sequencing and fluorescence in situ hybridization. Clinical follow-up data were obtained for a subset of the patients. No tumor recurrence or tumor-related death was observed in any of the patients. Immunohistochemistry and molecular analyses led to the reclassification of 3 tumors as ccRCC and TFE3 translocation carcinomas. The immunohistochemical profile of ccpRCC and RAT samples was very similar but not identical, differing from both ccRCC and pRCC. Especially, the parafibromin and hKIM-1 expression exhibited differences in ccpRCC/RAT compared with ccRCC and pRCC. Genetic analysis revealed VHL mutations in 2/27 (7%) and 1/7 (14%) ccpRCC and RAT samples, respectively. Fluorescence in situ hybridization analysis disclosed a 3p loss in 2/20 (10%) ccpRCC samples. ccpRCC and RAT have a specific morphologic and immunohistochemical profile, but they share similarities with the more aggressive renal tumors. On the basis of our results, we regard ccpRCC/RAT as a distinct entity of RCCs.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

Rare genomic structural variants in complex disease: lessons from the replication of associations with obesity.

The limited ability of common variants to account for the genetic contribution to complex disease has prompted searches for rare variants of large effect, to partly explain the 'missing heritability'. Analyses of genome-wide genotyping data have identified genomic structural variants (GSVs) as a source of such rare causal variants. Recent studies have reported multiple GSV loci associated with risk of obesity. We attempted to replicate these associations by similar analysis of two familial-obesity case-control cohorts and a population cohort, and detected GSVs at 11 out of 18 loci, at frequencies similar to those previously reported. Based on their reported frequencies and effect sizes (OR≥25), we had sufficient statistical power to detect the large majority (80%) of genuine associations at these loci. However, only one obesity association was replicated. Deletion of a 220 kb region on chromosome 16p11.2 has a carrier population frequency of 2×10(-4) (95% confidence interval [9.6×10(-5)-3.1×10(-4)]); accounts overall for 0.5% [0.19%-0.82%] of severe childhood obesity cases (P = 3.8×10(-10); odds ratio = 25.0 [9.9-60.6]); and results in a mean body mass index (BMI) increase of 5.8 kg.m(-2) [1.8-10.3] in adults from the general population. We also attempted replication using BMI as a quantitative trait in our population cohort; associations with BMI at or near nominal significance were detected at two further loci near KIF2B and within FOXP2, but these did not survive correction for multiple testing. These findings emphasise several issues of importance when conducting rare GSV association, including the need for careful cohort selection and replication strategy, accurate GSV identification, and appropriate correction for multiple testing and/or control of false discovery rate. Moreover, they highlight the potential difficulty in replicating rare CNV associations across different populations. Nevertheless, we show that such studies are potentially valuable for the identification of variants making an appreciable contribution to complex disease.

Physical activity attenuates the influence of FTO variants on obesity risk: a meta-analysis of 218,166 adults and 19,268 children.

BACKGROUND: The FTO gene harbors the strongest known susceptibility locus for obesity. While many individual studies have suggested that physical activity (PA) may attenuate the effect of FTO on obesity risk, other studies have not been able to confirm this interaction. To confirm or refute unambiguously whether PA attenuates the association of FTO with obesity risk, we meta-analyzed data from 45 studies of adults (n = 218,166) and nine studies of children and adolescents (n = 19,268). METHODS AND FINDINGS: All studies identified to have data on the FTO rs9939609 variant (or any proxy [r(2)>0.8]) and PA were invited to participate, regardless of ethnicity or age of the participants. PA was standardized by categorizing it into a dichotomous variable (physically inactive versus active) in each study. Overall, 25% of adults and 13% of children were categorized as inactive. Interaction analyses were performed within each study by including the FTO×PA interaction term in an additive model, adjusting for age and sex. Subsequently, random effects meta-analysis was used to pool the interaction terms. In adults, the minor (A-) allele of rs9939609 increased the odds of obesity by 1.23-fold/allele (95% CI 1.20-1.26), but PA attenuated this effect (p(interaction)  = 0.001). More specifically, the minor allele of rs9939609 increased the odds of obesity less in the physically active group (odds ratio  = 1.22/allele, 95% CI 1.19-1.25) than in the inactive group (odds ratio  = 1.30/allele, 95% CI 1.24-1.36). No such interaction was found in children and adolescents. CONCLUSIONS: The association of the FTO risk allele with the odds of obesity is attenuated by 27% in physically active adults, highlighting the importance of PA in particular in those genetically predisposed to obesity.

Response to antiretroviral treatment in HIV-1-infected individuals with allelic variants of the multidrug resistance transporter 1: a pharmacogenetics study.

BACKGROUND:HIV-1-infected patients vary considerably by their response to antiretroviral treatment, drug concentrations in plasma, toxic events, and rate of immune recovery. This variability could have a genetic basis. We did a pharmacogenetics study to analyse the association between response to antiretroviral treatment and allelic variants of several genes. METHODS:In 123 patients, we did PCR analyses of the gene for the multidrug-resistance transporter (MDR1), which codes for P-glycoprotein, of genes coding for isoenzymes of cytochrome P450, CYP3A4, CYP3A5, CYP2D6, and CYP2C19, and of the gene for the chemokine receptor CCR5. We measured concentrations in plasma of the antiretroviral agents efavirenz and nelfinavir by high-performance liquid-chromatography, and measured levels of P-glycoprotein expression, CD4-cell count, and HIV-1 viraemia. FINDINGS: Median drug concentrations in patients with the MDR1 3435 TT, CT, and CC genotypes were at the 30th, 50th, and 75th percentiles, respectively (p=0.0001). In patients with CYP2D6 extensive-metaboliser or poor-metaboliser alleles, median drug concentrations were at percentiles 45 and 62.5, respectively (p=0.04). Patients with the MDR1 TT genotype 6 months after starting treatment had a greater rise in CD4-cell count (257 cells/microL) than patients with the CT (165 cells/microL) and CC (121 cells/microL) genotype (p=0.0048), and the best recovery of naïve CD4-cells. INTERPRETATION:The polymorphism MDR1 3435 C/T predicts immune recovery after initiation of antiretroviral treatment. This finding suggests that P-glycoprotein has an important role in admittance of antiretroviral drugs to restricted compartments in vivo.

Genetic architecture of ambulatory blood pressure in the general population: insights from cardiovascular gene-centric array.

Genetic determinants of blood pressure are poorly defined. We undertook a large-scale, gene-centric analysis to identify loci and pathways associated with ambulatory systolic and diastolic blood pressure. We measured 24-hour ambulatory blood pressure in 2020 individuals from 520 white European nuclear families (the Genetic Regulation of Arterial Pressure of Humans in the Community Study) and genotyped their DNA using the Illumina HumanCVD BeadChip array, which contains ≈50 000 single nucleotide polymorphisms in >2000 cardiovascular candidate loci. We found a strong association between rs13306560 polymorphism in the promoter region of MTHFR and CLCN6 and mean 24-hour diastolic blood pressure; each minor allele copy of rs13306560 was associated with 2.6 mm Hg lower mean 24-hour diastolic blood pressure (P=1.2×10(-8)). rs13306560 was also associated with clinic diastolic blood pressure in a combined analysis of 8129 subjects from the Genetic Regulation of Arterial Pressure of Humans in the Community Study, the CoLaus Study, and the Silesian Cardiovascular Study (P=5.4×10(-6)). Additional analysis of associations between variants in gene ontology-defined pathways and mean 24-hour blood pressure in the Genetic Regulation of Arterial Pressure of Humans in the Community Study showed that cell survival control signaling cascades could play a role in blood pressure regulation. There was also a significant overrepresentation of rare variants (minor allele frequency: <0.05) among polymorphisms showing at least nominal association with mean 24-hour blood pressure indicating that a considerable proportion of its heritability may be explained by uncommon alleles. Through a large-scale gene-centric analysis of ambulatory blood pressure, we identified an association of a novel variant at the MTHFR/CLNC6 locus with diastolic blood pressure and provided new insights into the genetic architecture of blood pressure.

Orosomucoid system: 17 additional orosomucoid variants and proposal for a new nomenclature.

There are two forms of orosomucoid (ORM) in the sera of most individuals. They are encoded by two separate but closely linked loci, ORM1 and ORM2. A number of variants have been identified in various populations. Duplication and nonexpression are also observed in some populations. Thus, the ORM system is very complicated and its nomenclature is very confusing. In order to propose a new nomenclature, ORM variants detected by several laboratories have been compared and characterized by isoelectric focusing (IEF) followed by immunoprinting. A total of 57 different alleles including 17 new ones were identified. The 27 alleles were assigned to the ORM1 locus, and the others to the ORM2 locus. The designations ORM*F1, ORM1*F2, ORM1*S and ORM2*M were adopted for the four common alleles instead of ORM1*1, ORM1*3, ORM1*2 and ORM2*1 (ORM2*A), respectively. The variants were designated alpha numerically according to their relative mobilities after IEF in a pH gradient of 4.5-5.4 with Triton X-100 and glycerol. For the duplicated genes a prefix is added to a combined name of two alleles, e.g. ORM1*dB9S. Silent alleles were named ORM1*Q0 and ORM2*Q0 conventionally. In addition, the effects of diseases to ORM band patterns after IEF are also discussed.