PubMed HealthSearch

SEARCH · PubMed Health

Results for “positive selection”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Evaluating selection at intermediate scales within genes provides robust identification of genes under positive selection in M. tuberculosis clinical isolates.

Multiple studies have reported genes in the M. tuberculosis (Mtb) genome that are under diversifying selection, based on genetic variants among Mtb clinical isolates. These might reflect adaptions to selection pressures associated with modern clinical treatment of TB. Many, but not all, of these genes under selection are related to drug resistance. Most of these studies have evaluated selection at the gene-level. However, positive selection can be evaluated on different scales, including individual sites (codons) and local regions within an ORF. In this paper, we use GenomegaMap, a Bayesian method for calculating selection, to evaluate selection of genes in the Mtb genome at all three levels. We present evidence that the intermediate analysis (windows of codons) yields the most credible list of candidate genes under selection (excluding PPE and PE_PGRS genes, which are predicted less reliably due to frequent sequencing errors). A further advantage of this approach is that it identifies specific regions within proteins that are under selective pressure, which is useful for structural and functional interpretation. In an analysis of two separate collections of Mtb clinical isolates (from Moldova; and a globally-representative set), we observed 53 and 173 significant genes under selection, with 36% overlap. The lists of genes under selection include many drug-resistance genes, as well as other genes that have previously been reported to be under selection (resR, phoR). The specific regions under selection identified within drug-resistance genes are shown to correspond to protein structural features known to be involved in resistance, supporting accuracy of the method. Positive selection in several ESX-1-related genes was also observed, suggesting adaptation to immune pressure.

adaptation

Lineage-specific targets of positive selection in three leaf beetles correspond with defence capacity against their shared parasitoid wasp.

Parasitoid wasps are major causes of mortality of many species, making host immune defences a common target of adaptive evolution, though such targets outside model species are poorly understood. In this study, we used two tests of positive selection to compare across three closely related Galerucella leaf beetles that show substantial differences in their phenotypic response to the shared parasitoid wasp Asecodes parviclava, their main natural enemy. Using a codon-based test, which detects excess amino acid fixations per locus along each species' lineage, we found more evidence of positive selection on parasitoid-relevant immune genes in the species with the strongest immunocompetence (G. pusilla) compared with the species having weaker immunocompetence (G. tenella and G. calmariensis). Moreover, genes coding for the early phases in the immune response cascade were predominantly among the positively selected immune genes, providing targets for future functional genomic study to pin-point connections between genotypic and phenotypic differences in defences towards a parasitoid wasp. In contrast, genome-wide analyses of the haplotype frequency spectrum, which quantify selection over recent evolutionary time scales, revealed similar signatures of positive selection on immune genes across species. These results advance the field of host-parasitoid dynamics by providing novel insights into the tempo and mode of insect host evolutionary dynamics, and offering a framework for making genotype to phenotype connections for immunocompetence phenotypes.

Animals

Molecular dynamics simulations of positively selected codons in FcγRI reveal novel biochemical binding properties.

FcγRI is a high-affinity receptor for IgG, associated with autoimmune disease pathology and determines clinical responses to antibody-based immunotherapies. FcγRI has a complex evolutionary history that is not fully understood, and to address this we explored signatures of positive selection in the receptor's functional gene, FCGR1A, using codon-based selection tests on aligned 1-1 orthologous sequences from placental mammals (n = 32). Signatures of positive selection have occurred at several locations within the gene, with two sites (H148 (M2a ω 0.997 & M8 ω = 0.993)) and (W149 (M2a ω = 0.999 & M8 ω = 1.000)) exhibiting highest posterior probabilities, suggesting strong evidence of positive selection; these positions are known to form one of the FcγRI-IgG binding interfaces. We employed ancestral reconstruction to statistically infer prior codon sequences at these sites and identified ancestral H148P and W149R codons at different nodes in the phylogeny. Employing molecular dynamics simulations, we determined how evolutionary changes at these sites may have influenced the binding of FcγRI-IgG of modern-day Homo sapiens. Measuring RMSD, free energy, radius of gyration, hydrogen bond formation, and analyzing free energy landscapes, we demonstrate that structural instability between mutant structures vs the WT counterpart; however, overall binding potential increases at position 148, yet decreases at 149 in potential. H148P protonation at physiological pH remains similar, yet during acidotic calculations, protonation is likely reduced, with predicted reduction in affinity for IgG. While ancestral W149R substitutions demonstrate an implication for electron conjugation. Examining key sites at this binding FcγRI-IgG interface, our data demonstrate that these two codons have evolved in humans to be relatively insensitive to shifts in pH promoting a more stable interaction with the Fc portion of IgG during diseases that promote acidosis.

