Inference of selective forces on house mouse genomes during secondary contact in East Asia.

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca).

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

Soluble SLAMF7 is generated by alternative splicing in multiple myeloma cells.

Not available.

Phylogenetics and Population Genetics of the Asian House Shrew, Suncus murinus-S. montanus Species Complex, Inferred From Whole-Genome and Mitochondrial DNA Sequences, with Special Reference to the Ryukyu Archipelago, Japan.

The house shrew (Suncus murinus-S. montanus species complex) colonized regions across southern Asia and the Indian Ocean following human activity. The house shrew is distributed on islands of the Ryukyu Archipelago, the southernmost part of Japan, but the evolutionary history of the shrew on those islands and possible associations between these populations and humans remain unknown. In this study, we conducted phylogenetic and population genetic analyses based on both nuclear and mitochondrial genome sequences of house shrews. Phylogenetic analyses based on mitochondrial cytochrome b (cytb) sequences revealed that shrews from the Ryukyu Archipelago showed strong genetic affinity to Vietnamese and southern Chinese shrews. Demographic analyses of cytb sequences indicated a rapid population expansion event affecting the haplotype group in Vietnam, southern China, and the Ryukyu Archipelago 3300-7900 years ago. Furthermore, gene flow between Ryukyu (Yonaguni Island) and Taiwan and between Ryukyu and Vietnam inferred from f4 statistics of the nuclear genomes suggested repeated immigration to Ryukyu in recent years. The present study demonstrates that the Nagasaki population has a different origin from the Ryukyu population. These findings elucidate the complex pattern of genetic admixture in house shrews and provide insights into their evolutionary history.

House mouse Mus musculus dispersal in East Eurasia inferred from 98 newly determined complete mitochondrial genome sequences.

The Eurasian house mouse Mus musculus is useful for tracing prehistorical human movement related to the spread of farming. We determined whole mitochondrial DNA (mtDNA) sequences (ca. 16,000 bp) of 98 wild-derived individuals of two subspecies, M. m. musculus (MUS) and M. m. castaneus (CAS). We revealed directional dispersals reaching as far as the Japanese Archipelago from their homelands. Our phylogenetic analysis indicated that the eastward movement of MUS was characterised by five step-wise regional extension events: (1) broad spatial expansion into eastern Europe and the western part of western China, (2) dispersal to the eastern part of western China, (3) dispersal to northern China, (4) dispersal to the Korean Peninsula and (5) colonisation and expansion in the Japanese Archipelago. These events were estimated to have occurred during the last 2000-18,000 years. The dispersal of CAS was characterised by three events: initial divergences (ca. 7000-9000 years ago) of haplogroups in northernmost China and the eastern coast of India, followed by two population expansion events that likely originated from the Yangtze River basin to broad areas of South and Southeast Asia, including Sri Lanka, Bangladesh and Indonesia (ca. 4000-6000 years ago) and to Yunnan, southern China and the Japanese Archipelago (ca. 2000-3500). This study provides a solid framework for the spatiotemporal movement of the human-associated organisms in Holocene Eastern Eurasia using whole mtDNA sequences, reliable evolutionary rates and accurate branching patterns. The information obtained here contributes to the analysis of a variety of animals and plants associated with prehistoric human migration.

Conversion of human fibroblasts into multipotent cells by cell-penetrating peptides.

The potential application of human induced pluripotent stem cells (hiPSCs) brings great expectations to regenerative medicine. However, several safety concerns, such as oncogenic transformation, remain. A number of methods have been developed to produce hiPSCs with potentially reduced risks. Cell-penetrating peptides (CPPs) are expected to improve the efficiency of nonviral reprogramming by delivering biologically active molecules into cells. Here, we show that the transfection of CPPs alone into normal adult human fibroblasts generated embryonic body (EB)-like cell clusters in the absence of reprogramming factors. The CPP-generated cell clusters were positive for a set of multipotency markers and differentiated into endodermal, ectodermal, and mesodermal cells in vitro. These results suggest that CPPs converted normal human adult somatic cells into multipotent cells. Moreover, we show that CPPs dissociated histone deacetylase 1 and lysine-specific demethylase 1 from the promoter/enhancer regions of reprogramming factors to reactivate their expression. This is the first report of an easy and quick method for somatic cell reprogramming by CPPs and a novel mechanism of reprogramming. The potential application of CPP-generated multipotent cells resolves several concerns, especially safety issues, in regenerative medicine.

