A common allele of HLA is associated with asymptomatic SARS-CoV-2 infection.

Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic1-4. Although most global efforts have focused on severe illness in COVID-19, examining asymptomatic infection provides a unique opportunity to consider early immunological features that promote rapid viral clearance. Here, postulating that variation in the human leukocyte antigen (HLA) loci may underly processes mediating asymptomatic infection, we enrolled 29,947 individuals, for whom high-resolution HLA genotyping data were available, in a smartphone-based study designed to track COVID-19 symptoms and outcomes. Our discovery cohort (n = 1,428) comprised unvaccinated individuals who reported a positive test result for SARS-CoV-2. We tested for association of five HLA loci with disease course and identified a strong association between HLA-B*15:01 and asymptomatic infection, observed in two independent cohorts. Suggesting that this genetic association is due to pre-existing T cell immunity, we show that T cells from pre-pandemic samples from individuals carrying HLA-B*15:01 were reactive to the immunodominant SARS-CoV-2 S-derived peptide NQKLIANQF. The majority of the reactive T cells displayed a memory phenotype, were highly polyfunctional and were cross-reactive to a peptide derived from seasonal coronaviruses. The crystal structure of HLA-B*15:01-peptide complexes demonstrates that the peptides NQKLIANQF and NQKLIANAF (from OC43-CoV and HKU1-CoV) share a similar ability to be stabilized and presented by HLA-B*15:01. Finally, we show that the structural similarity of the peptides underpins T cell cross-reactivity of high-affinity public T cell receptors, providing the molecular basis for HLA-B*15:01-mediated pre-existing immunity.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region.

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

Comparison between qPCR and RNA-seq reveals challenges of quantifying HLA expression.

Human leukocyte antigen (HLA) class I and II loci are essential elements of innate and acquired immunity. Their functions include antigen presentation to T cells leading to cellular and humoral immune responses, and modulation of NK cells. Their exceptional influence on disease outcome has now been made clear by genome-wide association studies. The exons encoding the peptide-binding groove have been the main focus for determining HLA effects on disease susceptibility/pathogenesis. However, HLA expression levels have also been implicated in disease outcome, adding another dimension to the extreme diversity of HLA that impacts variability in immune responses across individuals. To estimate HLA expression, immunogenetic studies traditionally rely on quantitative PCR (qPCR). Adoption of alternative high-throughput technologies such as RNA-seq has been hampered by technical issues due to the extreme polymorphism at HLA genes. Recently, however, multiple bioinformatic methods have been developed to accurately estimate HLA expression from RNA-seq data. This opens an exciting opportunity to quantify HLA expression in large datasets but also brings questions on whether RNA-seq results are comparable to those by qPCR. In this study, we analyze three classes of expression data for HLA class I genes for a matched set of individuals: (a) RNA-seq, (b) qPCR, and (c) cell surface HLA-C expression. We observed a moderate correlation between expression estimates from qPCR and RNA-seq for HLA-A, -B, and -C (0.2 ≤ rho ≤ 0.53). We discuss technical and biological factors which need to be accounted for when comparing quantifications for different molecular phenotypes or using different techniques.

Uniparental markers reveal new insights on subcontinental ancestry and sex-biased admixture in Brazil.

The Brazilian population is a product of asymmetric admixture among European men and Amerindian and African women. However, Brazilian subcontinental ancestry is scarcely documented, especially regarding its African roots. Here, we aimed to unveil the uniparental continental and subcontinental contributions from distinct Brazilian regions, including South (n = 43), Southeast (n = 71), the poorly genetically characterized Central-Western region (n = 323), and a subset of unique Brazilian Amerindians (n = 24), in the context of their genome-wide ancestral contributions. The overwhelming majority of European Y haplogroups (85%) contrast sharply with the predominant African and Amerindian mtDNA haplogroups (73.2%) in admixed populations, whereas in Amerindians, non-Native haplogroups could only be detected through the paternal line. Our in-depth investigation of uniparental markers showed signals of an Andean and Central-Brazilian Amerindian maternal contribution to Southeastern and Central-Western Brazil (83.1 ± 2.1% and 56.9 ± 0.2%, respectively), the last having the highest paternal Amerindian ancestry yet described for an admixed Brazilian region (9.7%) and contrasting with higher Southern-Brazilian Amerindian contribution to Southern Brazil (59.6 ± 1%). Unlike the higher African Bantu contribution previously reported for the South and Southeast, a relevant Western African non-Bantu contribution was detected in those regions (85.7 ± 5% and 71.8 ± 10.8% respectively). In contrast, a higher Bantu contribution was described for the first time in the Central-West (64.8 ± 1.3% maternal and 86.9 ± 9.6% paternal). We observed sex-biased signatures consistent with the historically recorded Brazilian colonization and added new insights in the subcontinental maternal ancestry of Brazilians from regions never studied at this level.

