Derivation of Pluripotent Stem Cells with In Vivo Embryonic and Extraembryonic Potency.

Of all known cultured stem cell types, pluripotent stem cells (PSCs) sit atop the landscape of developmental potency and are characterized by their ability to generate all cell types of an adult organism. However, PSCs show limited contribution to the extraembryonic placental tissues in vivo. Here, we show that a chemical cocktail enables the derivation of stem cells with unique functional and molecular features from mice and humans, designated as extended pluripotent stem (EPS) cells, which are capable of chimerizing both embryonic and extraembryonic tissues. Notably, a single mouse EPS cell shows widespread chimeric contribution to both embryonic and extraembryonic lineages in vivo and permits generating single-EPS-cell-derived mice by tetraploid complementation. Furthermore, human EPS cells exhibit interspecies chimeric competency in mouse conceptuses. Our findings constitute a first step toward capturing pluripotent stem cells with extraembryonic developmental potentials in culture and open new avenues for basic and translational research. VIDEO ABSTRACT.

Transcription factor Hoxb5 reprograms B cells into functional T lymphocytes.

Deletion of master regulators of the B cell lineage reprograms B cells into T cells. Here we found that the transcription factor Hoxb5, which is expressed in uncommitted hematopoietic progenitor cells but is not present in cells committed to the B cell or T cell lineage, was able to reprogram pro-pre-B cells into functional early T cell lineage progenitors. This reprogramming started in the bone marrow and was completed in the thymus and gave rise to T lymphocytes with transcriptomes, hierarchical differentiation, tissue distribution and immunological functions that closely resembled those of their natural counterparts. Hoxb5 repressed B cell 'master genes', activated regulators of T cells and regulated crucial chromatin modifiers in pro-pre-B cells and ultimately drove the B cell fate-to-T cell fate conversion. Our results provide a de novo paradigm for the generation of functional T cells through reprogramming in vivo.

Publisher Correction: Transcription factor Hoxb5 reprograms B cells into functional T lymphocytes.

In the version of this article initially published, some identification of the supplementary information was incorrect. The items originally called Supplementary Tables 1, 2, 3, 4 and 5 should be Source Data Figures 1, 2, 4, 5 and 7, respectively; those originally called Supplementary Tables 6, 7 and 8 should be Supplementary Tables 1, 2 and 3, respectively; and those originally called Source Data Figures 1, 2, 4, 5 and 7 should be Supplementary Tables 4, 5, 6, 7 and 8, respectively. The errors have been corrected in the HTML version of the article.

Mitochondrial base editor induces substantial nuclear off-target mutations.

DddA-derived cytosine base editors (DdCBEs)-which are fusions of split DddA halves and transcription activator-like effector (TALE) array proteins from bacteria-enable targeted C•G-to-T•A conversions in mitochondrial DNA1. However, their genome-wide specificity is poorly understood. Here we show that the mitochondrial base editor induces extensive off-target editing in the nuclear genome. Genome-wide, unbiased analysis of its editome reveals hundreds of off-target sites that are TALE array sequence (TAS)-dependent or TAS-independent. TAS-dependent off-target sites in the nuclear DNA are often specified by only one of the two TALE repeats, challenging the principle that DdCBEs are guided by paired TALE proteins positioned in close proximity. TAS-independent off-target sites on nuclear DNA are frequently shared among DdCBEs with distinct TALE arrays. Notably, they co-localize strongly with binding sites for the transcription factor CTCF and are enriched in topologically associating domain boundaries. We engineered DdCBE to alleviate such off-target effects. Collectively, our results have implications for the use of DdCBEs in basic research and therapeutic applications, and suggest the need to thoroughly define and evaluate the off-target effects of base-editing tools.

Polycomb Group Proteins Regulate Chromatin Architecture in Mouse Oocytes and Early Embryos.