Ancient hybridization and admixture in macaques (genus Macaca) inferred from whole genome sequences.

The evolutionary history of the stump-tailed macaque (Macaca arctoides) and its genetic relationship to other macaques is a subject of continuing controversy. Here, we have reported the first genome sequences of two stump-tailed macaques and one Assamese macaque (M. assamensis). Additionally, we have investigated the genetic diversity between macaque species and analyzed ancient hybridization events. Genome-wide analyses demonstrated that the stump-tailed macaque is more closely related to sinica species than to fascicularis/mulatta species. This topology contradicts the mitochondrial sequence-based phylogeny that places the stump-tailed macaque into the fascicularis/mulatta group. However, our results further show that stump-tailed macaques have genetic backgrounds distinct from sinica species, and present evidence of gene flows with rhesus macaques. We suggest that an ancient introgression occurred after stump-tailed macaques diverged from sinica species. The distinct gene flow between proto-arctoides and proto-mulatta resulted in the transfer of rhesus macaque-type mitochondria into proto-arctoides. The rhesus macaque-type mitochondria remained in populations because of genetic drift during the bottleneck. The PSMC results and morphological and geographic evidence are consistent with the mitochondria capture pattern in the stump-tailed macaque. The molecular clock estimates suggest that the mitochondrial transference into stump-tailed macaques occurred 0.4-1.4 million years ago. Furthermore, we detected extensive admixtures between different macaque species, indicating that gene flow has played an important role in the evolutionary history of the genus Macaca.

Comprehensive phylogenomic analysis reveals a novel cluster of simian endogenous retroviral sequences in Colobinae monkeys.

Simian retrovirus (SRV) is a type-D betaretrovirus infectious to the Old World monkeys causing a variety of symptoms. SRVs are also present in the Old World monkey genomes as endogenous forms, which are referred to as Simian endogenous retroviruses (SERVs). Although many SERV sequences have been identified in Cercopithecinae genomes, with potential of encoding all functional genes, the distribution of SERVs in genomes and evolutionary relationship between exogeneous SRVs and SERVs remains unclear. In this study, we comprehensively investigated seven draft genome sequences of the Old World monkeys, and identified a novel cluster of SERVs in the two Rhinopithecus (R. roxellana and R. bieti) genomes, which belong to the Colobinae subfamily. The Rhinopithecus genomes harbored higher copy numbers of SERVs than the Cercopithecinae genomes. A reconstructed phylogenetic tree showed that the Colobinae SERVs formed a distinct cluster from SRVs and Cercopithecinae SERVs, and more closely related to exogenous SRVs than Cercopithecinae SERVs. Three radical amino acid substitutions specific to Cercopithecinae SERVs, which potentially affect the infectious ability of SERVs, were also identified in the proviral envelope protein. In addition, we found many integration events of SERVs were genus- or species-specific, suggesting the recent activity of SERVs in the Old World monkey genomes. The results suggest that SERVs in Cercopithecinae and Colobinae monkeys were endogenized after the divergence of subfamilies and do not transmit across subfamilies. Our findings also support the hypothesis that Colobinae SERVs are direct ancestors of SRV-6, which has a different origin from the other exogenous SRVs. These findings shed novel insight into the evolutionary history of SERVs, and may help to understand the process of endogenization of SRVs.

Adaptive divergence with gene flow in incipient speciation of Miscanthus floridulus/sinensis complex (Poaceae).

