Comparative genomic analysis of the MHC: the evolution of class I duplication blocks, diversity and complexity from shark to man

Universal signatures of transposable element compartmentalization across eukaryotic genomes

The evolutionary mechanisms that drive the emergence of genome architecture remain poorly understood but can now be assessed with unprecedented power due to the massive accumulation of genome assemblies spanning phylogenetic diversity1,2. Transposable elements (TEs) are a rich source of large-effect mutations since they directly and indirectly drive genomic structural variation and changes in gene expression3. Here, we demonstrate universal patterns of TE compartmentalization across eukaryotic genomes spanning ~1.7 billion years of evolution, in which TEs colocalize with gene families under strong predicted selective pressure for dynamic evolution and involved in specific functions. For non-pathogenic species these genes represent families involved in defense, sensory perception and environmental interaction, whereas for pathogenic species, TE-compartmentalized genes are highly enriched for pathogenic functions. Many TE-compartmentalized gene families display signatures of positive selection at the molecular level. Furthermore, TE-compartmentalized genes exhibit an excess of high-frequency alleles for polymorphic TE insertions in fruit fly populations. We postulate that these patterns reflect selection for adaptive TE insertions as well as TE-associated structural variants. This process may drive the emergence of a shared TE-compartmentalized genome architecture across diverse eukaryotic lineages.

Genetic Diversity and Sequence Conservation of Peptide-Binding Regions of MHC Class I Genes in Pig, Cattle, Chimpanzee, and Human

