In human systemic lupus erythematosus (SLE), diverse autoantibodies accumulate over years before disease manifestation. Unaffected relatives of SLE patients frequently share a sustained production of autoantibodies with indiscriminable specificity, usually without ever acquiring the disease. We studied relations of IgG autoantibody profiles and peripheral blood activated regulatory T-cells (aTregs), represented by CD4+CD25bright T-cells that were regularly 70–90% Foxp3+. We found consistent positive correlations of broad-range as well as specific SLE-associated IgG with aTreg frequencies within unaffected relatives, but not patients or unrelated controls. Our interpretation: unaffected relatives with shared genetic factors compensated pathogenic effects by aTregs engaged in parallel with the individual autoantibody production. To study this further, we applied a novel analytic approach named coreferentiality that tests the indirect relatedness of parameters in respect to multivariate phenotype data. Results show that independently of their direct correlation, aTreg frequencies and specific SLE-associated IgG were likely functionally related in unaffected relatives: they significantly parallelled each other in their relations to broad-range immunoblot autoantibody profiles. In unaffected relatives, we also found coreferential effects of genetic variation in the loci encoding IL-2 and CD25. A model of CD25 functional genetic effects constructed by coreferentiality maximization suggests that IL-2-CD25 interaction, likely stimulating aTregs in unaffected relatives, had an opposed effect in SLE patients, presumably triggering primarily T-effector cells in this group. Coreferentiality modeling as we do it here could also be useful in other contexts, particularly to explore combined functional genetic effects.

Selective IgA deficiency (IgAD; serum IgA<0.07 g/l) is the most common form of human primary immune deficiency, affecting approximately 1∶600 individuals in populations of Northern European ancestry. The polygenic nature of IgAD is underscored by the recent identification of several new risk genes in a genome-wide association study. Among the characterized susceptibility loci, the association with specific HLA haplotypes represents the major genetic risk factor for IgAD. Despite the robust association, the nature and location of the causal variants in the HLA region remains unknown. To better characterize the association signal in this region, we performed a high-density SNP mapping of the HLA locus and imputed the genotypes of common HLA-B, -DRB1, and -DQB1 alleles in a combined sample of 772 IgAD patients and 1,976 matched controls from 3 independent European populations. We confirmed the complex nature of the association with the HLA locus, which is the result of multiple effects spanning the entire HLA region. The primary association signal mapped to the HLA-DQB1*02 allele in the HLA Class II region (combined P = 7.69×10−57; OR = 2.80) resulting from the combined independent effects of the HLA-B*0801-DRB1*0301-DQB1*02 and -DRB1*0701-DQB1*02 haplotypes, while additional secondary signals were associated with the DRB1*0102 (combined P = 5.86×10−17; OR = 4.28) and the DRB1*1501 (combined P = 2.24×10−35; OR = 0.13) alleles. Despite the strong population-specific frequencies of HLA alleles, we found a remarkable conservation of these effects regardless of the ethnic background, which supports the use of large multi-ethnic populations to characterize shared genetic association signals in the HLA region. We also provide evidence for the location of association signals within the specific extended haplotypes, which will guide future sequencing studies aimed at characterizing the precise functional variants contributing to disease pathogenesis.

Author Summary

The human leukocyte antigen (HLA) locus is robustly associated with many immune-mediated conditions. However, identification of the genetic variants contributing to the disease pathophysiology has been greatly hampered by the extensive chromosomal conservation within this genomic region. To better understand the association of the HLA locus in selective IgA deficiency (IgAD), we used an extensive genotyping database from a recent genome-wide association study (GWAS) to generate a high-density SNP map of this region in a combined sample of >2,700 individuals from 3 independent European populations. In addition, we took advantage of recent methodological advances to impute the more common HLA-B, -DRB1, and -DQB1 alleles in all subjects. We confirmed the strong disease-association of the HLA locus and identified several different signals located in specific conserved HLA haplotypes contributing independent risk or protection for IgAD. Further analysis of the chromosomal sequences associated with the associated HLA alleles allowed us to refine the mapping of the susceptibility variants. These findings represent the most comprehensive high-density SNP mapping of the HLA locus in IgAD to date and provide important new information as to the location of the genetic variants contributing to this common immune deficiency.

Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 × 10−8): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P ≤ 1 × 10−5. A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 × 10−3) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.

