The major histocompatibility complex class I gene repertoire was investigated in a large panel of rhesus macaques of Chinese origin. As observed in Indian animals, subjects of Chinese derivation display Mamu-B gene copy number variation, and the sum of expressed genes varies among haplotypes. In addition, these genes display differential transcription levels. The majority of the Mamu-B alleles discovered during this investigation appear to be unique for the population studied. Only one particular Mamu-B haplotype is shared between Indian and Chinese animals, and it must have been present in the progenitor stock. Hence, the data highlight the fact that most allelic polymorphism, and most of the Mamu-B haplotypes themselves, are of relatively recent origin and were most likely generated after the separation of the Indian and Chinese rhesus macaque populations.

T cell receptor (TCR) nucleotide sequences are often generated during analyses of T cell responses to pathogens or autoantigens. The most important region of the TCR is the third complementarity-determining region (CDR3) whose nucleotide sequence is unique to each T cell clone. The CDR3 interacts with the peptide and thus is important for recognizing pathogen or autoantigen epitopes. While conventions exist for identifying the various TCR chains, there is a lack of a concise nomenclature that would identify both the amino acid translation and nucleotide sequence of the CDR3. This deficiency makes the comparison of published TCR genetic and proteomic information difficult. To enhance information sharing among different databases and to facilitate computational assessment of clonotypic T cell repertoires, we propose a clonotype nomenclature. The rules for generating a clonotype identifier are simple and easy to follow, and have a built-in error-checking system. The identifier includes the V and J region, the CDR3 length as well as its human or mouse origin. The framework of this naming system could also be expanded to the B cell receptor.

T cell epitopes derived from polymorphic proteins or from proteins encoded by alternative reading frames (ARFs) play an important role in (tumor) immunology. Identification of these peptides is successfully performed with mass spectrometry. In a mass spectrometry-based approach, the recorded tandem mass spectra are matched against hypothetical spectra generated from known protein sequence databases. Commonly used protein databases contain a minimal level of redundancy, and thus, are not suitable data sources for searching polymorphic T cell epitopes, either in normal or ARFs. At the same time, however, these databases contain much non-polymorphic sequence information, thereby complicating the matching of recorded and theoretical spectra, and increasing the potential for finding false positives. Therefore, we created a database with peptides from ARFs and peptide variation arising from single nucleotide polymorphisms (SNPs). It is based on the human mRNA sequences from the well-annotated reference sequence (RefSeq) database and associated variation information derived from the Single Nucleotide Polymorphism Database (dbSNP). In this process, we removed all non-polymorphic information. Investigation of the frequency of SNPs in the dbSNP revealed that many SNPs are non-polymorphic "SNPs". Therefore, we removed those from our dedicated database, and this resulted in a comprehensive high quality database, which we coined the Human Short Peptide Variation Database (HSPVdb). The value of our HSPVdb is shown by identification of the majority of published polymorphic SNP- and/or ARF-derived epitopes from a mass spectrometry-based proteomics workflow, and by a large variety of polymorphic peptides identified as potential T cell epitopes in the HLA-ligandome presented by the Epstein-Barr virus cells.

Stem cells are unique in that they possess both the capacity to self-renew and thereby maintain their original pool as well as the capacity to differentiate into mature cells. In the past number of years, transcriptional profiling of enriched stem cell populations has been extensively performed in an attempt to identify a universal stem cell gene expression signature. While stem-cell-specific transcripts were identified in each case, this approach has thus far been insufficient to identify a universal group of core "stemness" genes ultimately responsible for self-renewal and multipotency. Similarly, in the hematopoietic system, comparisons of transcriptional profiles between different hematopoietic cell stages have had limited success in revealing core genes ultimately responsible for the initiation of differentiation and lineage specification. Here, we propose that the combined use of transcriptional profiling and genetic linkage analysis, an approach called "genetical genomics", can be a valuable tool to assist in the identification of genes and gene networks that specify "stemness" and cell fate decisions. We review past studies of hematopoietic cells that utilized transcriptional profiling and/or genetic linkage analysis, and discuss several potential future applications of genetical genomics.