Receptors, IgG

Detecting Interspecific Positive Selection Using Convolutional Neural Networks.

Traditional statistical methods using maximum likelihood and Bayesian inference can detect positive selection from an interspecific phylogeny and a codon sequence alignment based on model assumptions, but they are prone to false positives due to alignment errors and can lack power. These problems are particularly pronounced when faced with high levels of indels and divergence. To address these issues, we trained and tested convolutional neural network models on simulated data and achieved higher accuracy in detecting selection across a specific range of phylogenetic scenarios and evolutionary modes. This advantage is particularly evident when performing inference on noisy data prone to misalignments. Our method shows some ability to account for these errors, where most statistical frameworks fail to do so in a tractable manner. We explore the generalizability of our convolutional neural network models to unseen evolutionary scenarios and identify future avenues to achieve broader utility. Once trained, our convolutional neural network model is faster at test time, making it a scalable alternative to traditional statistical methods for large-scale, multigene analyses. In addition to binary classification (inference of the presence or absence of positive selection during the evolution of the sequences), we use saliency maps to understand what the model learns and observe how this could be leveraged for sitewise inference of positive selection.

Neural Networks, Computer

Positive selection and relaxed purifying selection contribute to rapid evolution of sex-biased genes in green seaweed Ulva.

BACKGROUND: The evolution of differences in gamete size and number between sexes is a cornerstone of sexual selection theories. The green macroalga Ulva, with incipient anisogamy and parthenogenetic gametes, provides a unique system to investigate theoretical predictions regarding the evolutionary pressures that drive the transition from isogamy to anisogamy, particularly in relation to gamete size differentiation and sexual selection. Its minimal gamete dimorphism and facultative parthenogenesis enable a rare window into early evolutionary steps toward anisogamy. RESULTS: By analyzing the expression profiles of sex-biased genes (SBGs) during gametogenesis, we found that SBGs evolve faster than unbiased genes, driven by higher rates of non-synonymous substitution (dN), indicating that SBGs are under stronger selective pressures. Mating type minus-biased genes (mt-BGs) exhibit higher dN/dS values than mating type plus-biased genes (mt+BGs), suggesting stronger selective pressures on mt-BGs, although this difference was not statistically significant (P = 0.08). Using branch-site and RELAX models, we found positive selection and relaxed purifying selection acting on a significant proportion of SBGs, particularly those associated with flagella function. CONCLUSIONS: This study highlights the selective pressures shaping anisogamy and provides insights into the molecular mechanisms underlying its evolution. The faster evolution of SBGs, particularly mt-BGs, and the positive selection on genes associated with motility, such as those related to flagella function, suggest the importance of enhanced gamete motility in the transition to anisogamy. These findings contribute to our understanding of sexual selection and the evolutionary forces that drive the differentiation of gamete size and number between sexes.

Selection, Genetic

Pervasive positive selection on X-linked ampliconic genes in primates.

Mammalian sex chromosomes harbour ampliconic gene families, which are multi-copy genes with ≥97% sequence identity, predominantly expressed in testis tissue and essential for male fertility. The amplification of testis-specific genes is conserved across mammals, yet the specific gene families that expand show striking lineage-specific variation. Previous studies suggest a dynamic turnover with adaptive evolution for several of these families, but their analysis has been limited by the quality of reference genomes of repetitive regions. To characterise the molecular evolutionary processes of ampliconic gene families on both sex chromosomes, we analysed telomere-to-telomere genome assemblies from eight primate species spanning 25 million years of evolution. We identified 53 X-linked and 19 Y-linked ampliconic gene families with dynamic copy number variation. Gene conversion through palindromic pairing and tandem arrays maintained high sequence similarity despite accumulating mutations. X-linked families maintained conserved chromosomal positions despite copy number changes, whereas Y-linked families showed frequent positional turnover. Strikingly, multiple X-linked families (GAGE, SSX, CSAG, and VCX) showed pervasive positive selection across the primate phylogeny and multiple (MAGEB, CT45, HSFX) showed lineage specific positive selection. Y-linked families predominantly evolve under purifying selection. Examining intraspecific copy number variation of the X-linked ampliconic families in chimpanzees, humans, and gorillas, we found variation among individuals but clear differences between species, with the largest families varying the most. These patterns could suggest that sperm competition, meiotic drive, or dosage-dependent selection drive the rapid, lineage-specific evolution of testis-expressed ampliconic genes in primates.