In mammals, chromatin organization undergoes drastic reorganization during oocyte development. However, the dynamics of three-dimensional chromatin structure in this process is poorly characterized. Using low-input Hi-C (genome-wide chromatin conformation capture), we found that a unique chromatin organization gradually appears during mouse oocyte growth. Oocytes at late stages show self-interacting, cohesin-independent compartmental domains marked by H3K27me3, therefore termed Polycomb-associating domains (PADs). PADs and inter-PAD (iPAD) regions form compartment-like structures with strong inter-domain interactions among nearby PADs. PADs disassemble upon meiotic resumption from diplotene arrest but briefly reappear on the maternal genome after fertilization. Upon maternal depletion of Eed, PADs are largely intact in oocytes, but their reestablishment after fertilization is compromised. By contrast, depletion of Polycomb repressive complex 1 (PRC1) proteins attenuates PADs in oocytes, which is associated with substantial gene de-repression in PADs. These data reveal a critical role of Polycomb in regulating chromatin architecture during mammalian oocyte growth and early development.

Single-cell EpiChem jointly measures drug-chromatin binding and multimodal epigenome.

Studies of molecular and cellular functions of small-molecule inhibitors in cancer treatment, eliciting effects by targeting genome and epigenome associated proteins, requires measurement of drug-target engagement in single-cell resolution. Here we present EpiChem for in situ single-cell joint mapping of small molecules and multimodal epigenomic landscape. We demonstrate single-cell co-assays of three small molecules together with histone modifications, chromatin accessibility or target proteins in human colorectal cancer (CRC) organoids. Integrated multimodal analysis reveals diverse drug interactions in the context of chromatin states within heterogeneous CRC organoids. We further reveal drug genomic binding dynamics and adaptive epigenome across cell types after small-molecule drug treatment in CRC organoids. This method provides a unique tool to exploit the mechanisms of cell type-specific drug actions.

CoBATCH for High-Throughput Single-Cell Epigenomic Profiling.

An efficient, generalizable method for genome-wide mapping of single-cell histone modifications or chromatin-binding proteins is lacking. Here, we develop CoBATCH, combinatorial barcoding and targeted chromatin release, for single-cell profiling of genomic distribution of chromatin-binding proteins in cell culture and tissue. Protein A in fusion to Tn5 transposase is enriched through specific antibodies to genomic regions, and Tn5 generates indexed chromatin fragments ready for library preparation and sequencing. Importantly, this strategy enables not only low-input epigenomic profiling in intact tissues but also measures scalable up to tens of thousands of single cells per experiment under both native and cross-linked conditions. CoBATCH produces ∼12,000 reads/cell with extremely low background. Mapping of endothelial cell lineages from ten embryonic mouse organs through CoBATCH allows for efficient deciphering of epigenetic heterogeneity of cell populations and cis-regulatory mechanisms. Thus, obviating specialized devices, CoBATCH is broadly applicable and easily deployable for single-cell profiling of protein-DNA interactions.

Single-cell joint detection of chromatin occupancy and transcriptome enables higher-dimensional epigenomic reconstructions.

Deciphering mechanisms in cell-fate decisions requires single-cell holistic reconstructions of multidimensional epigenomic states in transcriptional regulation. Here we develop CoTECH, a combinatorial barcoding method allowing high-throughput single-cell joint detection of chromatin occupancy and transcriptome. We used CoTECH to examine bivalent histone marks (H3K4me3 and H3K27me3) with transcription from naive to primed mouse embryonic stem cells. We also derived concurrent bivalent marks in pseudosingle cells using transcriptome as an anchor for resolving pseudotemporal bivalency trajectories and disentangling a context-specific interplay between H3K4me3/H3K27me3 and transcription level. Next, we revealed the regulatory basis of endothelial-to-hematopoietic transition in two waves of hematopoietic cells and distinctive enhancer-gene-linking schemes guiding hemogenic endothelial cell emergence, indicating a unique epigenetic control of transcriptional regulation for hematopoietic stem cell priming. CoTECH provides an efficient framework for single-cell coassay of chromatin occupancy and transcription, thus enabling higher-dimensional epigenomic reconstructions.

Loss of sphingosine kinase 2 promotes the expansion of hematopoietic stem cells by improving their metabolic fitness.