High-throughput Interpretation of Killer-cell Immunoglobulin-like Receptor Short-read Sequencing Data with PING.

The killer-cell immunoglobulin-like receptor (KIR) complex on chromosome 19 encodes receptors that modulate the activity of natural killer cells, and variation in these genes has been linked to infectious and autoimmune disease, as well as having bearing on pregnancy and transplant outcomes. The medical relevance and high variability of KIR genes makes short-read sequencing an attractive technology for interrogating the region, providing a high-throughput, high-fidelity sequencing method that is cost-effective. However, because this gene complex is characterized by extensive nucleotide polymorphism, structural variation including gene fusions and deletions, and a high level of homology between genes, its interrogation at high resolution has been thwarted by bioinformatic challenges, with most studies limited to examining presence or absence of specific genes. Here, we present the PING (Pushing Immunogenetics to the Next Generation) pipeline, which incorporates empirical data, novel alignment strategies and a custom alignment processing workflow to enable high-throughput KIR sequence analysis from short-read data. PING provides KIR gene copy number classification functionality for all KIR genes through use of a comprehensive alignment reference. The gene copy number determined per individual enables an innovative genotype determination workflow using genotype-matched references. Together, these methods address the challenges imposed by the structural complexity and overall homology of the KIR complex. To determine copy number and genotype determination accuracy, we applied PING to European and African validation cohorts and a synthetic dataset. PING demonstrated exceptional copy number determination performance across all datasets and robust genotype determination performance. Finally, an investigation into discordant genotypes for the synthetic dataset provides insight into misaligned reads, advancing our understanding in interpretation of short-read sequencing data in complex genomic regions. PING promises to support a new era of studies of KIR polymorphism, delivering high-resolution KIR genotypes that are highly accurate, enabling high-quality, high-throughput KIR genotyping for disease and population studies.

Killer Cell Immunoglobulin-like Receptor Variants Are Associated with Protection from Symptoms Associated with More Severe Course in Parkinson Disease.

Immune dysfunction plays a role in the development of Parkinson disease (PD). NK cells regulate immune functions and are modulated by killer cell immunoglobulin-like receptors (KIR). KIR are expressed on the surface of NK cells and interact with HLA class I ligands on the surface of all nucleated cells. We investigated KIR-allelic polymorphism to interrogate the role of NK cells in PD. We sequenced KIR genes from 1314 PD patients and 1978 controls using next-generation methods and identified KIR genotypes using custom bioinformatics. We examined associations of KIR with PD susceptibility and disease features, including age at disease onset and clinical symptoms. We identified two KIR3DL1 alleles encoding highly expressed inhibitory receptors associated with protection from PD clinical features in the presence of their cognate ligand: KIR3DL1*015/HLA-Bw4 from rigidity (p c = 0.02, odds ratio [OR] = 0.39, 95% confidence interval [CI] 0.23-0.69) and KIR3DL1*002/HLA-Bw4i from gait difficulties (p c = 0.05, OR = 0.62, 95% CI 0.44-0.88), as well as composite symptoms associated with more severe disease. We also developed a KIR3DL1/HLA interaction strength metric and found that weak KIR3DL1/HLA interactions were associated with rigidity (pc = 0.05, OR = 9.73, 95% CI 2.13-172.5). Highly expressed KIR3DL1 variants protect against more debilitating symptoms of PD, strongly implying a role of NK cells in PD progression and manifestation.

