Mexican Biobank advances population and medical genomics of diverse ancestries.

Latin America continues to be severely underrepresented in genomics research, and fine-scale genetic histories and complex trait architectures remain hidden owing to insufficient data1. To fill this gap, the Mexican Biobank project genotyped 6,057 individuals from 898 rural and urban localities across all 32 states in Mexico at a resolution of 1.8 million genome-wide markers with linked complex trait and disease information creating a valuable nationwide genotype-phenotype database. Here, using ancestry deconvolution and inference of identity-by-descent segments, we inferred ancestral population sizes across Mesoamerican regions over time, unravelling Indigenous, colonial and postcolonial demographic dynamics2-6. We observed variation in runs of homozygosity among genomic regions with different ancestries reflecting distinct demographic histories and, in turn, different distributions of rare deleterious variants. We conducted genome-wide association studies (GWAS) for 22 complex traits and found that several traits are better predicted using the Mexican Biobank GWAS compared to the UK Biobank GWAS7,8. We identified genetic and environmental factors associating with trait variation, such as the length of the genome in runs of homozygosity as a predictor for body mass index, triglycerides, glucose and height. This study provides insights into the genetic histories of individuals in Mexico and dissects their complex trait architectures, both crucial for making precision and preventive medicine initiatives accessible worldwide.

A common NFKB1 variant detected through antibody analysis in UK Biobank predicts risk of infection and allergy.

Infectious agents contribute significantly to the global burden of diseases through both acute infection and their chronic sequelae. We leveraged the UK Biobank to identify genetic loci that influence humoral immune response to multiple infections. From 45 genome-wide association studies in 9,611 participants from UK Biobank, we identified NFKB1 as a locus associated with quantitative antibody responses to multiple pathogens, including those from the herpes, retro-, and polyoma-virus families. An insertion-deletion variant thought to affect NFKB1 expression (rs28362491), was mapped as the likely causal variant and could play a key role in regulation of the immune response. Using 121 infection- and inflammation-related traits in 487,297 UK Biobank participants, we show that the deletion allele was associated with an increased risk of infection from diverse pathogens but had a protective effect against allergic disease. We propose that altered expression of NFKB1, as a result of the deletion, modulates hematopoietic pathways and likely impacts cell survival, antibody production, and inflammation. Taken together, we show that disruptions to the tightly regulated immune processes may tip the balance between exacerbated immune responses and allergy, or increased risk of infection and impaired resolution of inflammation.

The UK Biobank resource with deep phenotyping and genomic data.

The UK Biobank project is a prospective cohort study with deep genetic and phenotypic data collected on approximately 500,000 individuals from across the United Kingdom, aged between 40 and 69 at recruitment. The open resource is unique in its size and scope. A rich variety of phenotypic and health-related information is available on each participant, including biological measurements, lifestyle indicators, biomarkers in blood and urine, and imaging of the body and brain. Follow-up information is provided by linking health and medical records. Genome-wide genotype data have been collected on all participants, providing many opportunities for the discovery of new genetic associations and the genetic bases of complex traits. Here we describe the centralized analysis of the genetic data, including genotype quality, properties of population structure and relatedness of the genetic data, and efficient phasing and genotype imputation that increases the number of testable variants to around 96 million. Classical allelic variation at 11 human leukocyte antigen genes was imputed, resulting in the recovery of signals with known associations between human leukocyte antigen alleles and many diseases.

Fine-mapping inflammatory bowel disease loci to single-variant resolution.

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Common allotypes of ER aminopeptidase 1 have substrate-dependent and highly variable enzymatic properties.