Hematopoietic stem cells (HSCs) have reduced capacities to properly maintain and replenish the hematopoietic system during myelosuppressive injury or aging. Expanding and rejuvenating HSCs for therapeutic purposes has been a long-sought goal with limited progress. Here, we show that the enzyme Sphk2 (sphingosine kinase 2), which generates the lipid metabolite sphingosine-1-phosphate, is highly expressed in HSCs. The deletion of Sphk2 markedly promotes self-renewal and increases the regenerative potential of HSCs. More importantly, Sphk2 deletion globally preserves the young HSC gene expression pattern, improves the function, and sustains the multilineage potential of HSCs during aging. Mechanistically, Sphk2 interacts with prolyl hydroxylase 2 and the Von Hippel-Lindau protein to facilitate HIF1α ubiquitination in the nucleus independent of the Sphk2 catalytic activity. Deletion of Sphk2 increases hypoxic responses by stabilizing the HIF1α protein to upregulate PDK3, a glycolysis checkpoint protein for HSC quiescence, which subsequently enhances the function of HSCs by improving their metabolic fitness; specifically, it enhances anaerobic glycolysis but suppresses mitochondrial oxidative phosphorylation and generation of reactive oxygen species. Overall, targeting Sphk2 to enhance the metabolic fitness of HSCs is a promising strategy to expand and rejuvenate functional HSCs.

GATA4 Regulates Developing Endocardium Through Interaction With ETS1.

BACKGROUND: The pioneer transcription factor (TF) GATA4 (GATA Binding Protein 4) is expressed in multiple cardiovascular lineages and is essential for heart development. GATA4 lineage-specific occupancy in the developing heart underlies its lineage specific activities. Here, we characterized GATA4 chromatin occupancy in cardiomyocyte and endocardial lineages, dissected mechanisms that control lineage specific occupancy, and analyzed GATA4 regulation of endocardial gene expression. METHODS: We mapped GATA4 chromatin occupancy in cardiomyocyte and endocardial cells of embryonic day 12.5 (E12.5) mouse heart using lineage specific, Cre-activated biotinylation of GATA4. Regulation of GATA4 pioneering activity was studied in cell lines stably overexpressing GATA4. GATA4 regulation of endocardial gene expression was analyzed using single cell RNA sequencing and luciferase reporter assays. RESULTS: Cardiomyocyte-selective and endothelial-selective GATA4 occupied genomic regions had features of lineage specific enhancers. Footprints within cardiomyocyte- and endothelial-selective GATA4 regions were enriched for NKX2-5 (NK2 homeobox 5) and ETS1 (ETS Proto-Oncogene 1) motifs, respectively, and both of these TFs interacted with GATA4 in co-immunoprecipitation assays. In stable NIH3T3 cell lines expressing GATA4 with or without NKX2-5 or ETS1, the partner TFs re-directed GATA4 pioneer binding and augmented its ability to open previously inaccessible regions, with ETS1 displaying greater potency as a pioneer partner than NKX2-5. Single-cell RNA sequencing of embryonic hearts with endothelial cell-specific Gata4 inactivation identified Gata4-regulated endocardial genes, which were adjacent to GATA4-bound, endothelial regions enriched for both GATA4 and ETS1 motifs. In reporter assays, GATA4 and ETS1 cooperatively stimulated endothelial cell enhancer activity. CONCLUSIONS: Lineage selective non-pioneer TFs NKX2-5 and ETS1 guide the activity of pioneer TF GATA4 to bind and open chromatin and create active enhancers and mechanistically link ETS1 interaction to GATA4 regulation of endocardial development.

Generation and validation of versatile inducible CRISPRi embryonic stem cell and mouse model.

Clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR-associated (Cas) 9 has been widely used far beyond genome editing. Fusions of deactivated Cas9 (dCas9) to transcription effectors enable interrogation of the epigenome and controlling of gene expression. However, the large transgene size of dCas9-fusion hinders its applications especially in somatic tissues. Here, we develop a robust CRISPR interference (CRISPRi) system by transgenic expression of doxycycline (Dox) inducible dCas9-KRAB in mouse embryonic stem cells (iKRAB ESC). After introduction of specific single-guide RNAs (sgRNAs), the induced dCas9-KRAB efficiently maintains gene inactivation, although it modestly down-regulates the expression of active genes. The proper timing of Dox addition during cell differentiation or reprogramming allows us to study or screen spatiotemporally activated promoters or enhancers and thereby the gene functions. Furthermore, taking the ESC for blastocyst injection, we generate an iKRAB knock-in (KI) mouse model that enables the shutdown of gene expression and loss-of-function (LOF) studies ex vivo and in vivo by a simple transduction of gRNAs. Thus, our inducible CRISPRi ESC line and KI mouse provide versatile and convenient platforms for functional interrogation and high-throughput screens of specific genes and potential regulatory elements in the setting of development or diseases.