Young incipient species provide ideal materials for untangling the process of ecological speciation in the presence of gene flow. The Miscanthus floridulus/sinensis complex exhibits diverse phenotypic and ecological differences despite recent divergence (approximately 1.59 million years ago). To elucidate the process of genetic differentiation during early stages of ecological speciation, we analyzed genomic divergence in the Miscanthus complex using 72 randomly selected genes from a newly assembled transcriptome. In this study, rampant gene flow was detected between species, estimated as M = 3.36 × 10(-9) to 1.20 × 10(-6) , resulting in contradicting phylogenies across loci. Nevertheless, beast analyses revealed the species identity and the effects of extrinsic cohesive forces that counteracted the non-stop introgression. As expected, early in speciation with gene flow, only 3-13 loci were highly diverged; two to five outliers (approximately 2.78-6.94% of the genome) were characterized by strong linkage disequilibrium, and asymmetrically distributed among ecotypes, indicating footprints of diversifying selection. In conclusion, ecological speciation of incipient species of Miscanthus probably followed the parapatric model, whereas allopatric speciation cannot be completely ruled out, especially between the geographically isolated northern and southern M. sinensis, for which no significant gene flow across oceanic barriers was detected. Divergence between local ecotypes in early-stage speciation began at a few genomic regions under the influence of natural selection and divergence hitchhiking that overcame gene flow.

Ecological genomics in Xanthomonas: the nature of genetic adaptation with homologous recombination and host shifts.

BACKGROUND: Comparative genomics provides insights into the diversification of bacterial species. Bacterial speciation usually takes place with lasting homologous recombination, which not only acts as a cohering force between diverging lineages but brings advantageous alleles favored by natural selection, and results in ecologically distinct species, e.g., frequent host shift in Xanthomonas pathogenic to various plants. RESULTS: Using whole-genome sequences, we examined the genetic divergence in Xanthomonas campestris that infected Brassicaceae, and X. citri, pathogenic to a wider host range. Genetic differentiation between two incipient races of X. citri pv. mangiferaeindicae was attributable to a DNA fragment introduced by phages. In contrast to most portions of the genome that had nearly equivalent levels of genetic divergence between subspecies as a result of the accumulation of point mutations, 10% of the core genome involving with homologous recombination contributed to the diversification in Xanthomonas, as revealed by the correlation between homologous recombination and genomic divergence. Interestingly, 179 genes were under positive selection; 98 (54.7%) of these genes were involved in homologous recombination, indicating that foreign genetic fragments may have caused the adaptive diversification, especially in lineages with nutritional transitions. Homologous recombination may have provided genetic materials for the natural selection, and host shifts likely triggered ecological adaptation in Xanthomonas. To a certain extent, we observed positive selection nevertheless contributed to ecological divergence beyond host shifting. CONCLUSION: Altogether, mediated with lasting gene flow, species formation in Xanthomonas was likely governed by natural selection that played a key role in helping the deviating populations to explore novel niches (hosts) or respond to environmental cues, subsequently triggering species diversification.

Whole-genome sequencing of tibetan macaque (Macaca Thibetana) provides new insight into the macaque evolutionary history.