Journal Article

Identifying loci under positive selection in complex population histories.

Detailed modeling of a species' history is of prime importance for understanding how natural selection operates over time. Most methods designed to detect positive selection along sequenced genomes, however, use simplified representations of past histories as null models of genetic drift. Here, we present the first method that can detect signatures of strong local adaptation across the genome using arbitrarily complex admixture graphs, which are typically used to describe the history of past divergence and admixture events among any number of populations. The method-called graph-aware retrieval of selective sweeps (GRoSS)-has good power to detect loci in the genome with strong evidence for past selective sweeps and can also identify which branch of the graph was most affected by the sweep. As evidence of its utility, we apply the method to bovine, codfish, and human population genomic data containing panels of multiple populations related in complex ways. We find new candidate genes for important adaptive functions, including immunity and metabolism in understudied human populations, as well as muscle mass, milk production, and tameness in specific bovine breeds. We are also able to pinpoint the emergence of large regions of differentiation owing to inversions in the history of Atlantic codfish.

Animals

Proteomic Analysis of Three Independent Series of Sequential Cystic Fibrosis Strains in an International Pseudomonas aeruginosa Reference Panel Indicates Positive Selection in Late Infection Strains.

Pseudomonas aeruginosa is a highly diverse, adaptable Gram-negative bacterium that thrives in many environments and is a frequent cause of chronic opportunistic infections in people with cystic fibrosis (CF). P. aeruginosa adapts over time of colonization to facilitate chronic infection, including loss of virulence factors; however, proteomic analyses of the adaptation to chronic infection have been limited. We previously assembled and characterized an international panel of P. aeruginosa strains from diverse clinical and geographical sources, including three sets of sequential CF isolates, enabling identification of conserved adaptations linked to chronic CF lung colonization. We compared the proteomes of eight strains (three early and five late infection) to assess whether common proteomic changes emerged during colonization across the sequential isolates. We identified 11 proteins showing increased abundance in late isolates in all three series, many of which are associated with virulence, regulation of virulence, or response to hypoxia. These include CF inhibitory factor repressor (CifR); WspR; 2 two-component response regulators (PA2572 and PA3702), and transcriptional regulator (PA2551). Moreover, we identified three proteins (PA2572, PA3819, and PA5028) that showed increased abundance in all five late isolates. The probability of this being random is 5.06 × 10-53 and therefore, strong evidence of positive selection. The increased abundance of PA2573 and PA3819 appears to improve the fitness of P. aeruginosa in response to antibiotics and oxidative stress. All three proteins share a tyrosine phosphorylation motif, suggesting a common regulatory mechanism. Overall, despite substantial diversity across P. aeruginosa, common adaptations occur in the CF lung.

Pseudomonas aeruginosa

Disruption of mitonuclear coadaptation and compensatory evolution after an extreme dietary shift in carnivorous butterflies.

Mitochondrial function depends on tight coordination between mitochondrial and nuclear genomes, which requires long-term coevolution to maintain mitonuclear coadaptation. While mitonuclear incompatibility is typically studied in the context of hybridization, other evolutionary scenarios that may disrupt coadaptation between the two genomes remain less explored. Here, we propose that extreme ecological niche shifts may disrupt mitonuclear coadaptation, which we test in carnivorous Miletinae butterflies with an extreme dietary transition. By generating high-quality genome assemblies, we found that Miletinae exhibit extensive chromosomal rearrangements. Comparative phylogenomic analyses revealed a striking asymmetric mitonuclear evolutionary response: Miletinae exhibit elevated mitochondrial nucleotide substitution rates compared to phytophagous relatives, whereas nuclear rates remain stable. This shift reverses the typical lepidopteran pattern where nuclear rates exceed mitochondrial rates. Interestingly, this mitochondrial acceleration is driven primarily by relaxed purifying selection rather than positive selection. To sustain mitochondrial function, the nuclear genome of Miletinae underwent pervasive, multilayered compensatory evolution. We detected strong signatures of positive selection and accelerated evolution in nuclear genes directly interacting with mitochondrial components across oxidative phosphorylation (OXPHOS) complexes, the mitochondrial translation, and replication and transcription machinery. Furthermore, this nuclear compensatory response extends to systems governing mitochondrial homeostasis, including protein quality control and RNA degradation and stabilization. Our results support a model in which extreme ecological transitions can disrupt ancestral mitonuclear coadaptation and promote the emergence of a new coadapted state through systemic nuclear compensation. This study broadens the conceptual framework of mitonuclear coevolution and highlights its role in facilitating evolutionary persistence after major ecological shifts.

