Comparative genomic analysis of mammalian NKG2D ligand family genes provides insights into their origin and evolution

Regulation of NKG2D Stress Ligands and Its Relevance in Cancer Progression

Under cellular distress, multiple facets of normal homeostatic signaling are altered or disrupted. In the context of the immune landscape, external and internal stressors normally promote the expression of natural killer group 2 member D (NKG2D) ligands that allow for the targeted recognition and killing of cells by NKG2D receptor-bearing effector populations. The presence or absence of NKG2D ligands can heavily influence disease progression and impact the accessibility of immunotherapy options. In cancer, tumor cells are known to have distinct regulatory mechanisms for NKG2D ligands that are directly associated with tumor progression and maintenance. Therefore, understanding the regulation of NKG2D ligands in cancer will allow for targeted therapeutic endeavors aimed at exploiting the stress response pathway. In this review, we summarize the current understanding of regulatory mechanisms controlling the induction and repression of NKG2D ligands in cancer. Additionally, we highlight current therapeutic endeavors targeting NKG2D ligand expression and offer our perspective on considerations to further enhance the field of NKG2D ligand biology.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Mechanisms of Natural Killer Cell Evasion Through Viral Adaptation

The continuous interactions between host and pathogens during their coevolution have shaped both the immune system and the countermeasures used by pathogens. Natural killer (NK) cells are innate lymphocytes that are considered central players in the antiviral response. Not only do they express a variety of inhibitory and activating receptors to discriminate and eliminate target cells but they can also produce immunoregulatory cytokines to alert the immune system. Reciprocally, several unrelated viruses including cytomegalovirus, human immunodeficiency virus, influenza virus, and dengue virus have evolved a multitude of mechanisms to evade NK cell function, such as the targeting of pathways for NK cell receptors and their ligands, apoptosis, and cytokine-mediated signaling. The studies discussed in this article provide further insights into the antiviral function of NK cells and the pathways involved, their constituent proteins, and ways in which they could be manipulated for host benefit.

Deletion of a Seminal Gene Cluster Reinforces a Crucial Role of SVS2 in Male Fertility

Multiple genes, whose functions or expression are overlapping, compensate for the loss of one gene. A gene cluster in the mouse genome encodes five seminal vesicle proteins (SVS2, SVS3, SVS4, SVS5, and SVS6). These proteins are produced by male rodents and function in formation of the copulatory plug following mating. SVS2 plays an essential role in the successful internal fertilization by protecting the sperm membrane against a uterine immune attack. We hypothesized that the four remaining seminal vesicle proteins (SVPs) of this gene cluster may partially/completely compensate for the deficiency of SVS2. For confirming our hypothesis, we generated mice lacking the entire SVP-encoding gene cluster and compared their fecundity with Svs2-deficient (Svs2^−/−) mice; that is, mice deficient in Svs2 alone. A single loxP site remained after the deletion of the Svs2 gene. Therefore, we inserted another loxP site by combining the CRISPR/Cas9 system with single-stranded oligodeoxynucleotides (ssODN). Male mice lacking the entire SVP-encoding gene cluster (Svs2–6^−/− mice) and thereby all five SVP proteins, generated by the deletion of 100kbp genomic DNA, showed low fecundity. However, the fecundity level was comparable with that from Svs2^−/− male mice. Our results demonstrate that SVS3, SVS4, SVS5, and SVS6 do not function in the protection of sperm against a uterine immune attack in the absence of SVS2. Thus, Svs2 is the critical gene in the SVP gene cluster.

Structure of MHC class I-like MILL2 reveals heparan-sulfate binding and interdomain flexibility

The MILL family, composed of MILL1 and MILL2, is a group of nonclassical MHC class I molecules that occur in some orders of mammals. It has been reported that mouse MILL2 is involved in wound healing; however, the molecular mechanisms remain unknown. Here, we determine the crystal structure of MILL2 at 2.15 Å resolution, revealing an organization similar to classical MHC class I. However, the α1-α2 domains are not tightly fixed on the α3-β₂m domains, indicating unusual interdomain flexibility. The groove between the two helices in the α1-α2 domains is too narrow to permit ligand binding. Notably, an unusual basic patch on the α3 domain is involved in the binding to heparan sulfate which is essential for MILL2 interactions with fibroblasts. These findings suggest that MILL2 has a unique structural architecture and physiological role, with binding to heparan sulfate proteoglycans on fibroblasts possibly regulating cellular recruitment in biological events.

