IgA measurements in over 12 000 Swedish twins reveal sex differential heritability and regulatory locus near CD30L

Genetic regulation of serum IgA levels and susceptibility to common immune, infectious, kidney, and cardio-metabolic traits

Immunoglobulin A (IgA) mediates mucosal responses to food antigens and the intestinal microbiome and is involved in susceptibility to mucosal pathogens, celiac disease, inflammatory bowel disease, and IgA nephropathy. We performed a genome-wide association study of serum IgA levels in 41,263 individuals of diverse ancestries and identified 20 genome-wide significant loci, including 9 known and 11 novel loci. Co-localization analyses with expression QTLs prioritized candidate genes for 14 of 20 significant loci. Most loci encoded genes that produced immune defects and IgA abnormalities when genetically manipulated in mice. We also observed positive genetic correlations of serum IgA levels with IgA nephropathy, type 2 diabetes, and body mass index, and negative correlations with celiac disease, inflammatory bowel disease, and several infections. Mendelian randomization supported elevated serum IgA as a causal factor in IgA nephropathy. African ancestry was consistently associated with higher serum IgA levels and greater frequency of IgA-increasing alleles compared to other ancestries. Our findings provide novel insights into the genetic regulation of IgA levels and its potential role in human disease.

A Comprehensive Evaluation of the Relationship Between Different IgG and IgA Anti-Modified Protein Autoantibodies in Rheumatoid Arthritis

Seropositive rheumatoid arthritis (RA) is characterized by the presence of rheumatoid factor (RF) and anti-citrullinated protein autoantibodies (ACPA) with different fine-specificities. Yet, other serum anti-modified protein autoantibodies (AMPA), e.g. anti-carbamylated (Carb), -acetylated (KAc), and malondialdehyde acetaldehyde (MAA) modified protein antibodies, have been described. In this comprehensive study, we analyze 30 different IgG and IgA AMPA reactivities to Cit, Carb, KAc, and MAA antigens detected by ELISA and autoantigen arrays in N=1985 newly diagnosed RA patients. Association with patient characteristics such as smoking and disease activity were explored. Carb and KAc reactivities by different assays were primarily seen in patients also positive for anti-citrulline reactivity. Modified vimentin (mod-Vim) peptides were used for direct comparison of different AMPA reactivities, revealing that IgA AMPA recognizing mod-Vim was mainly detected in subsets of patients with high IgG anti-Cit-Vim levels and a history of smoking. IgG reactivity to acetylation was mainly detected in a subset of patients with Cit and Carb reactivity. Anti-acetylated histone reactivity was RA-specific and associated with high anti-CCP2 IgG levels, multiple ACPA fine-specificities, and smoking status. This reactivity was also found to be present in CCP2+ RA-risk individuals without arthritis. Our data further demonstrate that IgG autoreactivity to MAA was increased in RA compared to controls with highest levels in CCP2+ RA, but was not RA-specific, and showed low correlation with other AMPA. Anti-MAA was instead associated with disease activity and was not significantly increased in CCP2+ individuals at risk of RA. Notably, RA patients could be subdivided into four different subsets based on their AMPA IgG and IgA reactivity profiles. Our serology results were complemented by screening of monoclonal antibodies derived from single B cells from RA patients for the same antigens as the RA cohort. Certain CCP2+ clones had Carb or Carb+KAc+ multireactivity, while such reactivities were not found in CCP2- clones. We conclude that autoantibodies exhibiting different patterns of ACPA fine-specificities as well as Carb and KAc reactivity are present in RA and may be derived from multireactive B-cell clones. Carb and KAc could be considered reactivities within the "Cit-umbrella" similar to ACPA fine-specificities, while MAA reactivity is distinctly different.

Multi-method genome- and epigenome-wide studies of inflammatory protein levels in healthy older adults

BACKGROUND

The molecular factors which control circulating levels of inflammatory proteins are not well understood. Furthermore, association studies between molecular probes and human traits are often performed by linear model-based methods which may fail to account for complex structure and interrelationships within molecular datasets.

METHODS

In this study, we perform genome- and epigenome-wide association studies (GWAS/EWAS) on the levels of 70 plasma-derived inflammatory protein biomarkers in healthy older adults (Lothian Birth Cohort 1936; n = 876; Olink® inflammation panel). We employ a Bayesian framework (BayesR+) which can account for issues pertaining to data structure and unknown confounding variables (with sensitivity analyses using ordinary least squares- (OLS) and mixed model-based approaches).

RESULTS