Population-specific diversity of the immunoglobulin constant heavy G chain (IGHG) genes.

Human immunoglobulin G (IgG) molecules, IgG1, IgG2 and IgG3, exhibit substantial inter-individual variation in their constant heavy chain regions, as discovered by serological methods. This polymorphism is encoded by the IGHG1, IGHG2, and IGHG3 genes and may influence antibody function. We sequenced the coding fragments of these genes in 95 European Americans, 94 African Americans, and 94 Black South Africans. Striking differences were observed between the population groups, including extremely low amino acid sequence variation in IGHG1 among South Africans, and higher IGHG2 and IGHG3 diversity in individuals of African descent compared to individuals of European descent. Molecular definition of the loci illustrates a greater level of allelic polymorphism than previously described, including the presence of common IGHG2 and IGHG3 variants that were indistinguishable serologically. Comparison of our data with the 1000 Genome Project sequences indicates overall agreement between the datasets, although some inaccuracies in the 1000 Genomes Project are likely. These data represent the most comprehensive analysis of IGHG polymorphisms across major populations, which can now be applied to deciphering their functional impact.

Variation in genes implicated in B-cell development and antibody production affects susceptibility to pemphigus.

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease characterized by the presence of pathogenic autoantibodies against desmoglein 1, a component of intercellular desmosome junctions. PF occurs sporadically across the globe and is endemic in some Brazilian regions. Because PF is a B-cell-mediated disease, we aimed to study the impact of variants within genes encoding molecules involved in the different steps of B-cell development and antibody production on the susceptibility of endemic PF. We analysed 3,336 single nucleotide polymorphisms (SNPs) from 167 candidate genes genotyped with Illumina microarray in a cohort of 227 PF patients and 193 controls. After quality control and exclusion of non-informative and redundant SNPs, 607 variants in 149 genes remained in the logistic regression analysis, in which sex and ancestry were included as covariates. Our results revealed 10 SNPs within or nearby 11 genes that were associated with susceptibility to endemic PF (OR >1.56; p < 0.005): rs6657275*G (TGFB2); rs1818545*A (RAG1/RAG2/IFTAP);rs10781530*A (PAXX), rs10870140*G and rs10781522*A (TRAF2); rs535068*A (TNFRSF1B); rs324011*A (STAT6);rs6432018*C (YWHAQ); rs17149161*C (YWHAG); and rs2070729*C (IRF1). Interestingly, these SNPs have been previously associated with differential gene expression, mostly in peripheral blood, in publicly available databases. For the first time, we show that polymorphisms in genes involved in B-cell development and antibody production confer differential susceptibility to endemic PF, and therefore are candidates for possible functional studies to understand immunoglobulin gene rearrangement and its impact on diseases.

Genetic association and differential expression of HLAComplexGroup lncRNAs in pemphigus.

BACKGROUND: Pemphigus is a group of bullous diseases characterized by acantholysis and skin blisters. As for other autoimmune diseases, the strongest genetic associations found so far for pemphigus foliaceus (PF) and vulgaris (PV) are with alleles of HLA genes. However, apart from protein-coding genes, the MHC region includes a set of poorly explored long non-coding RNA (lncRNA) genes, the HLA complex group (HCG). OBJECTIVES: To investigate if HCG lncRNA alleles are associated with pemphigus susceptibility. METHODS AND RESULTS: We analyzed SNPs in 13 HCG lncRNA genes, both in PV (Germany: 241 patients; 1,188 controls) and endemic PF (Brazil: 227 patients; 194 controls), applying multivariate logistic regression. We found 55 associations with PV (pcorr < 0.01) and nine with endemic PF (pcorr < 0.05), the majority located in TSBP1-AS1 (which includes HCG23) and HCG27 lncRNA genes, independently of HLA alleles previously associated with pemphigus. The association of TSBP1-AS1 rs3129949*A allele was further replicated in sporadic PF (p = 0.027, OR = 0.054; 75 patients and 150 controls, all from Germany). Next, we evaluated the expression levels of TSBP1-AS1, TSBP1, HCG23, and HCG27 in blood mononuclear cells of Brazilian patients and controls. HCG27 was upregulated in endemic PF (p = 0.035, log2 FC = 1.3), while TSBP1-AS1 was downregulated in PV (p = 0.029, log2 FC = -1.29). The same expression patterns were also seen in cultured keratinocytes stimulated with IgG antibodies from patients and controls from Germany. TSBP1 mRNA levels were also decreased in endemic PF blood cells (p = 0.042, log2 FC = -2.14). TSBP1-AS1 and HCG27 were also observed downregulated in CD19+ cells of endemic PF (p < 0.01, log2 FC = -0.226 and -0.46 respectively). CONCLUSIONS: HCG lncRNAs are associated with susceptibility to pemphigus, being TSBP1-AS1 and HCG27 also differentially expressed in distinct cell populations. These results suggest a role for HCG lncRNAs in pemphigus autoimmunity.

