There is growing evidence of shared risk alleles for complex traits (pleiotropy), including autoimmune and neuropsychiatric diseases. This might be due to sharing among all individuals (whole-group pleiotropy) or a subset of individuals in a genetically heterogeneous cohort (subgroup heterogeneity). Here we describe the use of a well-powered statistic, BUHMBOX, to distinguish between those two situations using genotype data. We observed a shared genetic basis for 11 autoimmune diseases and type 1 diabetes (T1D; P < 1 × 10(-4)) and for 11 autoimmune diseases and rheumatoid arthritis (RA; P < 1 × 10(-3)). This sharing was not explained by subgroup heterogeneity (corrected PBUHMBOX > 0.2; 6,670 T1D cases and 7,279 RA cases). Genetic sharing between seronegative and seropostive RA (P < 1 × 10(-9)) had significant evidence of subgroup heterogeneity, suggesting a subgroup of seropositive-like cases within seronegative cases (PBUHMBOX = 0.008; 2,406 seronegative RA cases). We also observed a shared genetic basis for major depressive disorder (MDD) and schizophrenia (P < 1 × 10(-4)) that was not explained by subgroup heterogeneity (PBUHMBOX = 0.28; 9,238 MDD cases).

Meta-analysis strategies have become critical to augment power of genome-wide association studies (GWAS). To reduce genotyping or sequencing cost, many studies today utilize shared controls, and these individuals can inadvertently overlap among multiple studies. If these overlapping individuals are not taken into account in meta-analysis, they can induce spurious associations. In this article, we propose a general framework for adjusting association statistics to account for overlapping subjects within a meta-analysis. The key idea of our method is to transform the covariance structure of the data, so it can be used in downstream analyses. As a result, the strategy is very flexible and allows a wide range of meta-analysis methods, such as the random effects model, to account for overlapping subjects. Using simulations and real datasets, we demonstrate that our method has utility in meta-analyses of GWAS, as well as in a multi-tissue mouse expression quantitative trait loci (eQTL) study where our method increases the number of discovered eQTL by up to 19% compared with existing methods.

Genome-wide strategies have driven the discovery of more than 300 susceptibility loci for autoimmune diseases. However, for almost all loci, understanding of the mechanisms leading to autoimmunity remains limited, and most variants that are likely to be causal are in non-coding regions of the genome. A critical next step will be to identify the in vivo and ex vivo immunophenotypes that are affected by risk variants. To do this, key cell types and cell states that are implicated in autoimmune diseases will need to be defined. Functional genomic annotations from these cell types and states can then be used to resolve candidate genes and causal variants. Together with longitudinal studies, this approach may yield pivotal insights into how autoimmunity is triggered.

There has been intense debate over the immunological basis of schizophrenia, and the potential utility of adjunct immunotherapies. The major histocompatibility complex is consistently the most powerful region of association in genome-wide association studies (GWASs) of schizophrenia and has been interpreted as strong genetic evidence supporting the immune hypothesis. However, global pathway analyses provide inconsistent evidence of immune involvement in schizophrenia, and it remains unclear whether genetic data support an immune etiology per se. Here we empirically test the hypothesis that variation in immune genes contributes to schizophrenia. We show that there is no enrichment of immune loci outside of the MHC region in the largest genetic study of schizophrenia conducted to date, in contrast to 5 diseases of known immune origin. Among 108 regions of the genome previously associated with schizophrenia, we identify 6 immune candidates (DPP4, HSPD1, EGR1, CLU, ESAM, NFATC3) encoding proteins with alternative, nonimmune roles in the brain. While our findings do not refute evidence that has accumulated in support of the immune hypothesis, they suggest that genetically mediated alterations in immune function may not play a major role in schizophrenia susceptibility. Instead, there may be a role for pleiotropic effects of a small number of immune genes that also regulate brain development and plasticity. Whether immune alterations drive schizophrenia progression is an important question to be addressed by future research, especially in light of the growing interest in applying immunotherapies in schizophrenia.

OBJECTIVES: Our objective was to estimate the risk of developing rheumatoid arthritis (RA) associated with a family history of non-RA arthritis-related diseases. This familial co-aggregation is of clinical interest since it is often encountered when assessing family history of RA specifically, but also informative on the genetic overlap between these diseases. Since anticitrullinated peptide antibodies/rheumatoid factor (RF)-positive and RF-negative RA have both specific and shared genetic factors, the familial co-aggregation was assessed separately for seropositive and seronegative disease. METHODS: Nested case-control study in prospectively recorded Swedish total population data. The Multi-Generation Register identified first-degree relatives. RA and arthritis-related diseases were ascertained through the nationwide patient register. RA serology was based on International Classification of Diseases tenth revision coded diagnoses, mainly reflecting RF. Familial risks were calculated using conditional logistic regression. Results were replicated using the Swedish rheumatology register. RESULTS: Familial co-aggregation was found between RA and every studied arthritis-related disease, but the magnitude varied widely, from juvenile idiopathic arthritis (JIA) (seropositive RA OR=3.98 (3.01 to 5.26); seronegative RA OR=5.70 (3.47 to 9.36)) to osteoarthritis (seropositive RA OR=1.03 (1.00 to 1.06); seronegative RA OR=1.05 (1.00 to 1.09)). The familial co-aggregation pattern of non-RA arthritis-related diseases was overall similar for seropositive and seronegative RA. Among those with family history of RA, relatives' other arthritis-related diseases conferred little or no additional risk. CONCLUSIONS: Although family history of several arthritis-related diseases may be useful to predict RA (eg, lupus and JIA), others (eg, osteoarthritis and arthralgia) are less useful. Seropositive and seronegative RA had rather similar familial co-aggregation patterns with arthritis-related diseases, suggesting that the two RA subsets are similar in the genetic factors that overlap with these diseases.