The MILL (MHC-I-like located near the leukocyte receptor complex) family is a group of related nonclassical MHC-I molecules. Here the authors present the crystal structure of MILL2, which reveals an unusual interdomain flexibility, and show that MILL2 binds heparan sulfate on the surface of fibroblasts through a basic patch.

Diversity of ULBP5 in Old-World monkeys (Cercopithecidae) and divergence of the ULBP gene family in primates

Non-human primates such as rhesus macaque and cynomolgus macaque are important animals for medical research. These species are classified as Old-World monkeys (Cercopithecidae), in which the immune-related genome structure is characterized by gene duplications. In the present study, we investigated polymorphisms in two genes for ULBP5 encoding ligands for NKG2D. We found 18 and 11 ULBP5.1 alleles and 11 and 13 ULBP5.2 alleles in rhesus macaques and cynomolgus macaques, respectively. In addition, phylogenetic analyses revealed that ULBP5.2 diverged from a branch of ULBP5.1. These data suggested that human ULBP genes diverged from an ancestral gene of ULBP2-ULBP5 and that ULBP6/RAET1L, specifically identified in human, diverged from an ancestral ULBP2 by a recent gene duplication after the diversification of homininae (human and other higher great apes), which were consistent with the findings in our previous analysis of ULBP2 genes in rhesus and cynomolgus macaques.

Anti-NKG2D mAb: A New Treatment for Crohn's Disease?

Crohn’s disease (CD) and ulcerative colitis (UC) are immunologically-mediated, debilitating conditions resulting from destructive inflammation of the gastrointestinal tract. The pathogenesis of IBD is incompletely understood, but is considered to be the result of an abnormal immune response with a wide range of cell types and proteins involved. Natural Killer Group 2D (NKG2D) is an activating receptor constitutively expressed on human Natural Killer (NK), γδ T, mucosal-associated invariant T (MAIT), CD56⁺ T, and CD8⁺ T cells. Activation of NKG2D triggers cellular proliferation, cytokine production, and target cell killing. Research into the NKG2D mechanism of action has primarily been focused on cancer and viral infections where cytotoxicity evasion is a concern. In human inflammatory bowel disease (IBD) this system is less characterized, but the ligands have been shown to be highly expressed during intestinal inflammation and the following receptor activation may contribute to tissue degeneration. A recent phase II clinical trial showed that an antibody against NKG2D induced clinical remission of CD in some patients, suggesting NKG2D and its ligands to be of importance in the pathogenesis of CD. This review will describe the receptor and its ligands in intestinal tissues and the clinical potential of blocking NKG2D in Crohn’s disease.

Comparative genomics of the NKG2D ligand gene family

NKG2D ligands (NKG2DLs) are a group of stress-inducible major histocompatibility complex (MHC) class I-like molecules that act as a danger signal alerting the immune system to the presence of abnormal cells. In mammals, two families of NKG2DL genes have been identified: the MIC gene family encoded in the MHC region and the ULBP gene family encoded outside the MHC region in most species. Some mammals have a third family of NKG2DL-like class I genes which we named MILL (MHC class I-like located near the leukocyte receptor complex). Despite the fact that MILL genes are more closely related to MIC genes than ULBP genes are to MIC genes, MILL molecules do not function as NKG2DLs, and their function remains unknown. With the progress of mammalian genome projects, information on the MIC, ULBP, and MILL gene families became available in many mammalian species. Here, we summarize such information and discuss the origin and evolution of the NKG2DL gene family from the viewpoint of host-pathogen coevolution.

Genetics, genomics, and evolutionary biology of NKG2D ligands

Human and mouse NKG2D ligands (NKG2DLs) are absent or only poorly expressed by most normal cells but are upregulated by cell stress, hence, alerting the immune system in case of malignancy or infection. Although these ligands are numerous and highly variable (at genetic, genomic, structural, and biochemical levels), they all belong to the major histocompatibility complex class I gene superfamily and bind to a single, invariant, receptor: NKG2D. NKG2D (CD314) is an activating receptor expressed on NK cells and subsets of T cells that have a key role in the recognition and lysis of infected and tumor cells. Here, we review the molecular diversity of NKG2DLs, discuss the increasing appreciation of their roles in a variety of medical conditions, and propose several explanations for the evolutionary force(s) that seem to drive the multiplicity and diversity of NKG2DLs while maintaining their interaction with a single invariant receptor.