An immunogenetic view of COVID-19.

Meeting the challenges brought by the COVID-19 pandemic requires an interdisciplinary approach. In this context, integrating knowledge of immune function with an understanding of how genetic variation influences the nature of immunity is a key challenge. Immunogenetics can help explain the heterogeneity of susceptibility and protection to the viral infection and disease progression. Here, we review the knowledge developed so far, discussing fundamental genes for triggering the innate and adaptive immune responses associated with a viral infection, especially with the SARS-CoV-2 mechanisms. We emphasize the role of the HLA and KIR genes, discussing what has been uncovered about their role in COVID-19 and addressing methodological challenges of studying these genes. Finally, we comment on questions that arise when studying admixed populations, highlighting the case of Brazil. We argue that the interplay between immunology and an understanding of genetic associations can provide an important contribution to our knowledge of COVID-19.

Characterization of serum cytokines and circulating microRNAs that are predicted to regulate inflammasome genes in cutaneous leishmaniasis patients.

Leishmaniasis is a neglected disease caused by an intracellular protozoan parasite of the genus Leishmania. Infection starts when this protozoan replicates in a phagolysosomal compartment in macrophages, after evading host immune responses. The balance of Th1 and Th2 immune responses is crucial in leishmaniasis because it will determine whether the infection will be under control or if clinical complications will occur. The inflammasome, which is activated during Leishmania infection, involves the action of caspase-1 and release of the proinflammatory cytokines interleukin-1ß and interleukin-18. Together, they contribute to the maintenance of an inflammatory response and pyroptosis. Here, we evaluated the serum levels of cytokines and the expression of circulating microRNAs related to inflammasome regulation in twenty-seven patients with cutaneous leishmaniasis in comparison to nine healthy individuals, in the context of the inflammasome activation. Evaluation of serum cytokines activation (IL-1ß, IL-2, IL-4, IL-6, IL-10, and IL-17) was performed by flow cytometry using CBA kits (cytometric beads array) while the expression of circulating microRNAs (miR-7, miR-133a, miR-146b, miR-155, miR-223, miR-328, and miR-342) in plasma was measured by quantitative polymerase chain reaction. Our results showed an increase of the expression of miR-7-5p (p < 10-5), miR-133a (p = 0.034), miR-146b (p = 0.003), miR-223-3p (p = 10-5), and miR-328-3p (p = 0.002), and cytokine levels for IL-1ß (p = 0.0005), IL-6 (p = 0.001), and IL-17 (p = 0.001) in patients with cutaneous leishmaniasis compared to the controls. These results suggest that microRNAs and cytokines can play an important role in regulating the human immune responses to Leishmania infection. Our findings may contribute to the understanding of the mechanisms of the gene regulation during the cutaneous leishmaniasis and to the identification of possible biomarkers of the infection.

A genetic variant in microRNA-146a is associated with sporadic breast cancer in a Southern Brazilian Population.

