Evidence for high bi-allelic expression of activating Ly49 receptors

Binary outcomes of enhancer activity underlie stable random monoallelic expression

Mitotically stable random monoallelic gene expression (RME) is documented for a small percentage of autosomal genes. We developed an in vivo genetic model to study the role of enhancers in RME using high-resolution single-cell analysis of natural killer (NK) cell receptor gene expression and enhancer deletions in the mouse germline. Enhancers of the RME NK receptor genes were accessible and enriched in H3K27ac on silent and active alleles alike in cells sorted according to allelic expression status, suggesting enhancer activation and gene expression status can be decoupled. In genes with multiple enhancers, enhancer deletion reduced gene expression frequency, in one instance converting the universally expressed gene encoding NKG2D into an RME gene, recapitulating all aspects of natural RME including mitotic stability of both the active and silent states. The results support the binary model of enhancer action, and suggest that RME is a consequence of general properties of gene regulation by enhancers rather than an RME-specific epigenetic program. Therefore, many and perhaps all genes may be subject to some degree of RME. Surprisingly, this was borne out by analysis of several genes that define different major hematopoietic lineages, that were previously thought to be universally expressed within those lineages: the genes encoding NKG2D, CD45, CD8α, and Thy-1. We propose that intrinsically probabilistic gene allele regulation is a general property of enhancer-controlled gene expression, with previously documented RME representing an extreme on a broad continuum.

Influence of Self-MHC Class I Recognition on the Dynamics of NK Cell Responses to Cytomegalovirus Infection

Although interactions between inhibitory Ly49 receptors and their self-MHC class I ligands in C57BL/6 mice are known to limit NK cell proliferation during mouse CMV (MCMV) infection, we created a 36-marker mass cytometry (CyTOF) panel to investigate how these inhibitory receptors impact the NK cell response to MCMV in other phenotypically measurable ways. More than two thirds of licensed NK cells (i.e., those expressing Ly49C, Ly49I, or both) in uninfected mice had already differentiated into NK cells with phenotypes indicative of Ag encounter (KLRG1+Ly6C-) or memory-like status (KLRG1+Ly6C+). These pre-existing KLRG1+Ly6C+ NK cells resembled known Ag-specific memory NK cell populations in being less responsive to IL-18 and IFN-α stimulation in vitro and by selecting for NK cell clones with elevated expression of a Ly49 receptor. During MCMV infection, the significant differences between licensed and unlicensed (Ly49C-Ly49I-) NK cells disappeared within both CMV-specific (Ly49H+) and nonspecific (Ly49H-) responses. This lack of heterogeneity carried into the memory phase, with only a difference in CD16 expression manifesting between licensed and unlicensed MCMV-specific memory NK cell populations. Our results suggest that restricting proliferation is the predominant effect licensing has on the NK cell population during MCMV infection, but the inhibitory Ly49-MHC interactions that take place ahead of infection contribute to their limited expansion by shrinking the pool of licensed NK cells capable of robustly responding to new challenges.

Evidence for Prescribed NK Cell Ly-49 Developmental Pathways in Mice

Previous studies of NK cell inhibitory Ly-49 genes showed their expression is stochastic. However, relatively few studies have examined the mechanisms governing acquisition of inhibitory receptors in conjunction with activating Ly-49 receptors and NK cell development. We hypothesized that the surface expression of activating Ly-49 receptors is nonrandom and is influenced by inhibitory Ly-49 receptors. We analyzed NK cell "clusters" defined by combinatorial expression of activating (Ly-49H and Ly-49D) and inhibitory (Ly-49I and Ly-49G2) receptors in C57BL/6 mice. Using the product rule to evaluate the interdependencies of the Ly-49 receptors, we found evidence for a tightly regulated expression at the immature NK cell stage, with the highest interdependencies between clusters that express at least one activating receptor. Further analysis demonstrated that certain NK clusters predominated at the immature (CD27⁺CD11b^-), transitional (CD27⁺CD11b⁺), and mature (CD27^-CD11b^-) NK cell stages. Using parallel in vitro culture and in vivo transplantation of sorted NK clusters, we discovered nonrandom expression of Ly-49 receptors, suggesting that prescribed pathways of NK cluster differentiation exist. Our data infer that surface expression of Ly-49I is an important step in NK cell maturation. Ki-67 expression and cell counts confirmed that immature NK cells proliferate more than mature NK cells. We found that MHC class I is particularly important for regulation of Ly-49D and Ly-49G2, even though no known MHC class I ligand for these receptors is present in B6 mice. Our data indicate that surface expression of both activating and inhibitory Ly-49 receptors on NK cell clusters occurs in a nonrandom process correlated to their maturation stage.

