Genome-wide association studies (GWAS) have detected many disease associations. However, the reported variants tend to explain small fractions of risk, and there are doubts about issues such as the ...portability of findings over different ethnic groups or the relative roles of rare versus common variants in the genetic architecture of complex disease. Studying the degree of sharing of disease-associated variants across populations can help in solving these issues. We present a comprehensive survey of GWAS replicability across 28 diseases. Most loci and SNPs discovered in Europeans for these conditions have been extensively replicated using peoples of European and East Asian ancestry, while the replication with individuals of African ancestry is much less common. We found a strong and significant correlation of Odds Ratios across Europeans and East Asians, indicating that underlying causal variants are common and shared between the two ancestries. Moreover, SNPs that failed to replicate in East Asians map into genomic regions where Linkage Disequilibrium patterns differ significantly between populations. Finally, we observed that GWAS with larger sample sizes have detected variants with weaker effects rather than with lower frequencies. Our results indicate that most GWAS results are due to common variants. In addition, the sharing of disease alleles and the high correlation in their effect sizes suggest that most of the underlying causal variants are shared between Europeans and East Asians and that they tend to map close to the associated marker SNPs.
Since the publication of the Wellcome Trust Case Control Consortium (WTCCC) landmark study a decade ago, genome-wide association studies (GWAS) have led to the discovery of thousands of risk variants ...involved in disease etiology. This success story has two angles that are often overlooked. First, GWAS findings are highly replicable. This is an unprecedented phenomenon in complex trait genetics, and indeed in many areas of science, which in past decades have been plagued by false positives. At a time of increasing concerns about the lack of reproducibility, we examine the biological and methodological reasons that account for the replicability of GWAS and identify the challenges ahead. In contrast to the exemplary success of disease gene discovery, at present GWAS findings are not useful for predicting phenotypes. We close with an overview of the prospects for individualized prediction of disease risk and its foreseeable impact in clinical practice.
The decade of GWAS constitutes a clear improvement in the recent history of reproducibility in medical genetic research. The extremely high rates of replication imply that, for the first time, findings can be trusted.
Large numbers of false positives and the tiny effect size of genetic risk variants induced a change in incentives during the GWAS era, priming the requirement for large sample sizes and a culture of data sharing.
There is increasing interest in and need of new methodologies to better understand the genetic architecture of complex traits.
It is necessary to keep fostering a culture of compulsory replication to maintain the current high reliability in findings.
Although the success of GWAS has not translated into an ability to predict phenotypes based on genetic markers, polygenic and transcriptional risk scores (PRSs and TRSs) hold potential for stratification according to risk.
Systemic Lupus Erythematosus (SLE) is characterized by B cells lacking IgD and CD27 (double negative; DN). We show that DN cell expansions reflected a subset of CXCR5- CD11c+ cells (DN2) representing ...pre-plasma cells (PC). DN2 cells predominated in African-American patients with active disease and nephritis, anti-Smith and anti-RNA autoantibodies. They expressed a T-bet transcriptional network; increased Toll-like receptor-7 (TLR7); lacked the negative TLR regulator TRAF5; and were hyper-responsive to TLR7. DN2 cells shared with activated naive cells (aNAV), phenotypic and functional features, and similar transcriptomes. Their PC differentiation and autoantibody production was driven by TLR7 in an interleukin-21 (IL-21)-mediated fashion. An in vivo developmental link between aNAV, DN2 cells, and PC was demonstrated by clonal sharing. This study defines a distinct differentiation fate of autoreactive naive B cells into PC precursors with hyper-responsiveness to innate stimuli, as well as establishes prominence of extra-follicular B cell activation in SLE, and identifies therapeutic targets.
Display omitted
•Autoreactive CD27- IgD- CXCR5- CD11c+ (DN2) B cells expand in lupus patients•DN2 cells derive from naive cells and are poised to generate plasmablasts•DN2 B cells are hyper-responsive to Toll-like receptor-7 signaling•The properties of SLE DN2 B cells stem from distinct transcriptional networks
The role of extrafollicular B cells in human systemic lupus is unknown. Jenks et al. define the main components of this pathway and its prominence in severe disease. Its activation is mediated by hyper-responsiveness to Toll-like receptor-7 and leads to the generation of autoreactive antibody-secreting plasmablasts.
Summary Background Stricturing and penetrating complications account for substantial morbidity and health-care costs in paediatric and adult onset Crohn's disease. Validated models to predict risk ...for complications are not available, and the effect of treatment on risk is unknown. Methods We did a prospective inception cohort study of paediatric patients with newly diagnosed Crohn's disease at 28 sites in the USA and Canada. Genotypes, antimicrobial serologies, ileal gene expression, and ileal, rectal, and faecal microbiota were assessed. A competing-risk model for disease complications was derived and validated in independent groups. Propensity-score matching tested the effect of anti-tumour necrosis factor α (TNFα) therapy exposure within 90 days of diagnosis on complication risk. Findings Between Nov 1, 2008, and June 30, 2012, we enrolled 913 patients, 78 (9%) of whom experienced Crohn's disease complications. The validated competing-risk model included age, race, disease location, and antimicrobial serologies and provided a sensitivity of 66% (95% CI 51–82) and specificity of 63% (55–71), with a negative predictive value of 95% (94–97). Patients who received early anti-TNFα therapy were less likely to have penetrating complications (hazard ratio HR 0·30, 95% CI 0·10–0·89; p=0·0296) but not stricturing complication (1·13, 0·51–2·51; 0·76) than were those who did not receive early anti-TNFα therapy. Ruminococcus was implicated in stricturing complications and Veillonella in penetrating complications. Ileal genes controlling extracellular matrix production were upregulated at diagnosis, and this gene signature was associated with stricturing in the risk model (HR 1·70, 95% CI 1·12–2·57; p=0·0120). When this gene signature was included, the model's specificity improved to 71%. Interpretation Our findings support the usefulness of risk stratification of paediatric patients with Crohn's disease at diagnosis, and selection of anti-TNFα therapy. Funding Crohn's and Colitis Foundation of America, Cincinnati Children's Hospital Research Foundation Digestive Health Center.
Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and ...trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.
The Romani, the largest European minority group with approximately 11 million people 1, constitute a mosaic of languages, religions, and lifestyles while sharing a distinct social heritage. ...Linguistic 2 and genetic 3–8 studies have located the Romani origins in the Indian subcontinent. However, a genome-wide perspective on Romani origins and population substructure, as well as a detailed reconstruction of their demographic history, has yet to be provided. Our analyses based on genome-wide data from 13 Romani groups collected across Europe suggest that the Romani diaspora constitutes a single initial founder population that originated in north/northwestern India ∼1.5 thousand years ago (kya). Our results further indicate that after a rapid migration with moderate gene flow from the Near or Middle East, the European spread of the Romani people was via the Balkans starting ∼0.9 kya. The strong population substructure and high levels of homozygosity we found in the European Romani are in line with genetic isolation as well as differential gene flow in time and space with non-Romani Europeans. Overall, our genome-wide study sheds new light on the origins and demographic history of European Romani.
► The Romani diaspora originated in north/northwest India around 1.5 kya ► The European spread of the Romani people was via the Balkans starting ∼0.9 kya ► Footprints of bottlenecks, endogamy, and admixture in the European Romani genomes
Type 2 diabetes (T2D) is a very common disease in humans. Here we conduct a meta-analysis of genome-wide association studies (GWAS) with ~16 million genetic variants in 62,892 T2D cases and 596,424 ...controls of European ancestry. We identify 139 common and 4 rare variants associated with T2D, 42 of which (39 common and 3 rare variants) are independent of the known variants. Integration of the gene expression data from blood (n = 14,115 and 2765) with the GWAS results identifies 33 putative functional genes for T2D, 3 of which were targeted by approved drugs. A further integration of DNA methylation (n = 1980) and epigenomic annotation data highlight 3 genes (CAMK1D, TP53INP1, and ATP5G1) with plausible regulatory mechanisms, whereby a genetic variant exerts an effect on T2D through epigenetic regulation of gene expression. Our study uncovers additional loci, proposes putative genetic regulatory mechanisms for T2D, and provides evidence of purifying selection for T2D-associated variants.
Epstein-Barr virus (EBV), human herpes virus 4, has been classically associated with infectious mononucleosis, multiple sclerosis and several types of cancers. Many of these diseases show marked ...geographical differences in prevalence, which points to underlying genetic and/or environmental factors. Those factors may include a different susceptibility to EBV infection and viral copy number among human populations. Since EBV is commonly used to transform B-cells into lymphoblastoid cell lines (LCLs) we hypothesize that differences in EBV copy number among individual LCLs may reflect differential susceptibility to EBV infection. To test this hypothesis, we retrieved whole-genome sequenced EBV-mapping reads from 1,753 LCL samples derived from 19 populations worldwide that were sequenced within the context of the 1000 Genomes Project. An in silico methodology was developed to estimate the number of EBV copy number in LCLs and validated these estimations by real-time PCR. After experimentally confirming that EBV relative copy number remains stable over cell passages, we performed a genome wide association analysis (GWAS) to try detecting genetic variants of the host that may be associated with EBV copy number. Our GWAS has yielded several genomic regions suggestively associated with the number of EBV genomes per cell in LCLs, unraveling promising candidate genes such as CAND1, a known inhibitor of EBV replication. While this GWAS does not unequivocally establish the degree to which genetic makeup of individuals determine viral levels within their derived LCLs, for which a larger sample size will be needed, it potentially highlighted human genes affecting EBV-related processes, which constitute interesting candidates to follow up in the context of EBV related pathologies.
Gene expression profiling can be used to uncover the mechanisms by which loci identified through genome-wide association studies (GWAS) contribute to pathology. Given that most GWAS hits are in ...putative regulatory regions and transcript abundance is physiologically closer to the phenotype of interest, we hypothesized that summation of risk-allele-associated gene expression, namely a transcriptional risk score (TRS), should provide accurate estimates of disease risk. We integrate summary-level GWAS and expression quantitative trait locus (eQTL) data with RNA-seq data from the RISK study, an inception cohort of pediatric Crohn's disease. We show that TRSs based on genes regulated by variants linked to inflammatory bowel disease (IBD) not only outperform genetic risk scores (GRSs) in distinguishing Crohn's disease from healthy samples, but also serve to identify patients who in time will progress to complicated disease. Our dissection of eQTL effects may be used to distinguish genes whose association with disease is through promotion versus protection, thereby linking statistical association to biological mechanism. The TRS approach constitutes a potential strategy for personalized medicine that enhances inference from static genotypic risk assessment.
In order to evaluate whether rare regulatory variants in the vicinity of promoters are likely to impact gene expression, we conducted a novel burden test for enrichment of rare variants at the ...extremes of expression. After sequencing 2-kb promoter regions of 472 genes in 410 healthy adults, we performed a quadratic regression of rare variant count on bins of peripheral blood transcript abundance from microarrays, summing over ranks of all genes. After adjusting for common eQTLs and the major axes of gene expression covariance, a highly significant excess of variants with minor allele frequency less than 0.05 at both high and low extremes across individuals was observed. Further enrichment was seen in sites annotated as potentially regulatory by RegulomeDB, but a deficit of effects was associated with known metabolic disease genes. The main result replicates in an independent sample of 75 individuals with RNA-seq and whole-genome sequence information. Three of four predicted large-effect sites were validated by CRISPR/Cas9 knockdown in K562 cells, but simulations indicate that effect sizes need not be unusually large to produce the observed burden. Unusually divergent low-frequency promoter haplotypes were observed at 31 loci, at least 9 of which appear to be derived from Neandertal admixture, but these were not associated with divergent gene expression in blood. The overall burden test results are consistent with rare and private regulatory variants driving high or low transcription at specific loci, potentially contributing to disease.