Proteolysis is the general term to describe the process of protein degradation into peptides. Proteasomes are the main actors in cellular proteolysis, and their activity can be measured in in vitro digestion experiments. However, in vivo proteolysis can be different than what is measured in these experiments if other proteases participate or if proteasomal activity is different in vivo. The in vivo proteolysis can be measured only indirectly, by the analysis of peptides presented on MHC-I molecules. MHC-I presented peptides are protected from further degradation, thus enabling an indirect view on the underlying in vivo proteolysis. The ligands presented on different MHC-I molecules enable different views on this process; in combination, they might give a complete picture. Based on in vitro proteasome-only digestions and MHC-I ligand data, different proteolysis predictors have been developed. With new in vitro digestion and MHC-I ligand data sets, we benchmarked how well these predictors capture in vitro proteasome-only activity and in vivo whole-cell proteolysis, respectively. Even though the in vitro proteasome digestion patterns were best captured by methods trained on such data (ProteaSMM and NetChop 20S), the in vivo whole-cell proteolysis was best predicted by a method trained on MHC-I ligand data (NetChop Cterm). Follow-up analysis showed that the likely source of this difference is the activity from proteases other than the proteasome, such as TPPII. This non-proteasomal in vivo activity is captured by NetChop Cterm and should be taken into account in MHC-I ligand predictions.

The present analysis made attempts to resolve discrepancies in descriptions of eutherian tumor necrosis factor ligand genes implicated in cell signalling pathways, as well as in major physiological and pathological processes. Among 455 potential coding sequences, the eutherian comparative genomic analysis protocol annotated 211 complete coding sequences using public genomic sequence assemblies. The most comprehensive third party data gene data set first described 8 superclusters of eutherian tumor necrosis factor ligand genes, including 18 major gene clusters. The integrated gene annotations, phylogenetic analysis, and protein molecular evolution analysis proposed new classification and nomenclature of eutherian tumor necrosis factor ligand genes, as new framework of future experiments.

Bony fish encode multiple multi-gene families of membrane receptors that are comprised of immunoglobulin (Ig) domains and are predicted to function in innate immunity. One of these families, the diverse immunoglobulin (Ig) domain-containing protein (DICP) genes, maps to three chromosomal loci in zebrafish. Most DICPs possess one or two Ig ectodomains and include membrane-bound and secreted forms. Membrane-bound DICPs include putative inhibitory and activating receptors. Recombinant DICP Ig domains bind lipids with varying specificity, a characteristic shared with mammalian CD300 and TREM family members. Numerous DICP transcripts amplified from different lines of zebrafish did not match the zebrafish reference genome sequence suggesting polymorphic and haplotypic variation. The expression of DICPs in three different lines of zebrafish has been characterized employing PCR-based strategies. Certain DICPs exhibit restricted expression in adult tissues whereas others are expressed ubiquitously. Transcripts of a subset of DICPs can be detected during embryonic development suggesting roles in embryonic immunity or other developmental processes. Transcripts representing 11 previously uncharacterized DICP sequences were identified. The assignment of two of these sequences to an unplaced genomic scaffold resulted in the identification of an alternative DICP haplotype that is linked to a MHC class I Z lineage haplotype on zebrafish chromosome 3. The linkage of DICP and MHC class I genes also is observable in the genomes of the related grass carp (Ctenopharyngodon idellus) and common carp (Cyprinus carpio) suggesting that this is a shared character with the last common Cyprinidae ancestor.

The main function of HLA class I molecules is to present pathogen-derived peptides to cytotoxic T lymphocytes. This function is assumed to drive the maintenance of an extraordinary amount of polymorphism at each HLA locus, providing an immune advantage to heterozygote individuals capable to present larger repertories of peptides than homozygotes. This seems contradictory, however, with a reduced diversity at individual HLA loci exhibited by some isolated populations. This study shows that the level of functional diversity predicted for the two HLA-A and HLA-B genes considered simultaneously is similar (almost invariant) between 46 human populations, even when a reduced diversity exists at each locus. We thus propose that HLA-A and HLA-B evolved through a model of joint divergent asymmetric selection conferring all populations an equivalent immune potential. The distinct pattern observed for HLA-C is explained by its functional evolution towards killer cell immunoglobulin-like receptor (KIR) activity regulation rather than peptide presentation.