Two Lineages of KLRA with Contrasting Transcription Patterns Have Been Conserved at a Single Locus during Ruminant Speciation

Cattle possess the most diverse repertoire of NK cell receptor genes among all mammals studied to date. Killer cell receptor genes encoded within the NK complex and killer cell Ig-like receptor genes encoded within the leukocyte receptor complex have both been expanded and diversified. Our previous studies identified two divergent and polymorphic KLRA alleles within the NK complex in the Holstein-Friesian breed of dairy cattle. By examining a much larger cohort and other ruminant species, we demonstrate the emergence and fixation of two KLRA allele lineages (KLRA*01 and -*02) at a single locus during ruminant speciation. Subsequent recombination events between these allele lineages have increased the frequency of KLRA*02 extracellular domains. KLRA*01 and KLRA*02 transcription levels contrasted in response to cytokine stimulation, whereas homozygous animals consistently transcribed higher levels of KLRA, regardless of the allele lineage. KLRA*02 mRNA levels were also generally higher than KLRA*01 Collectively, these data point toward alternative functional roles governed by KLRA genotype and allele lineage. On a background of high genetic diversity of NK cell receptor genes, this KLRA allele fixation points to fundamental and potentially differential function roles.

Blastocyst transfer in mice alters the placental transcriptome and growth

Assisted reproduction technologies (ART) are becoming increasingly common. Therefore, how these procedures influence gene regulation and feto-placental development are important to explore. Here, we assess the effects of blastocyst transfer on mouse placental growth and transcriptome. C57Bl/6 blastocysts were transferred into uteri of B6D2F1 pseudopregnant females and dissected at embryonic day 10.5 for analysis. Compared to non-transferred controls, placentas from transferred conceptuses weighed less even though the embryos were larger on average. This suggested a compensatory increase in placental efficiency. RNA-sequencing of whole male placentas revealed 543 differentially expressed genes (DEGs) after blastocyst transfer: 188 and 355 genes were down-regulated and up-regulated, respectively. DEGs were independently validated in male and female placentas. Bioinformatic analyses revealed that DEGs represented expression in all major placental cell types and included genes that are critical for placenta development and/or function. Furthermore, the direction of transcriptional change in response to blastocyst transfer implied an adaptive response to improve placental function to maintain fetal growth. Our analysis revealed that CpG methylation at regulatory regions of two DEGs was unchanged in female transferred placentas and that DEGs had fewer gene-associated CpG islands (within ~20 kb region) compared to the larger genome. These data suggested that altered methylation at proximal promoter regions might not lead to transcriptional disruption in transferred placentas. Genomic clustering of some DEGs warrants further investigation of long-range, cis-acting epigenetic mechanisms including histone modifications together with DNA methylation. We conclude that embryo transfer, a protocol required for ART, significantly impacts the placental transcriptome and growth.

Activating Receptor Signals Drive Receptor Diversity in Developing Natural Killer Cells

It has recently been appreciated that NK cells exhibit many features reminiscent of adaptive immune cells. Considerable heterogeneity exists with respect to the ligand specificity of individual NK cells and as such, a subset of NK cells can respond, expand, and differentiate into memory-like cells in a ligand-specific manner. MHC I-binding inhibitory receptors, including those belonging to the Ly49 and KIR families, are expressed in a variegated manner, which creates ligand-specific diversity within the NK cell pool. However, how NK cells determine which inhibitory receptors to express on their cell surface during a narrow window of development is largely unknown. In this manuscript, we demonstrate that signals from activating receptors are critical for induction of Ly49 and KIR receptors during NK cell development; activating receptor-derived signals increased the probability of the Ly49 bidirectional Pro1 promoter to transcribe in the forward versus the reverse direction, leading to stable expression of Ly49 receptors in mature NK cells. Our data support a model where the balance of activating and inhibitory receptor signaling in NK cells selects for the induction of appropriate inhibitory receptors during development, which NK cells use to create a diverse pool of ligand-specific NK cells.

Ly49C Impairs NK Cell Memory in Mouse Cytomegalovirus Infection