We identified 13 SNPs associated with 13 proteins (n = 1 SNP each) concordant across OLS and Bayesian methods. We identified 3 CpG sites spread across 3 proteins (n = 1 CpG each) that were concordant across OLS, mixed-model and Bayesian analyses. Tagged genetic variants accounted for up to 45% of variance in protein levels (for MCP2, 36% of variance alone attributable to 1 polymorphism). Methylation data accounted for up to 46% of variation in protein levels (for CXCL10). Up to 66% of variation in protein levels (for VEGFA) was explained using genetic and epigenetic data combined. We demonstrated putative causal relationships between CD6 and IL18R1 with inflammatory bowel disease and between IL12B and Crohn's disease.

CONCLUSIONS

Our data may aid understanding of the molecular regulation of the circulating inflammatory proteome as well as causal relationships between inflammatory mediators and disease.

Rheumatoid arthritis patients display B-cell dysregulation already in the naïve repertoire consistent with defects in B-cell tolerance

B cells are postulated to be central in seropositive rheumatoid arthritis (RA). Here, we use exploratory mass cytometry (n = 23) and next-generation sequencing (n = 19) to study B-cell repertoire shifts in RA patients. Expression of several B-cell markers were significantly different in ACPA⁺ RA compared to healthy controls, including an increase in HLA-DR across subsets, CD22 in clusters of IgM⁺ B cells and CD11c in IgA⁺ memory. Moreover, both IgA⁺ and IgG⁺ double negative (IgD⁻ CD27⁻) CD11c⁺ B cells were increased in ACPA⁺ RA, and there was a trend for elevation in a CXCR5/CCR6^high transitional B-cell cluster. In the RA BCR repertoire, there were significant differences in subclass distribution and, notably, the frequency of VH with low somatic hypermutation (SHM) was strikingly higher, especially in IgG1 (p < 0.0001). Furthermore, both ACPA⁺ and ACPA⁻ RA patients had significantly higher total serum IgA and IgM compared to controls, based on serology of larger cohorts (n = 3494 IgA; n = 397 IgM). The observed elevated Ig-levels, distortion in IgM⁺ B cells, increase in double negative B cells, change in B-cell markers, and elevation of unmutated IgG⁺ B cells suggests defects in B-cell tolerance in RA. This may represent an underlying cause of increased polyreactivity and autoimmunity in RA.

Human genetic variants and age are the strongest predictors of humoral immune responses to common pathogens and vaccines

BACKGROUND

Humoral immune responses to infectious agents or vaccination vary substantially among individuals, and many of the factors responsible for this variability remain to be defined. Current evidence suggests that human genetic variation influences (i) serum immunoglobulin levels, (ii) seroconversion rates, and (iii) intensity of antigen-specific immune responses. Here, we evaluated the impact of intrinsic (age and sex), environmental, and genetic factors on the variability of humoral response to common pathogens and vaccines.

METHODS

We characterized the serological response to 15 antigens from common human pathogens or vaccines, in an age- and sex-stratified cohort of 1000 healthy individuals (Milieu Intérieur cohort). Using clinical-grade serological assays, we measured total IgA, IgE, IgG, and IgM levels, as well as qualitative (serostatus) and quantitative IgG responses to cytomegalovirus, Epstein-Barr virus, herpes simplex virus 1 and 2, varicella zoster virus, Helicobacter pylori, Toxoplasma gondii, influenza A virus, measles, mumps, rubella, and hepatitis B virus. Following genome-wide genotyping of single nucleotide polymorphisms and imputation, we examined associations between ~ 5 million genetic variants and antibody responses using single marker and gene burden tests.

RESULTS

We identified age and sex as important determinants of humoral immunity, with older individuals and women having higher rates of seropositivity for most antigens. Genome-wide association studies revealed significant associations between variants in the human leukocyte antigen (HLA) class II region on chromosome 6 and anti-EBV and anti-rubella IgG levels. We used HLA imputation to fine map these associations to amino acid variants in the peptide-binding groove of HLA-DRβ1 and HLA-DPβ1, respectively. We also observed significant associations for total IgA levels with two loci on chromosome 2 and with specific KIR-HLA combinations.

CONCLUSIONS

Using extensive serological testing and genome-wide association analyses in a well-characterized cohort of healthy individuals, we demonstrated that age, sex, and specific human genetic variants contribute to inter-individual variability in humoral immunity. By highlighting genes and pathways implicated in the normal antibody response to frequently encountered antigens, these findings provide a basis to better understand disease pathogenesis.

TRIALS REGISTRATION

