Search

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., (© 2024. The Author(s).)

CD8+ T cells are the dominant lymphocyte population in multiple sclerosis (MS) lesions where they are highly clonally expanded. The clonal identity, function, and antigen specificity of CD8+ T cells in MS are not well understood. Here we report a comprehensive single-cell RNA-seq and T cell receptor (TCR)-seq analysis of the cerebrospinal fluid (CSF) and blood from a cohort of treatment-naïve MS patients and control participants. A small subset of highly expanded and activated CD8+ T cells were enriched in the CSF in MS that displayed high activation, cytotoxicity and tissue-homing transcriptional profiles. Using a combination of unbiased and targeted antigen discovery approaches, MS-derived CD8+ T cell clonotypes recognizing Epstein-Barr virus (EBV) antigens and multiple novel mimotopes were identified. These findings shed vital insight into the role of CD8+ T cells in MS and pave the way towards disease biomarkers and therapeutic targets.

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., Competing Interests: Competing interests. The authors declare no competing interests.

Although B cells are implicated in multiple sclerosis (MS) pathophysiology, a predictive or diagnostic autoantibody remains elusive. In this study, the Department of Defense Serum Repository (DoDSR), a cohort of over 10 million individuals, was used to generate whole-proteome autoantibody profiles of hundreds of patients with MS (PwMS) years before and subsequently after MS onset. This analysis defines a unique cluster in approximately 10% of PwMS who share an autoantibody signature against a common motif that has similarity with many human pathogens. These patients exhibit antibody reactivity years before developing MS symptoms and have higher levels of serum neurofilament light (sNfL) compared to other PwMS. Furthermore, this profile is preserved over time, providing molecular evidence for an immunologically active preclinical period years before clinical onset. This autoantibody reactivity was validated in samples from a separate incident MS cohort in both cerebrospinal fluid and serum, where it is highly specific for patients eventually diagnosed with MS. This signature is a starting point for further immunological characterization of this MS patient subset and may be clinically useful as an antigen-specific biomarker for high-risk patients with clinically or radiologically isolated neuroinflammatory syndromes., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

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., (© 2023 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Pemphigus foliaceus (PF) is an autoimmune skin blistering disease characterized by antidesmoglein-1 IgG production, with an endemic form (EPF) in Brazil. Genetic and epigenetic factors have been associated with EPF, but its etiology is still not fully understood. To evaluate the genetic association of histone (de)acetylation-related genes with EPF susceptibility, we evaluated 785 polymorphisms from 144 genes, for 227 EPF patients and 194 controls. Carriers of HDAC4_rs4852054*A were more susceptible (OR = 1.79, p = 0.0038), whereas those with GSE1_rs13339618*A (OR = 0.57, p = 0.0011) and homozygotes for PHF21A_rs4756055*A (OR = 0.39, p = 0.0006) were less susceptible to EPF. These variants were not associated with sporadic PF (SPF) in German samples of 75 SPF patients and 150 controls, possibly reflecting differences in SPF and EPF pathophysiology. We further evaluated the expression of histone (de)acetylation-related genes in CD4 + T lymphocytes, using RNAseq. In these cells, we found a higher expression of KAT2B , PHF20, and ZEB2 and lower expression of KAT14 and JAD1 in patients with active EPF without treatment compared to controls from endemic regions. The encoded proteins cause epigenetic modifications related to immune cell differentiation and cell death, possibly affecting the immune response in patients with PF.

Background: The HLA-DRB1 gene in the major histocompatibility complex (MHC) region in chromosome 6p21 is the strongest genetic factor identified as influencing multiple sclerosis (MS) susceptibility. DNA methylation changes associated with MS have been consistently detected at the MHC region. However, understanding the full scope of epigenetic regulations of the MHC remains incomplete, due in part to the limited coverage of this region by standard whole genome bisulfite sequencing or array-based methods., Methods: We developed and validated an MHC capture protocol coupled with bisulfite sequencing and conducted a comprehensive analysis of the MHC methylation landscape in blood samples from 147 treatment naïve MS study participants and 129 healthy controls., Results: We identified 132 differentially methylated region (DMRs) within MHC region associated with disease status. The DMRs overlapped with established MS risk loci. Integration of the MHC methylome with human leukocyte antigen ( HLA ) genetic data indicate that the methylation changes are significantly associated with HLA genotypes. Using DNA methylation quantitative trait loci (mQTL) mapping and the causal inference test (CIT), we identified 643 cis-mQTL-DMRs paired associations, including 71 DMRs possibly mediating causal relationships between 55 single nucleotide polymorphisms (SNPs) and MS risk., Results: The results describe MS-associated methylation changes in MHC region and highlight the association between HLA genotypes and methylation changes. Results from the mQTL and CIT analyses provide evidence linking MHC region variations, methylation changes, and disease risk for MS., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Ma, Augusto, Montero-Martin, Caillier, Osoegawa, Cree, Hauser, Didonna, Hollenbach, Norman, Fernandez-Vina and Oksenberg.)

Studies have demonstrated that at least 20% of individuals infected with SARS-CoV-2 remain asymptomatic 1-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., (© 2023. The Author(s).)