NK cells possess inhibitory receptors that are responsible for self-MHC class I recognition; beyond their inhibitory function, accumulating evidence indicates that such receptors confer NK cell functional competence through an unclear process termed "licensing." Ly49C is the main self-specific inhibitory Ly49 receptor in H-2(b) C57BL/6 (B6) mice. We used B6 Ly49C-transgenic and B6 β2 microglobulin (β2m)-knockout Ly49C-transgenic mice to investigate the impact of licensing through this inhibitory receptor in precursor and mature NK cells. We found that self-specific inhibitory receptors affected NK cell precursor survival and proliferation at particular developmental stages in an MHC class I-dependent manner. The presence of Ly49C impacted the NK cell repertoire in a β2m-dependent manner, with reduced Ly49A(+), Ly49G2(+), and Ly49D(+) subsets, an increased DNAM-1(+) subset, and higher NKG2D expression. Licensed NK cells displayed a skewed distribution of the maturation stages, which was characterized by differential CD27 and CD11b expression, toward the mature phenotypes. We found that Ly49C-mediated licensing induced a split effect on NK cell functions, with increased cytokine-production capabilities following engagement of various activating receptors while cytotoxicity remained unchanged. Analysis of licensed NK cell functions in vivo, in a system of mouse CMV infection, indicated that licensing did not play a major role in the NK cell antiviral response during acute infection, but it strongly impaired the generation and/or persistence of memory NK cells. This study unravels multifaceted effects of licensing on NK cell populations and their functions.

Selection, tuning, and adaptation in mouse NK cell education

Natural killer (NK) cells recognize transformed cells with an array of germline-encoded inhibitory and activating receptors. Inhibitory Ly49 receptors bind major histocompatibility complex class I (MHC-I) molecules, providing a mechanism by which NK cells maintain self-tolerance yet eliminate cells expressing reduced levels of MHC-I. Additionally, MHC-I molecules are required for NK cell education, a process in which NK cells acquire responsiveness. In this review, we discuss three facets of MHC class I-dependent education of mouse NK cells: skewing of the inhibitory receptor repertoire, induction of functional responsiveness, and tuning in response to changes in MHC-I expression. We discuss prevailing models for education such as licensing and disarming and propose a model for positive selection of 'useful' NK cell subsets. Furthermore, we argue that both repertoire skewing and functional NK cell education may be altered in mature NK cells subject to changes in MHC-I input and suggest that this process may provide increased dynamics to the NK cell system.

Nucleosome Presence at AML-1 Binding Sites Inversely Correlates with Ly49 Expression: Revelations from an Informatics Analysis of Nucleosomes and Immune Cell Transcription Factors

Beyond its role in genomic organization and compaction, the nucleosome is believed to participate in the regulation of gene transcription. Here, we report a computational method to evaluate the nucleosome sensitivity for a transcription factor over a given stretch of the genome. Sensitive factors are predicted to be those with binding sites preferentially contained within nucleosome boundaries and lacking 10 bp periodicity. Based on these criteria, the Acute Myeloid Leukemia-1a (AML-1a) transcription factor, a regulator of immune gene expression, was identified as potentially sensitive to nucleosomal regulation within the mouse Ly49 gene family. This result was confirmed in RMA, a cell line with natural expression of Ly49, using MNase-Seq to generate a nucleosome map of chromosome 6, where the Ly49 gene family is located. Analysis of this map revealed a specific depletion of nucleosomes at AML-1a binding sites in the expressed Ly49A when compared to the other, silent Ly49 genes. Our data suggest that nucleosome-based regulation contributes to the expression of Ly49 genes, and we propose that this method of predicting nucleosome sensitivity could aid in dissecting the regulatory role of nucleosomes in general.

The nucleosome—a large protein complex with DNA wound around it—is the fundamental unit of genomic organization in the eukaryotic cell. More than just a DNA organizer, however, nucleosomes may control gene expression by interfering with the cell’s ability to access the wound-up DNA, as shown by recent research. In this report, we demonstrate a computational method for predicting which elements of the genome are sensitive to regulation by nucleosomes. As a proof-of-concept, we identify AML-1a binding sites—important sequences in DNA regulation—as being specifically nucleosome sensitive. We then show that AML-1a sites are specifically depleted of nucleosomes when a gene is expressed, indicating the ability for nucleosomes to suppress the expression of that gene. This finding confirms that nucleosomes are likely involved in genome regulation, and provides a method for predicting which areas of the genome are probably affected most by nucleosomes. This paper also highlights the usefulness of the Ly49 gene family in testing computer-derived genomic predictions, and is of interest to anyone studying how gene expression is regulated from cell to cell.

The distal upstream promoter in Ly49 genes, Pro1, is active in mature NK cells and T cells, does not require TATA boxes, and displays enhancer activity