Polymorphic variation of immune system proteins can drive variability of individual immune responses. Endoplasmic reticulum aminopeptidase 1 (ERAP1) generates antigenic peptides for presentation by major histocompatibility complex class I molecules. Coding SNPs in ERAP1 have been associated with predisposition to inflammatory rheumatic disease and shown to affect functional properties of the enzyme, but the interplay between combinations of these SNPs as they exist in allotypes has not been thoroughly explored. We used phased genotype data to estimate ERAP1 allotype frequency in 2504 individuals across five major human populations, generated highly pure recombinant enzymes corresponding to the ten most common ERAP1 allotypes, and systematically characterized their in vitro enzymatic properties. We find that ERAP1 allotypes possess a wide range of enzymatic activities, up to 60-fold, whose ranking is substrate dependent. Strikingly, allotype 10, previously associated with Behçet's disease, is consistently a low-activity outlier, suggesting that a significant percentage of individuals carry a subactive ERAP1 gene. Enzymatic analysis revealed that ERAP1 allotypes can differ in both catalytic efficiency and substrate affinity, differences that can change intermediate accumulation in multistep trimming reactions. Alterations in efficacy of an allosteric inhibitor that targets the regulatory site suggest that allotypic variation influences the communication between the regulatory and the active site. Our work defines the wide landscape of ERAP1 activity in human populations and demonstrates how common allotypes can induce substrate-dependent variability in antigen processing, thus contributing, in synergy with major histocompatibility complex haplotypes, to immune response variability and predisposition to chronic inflammatory conditions.

Correction: Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

[This corrects the article DOI: 10.1371/journal.pgen.1006866.].

ERAP1 association with ankylosing spondylitis is attributable to common genotypes rather than rare haplotype combinations.

We investigated the proposal that ankylosing spondylitis (AS) is associated with unusual ERAP1 genotypes. ERAP1 haplotypes were constructed for 213 AS cases and 46 rheumatoid arthritis controls using family data. Haplotypes were generated from five common ERAP1 single nucleotide polymorphisms (SNPs)-rs2287987 (M349V), rs30187 (K528R), rs10050860 (D575N), rs17482078 (R725Q), and rs27044 (Q730E). Haplotype frequencies were compared using Fisher's exact test. ERAP1 haplotypes imputed from the International Genetics of AS Consortium (IGAS) Immunochip study were also studied. In the family study, we identified only four common ERAP1 haplotypes ("VRNQE," "MKDRQ," "MRDRE," and "MKDRE") in both AS cases and controls apart from two rare (<0.5%) previously unreported haplotypes. There were no examples of the unusual ERAP1 haplotype combination ("*001/*005") previously reported by others in 53% of AS cases. As expected, K528-bearing haplotypes were increased in the AS family study (AS 43% vs. control 35%), due particularly to an increase in the MKDRQ haplotype (AS 35% vs. control 25%, P = 0.01). This trend was replicated in the imputed Immunochip data for the two K528-bearing haplotypes MKDRQ (AS 33% vs. controls 27%, P = 1.2 × 10-24) and MKDRE (AS 8% vs. controls 7%, P = 0.004). The ERAP1 association with AS is therefore predominantly attributable to common ERAP1 haplotypes and haplotype combinations.

Defining the genetic susceptibility to cervical neoplasia-A genome-wide association study.

A small percentage of women with cervical HPV infection progress to cervical neoplasia, and the risk factors determining progression are incompletely understood. We sought to define the genetic loci involved in cervical neoplasia and to assess its heritability using unbiased unrelated case/control statistical approaches. We demonstrated strong association of cervical neoplasia with risk and protective HLA haplotypes that are determined by the amino-acids carried at positions 13 and 71 in pocket 4 of HLA-DRB1 and position 156 in HLA-B. Furthermore, 36% (standard error 2.4%) of liability of HPV-associated cervical pre-cancer and cancer is determined by common genetic variants. Women in the highest 10% of genetic risk scores have approximately >7.1% risk, and those in the highest 5% have approximately >21.6% risk, of developing cervical neoplasia. Future studies should examine genetic risk prediction in assessing the risk of cervical neoplasia further, in combination with other screening methods.

Genetic insights into common pathways and complex relationships among immune-mediated diseases.

Shared aetiopathogenic factors among immune-mediated diseases have long been suggested by their co-familiality and co-occurrence, and molecular support has been provided by analysis of human leukocyte antigen (HLA) haplotypes and genome-wide association studies. The interrelationships can now be better appreciated following the genotyping of large immune disease sample sets on a shared SNP array: the 'Immunochip'. Here, we systematically analyse loci shared among major immune-mediated diseases. This reveals that several diseases share multiple susceptibility loci, but there are many nuances. The most associated variant at a given locus frequently differs and, even when shared, the same allele often has opposite associations. Interestingly, risk alleles conferring the largest effect sizes are usually disease-specific. These factors help to explain why early evidence of extensive 'sharing' is not always reflected in epidemiological overlap.