Macaques are the most widely distributed nonhuman primates and used as animal models in biomedical research. The availability of full-genome sequences from them would be essential to both biomedical and primate evolutionary studies. Previous studies have reported whole-genome sequences from rhesus macaque (Macaca mulatta) and cynomolgus macaque (M. fascicularis, CE), both of which belong to the fascicularis group. Here, we present a 37-fold coverage genome sequence of the Tibetan macaque (M. thibetana; TM). TM is an endemic species to China belonging to the sinica group. On the basis of mapping to the rhesus macaque genome, we identified approximately 11.9 million single-nucleotide variants), of which 3.9 million were TM specific, as assessed by comparison two Chinese rhesus macaques (CR) and two CE genomes. Some genes carried TM-specific homozygous nonsynonymous variants (TSHNVs), which were scored as deleterious in human by both PolyPhen-2 and SIFT (Sorting Tolerant From Intolerant) and were enriched in the eye disease genes. In total, 273 immune response and disease-related genes carried at least one TSHNV. The heterozygosity rates of two CRs (0.002617 and 0.002612) and two CEs (0.003004 and 0.003179) were approximately three times higher than that of TM (0.000898). Polymerase chain reaction resequencing of 18 TM individuals showed that 29 TSHNVs exhibited high allele frequencies, thus confirming their low heterozygosity. Genome-wide genetic divergence analysis demonstrated that TM was more closely related to CR than to CE. We further detected unusual low divergence regions between TM and CR. In addition, after applying statistical criteria to detect putative introgression regions (PIRs) in the TM genome, up to 239,620 kb PIRs (8.84% of the genome) were identified. Given that TM and CR have overlapping geographical distributions, had the same refuge during the Middle Pleistocene, and show similar mating behaviors, it is highly likely that there was an ancient introgression event between them. Moreover, demographic inferences revealed that TM exhibited a similar demographic history as other macaques until 0.5 Ma, but then it maintained a lower effective population size until present time. Our study has provided new insight into the macaque evolutionary history, confirming hybridization events between macaque species groups based on genome-wide data.

Complex patterns of cis-regulatory polymorphisms in ebony underlie standing pigmentation variation in Drosophila melanogaster.

Pigmentation traits in adult Drosophila melanogaster were used in this study to investigate how phenotypic variations in continuous ecological traits can be maintained in a natural population. First, pigmentation variation in the adult female was measured at seven different body positions in 20 strains from the Drosophila melanogaster Genetic Reference Panel (DGRP) originating from a natural population in North Carolina. Next, to assess the contributions of cis-regulatory polymorphisms of the genes involved in the melanin biosynthesis pathway, allele-specific expression levels of four genes were quantified by amplicon sequencing using a 454 GS Junior. Among those genes, ebony was significantly associated with pigmentation intensity of the thoracic segment. Detailed sequence analysis of the gene regulatory regions of this gene indicated that many different functional cis-regulatory alleles are segregating in the population and that variations outside the core enhancer element could potentially play important roles in the regulation of gene expression. In addition, a slight enrichment of distantly associated SNP pairs was observed in the ~10 kb cis-regulatory region of ebony, which suggested the presence of interacting elements scattered across the region. In contrast, sequence analysis in the core cis-regulatory region of tan indicated that SNPs within the region are significantly associated with allele-specific expression level of this gene. Collectively, the data suggest that the underlying genetic differences in the cis-regulatory regions that control intraspecific pigmentation variation can be more complex than those of interspecific pigmentation trait differences, where causal genetic changes are typically confined to modular enhancer elements.

CYP2D44 polymorphisms in cynomolgus and rhesus macaques.

Macaques, including cynomolgus and rhesus macaques, are important animal species used in drug metabolism studies. CYP2D44 is expressed in cynomolgus macaque liver and encodes a functional drug metabolizing enzyme, metabolizing typical human CYP2D substrates such as bufuralol and dextromethorphan. CYP2D44 is highly homologous to human CYP2D6 that is known to be polymorphic with a large inter-individual variation in metabolic activities, however, genetic polymorphisms have not been investigated in macaque CYP2D44. In the present study, screening of 78 cynomolgus and 40 rhesus macaques found a total of 67 variants, including 64 non-synonymous variants, 1 nonsense mutation, and 2 frameshift mutations, and 1 gene conversion, of which 14, 19, and 15 variants were unique to Indochinese cynomolgus macaques, Indonesian cynomolgus macaques, and Chinese rhesus macaques, respectively. Eleven of the 64 non-synonymous variants were located in substrate recognition sites, the regions important for protein function. By site-directed mutagenesis and metabolic assays, S175N, V185L, A235G, R242G, R245K, and N337D showed substantially decreased activity in bufuralol 1'-hydroxylation as compared with wild-type proteins. Moreover, two null alleles (c.128T>del and c.664G>T) were found in Indonesian cynomolgus macaques, but not in Indochinese cynomolgus macaques or Chinese rhesus macaques. These results suggest that genetic polymorphisms might account for the variability of CYP2D44-dependent metabolism in macaques.