Missing self recognition of MHC class I molecules is mediated in murine species primarily through the stochastic expression of CD94/NKG2 and Ly49 receptors on NK cells. Previous studies have suggested that the stochastic expression of Ly49 receptors is achieved through the use of an alternate upstream promoter, designated Pro1, that is active only in immature NK cells and operates via the mutually exclusive binding of transcription initiation complexes to closely opposed forward and reverse TATA boxes, with forward transcription being transiently required to activate the downstream promoters, Pro2/Pro3, that are subsequently responsible for transcription in mature NK cells. In this study, we report that Pro1 transcripts are not restricted to immature NK cells but are also found in mature NK cells and T cells, and that Pro1 fragments display strong promoter activity in mature NK cell and T cell lines as well as in immature NK cells. However, the strength of promoter activity in vitro does not correlate well with Ly49 expression in vivo and forward promoter activity is generally weak or undetectable, suggesting that components outside of Pro1 are required for efficient forward transcription. Indeed, conserved sequences immediately upstream and downstream of the core Pro1 region were found to inhibit or enhance promoter activity. Most surprisingly, promoter activity does not require either the forward or reverse TATA boxes, but is instead dependent on residues in the largely invariant central region of Pro1. Importantly, Pro1 displays strong enhancer activity, suggesting that this may be its principal function in vivo.

Epigenetic interplay between mouse endogenous retroviruses and host genes

Background

Transposable elements are often the targets of repressive epigenetic modifications such as DNA methylation that, in theory, have the potential to spread toward nearby genes and induce epigenetic silencing. To better understand the role of DNA methylation in the relationship between transposable elements and genes, we assessed the methylation state of mouse endogenous retroviruses (ERVs) located near genes.

Results

We found that ERVs of the ETn/MusD family show decreased DNA methylation when near transcription start sites in tissues where the nearby gene is expressed. ERVs belonging to the IAP family, however, are generally heavily methylated, regardless of the genomic environment and the tissue studied. Furthermore, we found full-length ETn and IAP copies that display differential DNA methylation between their two long terminal repeats (LTRs), suggesting that the environment surrounding gene promoters can prevent methylation of the nearby LTR. Spreading from methylated ERV copies to nearby genes was rarely observed, with the regions between the ERVs and genes apparently acting as a boundary, enriched in H3K4me3 and CTCF, which possibly protects the unmethylated gene promoter. Furthermore, the flanking regions of unmethylated ERV copies harbor H3K4me3, consistent with spreading of euchromatin from the host gene toward ERV insertions.

Conclusions

We have shown that spreading of DNA methylation from ERV copies toward active gene promoters is rare. We provide evidence that genes can be protected from ERV-induced heterochromatin spreading by either blocking the invasion of repressive marks or by spreading euchromatin toward the ERV copy.

Models and methods for analysis of lymphocyte repertoire generation, development, selection and evolution

T and B cell receptor repertoires are diversified by variable region gene rearrangement and selected based on functionality and lack of self-reactivity. Repertoires can also be defined based on phenotype and function rather than receptor specificity - such as the diversity of T helper cell subsets. Natural killer (NK) cell repertoires, in which each cell expresses a randomly chosen subset of its inhibitory receptor genes, and is educated based on self-MHC recognition by yet unknown mechanisms, are also phenotypic repertoires. Studying the generation, development and selection of lymphocyte repertoires, and their functions during immune responses, is essential for understanding the function of the immune system in healthy individuals and in immune deficient, autoimmune or cancer patients. The study of lymphocyte repertoires will enable clinical immunologists to develop better therapeutic monoclonal antibodies, vaccines, transplantation donor-recipient matching protocols, and other immune intervention strategies. The recent development of high-throughput methods for repertoire data collection - from multicolor flow cytometry through single-cell imaging to deep sequencing - presents us now, for the first time, with the ability to analyze and compare large samples of lymphocyte repertoires in health, aging and disease. The exponential growth of these datasets, however, challenges the theoretical immunology community to develop methods for data organization and analysis. Furthermore, the need to test hypotheses regarding immune function, and generate predictions regarding the outcomes of medical interventions, necessitates the development of complex mathematical and computational models, covering processes on multiple scales, from the genetic and molecular to the cellular and system scales.

Comprehensive analysis of transcript start sites in ly49 genes reveals an unexpected relationship with gene function and a lack of upstream promoters