As the only native insular Newfoundland canid between the extinction of the wolf in the 1930s and the recent arrival of coyotes, the red fox (Vulpes vulpes deletrix Bangs 1898) poses interesting questions about genetic distinctiveness and the post-glacial colonization history of the island's depauperate mammalian fauna. Here, we characterized genetic variability at the major histocompatibility complex (MHC) class II DR β1 domain (DRB1) locus in 28 red foxes from six sampling localities island-wide and compared it with mitochondrial control region (CR) diversity and DRB1 diversity in other canids. Our goals were to describe novel DRB1 alleles in a new canid population and to make inferences about the role of selection in maintaining their diversity. As in numerous studies of vertebrates, we found an order-of-magnitude higher nucleotide diversity at the DRB1 locus compared with the CR and significantly positive nonsynonymous-to-synonymous substitution ratios, indicative of selection in the distant past. Although the evidence is weaker, the Ewens-Watterson test of neutrality and the geographical distribution of variation compared with the CR suggest a role for selection over the evolutionary timescale of populations. We report the first genetic data from the DRB1 locus in the red fox and establish baseline information regarding immunogenetic variation in this island canid population which should inform continued investigations of population demography, adaptive genetic diversity, and wildlife disease in red foxes and related species.

Killer cell immunoglobulin-like receptors (KIR) are the most polymorphic receptors of natural killer (NK) cells. Their activity diversifies the functions of NK cells in the antiviral immune response, so the presence of certain KIR may affect transmission of HIV-1. The aim of the study was to evaluate the influence of KIR genes on the susceptibility to HIV-1 infection in the Polish population depending on the route of exposure. We determined the frequencies of activating (2DS1, 2DS2, 2DS3, 2DS4f, 2DS4del, 2DS5, 3DS1) and inhibitory (2DL1, 2DL2, 2DL3, 2DL5, 3DL1) KIRs in HIV-1-positive patients (n = 459), individuals exposed to HIV-1 but uninfected (EU, n = 118) and in uninfected, healthy blood donors (BD, n = 98). Analysis was performed using stepwise logistic regression. Apart from KIRs, CCR5-∆32, and CCR2-64I, alleles were also analyzed, as we knew or suspected that these features could affect susceptibility to HIV infection. The regression confirmed the protective effect of CCR5-∆32 (OR = 0.25, p = 0.006) and CCR2-64I (OR = 0.59, p = 0.032) against HIV infection. Among KIR genes, 2DL3 was found to be a protective factor (OR = 0.30, p = 0.015). A similar effect was seen for 3DS1 but only in intravenous drug users (IDUs) (OR = 0.30, p = 0.019), not in sexually exposed people. 2DL5 was found to be a factor facilitating HIV infection (OR = 2.13, p = 0.013). A similar effect was observed for 2DL2 but only in females (OR = 2.15, p = 0.040), and 2DS1 in IDUs (OR = 3.03, p = 0.022). Our results suggest a beneficial role of KIR3DS1 and 2DL3 supporting resistance to HIV infection and a harmful effect of 2DS1, 2DL5, and 2DL2 genes promoting HIV acquisition.