Summary

Luminescent bacteria (γ-Proteobacteria: Vibrionaceae) are found in complex bilobed light organs of both sepiolid and loliginid squids (Mollusca: Cephalopoda). Despite the existence of multiple strain colonization between Vibrio bacteria and loliginid squids, specificity at the genus level still exists and may influence interactions between symbiotic and free-living stages of the symbiont. The environmentally transmitted behaviour of Vibrio symbionts bestows a certain degree of recognition that exists prior and subsequent to the colonization process. Therefore, we identified bacterial genes required for successful colonization of loliginid light organs by examining transcripts solely expressed in either the light organ or free-living stages. Selective capture of transcribed sequences (SCOTS) was used to differentiate genes expressed by the same bacterium when thriving in two different environments (i.e. loliginid light organs and seawater). Genes specific for squid light organs included vulnibactin synthetase, outer membrane protein W and dihydroxy dehydratase, which have been associated with the maintenance of bacterial host associations in other systems. In contrast, genes that were solely expressed in the free-living condition consisted of transcripts recognized as important factors for bacterial survival in the environment. These transcripts included genes for methyl accepting chemotaxis proteins, arginine decarboxylase and chitinase. These results provide valuable information regarding mechanisms determining specificity, establishment, and maintenance of bacteria–squid associations.

Background

CD4+CD25+ regulatory T cells play an essential role in maintaining immune homeostasis and preventing autoimmunity. Therefore, defects in Treg development, maintenance or function have been associated with several human autoimmune diseases including Systemic Lupus Erythematosus (SLE), a systemic autoimmune disease characterized by loss of tolerance to nuclear components and significantly more frequent in females.

Results

To investigate the involvement of Treg in SLE pathogenesis, we determined the frequency of CD4+CD25+CD45RO+ T cells, which encompass the majority of Treg activity, in the PBMC of 148 SLE patients (76 patients were part of 54 families), 166 relatives and 117 controls. SLE patients and their relatives were recruited in several Portuguese hospitals and through the Portuguese Lupus Association. Control individuals were blood donors recruited from several regional blood donor centers. Treg frequency was significantly lower in SLE patients than healthy controls (z = -6.161, P < 0.00001) and intermediate in the relatives' group. Remarkably, this T cell subset was also lower in females, most strikingly in the control population (z = 4.121, P < 0.001). We further ascertained that the decreased frequency of Treg in SLE patients resulted from the specific reduction of bona fide FOXP3+CD4+CD25+ Treg. Treg frequency was negatively correlated with SLE activity index (SLEDAI) and titers of serum anti-dsDNA antibodies. Both Treg frequency and disease activity were modulated by IVIg treatment in a documented SLE case. The segregation of Treg frequency within the SLE families was indicative of a genetic trait. Candidate gene analysis revealed that specific variants of CTLA4 and TGFβ were associated with the decreased frequency of Treg in PBMC, while FOXP3 gene variants were associated with affection status, but not with Treg frequency.

Conclusion

SLE patients have impaired Treg production or maintenance, a trait strongly associated with SLE disease activity and autoantibody titers, and possibly resulting from the inability to convert FOXP3+CD25- into FOXP3+CD25+ T cells. Treg frequency is highly heritable within SLE families, with specific variants of the CTLA4 and TGFβ genes contributing to this trait, while FOXP3 contributes to SLE through mechanisms not involving a modulation of Treg frequency. These findings establish that the genetic components in SLE pathogenesis include genes related to Treg generation or maintenance.