Systematic classification of vertebrate chemokines based on conserved synteny and evolutionary history.

The genes involved in host defences are known to undergo rapid evolution. Therefore, it is often difficult to assign orthologs in multigene families among various vertebrate species. Chemokines are a large family of small cytokines that orchestrate cell migration in health and disease. Herein, we have surveyed the genomes of 18 representative vertebrate species for chemokine genes and identified a total of 553 genes. We have determined their orthologous relationships and classified them in accordance with the current systematic chemokine nomenclature system. Our study reveals an interesting evolutionary history that gave origin and diversification to the vertebrate chemokine superfamily.

Footprints of natural and artificial selection for photoperiod pathway genes in Oryza.

Asian rice, Oryza sativa, consists of two major subspecies, indica and japonica, which are physiologically differentiated and adapted to different latitudes. Genes for photoperiod sensitivity are likely targets of selection along latitude. We examined the footprints of natural and artificial selections for four major genes of the photoperiod pathway, namely PHYTOCHROME B (PhyB), HEADING DATE 1 (Hd1), HEADING DATE 3a (Hd3a), and EARLY HEADING DATE 1 (Ehd1), by investigation of the patterns of nucleotide polymorphisms in cultivated and wild rice. Geographical subdivision between tropical and subtropical O. rufipogon was found for all of the photoperiod genes in plants divided by the Tropic of Cancer (TOC). All of these genes, except for PhyB, were characterized by the existence of clades that split a long time ago and that corresponded to latitudinal subdivisions, and revealed a likely diversifying selection. Ssp. indica showed close affinity to tropical O. rufipogon for all genes, while ssp. japonica, which has a much wider range of distribution, displayed complex patterns of differentiation from O. rufipogon, which reflected various agricultural needs in relation to crop yield. In japonica, all genes, except Hd3a, were genetically differentiated at the TOC, while geographical subdivision occurred at 31°N in Hd3a, probably the result of varying photoperiods. Many other features of the photoperiod genes revealed domestication signatures, which included high linkage disequilibrium (LD) within genes, the occurrence of frequent and recurrent non-functional Hd1 mutants in cultivated rice, crossovers between subtropical and tropical alleles of Hd1, and significant LD between Hd1 and Hd3a in japonica and indica.

Mitochondrial-nuclear interactions and accelerated compensatory evolution: evidence from the primate cytochrome C oxidase complex.

Accelerated rates of mitochondrial protein evolution have been proposed to reflect Darwinian coadaptation for efficient energy production for mammalian flight and brain activity. However, several features of mammalian mtDNA (absence of recombination, small effective population size, and high mutation rate) promote genome degradation through the accumulation of weakly deleterious mutations. Here, we present evidence for "compensatory" adaptive substitutions in nuclear DNA- (nDNA) encoded mitochondrial proteins to prevent fitness decline in primate mitochondrial protein complexes. We show that high mutation rate and small effective population size, key features of primate mitochondrial genomes, can accelerate compensatory adaptive evolution in nDNA-encoded genes. We combine phylogenetic information and the 3D structure of the cytochrome c oxidase (COX) complex to test for accelerated compensatory changes among interacting sites. Physical interactions among mtDNA- and nDNA-encoded components are critical in COX evolution; amino acids in close physical proximity in the 3D structure show a strong tendency for correlated evolution among lineages. Only nuclear-encoded components of COX show evidence for positive selection and adaptive nDNA-encoded changes tend to follow mtDNA-encoded amino acid changes at nearby sites in the 3D structure. This bias in the temporal order of substitutions supports compensatory weak selection as a major factor in accelerated primate COX evolution.