MicroRNAs (miRNAs) play an essential role in gene expression and affect the development of tumours, including breast cancer (BC). Polymorphisms in miRNA genes can affect the interaction of miRNAs with their target messenger RNA by interfering, creating or disrupting target sites. The single nucleotide polymorphism (SNP) rs2910164, located in the seed region of miR146a, was shown to be associated with BC among different populations. In the present study, we investigated whether rs2910164 is associated with BC in 326 patients and 411 controls from a Brazilian population of predominantly European ancestry. The presence of the allele rs2910164*C was associated with an increased risk of BC (OR=1.4, 95% CI=1.03-1.85, p = 0.03). We also analysed publicly available RNA-seq data to evaluate if miR146a is differentially expressed in different subtypes of BC. Genotyping was performed by polymerase chain reaction with sequence-specific primers (PCR-SSP). By leveraging public data from TCGA database, we analysed 461 patients and found that miR146a is significantly more expressed in BC than in non-tumor tissue (1.47 fold, p = 0.02) and is expressed to a greater degree in aggressive BC subtypes.

Screening the full leucocyte receptor complex genomic region revealed associations with pemphigus that might be explained by gene regulation.

Pemphigus foliaceus (PF) is a blistering autoimmune skin disease rare in most of the world but endemic in certain regions of Brazil. PF is characterized by the detachment of epidermal cells and the presence of autoantibodies against desmoglein 1. In previous studies, we have shown that genetic polymorphisms and variable expression levels of certain leucocyte receptor complex (LRC) genes were associated with PF. However, the role of the LRC on PF susceptibility remained to be investigated. Here, we analysed 527 tag single nucleotide polymorphisms (SNPs) distributed within the 1·5 Mb LRC. After quality control, a total of 176 SNPs were analysed in 229 patients with PF and 194 controls. Three SNPs were associated with differential susceptibility to PF. The intergenic variant rs465169 [odds ratio (OR) = 1·50; P = 0·004] is located in a region that might regulate several immune-related genes, including VSTM1, LILRB1/2, LAIR1/2, LILRA3/4 and LENG8. The rs35336528 (OR = 3·44; P = 0·009) and rs1865097 (OR = 0·57; P = 0·005) SNPs in LENG8 and FCAR genes, respectively, were also associated with PF. Moreover, we found four haplotypes with SNPs within the KIR3DL2/3, LAIR2 and LILRB1 genes associated with PF (P < 0·05), which corroborate previously reported associations. Thus, our results confirm the importance of the LRC for differential susceptibility to PF and reveal new markers that might influence expression levels of several LRC genes, as well as candidates for further functional studies.

Sequence and Phylogenetic Analysis of the Untranslated Promoter Regions for HLA Class I Genes.

Polymorphisms located within the MHC have been linked to many disease outcomes by mechanisms not yet fully understood in most cases. Variants located within untranslated regions of HLA genes are involved in allele-specific expression and may therefore underlie some of these disease associations. We determined sequences extending nearly 2 kb upstream of the transcription start site for 68 alleles from 57 major lineages of classical HLA class I genes. The nucleotide diversity within this promoter segment roughly follows that seen within the coding regions, with HLA-B showing the highest (∼1.9%), followed by HLA-A (∼1.8%), and HLA-C showing the lowest diversity (∼0.9%). Despite its greater diversity, HLA-B mRNA expression levels determined in 178 European Americans do not vary in an allele- or lineage-specific manner, unlike the differential expression levels of HLA-A or HLA-C reported previously. Close proximity of promoter sequences in phylogenetic trees is roughly reflected by similarity of expression pattern for most HLA-A and -C loci. Although promoter sequence divergence might impact promoter activity, we observed no clear link between the phylogenetic structures as represented by pairwise nucleotide differences in the promoter regions with estimated differences in mRNA expression levels for the classical class I loci. Further, no pair of class I loci showed coordinated expression levels, suggesting that distinct mechanisms across loci determine their expression level under nonstimulated conditions. These data serve as a foundation for more in-depth analysis of the functional consequences of promoter region variation within the classical HLA class I loci.

Epigenetic regulation of differential HLA-A allelic expression levels.