Comprehensive analysis of the transcription start sites of the Ly49 genes of C57BL/6 mice using the oligo-capping 5′-RACE technique revealed that the genes encoding the “missing self” inhibitory receptors, Ly49A, C, G, and I, were transcribed from multiple broad regions in exon 1, in the intron1/exon2 region, and upstream of exon -1b. Ly49E was also transcribed in this manner, and uniquely showed a transcriptional shift from exon1 to exon 2 when NK cells were activated in vitro with IL2. Remarkably, a large proportion of Ly49E transcripts was then initiated from downstream of the translational start codon. By contrast, the genes encoding Ly49B and Q in myeloid cells, the activating Ly49D and H receptors in NK cells, and Ly49F in activated T cells, were predominantly transcribed from a conserved site in a pyrimidine-rich region upstream of exon 1. An ∼200 bp fragment from upstream of the Ly49B start site displayed tissue-specific promoter activity in dendritic cell lines, but the corresponding upstream fragments from all other Ly49 genes lacked detectable tissue-specific promoter activity. In particular, none displayed any significant activity in a newly developed adult NK cell line that expressed multiple Ly49 receptors. Similarly, no promoter activity could be found in fragments upstream of intron1/exon2. Collectively, these findings reveal a previously unrecognized relationship between the pattern of transcription and the expression/function of Ly49 receptors, and indicate that transcription of the Ly49 genes expressed in lymphoid cells is achieved in a manner that does not require classical upstream promoters.

Mice lacking Ly49E show normal NK cell development and provide evidence for probabilistic expression of Ly49E in NK cells and T cells

Ly49E is an unusual member of the Ly49 family that is expressed on fetal NK cells, epithelial T cells, and NKT cells, but not on resting adult NK cells. Ly49E(bgeo/bgeo) mice in which the Ly49E gene was disrupted by inserting a β-geo transgene were healthy, fertile, and had normal numbers of NK and T cells in all organs examined. Their NK cells displayed normal expression of Ly49 and other NK cell receptors, killed tumor and MHC class I-deficient cells efficiently, and produced normal levels of IFN-γ. In heterozygous Ly49E(+/bgeo) mice, the proportion of epidermal T cells, NKT cells, and IL-2-activated NK cells that expressed Ly49E was about half that found in wild-type mice. Surprisingly, although splenic T cells rarely expressed Ly49E, IL-2-activated splenic T cells from Ly49E(bgeo/bgeo) mice were as resistant to growth in G418 as NK cells and expressed similar levels of β-geo transcripts, suggesting that disruption of the Ly49E locus had increased its expression in these cells to the same level as that in NK cells. Importantly, however, the proportion of G418-resistant heterozygous Ly49E(+/bgeo) cells that expressed Ly49E from the wild-type allele was similar to that observed in control cells. Collectively, these findings demonstrate that Ly49E is not required for the development or homeostasis of NK and T cell populations or for the acquisition of functional competence in NK cells and provide compelling evidence that Ly49E is expressed in a probabilistic manner in adult NK cells and T cells.

The requirement for NKG2D in NK cell-mediated rejection of parental bone marrow grafts is determined by MHC class I expressed by the graft recipient

Natural killer (NK) cells provide a unique barrier to semiallogeneic bone marrow (BM) transplantation. In the setting where the parents donate to the F1 offspring, rejection of parental bone marrow occurs. This "hybrid resistance" is completely NK cell dependent, as T cells in the F1 recipient tolerate parental grafts. Previously, we demonstrated that rejection of BALB/c parental BM by (BALB/c × C57BL/6) F1-recipient NK cells is dependent on the NKG2D-activating receptor, whereas rejection of parental C57BL/6 BM does not require NKG2D. BALB/c and B6 mice possess different NKG2D ligand genes and express these ligands differently on reconstituting BM cells. Herein, we show that the requirement for NKG2D in rejection depends on the major histocompatibility complex haplotype of donor cells and not the differences in the expression of NKG2D ligands. NKG2D stimulation of NK cell-mediated rejection was required to overcome inhibition induced by H-2D(d) when it engaged an inhibitory Ly49 receptor, whereas rejection of parental BM expressing the ligand, H-2K(b), did not require NKG2D. Thus, interactions between the inhibitory receptors on F1 NK cells and parental major histocompatibility complex class I ligands determine whether activation via NKG2D is required to achieve the threshold for rejection of parental BM grafts.

Inhibitory Ly49 receptors on mouse natural killer cells

The Ly49 receptors, which are expressed in a stochastic manner on subsets of murine natural killer (NK) cells, T cells, and other cells, are encoded by the Klra gene family and include receptors with either inhibitory or activating function. All of the inhibitory Ly49 receptors are characterized by an immunoreceptor tyrosine-based inhibitory motif in their cytoplasmic domain, which upon phosphorylation recruits tyrosine or lipid phosphatases to dampen signals transmitted through other activating receptors. Most of the inhibitory Ly49 receptors recognize polymorphic epitopes on major histocompatibility complex (MHC) class I proteins as ligands. Here, we review the polymorphism, ligand specificity, and signaling capacity of the inhibitory Ly49 receptors and discuss how these molecules regulate NK cell development and function.