OBJECTIVES: A recent study identified 16 genetic variants associated with N-glycosylation of human IgG. Several of the genomic regions where these single nucleotide polymorphisms (SNPs) reside have also been associated with autoimmune disease (AID) susceptibility, suggesting there may be pleiotropy (genetic sharing) between loci controlling both N-glycosylation and AIDs. We investigated this by testing variants associated with levels of IgG N-glycosylation for association with rheumatoid arthritis (RA) susceptibility using a Mendelian randomisation study, and testing a subset of these variants in a less well-powered study of treatment response and severity. METHODS: SNPs showing association with IgG N-glycosylation were analysed for association with RA susceptibility in 14 361 RA cases and 43 923 controls. Five SNPs were tested for association with response to anti-tumour necrosis factor (TNF) therapy in 1081 RA patient samples and for association with radiological disease severity in 342 patients. RESULTS: Only one SNP (rs9296009) associated with N-glycosylation showed an association (p=6.92×10(-266)) with RA susceptibility, although this was due to linkage disequilibrium with causal human leukocyte antigen (HLA) variants. Four regions of the genome harboured SNPs associated with both traits (shared loci); although statistical analysis indicated that the associations observed for the two traits are independent. No SNPs showed association with response to anti-TNF therapy. One SNP rs12342831 was modestly associated with Larsen score (p=0.05). CONCLUSIONS: In a large, well-powered cohort of RA patients, we show SNPs driving levels of N-glycosylation have no association with RA susceptibility, indicating colocalisation of associated SNPs are not necessarily indicative of a shared genetic background or a role for glycosylation in disease susceptibility.

Hundreds of genomic loci have been associated with a significant number of immune-mediated diseases, and a large proportion of these associated loci are shared among traits. Both the molecular mechanisms by which these loci confer disease susceptibility and the extent to which shared loci are implicated in a common pathogenesis are unknown. We therefore sought to dissect the functional components at loci shared between two autoimmune diseases: coeliac disease (CeD) and rheumatoid arthritis (RA). We used a cohort of 12 381 CeD cases and 7827 controls, and another cohort of 13 819 RA cases and 12 897 controls, all genotyped with the Immunochip platform. In the joint analysis, we replicated 19 previously identified loci shared by CeD and RA and discovered five new non-HLA loci shared by CeD and RA. Our fine-mapping results indicate that in nine of 24 shared loci the associated variants are distinct in the two diseases. Using cell-type-specific histone markers, we observed that loci which pointed to the same variants in both diseases were enriched for marks of promoters active in CD14+ and CD34+ immune cells (P < 0.001), while loci pointing to distinct variants in one of the two diseases showed enrichment for marks of more specialized cell types, like CD4+ regulatory T cells in CeD (P < 0.0001) compared with Th17 and CD15+ in RA (P = 0.0029).

Although genome-wide association studies have identified over 100 risk loci that explain ∼33% of familial risk for prostate cancer (PrCa), their functional effects on risk remain largely unknown. Here we use genotype data from 59,089 men of European and African American ancestries combined with cell-type-specific epigenetic data to build a genomic atlas of single-nucleotide polymorphism (SNP) heritability in PrCa. We find significant differences in heritability between variants in prostate-relevant epigenetic marks defined in normal versus tumour tissue as well as between tissue and cell lines. The majority of SNP heritability lies in regions marked by H3k27 acetylation in prostate adenoc7arcinoma cell line (LNCaP) or by DNaseI hypersensitive sites in cancer cell lines. We find a high degree of similarity between European and African American ancestries suggesting a similar genetic architecture from common variation underlying PrCa risk. Our findings showcase the power of integrating functional annotation with genetic data to understand the genetic basis of PrCa.

BACKGROUND: Myocardial mass is a key determinant of cardiac muscle function and hypertrophy. Myocardial depolarization leading to cardiac muscle contraction is reflected by the amplitude and duration of the QRS complex on the electrocardiogram (ECG). Abnormal QRS amplitude or duration reflect changes in myocardial mass and conduction, and are associated with increased risk of heart failure and death. OBJECTIVES: This meta-analysis sought to gain insights into the genetic determinants of myocardial mass. METHODS: We carried out a genome-wide association meta-analysis of 4 QRS traits in up to 73,518 individuals of European ancestry, followed by extensive biological and functional assessment. RESULTS: We identified 52 genomic loci, of which 32 are novel, that are reliably associated with 1 or more QRS phenotypes at p < 1 × 10(-8). These loci are enriched in regions of open chromatin, histone modifications, and transcription factor binding, suggesting that they represent regions of the genome that are actively transcribed in the human heart. Pathway analyses provided evidence that these loci play a role in cardiac hypertrophy. We further highlighted 67 candidate genes at the identified loci that are preferentially expressed in cardiac tissue and associated with cardiac abnormalities in Drosophila melanogaster and Mus musculus. We validated the regulatory function of a novel variant in the SCN5A/SCN10A locus in vitro and in vivo. CONCLUSIONS: Taken together, our findings provide new insights into genes and biological pathways controlling myocardial mass and may help identify novel therapeutic targets.