Genealogical characterization of regional populations and dorsal coat color variation in the house mouse Mus musculus from Asia based on haplotype structure analysis of a gene-rich region harboring Mc1r.

We analyzed 196 haplotype sequences from a gene-rich region (250 kb) that includes Mc1r, a gene involved in coat color regulation, to gain insight into the evolution of coat color variation in subspecies of the house mouse Mus musculus. Phylogenetic networks revealed haplotype groups from the major subspecies of M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS). Using haplotype sequences assigned to each of CAS and MUS through phylogenetic analysis, we proposed migration routes associated with prehistoric humans from west to east across Eurasia. Comparing nucleotide diversity among subspecies-specific haplotypes in different geographic areas showed a marked reduction during migration, particularly in MUS-derived haplotypes from Korea and Japan, suggesting intensive population bottlenecks during migration. We found that a C>T polymorphism at site 302 (c.302C>T) in the Mc1r coding region correlated with a darkening of dorsal fur color in both CAS and MUS. However, C/C homozygous mice in MUS showed marked variation in lightness, indicating the possibility of another genetic determinant that affects the lightness of dorsal fur color. Detailed sequence comparisons of haplotypes revealed that short fragments assigned to DOM were embedded in CAS-assigned fragments, indicating ancient introgression between subspecies. The estimated age of c.302C>T also supports the hypothesis that genetic interaction between subspecies occurred in ancient times. This suggests that the genome of M. musculus evolved through gene flow between subspecies over an extended period before the movement of the species in conjunction with prehistoric humans.

c-FOS is an integral component of the IKZF1 transactivator complex and mediates lenalidomide resistance in multiple myeloma.

BACKGROUND: The immunomodulatory drug lenalidomide, which is now widely used for the treatment of multiple myeloma (MM), exerts pharmacological action through the ubiquitin-dependent degradation of IKZF1 and subsequent down-regulation of interferon regulatory factor 4 (IRF4), a critical factor for the survival of MM cells. IKZF1 acts principally as a tumour suppressor via transcriptional repression of oncogenes in normal lymphoid lineages. In contrast, IKZF1 activates IRF4 and other oncogenes in MM cells, suggesting the involvement of unknown co-factors in switching the IKZF1 complex from a transcriptional repressor to an activator. The transactivating components of the IKZF1 complex might promote lenalidomide resistance by residing on regulatory regions of the IRF4 gene to maintain its transcription after IKZF1 degradation. METHODS: To identify unknown components of the IKZF1 complex, we analyzed the genome-wide binding of IKZF1 in MM cells using chromatin immunoprecipitation-sequencing (ChIP-seq) and screened for the co-occupancy of IKZF1 with other DNA-binding factors on the myeloma genome using the ChIP-Atlas platform. RESULTS: We found that c-FOS, a member of the activator protein-1 (AP-1) family, is an integral component of the IKZF1 complex and is primarily responsible for the activator function of the complex in MM cells. The genome-wide screening revealed the co-occupancy of c-FOS with IKZF1 on the regulatory regions of IKZF1-target genes, including IRF4 and SLAMF7, in MM cells but not normal bone marrow progenitors, pre-B cells or mature T-lymphocytes. c-FOS and IKZF1 bound to the same consensus sequence as the IKZF1 complex through direct protein-protein interactions. The complex also includes c-JUN and IKZF3 but not IRF4. Treatment of MM cells with short-hairpin RNA against FOS or a selective AP-1 inhibitor significantly enhanced the anti-MM activity of lenalidomide in vitro and in two murine MM models. Furthermore, an AP-1 inhibitor mitigated the lenalidomide resistance of MM cells. CONCLUSIONS: C-FOS determines lenalidomide sensitivity and mediates drug resistance in MM cells as a co-factor of IKZF1 and thus, could be a novel therapeutic target for further improvement of the prognosis of MM patients.