Marsupials and monotremes possess a novel family of MHC class I genes that is lost from the eutherian lineage

Background

Major histocompatibility complex (MHC) class I genes are found in the genomes of all jawed vertebrates. The evolution of this gene family is closely tied to the evolution of the vertebrate genome. Family members are frequently found in four paralogous regions, which were formed in two rounds of genome duplication in the early vertebrates, but in some species class Is have been subject to additional duplication or translocation, creating additional clusters. The gene family is traditionally grouped into two subtypes: classical MHC class I genes that are usually MHC-linked, highly polymorphic, expressed in a broad range of tissues and present endogenously-derived peptides to cytotoxic T-cells; and non-classical MHC class I genes generally have lower polymorphism, may have tissue-specific expression and have evolved to perform immune-related or non-immune functions. As immune genes can evolve rapidly and are subject to different selection pressure, we hypothesised that there may be divergent, as yet unannotated or uncharacterised class I genes.

Results

Application of a novel method of sensitive genome searching of available vertebrate genome sequences revealed a new, extensive sub-family of divergent MHC class I genes, denoted as UT, which has not previously been characterized. These class I genes are found in both American and Australian marsupials, and in monotremes, at an evolutionary chromosomal breakpoint, but are not present in non-mammalian genomes and have been lost from the eutherian lineage. We show that UT family members are expressed in the thymus of the gray short-tailed opossum and in other immune tissues of several Australian marsupials. Structural homology modelling shows that the proteins encoded by this family are predicted to have an open, though short, antigen-binding groove.

Conclusions

We have identified a novel sub-family of putatively non-classical MHC class I genes that are specific to marsupials and monotremes. This family was present in the ancestral mammal and is found in extant marsupials and monotremes, but has been lost from the eutherian lineage. The function of this family is as yet unknown, however, their predicted structure may be consistent with presentation of antigens to T-cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1745-4) contains supplementary material, which is available to authorized users.

Divergence and diversity of ULBP2 genes in rhesus and cynomolgus macaques

Non-human primates such as rhesus macaque and cynomolgus macaque are important animals for medical research fields and they are classified as Old World monkey, in which genome structure is characterized by gene duplications. In the present study, we investigated polymorphisms in two genes for ULBP2 molecules that are ligands for NKG2D. A total of 15 and 11 ULBP2.1 alleles and 11 and 10 ULBP2.2 alleles were identified in rhesus macaques and cynomolgus macaques, respectively. Nucleotide sequences of exons for extra cellular domain were highly polymorphic and more than 70 % were non-synonymous variations in both ULBP2.1 and ULBP2.2. In addition, phylogenetic analyses revealed that the ULBP2.2 was diverged from a branch of ULBP2.1 along with ULBP2s of higher primates. Moreover, when 3D structural models were constructed for the rhesus ULBP2 molecules, residues at presumed contact sites with NKG2D were polymorphic in ULBP2.1 and ULBP2.2 in the rhesus macaque and cynomolgus macaque, respectively. These observations suggest that amino acid replacements at the interaction sites with NKG2D might shape a specific nature of ULBP2 molecules in the Old World monkeys.

Genome of the Chinese tree shrew

Chinese tree shrews (Tupaia belangeri chinensis) possess many features valuable in animals used as experimental models in biomedical research. Currently, there are numerous attempts to employ tree shrews as models for a variety of human disorders: depression, myopia, hepatitis B and C virus infections, and hepatocellular carcinoma, to name a few. Here we present a publicly available annotated genome sequence for the Chinese tree shrew. Phylogenomic analysis of the tree shrew and other mammalians highly support its close affinity to primates. By characterizing key factors and signalling pathways in nervous and immune systems, we demonstrate that tree shrews possess both shared common and unique features, and provide a genetic basis for the use of this animal as a potential model for biomedical research.

Structural basis of mouse cytomegalovirus m152/gp40 interaction with RAE1γ reveals a paradigm for MHC/MHC interaction in immune evasion