The human major histocompatibility complex (MHC) is contained within about 4 Mb on the short arm of chromosome 6 and is recognised as the most variable region in the human genome. The primary aim of the MHC Haplotype Project was to provide a comprehensively annotated reference sequence of a single, human leukocyte antigen-homozygous MHC haplotype and to use it as a basis against which variations could be assessed from seven other similarly homozygous cell lines, representative of the most common MHC haplotypes in the European population. Comparison of the haplotype sequences, including four haplotypes not previously analysed, resulted in the identification of >44,000 variations, both substitutions and indels (insertions and deletions), which have been submitted to the dbSNP database. The gene annotation uncovered haplotype-specific differences and confirmed the presence of more than 300 loci, including over 160 protein-coding genes. Combined analysis of the variation and annotation datasets revealed 122 gene loci with coding substitutions of which 97 were non-synonymous. The haplotype (A3-B7-DR15; PGF cell line) designated as the new MHC reference sequence, has been incorporated into the human genome assembly (NCBI35 and subsequent builds), and constitutes the largest single-haplotype sequence of the human genome to date. The extensive variation and annotation data derived from the analysis of seven further haplotypes have been made publicly available and provide a framework and resource for future association studies of all MHC-associated diseases and transplant medicine.

Comparisons of MHC gene content and diversity among closely related species can provide insights into the evolutionary mechanisms shaping immune system variation. After chimpanzees and bonobos, gorillas are humans' closest living relatives; but in contrast, relatively little is known about the structure and variation of gorilla MHC class I genes (Gogo). Here, we combined long-range amplifications and long-read sequencing technology to analyze full-length MHC class I genes in 35 gorillas. We obtained 50 full-length genomic sequences corresponding to 15 Gogo-A alleles, 4 Gogo-Oko alleles, 21 Gogo-B alleles, and 10 Gogo-C alleles including 19 novel coding region sequences. We identified two previously undetected MHC class I genes related to Gogo-A and Gogo-B, respectively, thereby illustrating the potential of this approach for efficient and highly accurate MHC genotyping. Consistent with their phylogenetic position within the hominid family, individual gorilla MHC haplotypes share characteristics with humans and chimpanzees as well as orangutans suggesting a complex history of the MHC class I genes in humans and the great apes. However, the overall MHC class I diversity appears to be low further supporting the hypothesis that gorillas might have experienced a reduction of their MHC repertoire.

HLA-G is known for its strictly restricted tissue distribution. HLA-G expression could be detected in immune privileged organs and many tumor entities such as leukemia, multiple myeloma, and non-Hodgkin and Hodgkin's lymphoma. This functional variability from mediation of immune tolerance to facilitation of tumor immune evasion strategies might translate to a differential NK cell inhibition between immune-privileged organs and tumor cells. The biophysical invariability of the HLA-G heavy chain and its contrary diversity in immunity implicates a strong influence of the bound peptides on the pHLA-G structure. The aim was to determine if HLA-G displays a tissue-specific peptide repertoire. Therefore, using soluble sHLA-G technology, we analyzed the K562 and HDLM-2 peptide repertoires. Although both cell lines possess a comparable proteome and recruit HLA-G-restricted peptides through the same peptide-loading pathway, the peptide features appear to be cell specific. HDLM-2 derived HLA-G peptides are anchored by an Arg at p1 and K562-derived peptides are anchored by a Lys. At p2, no anchor motif could be determined while peptides were anchored at pΩ with a Leu and showed an auxiliary anchor motif Pro at p3. To appreciate if the peptide anchor alterations are due to a cell-specific differential peptidome, we performed analysis of peptide availability within the different cell types. Yet, the comparison of the cell-specific proteome and HLA-G-restricted ligandome clearly demonstrates a tissue-specific peptide selection by HLA-G molecules. This exclusive and unexpected observation suggests an exquisite immune function of HLA-G.

The IPD-MHC Database represents the official repository for non-human major histocompatibility complex (MHC) sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. IPD-MHC gathers allelic MHC class I and class II sequences from classical and non-classical MHC loci from various non-human animals including pets, farmed and experimental model animals. So far, Atlantic salmon and rainbow trout are the only teleost fish species with MHC class I and class II sequences present. For the remaining teleost or ray-finned species, data on alleles originating from given classical locus is scarce hampering their inclusion in the database. However, a fast expansion of sequenced genomes opens for identification of classical loci where high-throughput sequencing (HTS) will enable typing of allelic variants in a variety of new teleost or ray-finned species. HTS also opens for large-scale studies of salmonid MHC diversity challenging the current database nomenclature and analysis tools. Here we establish an Illumina approach to identify allelic MHC diversity in Atlantic salmon, using animals from an endangered wild population, and alter the salmonid MHC nomenclature to accommodate the expected sequence expansions.