Comparative analyses of MHC gene diversity and evolution across different species could offer valuable insights into the evolution of MHC genes. Intra- and inter-species sequence diversity and conservation of 12 classical major histocompatibility complex (MHC) class I genes from cattle, chimpanzees, pigs, and humans was analyzed using 20 representative allelic groups for each gene. The combined analysis of paralogous loci for each species revealed that intra-locus amino-acid sequence variations in the peptide-binding region (PBR) of MHC I genes did not differ significantly between species, ranging from 8.44% for SLA to 10.75% for BoLA class I genes. In contrast, intraspecies differences in the non-PBRs of these paralogous genes were more pronounced, varying from 4.59% for SLA to 16.89% for HLA. Interestingly, the Shannon diversity index and rate of nonsynonymous substitutions for PBR were significantly higher in SLA and BoLA than those in Patr and HLA. Analysis of peptide-binding pockets across all analyzed MHC class I genes of the four species indicated that pockets A and E showed the lowest and highest diversity, respectively. The estimated divergence times suggest that primate and artiodactyl MHC class I genes diverged 60.41 Mya, and BoLA and SLA genes diverged 35.34 Mya. These results offer new insights into the conservation and diversity of MHC class I genes in various mammalian species.

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (β2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The β2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Targeting phagocytosis to enhance antitumor immunity

Many patients with metastatic or treatment-resistant cancer have experienced improved outcomes after immunotherapy that targets adaptive immune checkpoints. However, innate immune checkpoints, which can hinder the detection and clearance of malignant cells, are also crucial in tumor-mediated immune escape and may also serve as targets in cancer immunotherapy. In this review, we discuss the current understanding of immune evasion by cancer cells via disruption of phagocytic clearance, and the potential effects of blocking phagocytosis checkpoints on the activation of antitumor immune responses. We propose that a more effective combination immunotherapy strategy could be to exploit tumor-intrinsic processes that inhibit key innate immune surveillance processes, such as phagocytosis, and incorporate both innate and adaptive immune responses for treating patients with cancer.

Non-random Codon Usage of Synonymous and Non-synonymous Mutations in the Human HLA-A Gene

The structure and function of human leucocyte antigen (HLA-A) is well known and is an extremely variable protein. From the public HLA-A database, we chose 26 high frequency HLA-A alleles (45% of sequenced alleles). Using five arbitrary references from these alleles, we analyzed synonymous mutations at the third codon position (sSNP₃) and non-synonymous mutations (NSM). Both mutation types showed non-random locations of 29 sSNP₃ codons and 71 NSM codons in the five reference lists. Most sSNP₃ codons show identical mutation types with many mutations resulting from cytosine deamination. We proposed 23 ancestral parents of sSNP₃ in five reference sequences using conserved parents in five unidirectional codons and 18 majority parents in reciprocal codons. These 23 proposed ancestral parents show exclusive codon usage of G₃ or C₃ parents located on both DNA strands that mutate to A₃ or T₃ variants mostly (76%) by cytosine deamination The sSNP₃ and NSM show clear separation of the two variant types with most sSNP₃ located in conserved areas in exons 2, 3 and 4, compared to most NSM appearing in two Variable Areas with no sSNP₃ in the latter parts of exons 2 (α1) and 3 (α2). The Variable Areas contain NSM (polymorphic) residues at the center of the groove that bind the foreign peptide. We find distinctly different mutation patterns in NSM codons from those of sSNP₃. Namely, G-C to A-T mutation frequency was much smaller, suggesting that evolutional pressures of deamination and other mechanisms applied to the two areas are significantly different.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

Human leukocyte antigen super-locus: nexus of genomic supergenes, SNPs, indels, transcripts, and haplotypes

The human Major Histocompatibility Complex (MHC) or Human Leukocyte Antigen (HLA) super-locus is a highly polymorphic genomic region that encodes more than 140 coding genes including the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. In recent years, MHC genotyping and haplotyping using Sanger sequencing and next-generation sequencing (NGS) methods have produced many hundreds of genomic sequences of the HLA super-locus for comparative studies of the genetic architecture and diversity between the same and different haplotypes. In this special issue on 'The Current Landscape of HLA Genomics and Genetics', we provide a short review of some of the recent analytical developments used to investigate the SNP polymorphisms, structural variants (indels), transcription and haplotypes of the HLA super-locus. This review highlights the importance of using reference cell-lines, population studies, and NGS methods to improve and update our understanding of the mechanisms, architectural structures and combinations of human MHC genomic alleles (SNPs and indels) that better define and characterise haplotypes and their association with various phenotypes and diseases.

Comparative Genomics of the Major Histocompatibility Complex (MHC) of Felids

This review summarizes the current knowledge on the major histocompatibility complex (MHC) of the family Felidae. This family comprises an important domestic species, the cat, as well as a variety of free-living felids, including several endangered species. As such, the Felidae have the potential to be an informative model for studying different aspects of the biological functions of MHC genes, such as their role in disease mechanisms and adaptation to different environments, as well as the importance of genetic diversity for conservation issues in free-ranging or captive populations. Despite this potential, the current knowledge on the MHC in the family as a whole is fragmentary and based mostly on studies of the domestic cat and selected species of big cats. The overall structure of the domestic cat MHC is similar to other mammalian MHCs following the general scheme “centromere-MHC class I-MHC class III-MHC class II” with some differences in the gene contents. An unambiguously defined orthologue of the non-classical class I HLA-E gene has not been identified so far and the class II DQ and DP genes are missing or pseudogenized, respectively. A comparison with available genomes of other felids showed a generally high level of structural and sequence conservation of the MHC region. Very little and fragmentary information on in vitro and/or in vivo biological functions of felid MHC genes is available. So far, no association studies have indicated effects of MHC genetic diversity on a particular disease. No information is available on the role of MHC class I molecules in interactions with Natural Killer (NK) cell receptors or on the putative evolutionary interactions (co-evolution) of the underlying genes. A comparison of complex genomic regions encoding NK cell receptors (the Leukocyte Receptor Complex, LRC and the Natural Killer Cell Complex, NKC) in the available felid genomes showed a higher variability in the NKC compared to the LRC and the MHC regions. Studies of the genetic diversity of domestic cat populations and/or specific breeds have focused mainly on DRB genes. Not surprisingly, higher levels of MHC diversity were observed in stray cats compared to pure breeds, as evaluated by DRB sequencing as well as by MHC-linked microsatellite typing. Immunogenetic analysis in wild felids has only been performed on MHC class I and II loci in tigers, Namibian leopards and cheetahs. This information is important as part of current conservation tasks to assess the adaptive potential of endangered wild species at the human-wildlife interface, which will be essential for preserving biodiversity in a functional ecosystem.

A short HLA-DRA isoform binds the HLA-DR2 heterodimer on the outer domain of the peptide-binding site

The human leukocyte antigen (HLA) locus encodes a large group of proteins governing adaptive and innate immune responses. Among them, HLA class II proteins form α/β heterodimers on the membrane of professional antigen-presenting cells (APCs), where they display both, self and pathogen-derived exogenous antigens to CD4⁺ T lymphocytes. We have previously shown that a shorter HLA-DRA isoform (sHLA-DRA) lacking 25 amino acids can be presented onto the cell membrane via binding to canonical HLA-DR2 heterodimers. Here, we employed atomistic molecular dynamics simulations to decipher the binding position of sHLA-DRA and its structural impact on functional regions of the HLA-DR2 molecule. We show that a loop region exposed only in the short isoform (residues R69 to G83) is responsible for binding HLA-DR2 on the outer domain of the peptide-binding site, and experimentally validated the critical role of F76 in mediating such interaction. Additionally, sHLA-DRA allosterically modifies the peptide-binding pocket conformation. In summary, this study unravels key molecular mechanisms underlying sHLA-DRA function, providing important insights into the role of full-length proteins in structural modulation of HLA class II receptors.

Avian blood parasite richness decreases with major histocompatibility complex class I loci number

Major histocompatibility complex (MHC) genes are among the most polymorphic in the vertebrate genome. The high allele diversity is believed to be maintained primarily by sexual and pathogen-mediated balancing selection. The number of MHC loci also varies greatly across vertebrates, most notably across birds. MHC proteins play key roles in presenting antigens on the cell surface for recognition by T cells, with class I proteins specifically targeting intracellular pathogens. Here, we explore the hypothesis that MHC class I diversity (measured as loci number) coevolves with haemosporidian parasite burden of the host. Using data on 54 bird species, we demonstrate that high-MHC class I diversity is associated with significantly lower richness of Plasmodium, Haemoproteus as well as overall haemosporidian parasite lineages, the former thus indicating more efficient protection against intracellular pathogens. Nonetheless, the latter associations were only detected when MHC diversity was assessed using cloning and not 454 pyrosequencing-based studies, nor across all genotyping methods combined. Our results indicate that high-MHC class I diversity might play a key role in providing qualitative resistance against diverse haemosporidian parasites in birds, but further clarification is needed for the origin of contrasting results when using different genotyping methods for MHC loci quantification.

Mathematical modeling and stochastic simulations suggest that low-affinity peptides can bisect MHC1-mediated export of high-affinity peptides into "early"- and "late"-phases

The peptide loading complex (PLC) is a multi-protein complex of the endoplasmic reticulum (ER) which optimizes major histocompatibility I (MHC1)-mediated export of intracellular high-affinity peptides. Whilst, the molecular biology of MHC1-mediated export is well supported by empirical data, the stoichiometry, kinetics and spatio-temporal profile of the participating molecular entities are a matter of considerable debate. Here, a low-affinity peptide-driven (LAPD)-model of MHC1-mediated high-affinity peptide export is formulated, implemented, analyzed and simulated. The model is parameterized in terms of the contribution of the shunt reaction to the concentration of exportable MHC1. Theoretical analyses and simulation studies of the model suggest that low-affinity peptides can bisect MHC1-mediated export of high-affinity peptides into time-dependent distinct “early”- and “late”-phases. The net exportable MHC1 (eM1β(t)) is a function of the retrograde (rM1β(t))- and anterograde (aM1β(t))-derived fractions. The “early”-phase is dominated by the contribution of the retrograde/recyclable (rM1β≈61%,aM1β≈39%) pathway to exportable MHC1, is characterized by Tapasin-mediated peptide-editing and is ATP-independent. The “late”-phase on the other hand, is characterized by de novo PLC-assembly, rapid disassembly and a significant contribution of the anterograde pathway to exportable MHC1 (rM1β≈21%,aM1β≈79%). The shunt reaction is rate limiting and may integrate peptide translocation with PLC-assembly/disassembly thereby, regulating peptide export under physiological and pathological (viral infections, dysplastic alterations) conditions.

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Decellularized xenogeneic scaffolds in transplantation and tissue engineering: Immunogenicity versus positive cell stimulation

Seriously compromised function of some organs can only be restored by transplantation. Due to the shortage of human donors, the need to find another source of organs is of primary importance. Decellularized scaffolds of non-human origin are being studied as highly potential biomaterials for tissue engineering. Their biological nature and thus the ability to provide a naturally-derived environment for human cells to adhere and grow highlights their great advantage in comparison to synthetic scaffolds. Nevertheless, since every biomaterial implanted in the body generates immune reaction, studying the interaction of the scaffold with the surrounding tissues is necessary. This review aims to summarize current knowledge on the immunogenicity of semi-xenografts involved in transplantation. Moreover, positive aspects of the interaction between xenogeneic scaffold and human cells are discussed, focusing on specific roles of proteins associated with extracellular matrix in cell adhesion and signalling.

Molecular Characterization of MHC Class I Genes in Four Species of the Turdidae Family to Assess Genetic Diversity and Selection

In vertebrate animals, the molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immunity. MHC class I deals with intracellular pathogens (virus) in birds. MHC class I diversity depends on the consequence of local and global environment selective pressure and gene flow. Here, we evaluated the MHC class I gene in four species of the Turdidae family from a broad geographical area of northeast China. We isolated 77 MHC class I sequences, including 47 putatively functional sequences and 30 pseudosequences from 80 individuals. Using the method based on analysis of cloned amplicons (n = 25) for each species, we found two and seven MHC I sequences per individual indicating more than one MHC I locus identified in all sampled species. Results revealed an overall elevated genetic diversity at MHC class I, evidence of different selection patterns among the domains of PBR and non-PBR. Alleles are found to be divergent with overall polymorphic sites per species ranging between 58 and 70 (out of 291 sites). Moreover, transspecies alleles were evident due to convergent evolution or recent speciation for the genus. Phylogenetic relationships among MHC I show an intermingling of alleles clustering among the Turdidae family rather than between other passerines. Pronounced MHC I gene diversity is essential for the existence of species. Our study signifies a valuable tool for the characterization of evolutionary relevant difference across a population of birds with high conservational concerns.

Study on the roles of melatonin in regulating dermal fibroblast growth in Liaoning cashmere goats by transcriptome sequencing

In this study, the genes related to the Downy growth of Liaoning cashmere goats were screened for their expression with simultaneous melatonin administration, so as to investigate the effects of target genes on the proliferation of skin fibroblasts in this animal species. Genes related to the villus growth of skin fibroblasts were screened by in vitro transcriptome sequencing and verified by qPCR. In addition, gene overexpression and interference were used to study the effects of target genes on the proliferation of skin fibroblasts. Groups treated with M1_24H, M2_24H and M2_72H exhibited significant differences compared with the control group. Among them, the differentially expressed transcripts in the M2_72H group were significantly enriched in the TNF and NOD-like receptor signaling pathways, which are associated with the villus. In addition, eight differentially expressed genes were screened from the TNF and the NOD-like receptor signaling pathways. Verification by qPCR showed that the expression of TNF-α, IL-6, TNFAIP3, PYCARD and NFKBIA genes were significantly upregulated, which was consistent with the sequencing results. Melatonin treatments can significantly lead to an increase in the expression of IL-6 and TNF-α genes. Besides, melatonin treatments can affect cashmere growth in Liaoning cashmere goats by regulating several signaling pathways, including TNF, NOD-like receptor and NF-κB.

Multi-tissue transcriptomic analysis reveals that L-methionine supplementation maintains the physiological homeostasis of broiler chickens than D-methionine under acute heat stress

The objective of this study was to compare the effects of supplementation with two methionine isoforms, L-methionine (L-Met) or D-methionine (D-Met), on transcriptome expression in broiler chickens under acute heat stress. A total of 240 one-day-old chicks were randomly assigned to one of four treatments in a 2 × 2 factorial arrangement: thermo-neutral vs. acute heat-stress and L-Met vs. D-Met supplementation. On day 14, the heat-stressed group was exposed to 32°C for 5 h, while the others remained at 25°C. Six chicks were randomly selected per treatment and total RNA was isolated from whole blood, ileum, and liver tissues. Two RNA samples from each tissue of each treatment group were randomly selected and pooled in equal amounts. A total of 1.87 billion raw reads obtained from 36 samples (four treatments × three tissues × three composited replicates) were mapped to the reference genome build (Gallus_gallus-5.0) and used to identify differentially expressed genes (DEGs) using DESeq2. Functional enrichment of DEGs was tested using DAVID. Comparing the two isoforms of supplemented methionine, two, three, and ten genes were differentially expressed (> 1 or < -1 log₂ fold change) in whole blood, ileum, and liver, respectively. A total of 38, 71, and 16 genes were differentially expressed in response to the interaction between heat stress and Met isoforms in the blood, ileum, and liver, respectively. Three-tissue-specific DEGs were functionally enriched for regulation of cholesterol homeostasis and metabolism, glucose metabolism, and vascular patterning. Chicks fed with L-Met had lower immune (e.g., IL4I1 and SERPINI1) and intestinal angiogenic responses (e.g., FLT1 and FGD5), and stable glucose and lipid metabolism (e.g., PCK1 and LDLR) under heat stress conditions. In conclusion, unlike D-Met, L-Met supplementation seems to help maintain physiological homeostasis and enhances cellular defense systems against external stresses like high environmental temperature.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

A splice acceptor variant in HLA-DRA affects the conformation and cellular localization of the class II DR alpha-chain

Class II human leucocyte antigen (HLA) proteins are involved in the immune response by presenting pathogen-derived peptides to CD4⁺ T lymphocytes. At the molecular level, they are constituted by α/β-heterodimers on the surface of professional antigen-presenting cells. Here, we report that the acceptor variant (rs8084) in the HLA-DRA gene mediates the transcription of an alternative version of the α-chain lacking 25 amino acids in its extracellular domain. Molecular dynamics simulations suggest this isoform undergoes structural refolding which in turn affects its stability and cellular trafficking. The short HLA-DRA isoform cannot reach the cell surface, although it is still able to bind the corresponding β-chain. Conversely, it remains entrapped within the endoplasmic reticulum where it is targeted for degradation. Furthermore, we demonstrate that the short isoform can be transported to the cell membrane via interactions with the peptide-binding site of canonical HLA heterodimers. Altogether, our findings indicate that short HLA-DRA functions as a novel intact antigen for class II HLA molecules.

Marsupial chromosomics: bridging the gap between genomes and chromosomes

Marsupials have unique features that make them particularly interesting to study, and sequencing of marsupial genomes is helping to understand their evolution. A decade ago, it was a huge feat to sequence the first marsupial genome. Now, the advances in sequencing technology have made the sequencing of many more marsupial genomes possible. However, the DNA sequence is only one component of the structures it is packaged into: chromosomes. Knowing the arrangement of the DNA sequence on each chromosome is essential for a genome assembly to be used to its full potential. The importance of combining sequence information with cytogenetics has previously been demonstrated for rapidly evolving regions of the genome, such as the sex chromosomes, as well as for reconstructing the ancestral marsupial karyotype and understanding the chromosome rearrangements involved in the Tasmanian devil facial tumour disease. Despite the recent advances in sequencing technology assisting in genome assembly, physical anchoring of the sequence to chromosomes is required to achieve a chromosome-level assembly. Once chromosome-level assemblies are achieved for more marsupials, we will be able to investigate changes in the packaging and interactions between chromosomes to gain an understanding of the role genome architecture has played during marsupial evolution.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Transcriptomic analysis of chicken immune response to infection of different doses of Newcastle disease vaccine

Newcastle disease virus (NDV) is a contagious poultry paramyxovirus, leading to substantial economic losses to the poultry industry. Here, RNA-seq was carried out to investigate the altered expression of immune-related genes in chicken thymus within 96 h in response to NDV infection. In NDV-infected chicken thymus tissues, comparative transcriptome analysis revealed 1386 differentially expressed genes (DEGs) at 24 h with 989 up- and 397 down-regulated genes, 728 DEGs at 48 h with 567 up- and 161 down-regulated genes, 1514 DEGs at 72 h with 1016 up- and 498 down-regulated genes, and 1196 DEGs at 96 h with 522 up- and 674 down-regulated genes, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that these candidate targets mainly participate in biological processes or biochemical, metabolic and signal transduction processes. Notably, there is large enrichment in biological processes, cell components and metabolic processes, which may be related to NDV pathogenicity. In addition, the expression of five immune-related DEGs identified by RNA-seq was validated by quantitative real-time polymerase chain reaction (qRT-PCR). Our results indicated that the expression levels of AvBD5, IL16, IL22 and IL18R1 were obviously up-regulated, and Il-18 expression was also changed, but not significantly, which play key roles in the defense against NDV. Overall, we identified several candidate targets that may be involved in the regulation of NDV infection, which provide new insights into the complicated regulatory mechanisms of virus-host interactions, and explore new strategies for protecting chickens against the virus.

The major histocompatibility complex of Old World camelids: Class I and class I-related genes

The genomic structure of the Major Histocompatibility Complex (MHC) region and variation in selected MHC class I related genes in Old World camels, Camelus bactrianus and Camelus dromedaries were studied. The overall genomic organization of the camel MHC region follows a general pattern observed in other mammalian species and individual MHC loci appear to be well conserved. Selected MHC class I genes B-67 and BL3-7 exhibited unexpectedly low variability, even when compared to other camel MHC class I related genes MR1 and MICA. Interspecific SNP and allele sharing are relatively common, and frequencies of heterozygotes are usually low. Such a low variation in a genomic region generally considered as one of the most polymorphic in vertebrate genomes is unusual. Evolutionary relationships between MHC class I related genes and their counterparts from other species seem to be rather complex. Often, they do not follow the general evolutionary history of the species concerned. Close evolutionary relationships of individual MHC class I loci between camels, humans and dogs were observed. Based on the results of this study and on our data on MHC class II genes, the extent and the pattern of polymorphism of the MHC region of Old World camelids differed from most mammalian groups studied so far. Camels thus seem to be an important model for our understanding of the role of genetic diversity in immune functions, especially in the context of unique features of their immunoglobulin and T-cell receptor genes.

Comprehensive meta-analysis reveals an association of the HLA-DRB1*1602 allele with autoimmune diseases mediated predominantly by autoantibodies

The human leukocytes antigen (HLA)-DRB1*16:02 allele has been suggested to be associated with many autoimmune diseases. However, a validation of the results of the different studies by a comprehensive analysis of the corresponding meta data is lacking. In this study, we performed a meta-analysis of the association between HLA-DRB1*16:02 allele with various autoimmune disorders. Our analysis shows that HLA-DRB1*16:02 allele was associated with systemic lupus erythematosus, anti-N-Methyl-d-Aspartate receptor (NMDAR) encephalitis, Graves' disease, myasthenia gravis, neuromyelitis optica and antibody-associated systemic vasculitis with microscopic polyangiitis (AASV-MPA). However, no such association was found for multiple sclerosis, autoimmune hepatitis type 1, rheumatoid arthritis, type 1 diabetes and Vogt-Koyanagi-Harada syndrome. Re-analysis of the studies after their categorization into autoantibody-dependent and T cell-dependent autoimmune diseases revealed that the HLA-DRB1*16:02 allele was strongly associated with disorder predominantly mediated by autoantibodies (OR = 1.93; 95% CI = 1.63-2.28, P = 1.95 × 10^-14) but not with those predominantly mediated by T cells (OR = 1.08; 95% CI = 0.87-1.34, P = .474). In addition, amino acid sequence alignment of common HLA-DRB1 subtypes demonstrated that HLA-DRB1*16:02 carries a unique motif of amino acid residues at position 67-74 which encodes the third hypervariable region. Taken together, the distinct pattern of disease association and the unique amino acid sequence of the third hypervariable region of the HLA-DRB1 provide some hints on how HLA-DRB1*16:02 is involved in the pathogenesis of autoimmune diseases.

Deiminated proteins and extracellular vesicles as novel biomarkers in pinnipeds: Grey seal (Halichoerus gryptus) and harbour seal (Phoca vitulina)

Peptidylarginine deiminases (PADs) are phylogenetically conserved calcium-dependent enzymes which post-translationally convert arginine into citrulline in target proteins in an irreversible manner, leading to functional and structural changes in target proteins. Protein deimination can cause the generation of neo-epitopes, affect gene regulation and also allow for protein moonlighting and therefore facilitate multifaceted functions of the same protein. PADs are furthermore a key regulator of cellular release of extracellular vesicle (EVs), which are found in most body fluids and participate in cellular communication via transfer of cargo proteins and genetic material. In this study, post-translationally deiminated proteins and EVs were assessed in sera of two seal species, grey seal and harbour seal. We report a poly-dispersed population of serum-EVs, which were positive for phylogenetically conserved EV-specific markers and characterised by transmission electron microscopy. A number of deiminated proteins critical for immune and metabolic functions were identified in the seal sera and varied somewhat between the two species under study, while some targets were in common. EV profiles of the seal sera further revealed that key microRNAs for inflammation, immunity and hypoxia also vary between the two species. Protein deimination and EVs profiles may be useful biomarkers for assessing health status of sea mammals, which face environmental challenges, including opportunistic infection, pollution and shifting habitat due to global warming.

SLA-1 Genetic Diversity in Pigs: Extensive Analysis of Copy Number Variation, Heterozygosity, Expression, and Breed Specificity

Swine leukocyte antigens play indispensable roles in immune responses by recognizing a large number of foreign antigens and thus, their genetic diversity plays a critical role in their functions. In this study, we developed a new high-resolution typing method for pig SLA-1 and successfully typed 307 individuals from diverse genetic backgrounds including 11 pure breeds, 1 cross bred, and 12 cell lines. We identified a total of 52 alleles including 18 novel alleles and 9 SLA-1 duplication haplotypes, including 4 new haplotypes. We observed significant differences in the distribution of SLA-1 alleles among the different pig breeds, including the breed specific alleles. SLA-1 duplication was observed in 33% of the chromosomes and was especially high in the biomedical model breeds such as SNU (100%) and NIH (76%) miniature pigs. Our analysis showed that SLA-1 duplication is associated with the increased level of SLA-1 mRNA expression in porcine cells compared to that of the single copy haplotype. Therefore, we provide here the results of the most extensive genetic analysis on pig SLA-1.

The Major Histocompatibility Complex of Old World Camels-A Synopsis

This study brings new information on major histocompatibility complex (MHC) class III sub-region genes in Old World camels and integrates current knowledge of the MHC region into a comprehensive overview for Old World camels. Out of the MHC class III genes characterized, TNFA and the LY6 gene family showed high levels of conservation, characteristic for MHC class III loci in general. For comparison, an MHC class II gene TAP1, not coding for antigen presenting molecules but functionally related to MHC antigen presenting functions was studied. TAP1 had many SNPs, even higher than the MHC class I and II genes encoding antigen presenting molecules. Based on this knowledge and using new camel genomic resources, we constructed an improved genomic map of the entire MHC region of Old World camels. The MHC class III sub-region shows a standard organization similar to that of pig or cattle. The overall genomic structure of the camel MHC is more similar to pig MHC than to cattle MHC. This conclusion is supported by differences in the organization of the MHC class II sub-region, absence of functional DY genes, different organization of MIC genes in the MHC class I sub-region, and generally closer evolutionary relationships of camel and porcine MHC gene sequences analyzed so far.

Origin and evolution of the major histocompatibility complex class I region in eutherian mammals

Major histocompatibility complex (MHC) genes in vertebrates are vital in defending against pathogenic infections. To gain new insights into the evolution of MHC Class I (MHCI) genes and test competing hypotheses on the origin of the MHCI region in eutherian mammals, we studied available genome assemblies of nine species in Afrotheria, Xenarthra, and Laurasiatheria, and successfully characterized the MHCI region in six species. The following numbers of putatively functional genes were detected: in the elephant, four, one, and eight in the extended class I region, and κ and β duplication blocks, respectively; in the tenrec, one in the κ duplication block; and in the four bat species, one or two in the β duplication block. Our results indicate that MHCI genes in the κ and β duplication blocks may have originated in the common ancestor of eutherian mammals. In the elephant, tenrec, and all four bats, some MHCI genes occurred outside the MHCI region, suggesting that eutherians may have a more complex MHCI genomic organization than previously thought. Bat-specific three- or five-amino-acid insertions were detected in the MHCI α1 domain in all four bats studied, suggesting that pathogen defense in bats relies on MHCIs having a wider peptide-binding groove, as previously assayed by a bat MHCI gene with a three-amino-acid insertion showing a larger peptide repertoire than in other mammals. Our study adds to knowledge on the diversity of eutherian MHCI genes, which may have been shaped in a taxon-specific manner.

Long Noncoding RNA HCP5, a Hybrid HLA Class I Endogenous Retroviral Gene: Structure, Expression, and Disease Associations

The HCP5 RNA gene (NCBI ID: 10866) is located centromeric of the HLA-B gene and between the MICA and MICB genes within the major histocompatibility complex (MHC) class I region. It is a human species-specific gene that codes for a long noncoding RNA (lncRNA), composed mostly of an ancient ancestral endogenous antisense 3′ long terminal repeat (LTR, and part of the internal pol antisense sequence of endogenous retrovirus (ERV) type 16 linked to a human leukocyte antigen (HLA) class I promoter and leader sequence at the 5′-end. Since its discovery in 1993, many disease association and gene expression studies have shown that HCP5 is a regulatory lncRNA involved in adaptive and innate immune responses and associated with the promotion of some autoimmune diseases and cancers. The gene sequence acts as a genomic anchor point for binding transcription factors, enhancers, and chromatin remodeling enzymes in the regulation of transcription and chromatin folding. The HCP5 antisense retroviral transcript also interacts with regulatory microRNA and immune and cellular checkpoints in cancers suggesting its potential as a drug target for novel antitumor therapeutics.

MHC class I polymorphic Alu insertion (POALIN) allele and haplotype frequencies in the Arabs of the United Arab Emirates and other world populations

Polymorphic Alu insertions (POALINs) are found throughout the human genome and have been used in various studies to infer geographic origin of human populations. The main aim of this study was to determine the allele and haplotype frequencies of five POALINs, AluHF, AluHG, AluHJ, AluTF and AluMICB, within the major histocompatibility complex (MHC) class I region of 95 UAE Arabs, and correlate their frequencies to those of the HLA-A, HLA-C and HLA-B class I allele lineages. Evolutionary relationships between the POALINs of the Arabs and those previously studied in populations of African, Asian and European descent were compared. At each of the five Alu loci (AluHF, AluHG, AluHJ, AluTF and AluMICB), Alu insertion was designated as Alu(locus)*02 and absence was Alu(locus)*01. The AluHG insertion (AluHG*02) had the highest frequency (0.332), followed by AluHF*02 (0.300), AluHJ*02 (0.263), AluMICB*02 (0.111) and AluTF*02 (0.058). Of the 270 Alu-HLA haplotypes pairs in the UAE Arabs, 110 had no Alu insertion, and 54 had an Alu insertion at >50% per haplotype. An Alu insertion >75% per haplotype was found between AluMICB*02 and HLA-B*14, HLA-B*22, HLA-B*44, HLA-B*55, HLA-B*57 and HLA-B*73, and with HLA-C*01 and HLA-C*18; AluHJ*02 with HLA-A*01, HLA-A*19, HLA-A*24 and HLA-A*32; AluHG*02 with HLA-A*02 and HLA-B*18; and AluHF*02 with HLA-A*10. The genotyped allele and haplotype frequencies of the MHC POALINs in UAE Arabs were compared with the results of 30 previously published Asian, European, American and African populations. Phylogenetic and multidimensional scaling (MDS) analysis of the relative MHC POALINs allele and haplotype frequencies revealed that the UAE Arabs have a similar lineage to Caucasians and the most distant genetic relationship to the Waorani native American population of Ecuador. The structure of both the phylogenetic tree and the MDS analysis supports the Out of Africa theory of human evolution. The nature of the clusters suggests the Arabian Middle East represents a crossroads from which human populations migrated towards Asia in the east and Europe to the north-west.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Novel nucleotide variants in SLA-DOB and CD4 are associated with immune traits in pregnant sows

Reinforcing the immunity of pregnant sows can not only improve their own health condition but also increase the survival rate and healthy status of their piglets. This study aims to find single-nucleotide polymorphism (SNP) and molecular markers that are associated with the immune traits of pregnant sows. SLA-DOB and CD4 were selected as candidate genes, and blood samples were randomly collected from pregnant Landrace sows and used to detect T-lymphocyte subsets, interferon alpha, interleukin 6, Toll-like receptor 3, serum antibody immunoglobulin G, and porcine reproductive and respiratory syndrome virus-specific antibody. Then, association analyses were conducted for the polymorphic sites of candidate genes with immune traits. We found 12 mutations in the two genes and conducted an association study with eight of them. Our results indicated that among the eight mutations, SNP1, SNP2, and SNP3 of the SLA-DOB gene and Ins9, SNP10, and SNP11 in the CD4 gene are newly discovered mutations. Except for SNP1, SNP3, and SNP11, the other five SNPs are associated with at least one immune trait tested. Especially, SNP2 and Ins9 are significantly associated with at least one of the T-lymphocyte subgroups and at least one antibody. These novel mutations have potential important effects on the polymorphic loci of the above immune traits in pregnant sows. The results suggest that the SLA-DOB and CD4 genes and their genetic mutations can be considered as important candidate genes and mutations for the immunity of pregnant sows.

Specific MHC class I supertype associated with parasite infection and color morph in a wild lizard population

The major histocompatibility complex (MHC) is a large gene family that plays a central role in the immune system of all jawed vertebrates. Nonavian reptiles are underrepresented within the MHC literature and little is understood regarding the mechanisms maintaining MHC diversity in this vertebrate group. Here, we examined the relative roles of parasite-mediated selection and sexual selection in maintaining MHC class I diversity of a color polymorphic lizard. We discovered evidence for parasite-mediated selection acting via rare-allele advantage or fluctuating selection as ectoparasite load was significantly lower in the presence of a specific MHC supertype (functional clustering of alleles): supertype four. Based on comparisons between ectoparasite prevalence and load, and assessment of the impact of ectoparasite load on host fitness, we suggest that supertype four confers quantitative resistance to ticks or an intracellular tickborne parasite. We found no evidence for MHC-associated mating in terms of pair genetic distance, number of alleles, or specific supertypes. An association was uncovered between supertype four and male throat color morph. However, it is unlikely that male throat coloration acts as a signal of MHC genotype to conspecifics because we found no evidence to suggest that male throat coloration predicts male mating status. Overall, our results suggest that parasite-mediated selection plays a role in maintaining MHC diversity in this population via rare-allele advantage and/or fluctuating selection. Further work is required to determine whether sexual selection also plays a role in maintaining MHC diversity in agamid lizards.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title "Evolution of the immune system influences speciation rates in teleost fish", is unsubstantiated, and that their "pinpointing" of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

Positive selection in coding regions and motif duplication in regulatory regions of bottlenose dolphin MHC class II genes

The vertebrate immune response is mediated through highly adaptive, quickly evolving cell surface receptors, the major histocompatibility complex (MHC). MHC molecules bind and present a diverse array of pathogenic molecules and trigger a cascade of defenses. Use of MHC variation as a marker for population health has also evolved quickly following advances in sequencing methods. We applied a combination of traditional and next generation sequencing methodology to characterize coding (peptide binding region) and regulatory (proximal promoter) sequence variation in MHC Class II DQA and DQB genes between estuarine and coastal populations of the bottlenose dolphin, Tursiops truncatus, an apex predator whose health status is indicative of anthropogenic impacts on the ecosystem. The coding regions had 10 alleles each at DQA and DQB; the promoters had 6 and 7 alleles at DQA and DQB, respectively with variation within key regulatory motifs. Positive selection was observed for the coding regions of both genes while both coding and promoter regions exhibited geographic differences in allele composition that likely indicates diversifying selection across habitats. Most notable was the discovery of a complete duplication of a 14-bp T-box motif in the DQA promoter. Four class II promoter regions (DQA, DQB, DRA, DRB) were characterized in species from four cetacean families (Delphinidae, Monodontidae, Lipotidae, and Physeteridae) and revealed substantial promoter structural diversity across this order. Peptide binding regions may not be the only source of adaptive potential within cetacean MHC for responding to pathogenic threats. These findings are the first analysis of cetacean MHC regulatory motifs, which may divulge unique immunogenetic strategies among cetaceans and reveal how MHC transcriptional control continues to evolve. The combined MHC regulatory and coding data provide new genetic context for distinct vulnerability profiles between coastal and estuarine populations, which are key concerns for health and risk management.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title “Evolution of the immune system influences speciation rates in teleost fish”, is unsubstantiated, and that their “pinpointing” of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

The Egyptian Rousette Genome Reveals Unexpected Features of Bat Antiviral Immunity

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.

Swinepox virus vector-based vaccines: attenuation and biosafety assessments following subcutaneous prick inoculation

Swinepox virus (SPV) has several advantages as a potential clinical vector for a live vector vaccine. In this study, to obtain a safer and more efficient SPV vector, three SPV mutants, Δ003, Δ010, and ΔTK were successfully constructed. A virus replication experiment showed that these SPV mutants had lower replication abilities compared to wtSPV in 10 different host-derived cell lines. Animal experiments with mouse and rabbit models demonstrate that these three mutants and wtSPV did not cause any clinical signs of dermatitis. No fatalities were observed during a peritoneal challenge assay with these mutants and wtSPV in a mouse model. Additionally, the three mutants and wtSPV were not infectious at 60 h after vaccination in rabbit models. Furthermore, we evaluated biosafety, immunogenicity and effectiveness of the three mutants in 65 1-month-old piglets. The results show that there were no clinical signs of dermatitis in the Δ003 and ΔTK vaccination groups. However, mild signs were observed in the Δ010 vaccination groups when virus titres were high, and apparent clinical signs were observed at the sites of inoculation. Samples from all experimental pig groups were assessed by qPCR, and no SPV genomic DNA was found in five organs, faeces or blood. This suggests that the infectious abilities of wtSPV and the SPV mutants were poor and limited. In summary, this study indicates that two mutants of SPV, Δ003 and ΔTK, may be promising candidates for an attenuated viral vector in veterinary medicine.

AmpliSAS and AmpliHLA: Web Server Tools for MHC Typing of Non-Model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are web server tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de-novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human-specific version, it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Here we describe four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Inferring the evolution of the major histocompatibility complex of wild pigs and peccaries using hybridisation DNA capture-based sequencing

The major histocompatibility complex (MHC) is a key genomic model region for understanding the evolution of gene families and the co-evolution between host and pathogen. To date, MHC studies have mostly focused on species from major vertebrate lineages. The evolution of MHC classical (Ia) and non-classical (Ib) genes in pigs has attracted interest because of their antigen presentation roles as part of the adaptive immune system. The pig family Suidae comprises over 18 extant species (mostly wild), but only the domestic pig has been extensively sequenced and annotated. To address this, we used a DNA-capture approach, with probes designed from the domestic pig genome, to generate MHC data for 11 wild species of pigs and their closest living family, Tayassuidae. The approach showed good efficiency for wild pigs (~80% reads mapped, ~87× coverage), compared to tayassuids (~12% reads mapped, ~4× coverage). We retrieved 145 MHC loci across both families. Phylogenetic analyses show that the class Ia and Ib genes underwent multiple duplications and diversifications before suids and tayassuids diverged from their common ancestor. The histocompatibility genes mostly form orthologous groups and there is genetic differentiation for most of these genes between Eurasian and sub-Saharan African wild pigs. Tests of selection showed that the peptide-binding region of class Ib genes was under positive selection. These findings contribute to better understanding of the evolutionary history of the MHC, specifically, the class I genes, and provide useful data for investigating the immune response of wild populations against pathogens.

Discovery of gorilla MHC-C expressing C1 ligand for KIR

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.

Review of Current Conservation Genetic Analyses of Northeast Pacific Sharks

Conservation genetics is an applied science that utilizes molecular tools to help solve problems in species conservation and management. It is an interdisciplinary specialty in which scientists apply the study of genetics in conjunction with traditional ecological fieldwork and other techniques to explore molecular variation, population boundaries, and evolutionary relationships with the goal of enabling resource managers to better protect biodiversity and identify unique populations. Several shark species in the northeast Pacific (NEP) have been studied using conservation genetics techniques, which are discussed here. The primary methods employed to study population genetics of sharks have historically been nuclear microsatellites and mitochondrial (mt) DNA. These markers have been used to assess genetic diversity, mating systems, parentage, relatedness, and genetically distinct populations to inform management decisions. Novel approaches in conservation genetics, including next-generation DNA and RNA sequencing, environmental DNA (eDNA), and epigenetics are just beginning to be applied to elasmobranch evolution, physiology, and ecology. Here, we review the methods and results of past studies, explore future directions for shark conservation genetics, and discuss the implications of molecular research and techniques for the long-term management of shark populations in the NEP.

Limited MHC class I intron 2 repertoire variation in bonobos

Common chimpanzees (Pan troglodytes) experienced a selective sweep, probably caused by a SIV-like virus, which targeted their MHC class I repertoire. Based on MHC class I intron 2 data analyses, this selective sweep took place about 2-3 million years ago. As a consequence, common chimpanzees have a skewed MHC class I repertoire that is enriched for allotypes that are able to recognise conserved regions of the SIV proteome. The bonobo (Pan paniscus) shared an ancestor with common chimpanzees approximately 1.5 to 2 million years ago. To investigate whether the signature of this selective sweep is also detectable in bonobos, the MHC class I gene repertoire of two bonobo panels comprising in total 29 animals was investigated by Sanger sequencing. We identified 14 Papa-A, 20 Papa-B and 11 Papa-C alleles, of which eight, five and eight alleles, respectively, have not been reported previously. Within this pool of MHC class I variation, we recovered only 2 Papa-A, 3 Papa-B and 6 Papa-C intron 2 sequences. As compared to humans, bonobos appear to have an even more diminished MHC class I intron 2 lineage repertoire than common chimpanzees. This supports the notion that the selective sweep may have predated the speciation of common chimpanzees and bonobos. The further reduction of the MHC class I intron 2 lineage repertoire observed in bonobos as compared to the common chimpanzee may be explained by a founding effect or other subsequent selective processes.

MHC class I diversity in chimpanzees and bonobos

Major histocompatibility complex (MHC) class I genes are critically involved in the defense against intracellular pathogens. MHC diversity comparisons among samples of closely related taxa may reveal traces of past or ongoing selective processes. The bonobo and chimpanzee are the closest living evolutionary relatives of humans and last shared a common ancestor some 1 mya. However, little is known concerning MHC class I diversity in bonobos or in central chimpanzees, the most numerous and genetically diverse chimpanzee subspecies. Here, we used a long-read sequencing technology (PacBio) to sequence the classical MHC class I genes A, B, C, and A-like in 20 and 30 wild-born bonobos and chimpanzees, respectively, with a main focus on central chimpanzees to assess and compare diversity in those two species. We describe in total 21 and 42 novel coding region sequences for the two species, respectively. In addition, we found evidence for a reduced MHC class I diversity in bonobos as compared to central chimpanzees as well as to western chimpanzees and humans. The reduced bonobo MHC class I diversity may be the result of a selective process in their evolutionary past since their split from chimpanzees.

NOG-hIL-4-Tg, a new humanized mouse model for producing tumor antigen-specific IgG antibody by peptide vaccination

Immunodeficient mice transplanted with human peripheral blood mononuclear cells (PBMCs) are promising tools to evaluate human immune responses to vaccines. However, these mice usually develop severe graft-versus-host disease (GVHD), which makes estimation of antigen-specific IgG production after antigen immunization difficult. To evaluate antigen-specific IgG responses in PBMC-transplanted immunodeficient mice, we developed a novel NOD/Shi-scid-IL2rγ^null (NOG) mouse strain that systemically expresses the human IL-4 gene (NOG-hIL-4-Tg). After human PBMC transplantation, GVHD symptoms were significantly suppressed in NOG-hIL-4-Tg compared to conventional NOG mice. In kinetic analyses of human leukocytes, long-term engraftment of human T cells has been observed in peripheral blood of NOG-hIL-4-Tg, followed by dominant CD4+ T rather than CD8+ T cell proliferation. Furthermore, these CD4+ T cells shifted to type 2 helper (Th2) cells, resulting in long-term suppression of GVHD. Most of the human B cells detected in the transplanted mice had a plasmablast phenotype. Vaccination with HER2 multiple antigen peptide (CH401MAP) or keyhole limpet hemocyanin (KLH) successfully induced antigen-specific IgG production in PBMC-transplanted NOG-hIL-4-Tg. The HLA haplotype of donor PBMCs might not be relevant to the antibody secretion ability after immunization. These results suggest that the human PBMC-transplanted NOG-hIL-4-Tg mouse is an effective tool to evaluate the production of antigen-specific IgG antibodies.