MHC class I expression levels influence the strength of immune responses and represent another variable in determining outcome to disease beyond peptide binding alone. Identification of the HLA loci that vary in allelic expression levels and delineating the mechanism responsible for expression variation may provide the opportunity to modify their expression therapeutically. We have examined the expression levels of allelic lineages at the HLA-A locus in a sample of 216 European Americans using a real-time polymerase chain reaction assay, which amplifies all HLA-A lineages specifically with equal efficiency, and observed a gradient of expression that associates with HLA-A allelic lineage (R = 0.6, P = 5 × 10(-25)). DNA methylation of the HLA-A gene appears to contribute to the variation in HLA-A mRNA expression levels, as a significant inverse correlation was observed between HLA-A mRNA expression levels in untreated cells and the degree to which expression is increased after treatment of the cells with a DNA methyltransferase inhibitor (R = 0.6, P = 2.8 × 10(-6)). Further, deep-sequencing and immunoprecipitation assays revealed allelic lineage-specific methylation patterns within the HLA-A promoter region where increased DNA methylation levels correlated significantly with reduced HLA-A expression levels (R = 0.89, P = 3.7 × 10(-9)). These data demonstrate HLA-A allelic lineage-specific variation in expression levels, and DNA methylation as a likely factor in contributing to this variation.

Differential gene expression levels might explain association of LAIR2 polymorphisms with pemphigus.

The leukocyte-associated immunoglobulin-like receptor 1 (LAIR-1) is a collagen-binding inhibitory receptor important for the regulation of immune responses, expressed on the majority of peripheral blood mononuclear cells (PBMC). LAIR-2 is a soluble homolog that antagonizes LAIR-1 inhibitory function by binding the same ligands. We sought to investigate whether LAIR1 and LAIR2 single nucleotide polymorphisms (SNP) are associated with differential mRNA expression levels. We analyzed 14 SNPs of LAIR1 (6) and LAIR2 (8) by mass spectrometry-based genotyping and extracted mRNA from PBMC of 177 healthy subjects, followed by quantitative assays. Four SNPs of LAIR1 and two SNPs of LAIR2 mark differential mRNA levels in healthy individuals. To verify the biological relevance of these findings, we analyzed additional 515 individuals (282 patients and 233 controls) to check if LAIR1 and LAIR2 differential mRNA expression could be related to susceptibility to pemphigus foliaceus (PF), an autoimmune blistering skin disease endemic in Brazil. Two LAIR1 variants (rs56802430 G, OR = 1.52, p = 0.0329; rs11084332 C, OR = 0.57, p = 0.0022) and one LAIR2 (rs2287828 T+, OR = 1.9, p = 0.0097) contribute to differential susceptibility to PF. Furthermore, we demonstrate interactions among four LAIR2 SNPs (rs2042287, rs2287828, rs2277974 and rs114834145). A haplotype harboring these SNPs is strongly associated with higher LAIR2 mRNA levels (4.5-fold, p = 0.0069) and with higher susceptibility to PF (OR = 4.02, p = 0.0008). We suggest that LAIR1 and LAIR2 genetic variants are associated with regulation of gene expression and variable PF susceptibility, and show indirect association of LAIR2 differential mRNA expression with PF pathogenesis. Our data demonstrate how this relatively unknown disease can add invaluable knowledge regarding the role of LAIR1 and LAIR2 in immune responses.

Pemphigus is associated with KIR3DL2 expression levels and provides evidence that KIR3DL2 may bind HLA-A3 and A11 in vivo.