ClinicalTrials.gov , NCT01699893.

GWAS for serum galactose-deficient IgA1 implicates critical genes of the O-glycosylation pathway

Aberrant O-glycosylation of serum immunoglobulin A1 (IgA1) represents a heritable pathogenic defect in IgA nephropathy, the most common form of glomerulonephritis worldwide, but specific genetic factors involved in its determination are not known. We performed a quantitative GWAS for serum levels of galactose-deficient IgA1 (Gd-IgA1) in 2,633 subjects of European and East Asian ancestry and discovered two genome-wide significant loci, in C1GALT1 (rs13226913, P = 3.2 x 10⁻¹¹) and C1GALT1C1 (rs5910940, P = 2.7 x 10⁻⁸). These genes encode molecular partners essential for enzymatic O-glycosylation of IgA1. We demonstrated that these two loci explain approximately 7% of variability in circulating Gd-IgA1 in Europeans, but only 2% in East Asians. Notably, the Gd-IgA1-increasing allele of rs13226913 is common in Europeans, but rare in East Asians. Moreover, rs13226913 represents a strong cis-eQTL for C1GALT1 that encodes the key enzyme responsible for the transfer of galactose to O-linked glycans on IgA1. By in vitro siRNA knock-down studies, we confirmed that mRNA levels of both C1GALT1 and C1GALT1C1 determine the rate of secretion of Gd-IgA1 in IgA1-producing cells. Our findings provide novel insights into the genetic regulation of O-glycosylation and are relevant not only to IgA nephropathy, but also to other complex traits associated with O-glycosylation defects, including inflammatory bowel disease, hematologic disease, and cancer.

O-glycosylation is a common type of post-translational modification of proteins; specific abnormalities in the mechanism of O-glycosylation have been implicated in cancer, inflammatory and blood diseases. However, the molecular basis of abnormal O-glycosylation in these complex disorders is not known. We studied the genetic basis of defective O-glycosylation of serum immunoglobulin A1 (IgA1), that represents the key pathogenic defect in IgA nephropathy, the most common form of primary glomerulonephritis worldwide. We report our results of the first genome-wide association study for this trait using serum assays in 2,633 individuals of European and East-Asian ancestry. In our genome scan, we observed two significant signals with large effects, on chromosomes 7p21.3 and Xq24, jointly explaining about 7% of trait variability. These signals implicate two genes that encode molecular partners essential for enzymatic O-glycosylation of IgA1 and mucins, and represent potential new targets for therapy.

IgA1 Glycosylation Is Heritable in Healthy Twins

IgA nephropathy (IgAN) is the most common form of primary GN and an important cause of kidney failure. Characteristically, patients with IgAN have increased serum levels of undergalactosylated IgA1 (gd-IgA1). To assess the degree to which serum gd-IgA1 levels are genetically determined in healthy individuals, we determined serum IgA and gd-IgA1 levels by ELISA in a sample of 148 healthy female twins, including 27 monozygotic and 47 dizygotic pairs. Using the classic twin model, we found the heritability of serum gd-IgA1 and IgA levels to be 80% (95% confidence interval, 66% to 89%) and 46% (95% confidence interval, 15% to 69%), respectively. These data indicate that serum gd-IgA1 levels are highly heritable. Elucidating the genetic basis of this heritability will be important in understanding the pathogenesis of IgAN.

The higher frequency of IgA deficiency among Swedish twins is not explained by HLA haplotypes

Serum immunoglobulin A (IgA) concentrations were determined in 12 600 adult Swedish twins, applying a high-throughput reverse-phase protein microarray technique. The prevalence of IgA deficiency (IgAD) was found to be 1:241 in monozygotic (MZ) twins and 1:198 in dizygotic (DZ) twins. Hence, the prevalence in twins is markedly elevated as compared with the normal Swedish adult population (1:600). The twins did not show a difference in the frequency of HLA haplotypes in comparison with almost 40 000 healthy Swedish controls. As expected, the risk-conveying HLA alleles A*01, B*08 and DRB1*01 were overrepresented among the IgAD twins and were also associated with significantly lower mean serum IgA concentrations in the twin cohort. In contrast, significantly higher mean IgA concentrations were found among individuals carrying the protective HLA alleles B*07 and DRB1*15. Exome sequencing data from two MZ twin pairs discordant for the deficiency showed no differences between the siblings. Model fitting analyses derived a heritability of 35% and indicate that genetic influences are modestly important for IgAD. The probandwise concordance rates for IgAD were found to be 31% for MZ and 13% for DZ twins.