Animals

Comparative Phylogenetics Reveal Clade-specific Drivers of Recombination Rate Evolution Across Vertebrates.

Meiotic recombination is an integral cellular process, required for the production of viable gametes. Recombination rate is a fundamental genomic parameter, modulating genomic responses to selection. Our increasingly detailed understanding of its molecular underpinnings raises the prospect that we can gain insight into trait divergence by examining the molecular evolution of recombination genes from a pathway perspective, as in mammals, where protein-coding changes in later stages of the recombination pathway are connected to divergence in intra-clade recombination rate. Here, we leverage increased availability of avian and teleost genomes to reconstruct the evolution of the recombination pathway across two additional vertebrate clades: birds, which have higher and more variable rates of recombination and similar divergence times to mammals, and teleost fish, which have much deeper divergence times. Rates of molecular evolution of recombination genes are highly correlated between vertebrate clades and significantly elevated compared to control panels, suggesting that they experience similar selective pressures. Avian recombination genes are significantly more likely to exhibit signatures of positive selection than other clades, unrestricted to later stages of the pathway. Signatures of positive selection in genes linked to recombination rate variation in mammalian populations and those with signatures of positive selection across the avian phylogeny are highly correlated. In contrast, teleost fish recombination genes have significantly less evidence of positive selection despite high intra-clade recombination rate variability. Gaining clade-specific understanding of patterns of variation in recombination genes can elucidate drivers of recombination rate and thus, factors influencing genetic diversity, selection efficacy, and species divergence.

Animals

Molecular evolution of the members of the Snq2/Pdr18 subfamily of Pdr transporters in the Hemiascomycete yeasts.

The transporters of the ATP-Binding Cassette (ABC) Superfamily involved in the Multidrug Resistance (MDR) phenomena are also known as ABC-Pleiotropic Drug Resistance (PDR) proteins. The homologs of the Saccharomyces cerevisiae SNQ2 and PDR18 genes were identified in 171 yeast genomes, representing 68 different hemiascomycetous species. All early-divergent yeast species analyzed in this work lack Snq2/Pdr18 homologs, suggesting that the origin of these ABC-PDR genes in hemiascomycete yeasts resulted from a horizontal transfer event. The evolutionary pathway of the Snq2/Pdr18 protein subfamily in pathogenic Candida species was also reconstructed, revealing a main gene lineage leading to the Candida albicans SNQ2 gene. The results indicate that, after the gene duplication event at the origin of the SNQ2/PDR18 paralogs, the PDR18 ortholog has been under strong diversifying selection and suggest that a small portion of the sequence of the SNQ2 ancestral ortholog might have been under mild positive selection. The results also showed that strong positive selection was exerted over one of the two paralogs generated by the Whole Genome Duplication (WGD) event, corresponding to the duplicate at the origin of a "short-lived" WGD sublineage.

Evolution, Molecular

Contrasting signals of selection at the EDAR gene in global and Latin American populations.

The EDAR gene is a classic target of positive selection in humans, mainly through the nonsynonymous variant rs3827760 (EDARV370A) associated with ectodermal traits. Using high-resolution data from the 1000 Genomes Project, we combined sliding-window F_ST, BayeScan, and extended haplotype homozygosity (EHH) analyses to examine global and Latin American patterns of differentiation. Globally, a strong signal of positive selection was confirmed at EDAR, dominated by the rs3827760 haplotype background and its extended linkage disequilibrium structure. In contrast, within Latin America, differentiation reflected admixture-driven haplotype persistence rather than contemporary selection. A genome-wide FST scan comparing individuals from the upper and lower quartiles of Native American ancestry showed that EDAR lies among the most highly differentiated regions in this contrast, consistent with ancestry-driven haplotype structure rather than post-admixture adaptive evolution. These results indicate that EDAR retains its evolutionary signature globally but not within recently admixed populations, where demographic history rather than selection shapes its genetic landscape.

Humans

Multiomic Analyses Reveal the Molecular Mechanisms of Arid Adaptation in a Desert Rodent Species.