The lncRNA Hand2os1/Uph locus orchestrates heart development through regulation of precise expression of Hand2.

Exploration and dissection of potential actions and effects of long noncoding RNA (lncRNA) in animals remain challenging. Here, using multiple knockout mouse models and single cell RNA sequencing, we demonstrate that the divergent lncRNA Hand2os1/Uph has a key complex modulatory effect on the expression of its neighboring gene HAND2 and subsequently on heart development and function. Short deletion of the Hand2os1 promoter in mouse diminishes Hand2os1 transcription to ∼8-32%, but fails to affect HAND2 expression and yields no discernable heart phenotypes. Interestingly, full-length deletion of Hand2os1 in mouse causes moderate yet prevalent upregulation of HAND2 in hundreds of cardiac cells, leading to profound biological consequences, including dysregulated cardiac gene programs, congenital heart defects and perinatal lethality. We propose that the Hand2os1 locus dampens HAND2 expression to restrain cardiomyocyte proliferation, thereby orchestrating a balanced development of cardiac cell lineages. This study highlights the regulatory complexity of the lncRNA Hand2os1 on HAND2 expression, emphasizing the need for complementary genetic and single cell approaches to delineate the function and primary molecular effects of an lncRNA in animals.

Coupling of COPII vesicle trafficking to nutrient availability by the IRE1α-XBP1s axis.

The cytoplasmic coat protein complex-II (COPII) is evolutionarily conserved machinery that is essential for efficient trafficking of protein and lipid cargos. How the COPII machinery is regulated to meet the metabolic demand in response to alterations of the nutritional state remains largely unexplored, however. Here, we show that dynamic changes of COPII vesicle trafficking parallel the activation of transcription factor X-box binding protein 1 (XBP1s), a critical transcription factor in handling cellular endoplasmic reticulum (ER) stress in both live cells and mouse livers upon physiological fluctuations of nutrient availability. Using live-cell imaging approaches, we demonstrate that XBP1s is sufficient to promote COPII-dependent trafficking, mediating the nutrient stimulatory effects. Chromatin immunoprecipitation (ChIP) coupled with high-throughput DNA sequencing (ChIP-seq) and RNA-sequencing analyses reveal that nutritional signals induce dynamic XBP1s occupancy of promoters of COPII traffic-related genes, thereby driving the COPII-mediated trafficking process. Liver-specific disruption of the inositol-requiring enzyme 1α (IRE1α)-XBP1s signaling branch results in diminished COPII vesicle trafficking. Reactivation of XBP1s in mice lacking hepatic IRE1α restores COPII-mediated lipoprotein secretion and reverses the fatty liver and hypolipidemia phenotypes. Thus, our results demonstrate a previously unappreciated mechanism in the metabolic control of liver protein and lipid trafficking: The IRE1α-XBP1s axis functions as a nutrient-sensing regulatory nexus that integrates nutritional states and the COPII vesicle trafficking.

Replication-Independent Histone Turnover Underlines the Epigenetic Homeostasis in Adult Heart.

RATIONALE: Replication-independent histone turnover has been linked to cis-regulatory chromatin domains in cultured cell lines, but its molecular underpinnings and functional relevance in adult mammalian tissues remain yet to be defined. OBJECTIVE: We investigated regulatory functions of replication-independent histone turnover in chromatin states of postmitotic cardiomyocytes from adult mouse heart. METHODS AND RESULTS: We used H2B-GFP (histone 2B-green fluorescent protein) fusion protein pulse-and-chase approaches to measure histone turnover rate in mouse cardiomyocytes. Surprisingly, we found that the short histone half-life (≈2 weeks) contrasted dramatically with the long lifetime of cardiomyocytes, and rapid histone turnover regions corresponded to cis-regulatory domains of heart genes. Interestingly, recruitment of chromatin modifiers, including Polycomb EED (embryonic ectoderm development), was positively correlated with histone turnover rate at enhancers. Mechanistically, through directly interacting with and engaging the BAF (BRG1 [Brahma-related gene-1]-associated factor) complex for nucleosome exchange for stereotyped histone modifications from the free histone pool, EED augmented histone turnover to restrain enhancer overactivation. CONCLUSIONS: We propose a model in which replication-independent histone turnover reinforces robustness of local chromatin states for adult tissue homeostasis.