Although HLA-A3 and A11 have been reported to be ligands for KIR3DL2, evidence for any in vivo relevance of this interaction is still missing. To explore the functional importance of KIR3DL2 allelic variation, we analyzed the autoimmune disease pemphigus foliaceus, previously associated (lower risk) with activating KIR genes. KIR3DL2*001 was increased in patients (odds ratio (OR) = 2.04; p = 0.007). The risk was higher for the presence of both KIR3DL2*001 and HLA-A3 or A11 (OR = 3.76, p = 0.013), providing the first evidence that HLA-A3 and A11 may interact with KIR3DL2 in vivo. The nonsynonymous single nucleotide polymorphism 1190T (rs3745902) was associated with protection (OR = 0.52, p = 0.018). This SNP results in a threonine-to-methionine substitution. Individuals who have methionine in this position exhibit a lower percentage of KIR3DL2-positive natural killer (NK) cells and also lower intensity of KIR3DL2 on expressing natural killer cells; additionally, we show that the expression of KIR3DL2 is independent of other killer cell immunoglobulin-like receptors. Pemphigus foliaceus is a very unique complex disease strongly associated with immune-related genes. It is the only autoimmune disease known to be endemic, showing a strong correlation with environmental factors. Our data demonstrate that this relatively unknown autoimmune disease may facilitate understanding of the molecular mechanisms of KIR3DL2 ligand recognition.

KIR and HLA under pressure: evidences of coevolution across worldwide populations.

KIR (killer cell immunoglobulin-like receptors) and HLA (human leukocyte antigens) are two distinct gene families with remarkable importance for human immune responses. The recognition of HLA molecules by activating and inhibitory KIR promotes a balance of signals that regulates NK cell function and is especially important for the innate defense against pathogens and early placentation. There is no documented gametic association between these two gene families and no evidence of common regulation. However, due to the critical role of KIR recognition for immunity and reproduction, the possibility of KIR-HLA combinations being under selective pressure is not surprising. In this manuscript, we first summarize the HLA-KIR system, the HLA molecules that are the putative ligands for KIR, and then we review the evidences that suggest these two gene families are coevolving as an integrated system.

High-throughput complement component 4 genomic sequence analysis with C4Investigator.

The complement component 4 gene loci, composed of the C4A and C4B genes and located on chromosome 6, encodes for complement component 4 (C4) proteins, a key intermediate in the classical and lectin pathways of the complement system. The complement system is an important modulator of immune system activity and is also involved in the clearance of immune complexes and cellular debris. C4A and C4B gene loci exhibit copy number variation, with each composite gene varying between 0 and 5 copies per haplotype. C4A and C4B genes also vary in size depending on the presence of the human endogenous retrovirus (HERV) in intron 9, denoted by C4(L) for long-form and C4(S) for short-form, which affects expression and is found in both C4A and C4B. Additionally, human blood group antigens Rodgers and Chido are located on the C4 protein, with the Rodger epitope generally found on C4A protein, and the Chido epitope generally found on C4B protein. C4A and C4B copy number variation has been implicated in numerous autoimmune and pathogenic diseases. Despite the central role of C4 in immune function and regulation, high-throughput genomic sequence analysis of C4A and C4B variants has been impeded by the high degree of sequence similarity and complex genetic variation exhibited by these genes. To investigate C4 variation using genomic sequencing data, we have developed a novel bioinformatic pipeline for comprehensive, high-throughput characterization of human C4A and C4B sequences from short-read sequencing data, named C4Investigator. Using paired-end targeted or whole genome sequence data as input, C4Investigator determines the overall gene copy numbers, as well as C4A, C4B, C4(Rodger), C4(Ch), C4(L), and C4(S). Additionally, C4Ivestigator reports the full overall C4A and C4B aligned sequence, enabling nucleotide level analysis. To demonstrate the utility of this workflow we have analyzed C4A and C4B variation in the 1000 Genomes Project Data set, showing that these genes are highly poly-allelic with many variants that have the potential to impact C4 protein function.

MHConstructor: A high-throughput, haplotype-informed solution to the MHC assembly challenge.

The extremely high levels of genetic polymorphism within the human major histocompatibility complex (MHC) limit the usefulness of reference-based alignment methods for sequence assembly. We incorporate a short read de novo assembly algorithm into a workflow for novel application to the MHC. MHConstructor is a containerized pipeline designed for high-throughput, haplotype-informed, reproducible assembly of both whole genome sequencing and target-capture short read data in large, population cohorts. To-date, no other self-contained tool exists for the generation of de novo MHC assemblies from short read data. MHConstructor facilitates wide-spread access to high quality, alignment-free MHC sequence analysis.