The major histocompatibility complex (MHC) encodes proteins that are central for antigen presentation and pathogen elimination. MHC class I (MHC-I) genes have attracted a great deal of interest among researchers in ecology and evolution and have been partly characterized in a wide range of bird species. So far, the main focus has been on species within the bird orders Galliformes and Passeriformes, while Charadriiformes remain vastly underrepresented with only two species studied to date. These two Charadriiformes species exhibit striking differences in MHC-I characteristics and MHC-I diversity. We therefore set out to study a third species within Charadriiformes, the Icelandic subspecies of black-tailed godwits (Limosa limosa islandica). This subspecies is normally confined to parasite-poor environments, and we hence expected low MHC diversity. MHC-I was partially characterized first using Sanger sequencing and then using high-throughput sequencing (MiSeq) in 84 individuals. We verified 47 nucleotide alleles in open reading frame with classical MHC-I characteristics, and each individual godwit had two to seven putatively classical MHC alleles. However, in contrast to previous MHC-I data within Charadriiformes, we did not find any evidence of alleles with low sequence diversity, believed to represent non-classical MHC genes. The diversity and divergence of the godwits MHC-I genes to a large extent fell between the previous estimates within Charadriiformes. However, the MHC genes of the migratory godwits had few sites subject to positive selection, and one possible explanation could be a low exposure to pathogens.

In comparison to humans and chimpanzees, gorillas show low diversity at MHC class I genes (Gogo), as reflected by an overall reduced level of allelic variation as well as the absence of a functionally important sequence motif that interacts with killer cell immunoglobulin-like receptors (KIR). Here, we use recently generated large-scale genomic sequence data for a reassessment of allelic diversity at Gogo-C, the gorilla orthologue of HLA-C. Through the combination of long-range amplifications and long-read sequencing technology, we obtained, among the 35 gorillas reanalyzed, three novel full-length genomic sequences including a coding region sequence that has not been previously described. The newly identified Gogo-C*03:01 allele has a divergent recombinant structure that sets it apart from other Gogo-C alleles. Domain-by-domain phylogenetic analysis shows that Gogo-C*03:01 has segments in common with Gogo-B*07, the additional B-like gene that is present on some gorilla MHC haplotypes. Identified in ~ 50% of the gorillas analyzed, the Gogo-C*03:01 allele exclusively encodes the C1 epitope among Gogo-C allotypes, indicating its important function in controlling natural killer cell (NK cell) responses via KIR. We further explored the hypothesis whether gorillas experienced a selective sweep which may have resulted in a general reduction of the gorilla MHC class I repertoire. Our results provide little support for a selective sweep but rather suggest that the overall low Gogo class I diversity can be best explained by drastic demographic changes gorillas experienced in the ancient and recent past.

Systemic lupus erythematosus (SLE) is a complex disorder. Genetic association studies of complex disorders suffer from the following three major issues: phenotypic heterogeneity, false positive (type I error), and false negative (type II error) results. Hence, genes with low to moderate effects are missed in standard analyses, especially after statistical corrections. OASIS is a novel linkage disequilibrium clustering algorithm that can potentially address false positives and negatives in genome-wide association studies (GWAS) of complex disorders such as SLE. OASIS was applied to two SLE dbGAP GWAS datasets (6077 subjects; ∼0.75 million single-nucleotide polymorphisms). OASIS identified three known SLE genes viz. IFIH1, TNIP1, and CD44, not previously reported using these GWAS datasets. In addition, 22 novel loci for SLE were identified and the 5 SLE genes previously reported using these datasets were verified. OASIS methodology was validated using single-variant replication and gene-based analysis with GATES. This led to the verification of 60% of OASIS loci. New SLE genes that OASIS identified and were further verified include TNFAIP6, DNAJB3, TTF1, GRIN2B, MON2, LATS2, SNX6, RBFOX1, NCOA3, and CHAF1B. This study presents the OASIS algorithm, software, and the meta-analyses of two publicly available SLE GWAS datasets along with the novel SLE genes. Hence, OASIS is a novel linkage disequilibrium clustering method that can be universally applied to existing GWAS datasets for the identification of new genes.