Single-Cell Transcriptomics Reveals Chemotaxis-Mediated Intraorgan Crosstalk During Cardiogenesis.

RATIONALE: We hypothesized that the differentiation processes of cardiac progenitor cell (CP) from first and second heart fields (FHF and SHF) may undergo the unique instructive gene regulatory networks or signaling pathways, and the precise SHF progression is contingent on the FHF signaling developmental cues. OBJECTIVE: We investigated how the intraorgan communications control sequential building of discrete anatomic regions of the heart at single-cell resolution. METHODS AND RESULTS: By single-cell transcriptomic analysis of Nkx2-5 (NK2 homeobox 5) and Isl1 (ISL LIM homeobox 1) lineages at embryonic day 7.75, embryonic day 8.25, embryonic day 8.75, and embryonic day 9.25, we present a panoramic view of distinct CP differentiation hierarchies. Computational identifications of FHF- and SHF-CP descendants revealed that SHF differentiation toward cardiomyocytes underwent numerous step-like transitions, whereas earlier FHF progressed toward cardiomyocytes in a wave-like manner. Importantly, single-cell pairing analysis demonstrated that SHF-CPs were attracted to and expanded FHF-populated heart tube region through interlineage communications mediated by the chemotactic guidance (MIF [macrophage migration inhibitory factor]-CXCR2 [C-X-C motif chemokine receptor 2]). This finding was verified by pharmacological blockade of this chemotaxis in embryos manifesting limited SHF cell migration and contribution to the growth of the outflow tract and right ventricle but undetectable effects on the left ventricle or heart tube initiation. Genetic loss-of-function assay of Cxcr2 showed that the expression domain of CXCR4 was expanded predominantly at SHF. Furthermore, double knockout of Cxcr2/Cxcr4 exhibited defective SHF development, corroborating the redundant function. Mechanistically, NKX2-5 directly bound the Cxcr2 and Cxcr4 genomic loci and activated their transcription in SHF. CONCLUSIONS: Collectively, we propose a model in which the chemotaxis-mediated intraorgan crosstalk spatiotemporally guides the successive process of positioning SHF-CP and promoting primary heart expansion and patterning upon FHF-derived heart tube initiation.

Human exonization through differential nucleosome occupancy.

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, but the roles of nucleosome occupancy in shaping changes through evolution remain to be addressed. Here we present high-resolution nucleosome occupancy profiles for multiple tissues derived from human, macaque, tree shrew, mouse, and pig. Genome-wide comparison reveals conserved nucleosome occupancy profiles across both different species and tissue types. Notably, we found significantly higher levels of nucleosome occupancy in exons than in introns, a pattern correlated with the different exon-intron GC content. We then determined whether this biased occupancy may play roles in the origination of new exons through evolution, rather than being a downstream effect of exonization, through a comparative approach to sequentially trace the order of the exonization and biased nucleosome binding. By identifying recently evolved exons in human but not in macaque using matched RNA sequencing, we found that higher exonic nucleosome occupancy also existed in macaque regions orthologous to these exons. Presumably, such biased nucleosome occupancy facilitates the origination of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content. These data thus support a model that sites bound by nucleosomes are more likely to evolve into exons, which we term the "nucleosome-first" model.

Single-Cell Transcriptomic Analysis of Cardiac Progenitor Differentiation.

PURPOSE OF REVIEW: Emerging single-cell RNA sequencing technologies hold great promises to boost our understanding of the heterogeneity and molecular regulation of diverse cell phenotypes during organ development. In this review, we aimed at summarizing recent advances in employing single-cell transcriptomic analysis to depict the landscape of embryonic heart development, in particular, focusing on cardiac progenitor (CP) differentiation. RECENT FINDINGS: Recent studies unbiasedly cataloged and characterized cardiac cell types in the spatial and temporal resolution during early heart development. Pseudo-time analysis revealed a temporal continuum of the differentiation progress from embryonic day (E) 6.5 to E9.5, implicating early cardiac lineage restriction during mouse gastrulation. First and second heart field (FHF and SHF) CPs adopted different differentiation strategies and underwent distinct transcriptional regulation. Collectively, the comprehensive molecular atlases yield a rich resource for identification of the key cardiac regulators and signaling molecules within the key cardiac gene regulatory network (GRN) governing cardiac cell fate determinations. This review offers insights into the exquisite process and its regulation of CP differentiation at single-cell resolution. As single-cell technologies continuously grow and evolve, computational integration of multimodal single-cell data with well-designed experimental validation promises to further delineate molecular basis in deploying cardiac progenitors of distinct sources with anatomical information.