Introduction: Primarily driven by autoreactive B cells, pemphigus foliaceus (PF) is an uncommon autoimmune blistering skin disease of sporadic occurrence worldwide. However, PF reaches a prevalence of 3% in the endemic areas of Brazil, the highest ever registered for any autoimmune disease, which indicates environmental factors influencing the immune response in susceptible individuals. We aimed to provide insights into the immune repertoire of patients with PF living in the endemic region of the disease, compared to healthy individuals from the endemic region and a non-endemic area., Methods: We characterized the B-cell repertoire in i) nontreated patients (n=5); ii) patients under immunosuppressive treatment (n=5); iii) patients in remission without treatment (n=6); and two control groups iv) from the endemic (n=6) and v) non-endemic areas in Brazil (n=4). We used total RNA extracted from peripheral blood mononuclear cells and performed a comprehensive characterization of the variable region of immunoglobulin heavy chain (IGH) in IgG and IgM using next-generation sequencing., Results: Compared to individuals from a different area, we observed remarkably lower clonotype diversity in the B-cell immune repertoire of patients and controls from the endemic area ( p < 0.02), suggesting that the immune repertoire in the endemic area is under geographically specific and intense environmental pressure. Moreover, we observed longer CDR3 sequences in patients, and we identified differential disease-specific usage of IGHV segments, including increased IGHV3-30 and decreased IGHV3-23 in patients with active disease ( p < 0.04). Finally, our robust network analysis discovered clusters of CDR3 sequences uniquely observed in patients with PF., Discussion: Our results indicate that environmental factors, in addition to disease state, impact the characteristics of the repertoire. Our findings can be applied to further investigation of the environmental factors that trigger pemphigus and expand the knowledge for identifying new targeted and more effective therapies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Calonga-Solís, Olbrich, Ott, Adelman Cipolla, Malheiros, Künstner, Farias, Camargo, Petzl-Erler, Busch, Fähnrich and Augusto.)

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., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Trichoderma spp. are usually considered safe and normally used as biocontrol and biofertilization. Safety for human health is evaluated by several tests that detect various effects such as allergenicity, toxicity, infectivity, and pathogenicity. However, they do not evaluate the effects of the agent upon the immune system. The aim of this study was to investigate the interaction between T. stromaticum spores and mammalian cells to assess the immunomodulatory potential of the spores of this fungus. First, mouse macrophage cell line J774 and human macrophages were exposed to fungal spores and analyzed for structural features, through scanning and transmission electron microscopy. Then, various analysis were performed in human macrophages as to their effect in some functional and molecular aspects of the immune system through immunocytochemistry, flow cytometry and gene expression assays. We demonstrated that T. stromaticum spores induces autophagy and autophagy-related genes (ATGs) and downmodulate inflammatory mediators, including ROS, NLRP3, the cytokines IL-1β, IL-18, IL-12 and IL-10, as well as TLR2, TLR4, miR-146b and miR-155, which may lead to an augmented susceptibility to pathogens. Our study shows the extension of damages the biofungicide Tricovab® can cause in the innate immune response. Further studies are necessary to elucidate other innate and adaptive immune responses and, consequently, the safety of this fungus when in contact with humans., Competing Interests: The authors declare no conflict of interest., (© 2022 The Authors. Published by Elsevier B.V.)

The extraordinary variation of the human leukocyte antigen (HLA) molecules is critical for diversifying antigen presentation to T cells. Coupled with the rise of novel strains and rapidly evolving immune evasion by SARS-CoV-2 proteins, HLA-mediated immunity in COVID-19 is critically important but far from being fully understood. A growing number of studies have found the association of HLA variants with different COVID-19 outcomes and that HLA genotypes associate with differential immune responses against SARS-CoV-2. Prediction studies have shown that mutations in multiple viral strains, most concentrated in the Spike protein, affect the affinity between these mutant peptides and HLA molecules. Understanding the impact of this variation on T-cell responses is critical for comprehending the immunogenic mechanisms in both natural immunity and vaccine development., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