HLA and KIR Associations of Cervical Neoplasia.

Background: Cervical cancer is the fourth most common cancer in women, and we recently reported human leukocyte antigen (HLA) alleles showing strong associations with cervical neoplasia risk and protection. HLA ligands are recognized by killer immunoglobulin-like receptors (KIRs) expressed on a range of immune cell subsets, governing their proinflammatory activity. We hypothesized that the inheritance of particular HLA-KIR combinations would increase cervical neoplasia risk. Methods: Here, we used HLA and KIR dosages imputed from single-nucleotide polymorphism genotype data from 2143 cervical neoplasia cases and 13858 healthy controls of European decent. Results: The following 4 novel HLA alleles were identified in association with cervical neoplasia, owing to their linkage disequilibrium with known cervical neoplasia-associated HLA-DRB1 alleles: HLA-DRB3*9901 (odds ratio [OR], 1.24; P = 2.49 × 10-9), HLA-DRB5*0101 (OR, 1.29; P = 2.26 × 10-8), HLA-DRB5*9901 (OR, 0.77; P = 1.90 × 10-9), and HLA-DRB3*0301 (OR, 0.63; P = 4.06 × 10-5). We also found that homozygosity of HLA-C1 group alleles is a protective factor for human papillomavirus type 16 (HPV16)-related cervical neoplasia (C1/C1; OR, 0.79; P = .005). This protective association was restricted to carriers of either KIR2DL2 (OR, 0.67; P = .00045) or KIR2DS2 (OR, 0.69; P = .0006). Conclusions: Our findings suggest that HLA-C1 group alleles play a role in protecting against HPV16-related cervical neoplasia, mainly through a KIR-mediated mechanism.

Evaluating the Performance of Fine-Mapping Strategies at Common Variant GWAS Loci.

The growing availability of high-quality genomic annotation has increased the potential for mechanistic insights when the specific variants driving common genome-wide association signals are accurately localized. A range of fine-mapping strategies have been advocated, and specific successes reported, but the overall performance of such approaches, in the face of the extensive linkage disequilibrium that characterizes the human genome, is not well understood. Using simulations based on sequence data from the 1000 Genomes Project, we quantify the extent to which fine-mapping, here conducted using an approximate Bayesian approach, can be expected to lead to useful improvements in causal variant localization. We show that resolution is highly variable between loci, and that performance is severely degraded as the statistical power to detect association is reduced. We confirm that, where causal variants are shared between ancestry groups, further improvements in performance can be obtained in a trans-ethnic fine-mapping design. Finally, using empirical data from a recently published genome-wide association study for ankylosing spondylitis, we provide empirical confirmation of the behaviour of the approximate Bayesian approach and demonstrate that seven of twenty-six loci can be fine-mapped to fewer than ten variants.

Genetic diagnostic profiling in axial spondyloarthritis: a real world study.

OBJECTIVES: Spondyloarthritis (SpA) is often diagnosed late in the course of the disease and improved methods for early diagnosis are required. We have tested the ability of genetic profiling to diagnose axial SpA (axSpA) as a whole group, or ankylosing spondylitis (AS) alone, in a cohort of chronic back pain patients. METHODS: 282 patients were recruited from centres in the United Kingdom, Germany, Taiwan, Canada, Columbia and Turkey as part of the ASAS classification criteria for axSpA study (ASAS cohort). Subjects were classified according to the ASAS axSpA criteria, and the modified New York Criteria for AS. Patients were genotyped for ~200,000 immune-mediated disease SNPs using the Illumina Immunochip. RESULTS: We first established the predictive accuracy of genetic data comparing 9,638 healthy controls and 4,428 AS cases from the homogenous International Genetics of AS (IGAS) Consortium Immunochip study which showed excellent predictive power (AUC=0.91). Genetic risk scores had lower predictive power (AUC=0.83) comparing ASAS cohort axSpA cases meeting the ASAS imaging criteria with IGAS controls. Comparing genetic risk scores showed moderate discriminatory capacity between IGAS AS and ASAS imaging positive cases (AUC 0.67±0.05), indicating that significant differences in genetic makeup exist between the cohorts. CONCLUSIONS: In a clinical setting of referred back pain patients suspected to have axial SpA we were unable to use genetic data to construct a predictive model better than that based on existing clinical data. Potential confounding factors include significant heterogeneity in the ASAS cohort, possibly reflecting the disease heterogeneity of axSpA, or differences between centres in ascertainment or classification performance.