Natural killer (NK) cells are activated by engagement of the NKG2D receptor with ligands on target cells stressed by infection or tumorigenesis. Several human and rodent cytomegalovirus (CMV) immunoevasins down-regulate surface expression of NKG2D ligands. The mouse CMV MHC class I (MHC-I)-like m152/gp40 glycoprotein down-regulates retinoic acid early inducible-1 (RAE1) NKG2D ligands as well as host MHC-I. Here we describe the crystal structure of an m152/RAE1γ complex and confirm the intermolecular contacts by mutagenesis. m152 interacts in a pincer-like manner with two sites on the α1 and α2 helices of RAE1 reminiscent of the NKG2D interaction with RAE1. This structure of an MHC-I-like immunoevasin/MHC-I-like ligand complex explains the binding specificity of m152 for RAE1 and allows modeling of the interaction of m152 with classical MHC-I and of related viral immunoevasins.

Exceptionally high conservation of the MHC class I-related gene, MR1, among mammals

The major histocompatibility complex (MHC) class I-related gene, MR1, is a non-classical MHC class IA gene and is encoded outside the MHC region. The MR1 is responsible for activation of mucosal-associated invariant T (MAIT) cells expressing semi-invariant T cell receptors in the presence of bacteria, but its ligand has not been identified. A unique characteristic of MR1 is its high evolutionary conservation of the α1 and α2 domains corresponding to the peptide-binding domains of classical MHC class I molecules, showing about 90 % amino acid identity between human and mouse. To clarify the evolutionary history of MR1 and identify more critically conserved residues for the function of MR1, we searched for the MR1 gene using jawed vertebrate genome databases and isolated the MR1 cDNA sequences of marsupials (opossum and wallaby). A comparative genomic analysis indicated that MR1 is only present in placental and marsupial mammals and that the gene organization around MR1 is well conserved among analyzed jawed vertebrates. Moreover, the α1 and α2 domains, especially in amino acid residues presumably shaping a ligand-binding groove, were also highly conserved between placental and marsupial MR1. These findings suggest that the MR1 gene might have been established at its present location in a common ancestor of placental and marsupial mammals and that the shape of the putative ligand-binding groove in MR1 has been maintained, probably for presenting highly conserved component(s) of microbes to MAIT cells.

Immunogenetics of the NKG2D ligand gene family

NKG2D ligands (NKG2DLs) are a group of major histocompatibility complex (MHC) class I-like molecules, the expression of which is induced by cellular stresses such as infection, tumorigenesis, heat shock, tissue damage, and DNA damage. They act as a molecular danger signal alerting the immune system for infected or neoplastic cells. Mammals have two families of NKG2DL genes: the MHC-encoded MIC gene family and the ULBP gene family encoded outside the MHC region in most mammals. Rodents such as mice and rats lack the MIC family of ligands. Interestingly, some mammals have NKG2DL-like molecules named MILL that are phylogenetically related to MIC, but do not function as NKG2DLs. In this paper, we review our current knowledge of the MIC, ULBP, and MILL gene families in representative mammalian species and discuss the origin and evolution of the NKG2DL gene family.

Involvement of an NKG2D ligand H60c in epidermal dendritic T cell-mediated wound repair

Dendritic epidermal T cells (DETCs) found in mouse skin are NKG2D-positive γδ T cells involved in immune surveillance and wound repair. It is assumed that the interaction of an NKG2D receptor on DETCs and an MHC class I-like NKG2D ligand on keratinocytes activates DETCs, which then secrete cytokines promoting wound repair. However, direct evidence that DETC activation through NKG2D signaling promotes wound repair is not available. In the present study, we generated mAbs for an NKG2D ligand H60c previously suggested to be expressed specifically on skin keratinocytes. Local administration of H60c-specific mAb inhibited activation of DETCs and significantly delayed wound repair. Likewise, administration of NKG2D-specific mAb impaired wound repair to a similar extent. The delay in wound closure resulting from the blockade of the NKG2D pathway was comparable to that observed in γδ T cell-deficient mice. These results indicate that H60c/NKG2D interactions play a critical role in wound repair. Reassessment of binding affinities showed that H60c monomers bind to NKG2D with affinity (K(d) = 26 ± 3.2 nM) comparable to those of other high-affinity NKG2D ligands. H60c is transcribed not only in skin but also in tissues such as tongue and female reproductive tract known to contain epithelium-resident γδ T cells expressing invariant TCRs, suggesting a more general role for H60c in the maintenance of epithelial integrity.