Organisms living in desert habitats face multiple simultaneous pressures, such as high temperatures and arid, and the population dynamics and community diversity of small rodents are strongly affected by climate extremes. However, the potential mechanisms by which desert rodents adapt to arid remain largely unexplored. Here, we assembled a 3.18 Gb genome, including 25,812 protein-encoding genes, for Orientallactaga sibirica, which is widely distributed across both arid and semihumid environments in Eurasia. Orientallactaga sibirica has longer ears and hind limbs to enhance heat dissipation, which may be related to the positively selected genes, such as Fgf10, Fgf11, Hoxc4, Hoxd1, and Bmp4. The renal transcriptome revealed increased fat and carbohydrate metabolism for metabolic water production in O. sibirica residing in arid habitats. Pathways such as material metabolism, oxidative stress response, osmoregulation, and water and salt reabsorption were enriched in candidate genes, such as Avp, Ang, and Ace, under positive selection in O. sibirica. Moreover, amino acid replacement was observed in the protein sequences of seven candidate genes, including Aldh7a1, Lnpep, Wnk4, C1qc, and Awat2, and these specific amino acid replacements of genes such as Umod and Scnn1a were related to unique osmoregulation, osmotic protection, and water retention compensation mechanisms. Water deprivation under laboratory conditions induced the upregulation of Umod and Aldh7a1 expression, further supporting the results observed in the wild population. These findings demonstrate that the positively selected genes related to limb development and specific amino acid replacements in the genes Umod and Scnn1a for unique osmoregulation in the renal vascular system may contribute to arid adaptation in the desert rodent species O. sibirica. This study provides novel insights into the adaptive evolution of desert small mammals and can serve as a reference for future research on renal damage-related diseases, such as human kidney stones and salt-sensitive hypertension.

Animals

Primulina pan-genome reveals differential gene retention following whole-genome duplications and provides insights into edaphic specialization.

Primulina, a genus of >200 species specialized to extreme soils, provides a model for edaphic adaptation. We assemble seven genomes and construct a pan-genome spanning nine species from karst, Danxia, and acidic soils. Comparative analyses reveal that karst-adapted species have smaller genomes. Two lineage-specific whole-genome duplications (WGDs) exhibit biased duplicate loss in large gene families but preferential retention of transcription factors, indicating combined adaptive and nonadaptive forces. Pan-genome analyses identify ion channel and transporter genes enriched in variant hotspots and under positive selection in karst lineages. Candidate genes for drought and salt stress tolerance include ABC transporters and ion channels. Notably, an ABC transporter shows positive selection in karst species and unique structural variation in non-karst species. Together, our findings show that genome downsizing, biased post-WGD retention, and evolution of ion-transport pathways shape adaptation to extreme soils. The Primulina pan-genome provides a resource for dissecting mechanisms underlying edaphic specialization.

Gene Duplication

The Demographic History of Populations and Genomic Imprinting have Shaped the Transposon Patterns in Arabidopsis lyrata.

Purifying selection is expected to prevent the accumulation of transposable elements (TEs) within their host, especially when located in and around genes and if affected by epigenetic silencing. However, positive selection may favor the spread of TEs, causing genomic imprinting under parental conflict, as genomic imprinting allows parent-specific influence over resource accumulation to the progeny. Concomitantly, the number and frequency of TE insertions in natural populations are conditioned by demographic events. In this study, we aimed to test how demography and selective forces interact to affect the accumulation of TEs around genes, depending on their epigenetic silencing, with a particular focus on imprinted genes. To this aim, we compared the frequency and distribution of TEs in Arabidopsis lyrata from Europe and North America. Generally, we found that TE insertions showed a lower frequency when they were inserted in or near genes, especially TEs targeted by epigenetic silencing, suggesting purifying selection at work. We also found that many TEs were lost or got fixed in North American populations during the colonization and the postglacial range expansion from refugia of the species in North America, as well as during the transition to selfing, suggesting a potential "TE load." Finally, we found that silenced TEs increased in frequency and even tended to reach fixation when they were linked to imprinted genes. We conclude that in A. lyrata, genomic imprinting has spread in natural populations through demographic events and positive selection acting on silenced TEs, potentially under a parental conflict scenario.

DNA Transposable Elements

Evolutionary Reorganization of Transcriptomic Architecture Across a UVB Tolerance Gradient in Fish.