The genetic association of RUNX3 with ankylosing spondylitis can be explained by allele-specific effects on IRF4 recruitment that alter gene expression.

OBJECTIVES: To identify the functional basis for the genetic association of single nucleotide polymorphisms (SNP), upstream of the RUNX3 promoter, with ankylosing spondylitis (AS). METHODS: We performed conditional analysis of genetic association data and used ENCODE data on chromatin remodelling and transcription factor (TF) binding sites to identify the primary AS-associated regulatory SNP in the RUNX3 region. The functional effects of this SNP were tested in luciferase reporter assays. Its effects on TF binding were investigated by electrophoretic mobility gel shift assays and chromatin immunoprecipitation. RUNX3 mRNA levels were compared in primary CD8+ T cells of AS risk and protective genotypes by real-time PCR. RESULTS: The association of the RUNX3 SNP rs4648889 with AS (p<7.6×10(-14)) was robust to conditioning on all other SNPs in this region. We identified a 2 kb putative regulatory element, upstream of RUNX3, containing rs4648889. In reporter gene constructs, the protective rs4648889 'G' allele increased luciferase activity ninefold but significantly less activity (4.3-fold) was seen with the AS risk 'A' allele (p≤0.01). The binding of Jurkat or CD8+ T-cell nuclear extracts to the risk allele was decreased and IRF4 recruitment was reduced. The AS-risk allele also affected H3K4Me1 histone methylation and associated with an allele-specific reduction in RUNX3 mRNA (p<0.05). CONCLUSION: We identified a regulatory region upstream of RUNX3 that is modulated by rs4648889. The risk allele decreases TF binding (including IRF4) and reduces reporter activity and RUNX3 expression. These findings may have important implications for understanding the role of T cells and other immune cells in AS.

An ankylosing spondylitis-associated genetic variant in the IL23R-IL12RB2 intergenic region modulates enhancer activity and is associated with increased Th1-cell differentiation.

OBJECTIVES: To explore the functional basis for the association between ankylosing spondylitis (AS) and single-nucleotide polymorphisms (SNPs) in the IL23R-IL12RB2 intergenic region. METHODS: We performed conditional analysis on genetic association data and used epigenetic data on chromatin remodelling and transcription factor (TF) binding to identify the primary AS-associated IL23R-IL12RB2 intergenic SNP. Functional effects were tested in luciferase reporter assays in HEK293T cells and allele-specific TF binding was investigated by electrophoretic mobility gel shift assays. IL23R and IL12RB2 mRNA levels in CD4+ T cells were compared between cases homozygous for the AS-risk 'A' allele and the protective 'G' allele. The proportions of interleukin (IL)-17A+ and interferon (IFN)-γ+ CD4+ T-cells were measured by fluorescence-activated cell sorting and compared between these AS-risk and protective genotypes. RESULTS: Conditional analysis identified rs11209032 as the probable causal SNP within a 1.14 kb putative enhancer between IL23R and IL12RB2. Reduced luciferase activity was seen for the risk allele (p<0.001) and reduced H3K4me1 methylation observed in CD4+ T-cells from 'A/A' homozygotes (p=0.02). The binding of nuclear extract to the risk allele was decreased ∼3.5-fold compared with the protective allele (p<0.001). The proportion of IFN-γ+ CD4+ T-cells was increased in 'A/A' homozygotes (p=0.004), but neither IL23R nor IL12RB2 mRNA was affected. CONCLUSIONS: The rs11209032 SNP downstream of IL23R forms part of an enhancer, allelic variation of which may influence Th1-cell numbers. Homozygosity for the risk 'A' allele is associated with more IFN-γ-secreting (Th1) cells. Further work is necessary to explain the mechanisms for these important observations.

Resequencing and fine-mapping of the chromosome 12q13-14 locus associated with multiple sclerosis refines the number of implicated genes.