When DNA profile comparisons between a crime scene trace and a reference sample generate correspondence, the match probability has to be estimated, so that evaluation of the strength of the forensic DNA evidence can be made. The random match probability estimations require information on allele frequencies and an adjustment factor, referred to as theta (θ) or Fst, a co-ancestry correction factor for subpopulation effects. The θ value has been standardized for urban and isolated populations, but inconsistencies have been reported when it is specifically calculated for smaller and isolated populations, including Amerindian populations. Notably, attempts to characterize forensic markers of these minor populations have been extensively limited and more conservative estimates of the correction factor may be generated for each of them. Therefore, we estimate allele frequencies of 21 autosomal STR markers used for forensic testing and calculated relevant forensic parameters for the set. In addition, we featured the possible structure of five Brazilian Amerindian populations that have been genetically isolated for centuries so we could obtain the appropriate θ value for them. The sample consisted of 319 individuals: (1) 121 Kaingang, from two communities: Ivaí (KIV=61) and Rio das Cobras (KRC=60); and (2) 198 Guaranis from three communities: Mbya from Rio das Cobras (GRC=51), Guarani Ñandeva (GND=71) and Guarani Kaiowá (GKW=76). Between Guarani populations low (Rst=0.0402, p < 10 -4 ) to high (Rst=0.1557, p < 10 -5 ) differentiation was found. Regarding Guarani and Kaingang populations, intermediate (Rst=0.0590, p < 10 -5 ) to high (Rst=0.1604, p < 10 -5 ) differentiation was found. The two Kaingang populations showed very low differentiation between them (Rst=0.0017, p = 0.27), which justifies the union of both genetic data for forensic databases and calculations. The combined power of discrimination (PD) and the combined power of exclusion (PE) were calculated for each population, demonstrating the usefulness of this set of markers in forensic and kinship analysis regarding these populations. Considering the demographic heterogeneity of Amerindian populations in general, the Fst mean value (0.03) was evaluated regarding 43 different indigenous populations from the Americas, including Guaranis and Kaingangs. This result confirms the adequacy of the standardized θ value for the forensic random match probability estimations involving Amerindian populations., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The long search for the environmental trigger of the endemic pemphigus foliaceus (EPF, fogo selvagem) has not yet resulted in any tangible findings. Here, we searched for genetic associations and the differential expression of host genes involved in early viral infections and innate antiviral defense. Genetic variants could alter the structure, expression sites, or levels of the gene products, impacting their functions. By analyzing 3063 variants of 166 candidate genes in 227 EPF patients and 194 controls, we found 12 variants within 11 genes associated with differential susceptibility (p < 0.005) to EPF. The products of genes TRIM5, TPCN2, EIF4E, EIF4E3, NUP37, NUP50, NUP88, TPR, USP15, IRF8, and JAK1 are involved in different mechanisms of viral control, for example, the regulation of viral entry into the host cell or recognition of viral nucleic acids and proteins. Only two of nine variants were also associated in an independent German cohort of sporadic PF (75 patients, 150 controls), aligning with our hypothesis that antiviral host genes play a major role in EPF due to a specific virus−human interaction in the endemic region. Moreover, CCL5, P4HB, and APOBEC3G mRNA levels were increased (p < 0.001) in CD4+ T lymphocytes of EPF patients. Because there is limited or no evidence that these genes are involved in autoimmunity, their crucial role in antiviral responses and the associations that we observed support the hypothesis of a viral trigger for EPF, presumably a still unnoticed flavivirus. This work opens new frontiers in searching for the trigger of EPF, with the potential to advance translational research that aims for disease prevention and treatment.

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., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