PRC2 directly methylates GATA4 and represses its transcriptional activity.

Polycomb-repressive complex 2 (PRC2) promotes tissue-specific differentiation by depositing trimethylated histone H3 Lys 27 (H3K27me3) epigenetic marks to silence ectopic gene expression programs. Here, we show that EZH2, the catalytic subunit of PRC2, is required for cardiac morphogenesis. Both in vitro and in fetal hearts, EZH2 interacted with cardiac transcription factor GATA4 and directly methylated it at Lys 299. PRC2 methylation of GATA4 attenuated its transcriptional activity by reducing its interaction with and acetylation by p300. Our results reveal a new mechanism of PRC2-mediated transcriptional repression in which PRC2 methylates a transcription factor to inhibit its transcriptional activity.

Divergent Requirements for EZH1 in Heart Development Versus Regeneration.

RATIONALE: Polycomb repressive complex 2 is a major epigenetic repressor that deposits methylation on histone H3 on lysine 27 (H3K27me) and controls differentiation and function of many cells, including cardiac myocytes. EZH1 and EZH2 are 2 alternative catalytic subunits with partial functional redundancy. The relative roles of EZH1 and EZH2 in heart development and regeneration are unknown. OBJECTIVE: We compared the roles of EZH1 versus EZH2 in heart development and neonatal heart regeneration. METHODS AND RESULTS: Heart development was normal in Ezh1-/- (Ezh1 knockout) and Ezh2f/f::cTNT-Cre (Ezh2 knockout) embryos. Ablation of both genes in Ezh1-/-::Ezh2f/f::cTNT-Cre embryos caused lethal heart malformations, including hypertrabeculation, compact myocardial hypoplasia, and ventricular septal defect. Epigenome and transcriptome profiling showed that derepressed genes were upregulated in a manner consistent with total EZH dose. In neonatal heart regeneration, Ezh1 was required, but Ezh2 was dispensable. This finding was further supported by rescue experiments: cardiac myocyte-restricted re-expression of EZH1 but not EZH2 restored neonatal heart regeneration in Ezh1 knockout. In myocardial infarction performed outside of the neonatal regenerative window, EZH1 but not EZH2 likewise improved heart function and stimulated cardiac myocyte proliferation. Mechanistically, EZH1 occupied and activated genes related to cardiac growth. CONCLUSIONS: Our work unravels divergent mechanisms of EZH1 in heart development and regeneration, which will empower efforts to overcome epigenetic barriers to heart regeneration.

Isoform Evolution in Primates through Independent Combination of Alternative RNA Processing Events.

Recent RNA-seq technology revealed thousands of splicing events that are under rapid evolution in primates, whereas the reliability of these events, as well as their combination on the isoform level, have not been adequately addressed due to its limited sequencing length. Here, we performed comparative transcriptome analyses in human and rhesus macaque cerebellum using single molecule long-read sequencing (Iso-seq) and matched RNA-seq. Besides 359 million RNA-seq reads, 4,165,527 Iso-seq reads were generated with a mean length of 14,875 bp, covering 11,466 human genes, and 10,159 macaque genes. With Iso-seq data, we substantially expanded the repertoire of alternative RNA processing events in primates, and found that intron retention and alternative polyadenylation are surprisingly more prevalent in primates than previously estimated. We then investigated the combinatorial mode of these alternative events at the whole-transcript level, and found that the combination of these events is largely independent along the transcript, leading to thousands of novel isoforms missed by current annotations. Notably, these novel isoforms are selectively constrained in general, and 1,119 isoforms have even higher expression than the previously annotated major isoforms in human, indicating that the complexity of the human transcriptome is still significantly underestimated. Comparative transcriptome analysis further revealed 502 genes encoding selectively constrained, lineage-specific isoforms in human but not in rhesus macaque, linking them to some lineage-specific functions. Overall, we propose that the independent combination of alternative RNA processing events has contributed to complex isoform evolution in primates, which provides a new foundation for the study of phenotypic difference among primates.