Covid-19 has caused worldwide devastation. IFIH1 is a pattern recognition receptor that senses coronavirus RNA and triggers interferon production as a first line of viral immune defense. The role of IFIH1 polymorphism, rs1990760 (C>T; aaA946T) in the epidemiology of viral infection is well studied, and the minor allele T resists viral infection. Knock-in mice with mutated IFIH1 protein (946T) for this allele have enhanced interferon production and protection from lethal viral infection. The minor allele frequency (Tmaf) varies widely from Africans (0.06 to 0.35) to Chinese (0.19 to 0.23) to Caucasians (0.56 to 0.69). During the initial days of infection when the social restrictions were not imposed, I show that the infection rate in Italy was lower as expected from its higher Tmaf (0.56) than that in China (Tmaf for southern China, 0.23). The infection rate in the USA and Spain was intermediate between those two countries despite higher Caucasian overall Tmaf (0.69), perhaps due to a more admixed African population in these countries. These analyses suggest that African-Americans and Chinese with low Tmaf of rs1990760 are more vulnerable to SARS-COV2 infection, apart from other genetic factors or socioeconomic conditions in these population. Taken together, an IFN-beta supplement might aid in preventing COVID-19 infection and help in development of herd immunity.

Males and females often exhibit differences in behaviour, life histories, and ecology, many of which are typically reflected in their brains. Neuronal protection and maintenance include complex processes led by the microglia, which also interacts with metabolites such as hormones or immune components. Despite increasing interest in sex-specific brain function in laboratory animals, the significance of sex-specific immune activation in the brain of wild animals along with the variables that could affect it is widely lacking. Here, we use the Kentish plover (Charadrius alexandrinus) to study sex differences in expression of immune genes in the brain of adult males and females, in two wild populations breeding in contrasting habitats: a coastal sea-level population and a high-altitude inland population in China. Our analysis yielded 379 genes associated with immune function. We show a significant male-biased immune gene upregulation. Immune gene expression in the brain did not differ in upregulation between the coastal and inland populations. We discuss the role of dosage compensation in our findings and their evolutionary significance mediated by sex-specific survival and neuronal deterioration. Similar expression profiles in the coastal and inland populations suggest comparable genetic control by the microglia and possible similarities in pathogen pressures between habitats. We call for further studies on gene expression of males and females in wild population to understand the implications of immune function for life-histories and demography in natural systems.

NOD-like receptors (NLRs) are sensors of pathogen-associated molecular patterns with critical roles in the control of immune responses and programmed cell death. Recent studies have revealed inter-species differences in mammalian innate immune genes and a particular degeneration of nucleic acid sensing pathways in pangolins, which are currently investigated as potential hosts for zoonotic pathogens. Here, we used comparative genomics to determine which NLR genes are conserved or lost in pangolins and related mammals. We show that NOD2, which is implicated in sensing bacterial muramyl dipeptide and viral RNA, is a pseudogene in pangolins, but not in any other mammalian species investigated. NLRC4 and NAIP are absent in pangolins and canine carnivorans, suggesting convergent loss of cytoplasmic sensing of bacterial flagellin in these taxa. Among NLR family pyrin domain containing proteins (NLRPs), skin barrier-related NLRP10 has been lost in pangolins after the evolutionary divergence from Carnivora. Strikingly, pangolins lack all NLRPs associated with reproduction (germ cells and embryonic development) in other mammals, i.e., NLRP2, 4, 5, 7, 8, 9, 11, 13, and 14. Taken together, our study shows a massive degeneration of NLR genes in pangolins and suggests that these endangered mammals may have unique adaptations of innate immunity and reproductive cell biology.