BACKGROUNDVaccine-elicited adaptive immunity is a prerequisite for control of SARS-CoV-2 infection. Multiple sclerosis (MS) disease-modifying therapies (DMTs) differentially target humoral and cellular immunity. A comprehensive comparison of the effects of MS DMTs on SARS-CoV-2 vaccine-specific immunity is needed, including quantitative and functional B and T cell responses.METHODSSpike-specific Ab and T cell responses were measured before and following SARS-CoV-2 vaccination in a cohort of 80 study participants, including healthy controls and patients with MS in 6 DMT groups: untreated and treated with glatiramer acetate (GA), dimethyl fumarate (DMF), natalizumab (NTZ), sphingosine-1-phosphate (S1P) receptor modulators, and anti-CD20 mAbs. Anti-spike-Ab responses were assessed by Luminex assay, VirScan, and pseudovirus neutralization. Spike-specific CD4+ and CD8+ T cell responses were characterized by activation-induced marker and cytokine expression and tetramer.RESULTSAnti-spike IgG levels were similar between healthy control participants and patients with untreated MS and those receiving GA, DMF, or NTZ but were reduced in anti-CD20 mAb- and S1P-treated patients. Anti-spike seropositivity in anti-CD20 mAb-treated patients was correlated with CD19+ B cell levels and inversely correlated with cumulative treatment duration. Spike epitope reactivity and pseudovirus neutralization were reduced in anti-CD20 mAb- and S1P-treated patients. Spike-specific CD4+ and CD8+ T cell reactivity remained robust across all groups, except in S1P-treated patients, in whom postvaccine CD4+ T cell responses were attenuated.CONCLUSIONThese findings from a large cohort of patients with MS exposed to a wide spectrum of MS immunotherapies have important implications for treatment-specific COVID-19 clinical guidelines.FUNDINGNIH grants 1K08NS107619, K08NS096117, R01AI159260, R01NS092835, R01AI131624, and R21NS108159; NMSS grants TA-1903-33713 and RG1701-26628; Westridge Foundation; Chan Zuckerberg Biohub; Maisin Foundation.

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., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Although the genus Trichoderma is widely used as a biocontrol agent in crops, little is known about its potential impact on the human immune system. In mice, our group has shown that exposition to T. asperelloides spores lead to reduced neutrophil counts in the peripheral blood and in the peritoneal cavity. In addition, T. stromaticum spores produced an inflammatory infiltrate on mice lungs, reducing the levels of IFN-γ and IL-10 cytokines, reactive oxygen species, and receptors of microbial patterns. Here we demonstrate that the interaction of human peripheral neutrophils with T. stromaticum spores also leads to a reduced release of neutrophil extracellular traps (NETs) after induction with the NET-inducer agent phorbol 12-myristate 13-acetate. This interaction also reduced the expression levels of multiple microRNAs, such as miR-221, miR-222, miR-223 and miR-27a, as well as genes related to NETs, such as ELANE, MPO and PADI4. Furthermore, T. stromaticum spores affected the expression of the genes SOCS3, TLR4, CSNK2A1, GSDMD, and NFFKBIA, related to the activation of inflammatory immune responses in neutrophils. Overall, our results suggest T. stromaticum as a potential NET inhibitor and as an immunomodulatory agent. Since this fungus is used as biocontrol in crops, our findings point to the importance of advancing our knowledge on the effects of this bioagent on the human immune system. Finally, the study of the active compounds produced by the fungus is also important for the prospection of new drugs that could be used to block the exacerbation of inflammatory immune responses present in several human diseases., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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., (© 2021. The Author(s).)

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 (p corr < 0.01) and nine with endemic PF (p corr < 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, log 2 FC = 1.3), while TSBP1-AS1 was downregulated in PV (p = 0.029, log 2 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, log 2 FC = -2.14). TSBP1-AS1 and HCG27 were also observed downregulated in CD19 + cells of endemic PF (p < 0.01, log 2 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., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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.

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., Competing Interests: The authors have declared that no competing interests exist.

Pemphigus foliaceus (PF) is an autoimmune blistering disease of the skin, clinically characterized by erosions and, histopathologically, by acantholysis. PF is endemic in the Brazilian Central-Western region. Numerous single nucleotide polymorphisms (SNPs) have been shown to affect the susceptibility for PF, including SNPs at long non-coding RNA (lncRNA) genes, which are known to participate in many physiological and pathogenic processes, such as autoimmunity. Here, we investigated whether the genetic variation of immune-related lncRNA genes affects the risk for endemic and sporadic forms of PF. We analysed 692 novel SNPs for PF from 135 immune-related lncRNA genes in 227 endemic PF patients and 194 controls. The SNPs were genotyped by Illumina microarray and analysed by applying logistic regression at additive model, with correction for sex and population structure. Six associated SNPs were also evaluated in an independent German cohort of 76 sporadic PF patients and 150 controls. Further, we measured the expression levels of two associated lncRNA genes (LINC-PINT and LY86-AS1) by quantitative PCR, stratified by genotypes, in peripheral blood mononuclear cells of healthy subjects. We found 27 SNPs in 11 lncRNA genes associated with endemic PF (p < .05 without overlapping with protein-coding genes). Among them, the LINC-PINT SNP rs10228040*A (OR = 1.47, p = .012) was also associated with increased susceptibility for sporadic PF (OR = 2.28, p = .002). Moreover, the A+ carriers of LY86-AS1*rs12192707 mark lowest LY86-AS1 RNA levels, which might be associated with a decreasing autoimmune response. Our results suggest a critical role of lncRNA variants in immunopathogenesis of both PF endemic and sporadic forms., (© 2021 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

The KIR (killer-cell immunoglobulin-like receptor ) region is characterized by structural variation and high sequence similarity among genes, imposing technical difficulties for analysis. We undertook the most comprehensive study to date of KIR genetic diversity in a large population sample, applying next-generation sequencing in 2,130 United States European-descendant individuals. Data were analyzed using our custom bioinformatics pipeline specifically designed to address technical obstacles in determining KIR genotypes. Precise gene copy number determination allowed us to identify a set of uncommon gene-content KIR haplotypes accounting for 5.2% of structural variation. In this cohort, KIR2DL4 is the framework gene that most varies in copy number (6.5% of all individuals). We identified phased high-resolution alleles in large multi-locus insertions and also likely founder haplotypes from which they were deleted. Additionally, we observed 250 alleles at 5-digit resolution, of which 90 have frequencies ≥1%. We found sequence patterns that were consistent with the presence of novel alleles in 398 (18.7%) individuals and contextualized multiple orphan dbSNPs within the KIR complex. We also identified a novel KIR2DL1 variant, Pro151Arg, and demonstrated by molecular dynamics that this substitution is predicted to affect interaction with HLA-C. No previous studies have fully explored the full range of structural and sequence variation of KIR as we present here. We demonstrate that pairing high-throughput sequencing with state-of-art computational tools in a large cohort permits exploration of all aspects of KIR variation including determination of population-level haplotype diversity, improving understanding of the KIR system, and providing an important reference for future studies., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Amorim, Augusto, Nemat-Gorgani, Montero-Martin, Marin, Shams, Dandekar, Caillier, Parham, Fernández-Viña, Oksenberg, Norman and Hollenbach.)

Behçet disease is a multi-system disease associated with human leukocyte antigen (HLA) class I polymorphism. High-resolution next-generation sequencing (NGS) with haplotype analysis has not been performed previously for this disease. Sixty Egyptian patients diagnosed according to the International Study Group (ISG) criteria for Behçet disease and 160 healthy geographic and ethnic-matched controls were genotyped for HLA class I loci (HLA-A, B, C). For HLA class II loci (DRB1, DRB3/4/5, DQA1, DQB1, DPA1, DPB1), 40 control samples were genotyped. High-resolution HLA genotyping was performed using NGS and the results were analyzed. Clinical manifestations were oral ulcers (100%), genital ulcers (100%), eye (55%) and neurological (28%) and vascular involvement (35%). HLA-B*51:08 [odds ratio (OR) = 19·75, 95% confidence interval (CI) = 6·5-79; P < 0·0001], HLA-B*15:03 (OR = 12·15, 95% CI = 3·7-50·7; P < 0·0001), HLA-C*16:02 (OR = 6·53, 95% CI = 3-14; P < 0·0001), HLA-A*68:02 (OR = 3·14, 95% CI = 1·1-8·9; P < 0·01) were found to be associated with Behçet disease, as were HLA-DRB1*13:01 and HLA-DQB1*06:03 (OR = 3·39, 95% CI = 0·9-18·9; P = 0·04 for both). By contrast, HLA-A*03:01 (OR = 0·13, 95% CI = 0-0·8; P = 0·01) and HLA-DPB1*17:01 were found to be protective (OR = 0·27, 95% CI = 0·06-1·03; P = 0·02). We identified strong linkage disequilibrium between HLA-B*51:08 and C*16:02 and A*02:01 in a haplotype associated with Behçet disease. HLA-B*51:08 was significantly associated with legal blindness (OR = 2·98, 95% CI = 1·06-8·3; P = 0·01). In Egyptian Behçet patients, HLA-B*51:08 is the most common susceptibility allele and holds poor prognosis for eye involvement., (© 2021 British Society for Immunology.)

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., (© 2020 John Wiley & Sons Ltd.)

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 ( p c = 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., (Copyright © 2020 by The American Association of Immunologists, Inc.)

Regulation of NK cell activity is mediated through killer-cell immunoglobulin-like receptors (KIR) ability to recognize human leukocyte antigen (HLA) class I molecules as ligands. Interaction of KIR and HLA is implicated in viral infections, autoimmunity, and reproduction and there is growing evidence of the coevolution of these two independently segregating gene families. By leveraging KIR and HLA-C data from 1000 Genomes consortium we observed that the KIR2DL1 variant rs2304224 * T is associated with lower expression of HLA-C in individuals carrying the ligand HLA-C2 ( p = 0.0059). Using flow cytometry, we demonstrated that this variant is also associated with higher expression of KIR2DL1 on the NK cell surface ( p = 0.0002). Next, we applied next generation sequencing to analyze KIR2DL1 sequence variation in 109 Euro and 75 Japanese descendants. Analyzing the extended haplotype homozygosity, we show signals of positive selection for rs4806553 * G and rs687000 * G , which are in linkage disequilibrium with rs2304224 * T . Our results suggest that lower expression of HLA-C2 ligands might be compensated for higher expression of the receptor KIR2DL1 and bring new insights into the coevolution of KIR and HLA ., (Copyright © 2020 Vargas, Dourado, Amorim, Ho, Calonga-Solís, Issler, Marin, Beltrame, Petzl-Erler, Hollenbach and Augusto.)

Objectives: Genes of host immunity play an important role in disease pathogenesis and are determinants of clinical courses of infections, including hepatitis B virus (HBV). Killer-cell immunoglobulin-like receptor (KIR), expressed on the surface of natural killer cells (NK), regulate NK cell cytotoxicity by interacting with human leukocyte antigen (HLA) class I molecules and are candidates for influencing the course of HBV. This study evaluated whether variations in KIR gene content and HLA-C ligands are associated with HBV and with the development of liver cirrhosis and hepatocellular carcinoma., Methods: A Vietnamese study cohort (HBV n = 511; controls n = 140) was genotyped using multiplex sequence-specific polymerase chain reaction (PCR-SSP) followed by melting curve analysis., Results: The presence of the functional allelic group of KIR2DS4 was associated with an increased risk of chronic HBV (OR = 1.86, p corr = 0.02), while KIR2DL2+HLA-C1 (OR = 0.62, p corr = 0.04) and KIR2DL3+HLA-C1 (OR = 0.48, p corr = 0.04) were associated with a decreased risk. The pair KIR2DL3+HLA-C1 was associated with liver cirrhosis (OR = 0.40, p corr = 0.01). The presence of five or more activating KIR variants was associated with hepatocellular carcinoma (OR = 0.53, p corr = 0.04)., Conclusions: KIR gene content variation and combinations KIR-HLA influence the outcome of HBV infection., (Copyright © 2020 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Biological control agents (BCA) are an alternative to chemical pesticides and an emerging strategy to safely eliminate plant pathogens. Trichoderma spp. are the most common fungi used as BCAs. They produce spores that are released into the air and can potentially interact with immune system of mammals. We previously showed that Trichoderma affects expression of genes encoding pattern recognition receptors (PRRs) and cytokines in mice. PRRs are involved in the recognition of microorganisms and can lead to pro-tumoral signaling. Here, we evaluated if mice injected with low doses of murine melanoma exhibited increased development of lung tumor when treated with conidia of T. stromaticum . Mice treated with T. stromaticum and inoculated with B16-F10 melanoma cells exhibited significant increase in tumor uptake ( p = 0.006) and increased number of visible nodules in the lungs ( p = 0.015). We also analyzed mRNA expression levels of genes encoding PRRs in lung of mice exposed to T. stromaticum and demonstrated that mice treated with T. stromaticum conidia exhibited lower expression levels of Clec7a and increased expression of Tlr4 (toll like receptor 4) compared to non-treated controls. The expression levels of Clec7a and Tlr2 were increased in mice treated with T. stromaticum and inoculated with murine melanoma compared to controls only inoculated with melanoma. Our results demonstrate that intranasal exposition to T. stromaticum increases tumor in the B16-F10 model, which may raise concerns regarding the safety of its use in agriculture., (Copyright © 2020 Santos, Costa, de Freitas, de Oliveira, Santos, Silva, Santos, Dias, de Carvalho, Augusto and dos Santos.)

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., Competing Interests: Declaration of competing interest The authors declare that they have no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

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.

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease that occurs sporadically across the globe and is endemic in Brazil. Keratinocyte adhesion loss (acantholysis) is associated with high levels of anti-desmoglein 1 IgG autoantibodies, but the role of cell death is poorly understood in PF. Current evidence disqualifies apoptosis as the major cell death mechanism and no other process has yet been investigated. To approach the role of variation in genes responsible for cell death pathways in pemphigus susceptibility, we systematically investigated the frequencies of 1,167 polymorphisms from genes encoding products of all 12 well-established cell death cascades (intrinsic and extrinsic apoptosis, necrosis, necroptosis, ferroptosis, pyroptosis, parthanatos, entotic, NETotic, lysosome-dependent, autophagy-dependent, and immunogenic). By multivariate logistic regression, we compared allelic and genotypic frequencies of 227 PF patients and 194 controls obtained by microarray hybridization. We found 10 variants associated with PF ( p < 0.005), belonging to six cell death pathways: apoptosis ( TNF, TRAF2, CD36 , and PAK2 ), immunogenic cell death ( EIF2AK3, CD47 , and SIRPA ), necroptosis ( TNF and TRAF2 ), necrosis ( RAPGEF3 ), parthanatos ( HK1 ), and pyroptosis ( PRKN ). Five polymorphisms were associated with susceptibility: TNF rs1800630 * A ( OR = 1.9, p = 0.0003), CD36 rs4112274 * T ( OR = 2.14, p = 0.0015), CD47 rs12695175 * G ( OR = 1.77, p = 0.0043), SIRPA rs6075340 * A/A ( OR = 2.75, p = 0.0009), and HK1 rs7072268 * T ( OR = 1.48, p = 0.0045). Other five variants were associated with protection: TRAF2 rs10781522 * G ( OR = 0.64, p = 0.0014), PAK2 rs9325377 * A/A ( OR = 0.48, p = 0.0023), EIF2AK3 rs10167879 * T ( OR = 0.48, p = 0.0007), RAPGEF3 rs10747521 * A/A ( OR = 0.42, p = 0.0040), and PRKN rs9355950 * C ( OR = 0.57, p = 0.0004). Through functional annotation, we found that all associated alleles, with the exception of PRKN rs9355950 * C , were previously associated with differential gene expression levels in healthy individuals (mostly in skin and peripheral blood). Further functional validation of these genetic associations may contribute to the understanding of PF etiology and to the development of new drugs and therapeutic regimens for the disease., (Copyright © 2019 Bumiller-Bini, Cipolla, Spadoni, Augusto, Petzl-Erler, Beltrame and Boldt.)

Background: Pemphigus foliaceus (PF) is an epidermal autoimmune disease, characterized by the presence of autoantibodies against the desmosomal protein desmoglein 1. Genetic and environmental factors contribute to PF, a complex disease that is endemic in Brazil and Colombia and neighbouring countries, and in Tunisia. Long noncoding RNAs (lncRNAs) may participate in gene regulation by interacting with DNA, proteins and other RNAs. Dysregulation of lncRNAs has recently been recognized as an important coplayer in the onset or progression of complex diseases. In addition, single-nucleotide polymorphisms (SNPs) located in lncRNA genes have been associated with differential risk to cancer, autoimmunity and infection., Objectives: Here, we aimed to investigate whether SNPs in lncRNA genes are associated with differential susceptibility to endemic PF., Materials and Methods: We integrated data from the lncRNA SNP database with genome-wide genotype data obtained for 229 patients and 6681 controls. We tested the association between endemic PF and 2080 SNPs located in lncRNAs applying logistic regression., Results: The most significantly associated SNP was rs7144332 (OR = 1·63, P = 2·8 × 10 -6 ), located in the lncRNA gene AL110292·1. Results for five other SNPs were suggestive of association (P < 0·001). In silico analysis indicated that five of the six SNPs impact transcription, three may influence lncRNA's secondary structure, and three may alter microRNA-lncRNA interactions., Conclusions: We showed, for the first time, that variation in lncRNA genes may influence pemphigus pathogenesis. Our findings highlight the importance of lncRNA variation in autoimmune and possibly other complex diseases and suggest polymorphisms for functional validation., (© 2019 British Association of Dermatologists.)

Even though immunoglobulins are critical for immune responses and human survival, the diversity of the immunoglobulin heavy chain gene ( IGH ) is poorly known and mostly characterized only by serological methods. Moreover, this genomic region is not well-covered in genomic databases and genome-wide association studies due to particularities that impose technical difficulties for its analysis. Therefore, the IGH gene has never been systematically sequenced across populations. Here, we deliver an unprecedented and comprehensive characterization of the diversity of the IGHG1, IGHG2 , and IGHG3 gene segments, which encode the constant region of the most abundant circulating immunoglobulins: IgG1, IgG2, and IgG3, respectively. We used Sanger sequencing to analyze 357 individuals from seven different Brazilian populations, including five Amerindian, one Japanese-descendant and one Euro-descendant population samples. We discovered 28 novel IGHG alleles and provided evidence that some of them may have been originated by gene conversion between common alleles of different gene segments. The rate of synonymous substitutions was significantly higher than the rate of the non-synonymous substitutions for IGHG1 and IGHG2 ( p = 0.01 and 0.03, respectively), consistent with purifying selection. Fay and Wu's test showed significant negative values for most populations ( p < 0.001), which indicates that positive selection in an adjacent position may be shaping IGHG variation by hitchhiking of variants in the vicinity, possibly the regions that encode the Ig variable regions. This study shows that the variation in the IGH gene is largely underestimated. Therefore, exploring its nucleotide diversity in populations may provide valuable information for comprehension of its evolution, its impact on diseases and vaccine research.

The killer-cell immunoglobulin-like receptor ( KIR ) region comprises a fast-evolving family of genes that encode receptors for natural killer (NK) cells and have crucial role in host defense. Evolution of KIR was examined in the context of the human genome. Gene-content diversity and single nucleotide polymorphisms (SNP) in the KIR genes and flanking regions were compared to >660,000 genome-wide SNPs in over 800 individuals from 52 populations of the human genome diversity panel (HGDP). KIR allelic diversity was further examined using next generation sequencing in a subset of 56 individuals. We identified the SNP rs587560 located in KIR3DL3 as a marker of KIR2DL2 and KIR2DL3 and, consequently, Cen A and Cen B haplotypes. We also show that combinations of two KIR2DL4 SNPs ( rs35656676 and rs592645 ) distinguish KIR3DL1 from KIR3DS1 and also define the major KIR3DL1 high- and low-expressing alleles lineages. Comparing the diversity of the SNPs within the KIR region to remainder of the genome, we observed a high diversity for the centromeric KIR region consistent with balancing selection ( p < 0.01); in contrast, centromeric KIR diversity is significantly reduced in East Asian populations ( p < 0.01), indicating purifying selection. By analyzing SNP haplotypes in a region spanning ~500 kb that includes the KIR cluster, we observed evidence of strong positive selection in Africa for high-expressing KIR3DL1 alleles, favored over the low-expressing alleles ( p < 0.01). In sharp contrast, the strong positive selection ( p < 0.01) that we also observed in the telomeric KIR region in Oceanic populations tracked with a high frequency of KIR3DS1 . In addition, we demonstrated that worldwide frequency of high-expression KIR3DL1 alleles was correlated with virus with virus (r = 0.64, p < 10 -6 ) and protozoa (r = 0.69, p < 10 -6 ) loads, which points to selection globally on KIR3DL1 high-expressing alleles attributable to pathogen exposure.

The HLA-G and MICA genes are stimulated under inflammatory conditions and code for soluble (sMICA and sHLA-G) or membrane-bound molecules that exhibit immunomodulatory properties. It is still unclear whether they would have a synergistic or antagonistic effect on the immunomodulation of the inflammatory response, such as in chronic kidney disease (CKD), contributing to a better prognosis after the kidney transplantation. In this study, we went from genetic to plasma analysis, first evaluating the polymorphism of MICA, NKG2D and HLA-G in a cohort from Southern Brazil, subdivided in a control group of individuals (n = 75), patients with CKD (n = 94), and kidney-transplant (KT) patients (n = 64). MICA, NKG2D and HLA-G genotyping was performed by polymerase chain reaction with specific oligonucleotide probes, Taqman and Sanger sequencing, respectively. Levels of soluble forms of MICA and HLA-G were measured in plasma with ELISA. Case-control analysis showed that the individuals with haplotype HLA-G*01:01/UTR-4 have a lower susceptibility to develop chronic kidney disease (OR = 0.480; p = 0.032). Concerning the group of kidney-transplant patients, the HLA-G genotypes +3010 GC (rs1710) and +3142 GC (rs1063320) were associated with higher risk for allograft rejection (OR = 5.357; p = 0.013 and OR = 5.357, p = 0.013, respectively). Nevertheless, the genotype +3010 GG (OR = 0.136; p = 0.041) was associated with kidney allograft acceptance, suggesting that it is a protection factor for rejection. In addition, the phenotypic analysis revealed higher levels of sHLA-G (p = 0.003) and sMICA (p < 0.001) in plasma were associated with the development of CKD. For patients who were already under chronic pathological stress and underwent a kidney transplant, a high sMICA (p = 0.001) in pre-transplant proved to favor immunomodulation and allograft acceptance. Even so, the association of genetic factors with differential levels of soluble molecules were not evidenced, we displayed a synergistic effect of sMICA and sHLA-G in response to inflammation. This increase was observed in CKD patients, that when undergo transplantation, had this previous amount of immunoregulatory molecules as a positive factor for the allograft acceptance., Competing Interests: The authors have declared that no competing interests exist.

The human leukocyte antigen (HLA) are the most polymorphic genes in the human genome. Because of their importance for antigen recognition, HLA molecules play a central role in host defense and graft rejection upon transplantation. The aim of this study was to characterize allelic diversity of the classical HLA genes HLA-A, -B, -C, -DRA, -DRB1, -DQA1, -DQB1, -DPA1, -DPB1, and the non-classical class I genes HLA-E, -F and -G at high-resolution for a population of predominantly European ancestry from Curitiba, Brazil. Genotyping of 108 individuals was performed by next-generation sequencing on the MiSeq platform and also by Sanger sequencing. The genotype distributions of all loci were in accordance with Hardy-Weinberg equilibrium (P > 0.05) and a total of 202 HLA variants at second field resolution were observed for the 12 loci. The strongest linkage disequilibrium (r 2 = 1.0, P < 10 -5 ) was observed for the following pairs of alleles: HLA-B*42:01:01 ~ HLA-DRB1*03:02:01; HLA-B*14:02:01 ~ HLA-C*08:02:01; B*42:01:01 ~ HLA-C*17:01:01; HLA-DRB1*03:01:01 ~ HLA-DQB1*02:01:01 ~ DRB1*03:01:01 ~ HLA-DQB1*02:01:01; DRB1*13:01:01~ HLA-DQB1*06:03:01 and HLA-DRB1*09:01:02 ~ HLA-DQA1*03:02. This is the first study to characterize all 12 HLA genes at high resolution in a single population. On the basis of the allelic frequencies of worldwide populations and principal component analysis, we confirmed the similarity of the study population to European and other Euro-descendant populations., (© 2019 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

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., (© 2018 John Wiley & Sons Ltd.)

We here describe the KIR gene-content genotypes in the highly admixed population from Ilhéus, a city located in Bahia State, Northeast Brazil. Bahia is the State with the largest contribution of African ethnicity in comparison to other Brazilian areas and this is the first report of a population from this region. A total of 32 KIR gene-content genotypes have been found, from which 18 were observed in single individuals. The distinct KIR gene-content distribution observed in Ilhéus, along with the elevated diversity of genotypes point to importance of describing KIR polymorphism in unique populations from Brazil., (Copyright © 2018. Published by Elsevier Inc.)

Killer cell immunoglobulin-like receptor (KIR) genes encode cell surface molecules that recognize HLA molecules and modulate the activity of natural killer (NK) cells. KIR genes exhibit presence and absence polymorphism, which generates a variety of gene-content haplotypes in worldwide populations. KIR gene-content variation is implicated in many diseases and is also important for placentation and transplantation. Because of the complexity of KIR polymorphism, variation in this family is still mostly studied at the gene-content level, even with the advent of next-generation sequencing (NGS) methods. Gene-content determination is generally expensive and/or time-consuming. To overcome these difficulties, we developed a method based on multiplex polymerase chain reaction with specific sequence primers (PCR-SSP) followed by melting curve analysis that allows cost-effective, precise and fast generation of results. Our method was 100% concordant with a gel-based method and 99.9% concordant with presence and absence determination by NGS. The limit of detection for accurate typing was 30 ng of DNA (0.42 μM) with 260/230 and 260/280 ratios as low as 0.19 and of 0.44. In addition, we developed a user-friendly Java-based computational application called killerPeak that interprets the raw data generated by Viia7 or QuantStudio 7 quantitative PCR machines and reliably exports the final genotyping results in spreadsheet file format. The combination of a reliable method that requires low amount of DNA with an automated interpretation of results allows scaling the KIR genotyping in large cohorts with reduced turnaround time., (© 2018 John Wiley & Sons A/S. Published by John Wiley & Sons Ltd.)

Fungi from the widely distributed genus Trichoderma are of great biotechnological interest, being currently used in a vast range of applications. Here, we report that high-molecular weight fraction (HWF) derived from Trichoderma asperelloides ethanolic extract exhibits antibiotic activity against staphylococcal biofilms. The antibacterial and anti-biofilm properties of T. asperelloides extracts were evaluated by well-established assays in Staphylococcus aureus ATCC strains (29213 and 6538) and in one clinical isolate from bovine mastitis. The HWF from T. asperelloides eradicated S. aureus by causing substantial matrix de-structuring and biomass reduction (p < 10-5) at concentrations as low as 2.3 μg mL-1. Additionally, we present ultra-structure analysis by the use of scanning electron microscopy as well as transmission microscopy, which showed that T. asperelloides killed cells through cell wall and membrane disturbance. Remarkably, the HWF from T. asperelloides killed S. aureus and eradicated its biofilms in a greater performance than gentamicin (p < 10-5), a known potent antibiotic against S. aureus. Our results indicate that extract from T. asperelloides may represent a promising candidate for the development of new antibiotics against gram-positive bacteria., Competing Interests: The authors have declared that no competing interests exist.

Skin blisters of pemphigus foliaceus (PF) present concomitant deposition of autoantibodies and components of the complement system (CS), whose gene polymorphisms are associated with susceptibility to different autoimmune diseases. To investigate these in PF, we evaluated 992 single-nucleotide polymorphisms (SNPs) of 44 CS genes, genotyped through microarray hybridization in 229 PF patients and 194 controls. After excluding SNPs with minor allele frequency <1%, out of Hardy-Weinberg equilibrium in controls or in strong linkage disequilibrium ( r 2  ≥ 0.8), 201 SNPs remained for logistic regression. Polymorphisms of 11 genes were associated with PF. MASP1 encodes a crucial serine protease of the lectin pathway (rs13094773: OR = 0.5, p  = 0.0316; rs850309: OR = 0.23, p  = 0.03; rs3864098: OR = 1.53, p  = 0.0383; rs698104: OR = 1.52, p  = 0.0424; rs72549154: OR = 0.55, p  = 0.0453). C9 (rs187875: OR = 1.46, p  = 0.0189; rs700218: OR = 0.12, p  = 0.0471) and C8A (rs11206934: OR = 4.02, p  = 0.0323) encode proteins of the membrane attack complex (MAC) and C5AR1 (rs10404456: OR = 1.43, p  = 0.0155), a potent anaphylatoxin-receptor. Two encode complement regulators: MAC-blocking CD59 (rs1047581: OR = 0.62, p  = 0.0152) and alternative pathway-blocking CFH (rs34388368: OR = 2.57, p  = 0.0195). One encodes opsonin: C3 (rs4807895: OR = 2.52, p  = 0.0239), whereas four encode receptors for C3 fragments: CR1 (haplotype with rs6656401: OR = 1.37, p  = 0.0382), CR2 (rs2182911: OR = 0.23, p  = 0.0263), ITGAM (CR3, rs12928810: OR = 0.66, p  = 0.0435), and ITGAX (CR4, rs11574637: OR = 0.63, p  = 0.0056). Associations reinforced former findings, regarding differential gene expression, serum levels, C3, and MAC deposition on lesions. Deregulation of previously barely noticed processes, e.g., the lectin and alternative pathways and opsonization-mediated phagocytosis, also modulate PF susceptibility. The results open new crucial avenues for understanding disease etiology and may improve PF treatment through additional therapeutic targets.