Environmental stressors such as ultraviolet radiation impose strong selective pressures on organisms, yet how adaptation to such stressors shapes transcriptomic responses at the network level remains poorly understood. Although stratospheric ozone is recovering globally, substantial regional variation in UV exposure persists, particularly in high-altitude environments where extreme UV levels can occur. Here, we compared three fish models representing distinct biological responses to UVB exposure: wild-type zebrafish (Danio rerio), a melanin-deficient zebrafish mutant (nacre) lacking a major protective mechanism against UVB damage, and the high-altitude Andean killifish Orestias ascotanensis, a species naturally exposed to extreme UVB radiation. Together, these models define a gradient spanning physiological protection, impaired protection, and evolutionary adaptation to UVB stress. Using RNA-seq and protein-protein interaction networks, we show that transcriptomic responses differ markedly across this gradient. Wild-type and nacre zebrafish exhibited relatively limited transcriptomic changes (∼2%-2.4% of genes changing), whereas O. ascotanensis displayed a large-scale and highly coordinated response (∼21.6% of genes changing) characterized by functionally specialized networks enriched in DNA repair pathways. These differences involved not only transcriptomic magnitude but also marked reorganization of transcriptomic architecture. Integration with positive selection analyses revealed that positively selected genes were concentrated within highly interconnected regions of transcriptomic networks, consistent with adaptation involving network reorganization. Furthermore, ortholog-based analyses suggest that adaptive responses involve differential reorganization of a conserved functional background. Together, our results support a model in which adaptation to environmental stress is associated with the reorganization of conserved transcriptomic networks across physiological and evolutionary contexts, providing a systems-level perspective on the molecular basis of adaptation.

UVB radiation

Adaptive Evolution of the PFK Gene Family in Chinese Longsnout Catfish, Leiocassis longirostris.

The Chinese longsnout catfish is a typical carnivorous fish with a relatively weak ability to utilize glucose. However, the genomic basis for its glucose metabolic adaptation remains unclear. In this study, we used comparative genomics methods to systematically analyze the evolutionary characteristics of glucose metabolism-related genes in the Chinese longsnout catfish, focusing on gene family evolution, patterns of expansion and contraction, and selective pressures. The results indicate that glucose metabolism-related genes have undergone significant reshaping during evolution. Genes involved in glucose digestion, absorption, and insulin signaling pathways demonstrate a tendency toward contraction, while those associated with protein and lipid metabolism exhibit expansion. This pattern is consistent with the species' long-term adaptation to a high-protein, high-fat diet. Comparative analysis further revealed that, compared to fish with different dietary habits, certain key genes involved in glycolysis in the Chinese longsnout catfish exhibit a reduction in copy number. Molecular evolutionary analysis showed that key genes involved in glycolysis and gluconeogenesis (including hexokinase 2 (hk2), phosphofructokinase, muscle/platelet (pfkm/p)) exhibit signs of accelerated evolution or positive selection. Notably, the PFK gene family exhibits complex evolutionary characteristics resulting from the combined effects of gene contraction, rapid evolution, and positive selection. In summary, this study reveals the genomic evolutionary basis for the glucose metabolic adaptation of the Chinese longsnout catfish and identifies the PFK gene family as a key candidate for elucidating its unique glucose metabolic characteristics.

Leiocassis longirostris

Jarid2 is induced by TCR signalling and controls iNKT cell maturation.

Jarid2 is a reported component of three lysine methyltransferase complexes, polycomb repressive complex 2 (PRC2) that methylates histone 3 lysine 27 (H3K27), and GLP-G9a and SETDB1 complexes that methylate H3K9. Here we show that Jarid2 is upregulated upon TCR stimulation and during positive selection in the thymus. Mice lacking Jarid2 in T cells display an increase in the frequency of IL-4-producing promyelocytic leukemia zinc finger (PLZF)(hi) immature invariant natural killer T (iNKT) cells and innate-like CD8(+) cells; Itk-deficient mice, which have a similar increase of innate-like CD8(+) cells, show blunted upregulation of Jarid2 during positive selection. Jarid2 binds to the Zbtb16 locus, which encodes PLZF, and thymocytes lacking Jarid2 show increased PLZF and decreased H3K9me3 levels. Jarid2-deficient iNKT cells perturb Th17 differentiation, leading to reduced Th17-driven autoimmune pathology. Our results establish Jarid2 as a novel player in iNKT cell maturation that regulates PLZF expression by modulating H3K9 methylation.

Animals