Multiple sclerosis (MS) is a common chronic inflammatory disease of the central nervous system. Susceptibility to the disease is affected by both environmental and genetic factors. Genetic factors include haplotypes in the histocompatibility complex (MHC) and over 50 non-MHC loci reported by genome-wide association studies. Amongst these, we previously reported polymorphisms in chromosome 12q13-14 with a protective effect in individuals of European descent. This locus spans 288 kb and contains 17 genes, including several candidate genes which have potentially significant pathogenic and therapeutic implications. In this study, we aimed to fine-map this locus. We have implemented a two-phase study: a variant discovery phase where we have used next-generation sequencing and two target-enrichment strategies [long-range polymerase chain reaction (PCR) and Nimblegen's solution phase hybridization capture] in pools of 25 samples; and a genotyping phase where we genotyped 712 variants in 3577 healthy controls and 3269 MS patients. This study confirmed the association (rs2069502, P = 9.9 × 10(-11), OR = 0.787) and narrowed down the locus of association to an 86.5 kb region. Although the study was unable to pinpoint the key-associated variant, we have identified a 42 (genotyped and imputed) single-nucleotide polymorphism haplotype block likely to harbour the causal variant. No evidence of association at previously reported low-frequency variants in CYP27B1 was observed. As part of the study we compared variant discovery performance using two target-enrichment strategies. We concluded that our pools enriched with Nimblegen's solution phase hybridization capture had better sensitivity to detect true variants than the pools enriched with long-range PCR, whilst specificity was better in the long-range PCR-enriched pools compared with solution phase hybridization capture enriched pools; this result has important implications for the design of future fine-mapping studies.

Brief report: high-throughput sequencing of IL23R reveals a low-frequency, nonsynonymous single-nucleotide polymorphism that is associated with ankylosing spondylitis in a Han Chinese population.

OBJECTIVE: Ankylosing spondylitis (AS) is a highly heritable common inflammatory arthritis that targets the spine and sacroiliac joints of the pelvis, causing pain and stiffness and leading eventually to joint fusion. Although previous studies have shown a strong association of IL23R with AS in white Europeans, similar studies in East Asian populations have shown no association with common variants of IL23R, suggesting either that IL23R variants have no role or that rare genetic variants contribute. The present study was undertaken to screen IL23R to identify rare variants associated with AS in Han Chinese. METHODS: A 170-kb region containing IL23R and its flanking regions was sequenced in 50 patients with AS and 50 ethnically matched healthy control subjects from a Han Chinese population. In addition, the 30-kb region of peak association in white Europeans was sequenced in 650 patients with AS and 1,300 healthy controls. Validation genotyping was undertaken in 846 patients with AS and 1,308 healthy controls. RESULTS: We identified 1,047 variants, of which 729 were not found in the dbSNP genomic build 130. Several potentially functional rare variants in IL23R were identified, including one nonsynonomous single-nucleotide polymorphism (nsSNP), Gly(149) Arg (position 67421184 GA on chromosome 1). Validation genotyping showed that the Gly(149) Arg variant was associated with AS (odds ratio 0.61, P = 0.0054). CONCLUSION: This is the first study to implicate rare IL23R variants in the pathogenesis of AS. The results identified a low-frequency nsSNP with predicted loss-of-function effects that was protectively associated with AS in Han Chinese, suggesting that decreased function of the interleukin-23 (IL-23) receptor protects against AS. These findings further support the notion that IL-23 signaling has an important role in the pathogenesis of AS.

CoPheScan: phenome-wide association studies accounting for linkage disequilibrium.

Phenome-wide association studies (PheWAS) facilitate the discovery of associations between a single genetic variant with multiple phenotypes. For variants which impact a specific protein, this can help identify additional therapeutic indications or on-target side effects of intervening on that protein. However, PheWAS is restricted by an inability to distinguish confounding due to linkage disequilibrium (LD) from true pleiotropy. Here we describe CoPheScan (Coloc adapted Phenome-wide Scan), a Bayesian approach that enables an intuitive and systematic exploration of causal associations while simultaneously addressing LD confounding. We demonstrate its performance through simulation, showing considerably better control of false positive rates than a conventional approach not accounting for LD. We used CoPheScan to perform PheWAS of protein-truncating variants and fine-mapped variants from disease and pQTL studies, in 2275 disease phenotypes from the UK Biobank. Our results identify the complexity of known pleiotropic genes such as APOE, and suggest a new causal role for TGM3 in skin cancer.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens.

BACKGROUND: Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions. RESULTS: To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci. CONCLUSIONS: We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.