Sarkar et al.

SUPPLEMENTAL INFORMATION

Figure S1

A

10 20 30 40 50 MADKRKLQGE IDRCLKKVSE GVEQFEDIWQ KLHNAANANQ KEKYEADLKK 60 70 80 90 100 EIKKLQRLRD QIKTWVASNE IKDKRQLIEN RKLIETQMER FKVVERETKT Predicted 110 120 130 140 150 coiled-coil KAYSKEGLGL AQKVDPAQKE KEEVGQWLTN TIDTLNMQVD QFESEVESLS domain 160 170 180 190 200 (1-241) VQTRKKKGDK DKQDRIEGLK RHIEKHRYHV RMLETILRML DNDSILVDAI 210 220 230 240 250 RKIKDDVEYY VDSSQDPDFE ENEFLYDDLD LEDIPQALVA TSPPSHSHME 260 270 280 290 300 DEIFNQSSST PTSTTSSSPI PPSPANCTTE NSEDDKKRGR ST* DSEVSQSP * 310 320 330 340 350 AKNGSKPVHS NQHPQSPAVP PTYPSGPPPT TSALSSTPGN NGASTPAAPT 360 370 380 390 400 SALGPKASPA PSHNSGTPAP YAQAVAPPNA SGPSNAQPRP PSAQPSGGSG 410 420 430 440 450 Intrinsically GGSGGSSSNS NSGTGGGAGK QNGATSYSSV VADSPAEVTL SSSGGSSASS disordered region 460 470 480 490 500 (242-605) QALGPTSGPH NPAPSTSKES STAAPSGAGN VASGSGNNSG GPSLLVPLPV 510 520 530 540 550 NPPSSPTPSF SEAKAAGTLL NGPPQFSTTP EIKAPEPLSS LKSMAERAAI 560 570 580 590 600 SSGIEDPVPT LHLTDRDIIL SSTSAPPTSS QPPLQLSEVN IPLSLGVCPL 610 620 630 640 650 GPVSLTKEQL YQQAMEEAAW HHMPHPSDSE RIRQYLPRNP CPTPPYHHQM 660 670 680 690 700 NAR/CS PPPHSDTVEF YQRLSTETLF FIFYYLEGTK AQYLAAKALK KQSWRFHTKY 710 720 730 740 750 NOT box (654-751) MMWFQRHEEP KTITDEFEQG TYIYFDYEKW GQRKKEGFTF EYRYLEDRDLQ

B

2 Supplementary Figure 1: A. Sequence of the mouse CNOT3 (UNIPROT Q8K0V4). Vertical blue lines delineate the domains of the protein that are shown in Figure 1B. The green-shaded box indicates the the NOT box region, which is required for the interaction between CNOT3 and Aurora B (Figure 1C). Aurora B consensus phosphorylation sites are underlined. The red box and red text indicate the putative nuclear localization sequence. * indicates residues 292 and 294, which were mutated in this study. B. Phosphorylation sites detected on the CNOT3 protein by in vivo proteomic discovery mass pectrometry. Data was obtained from PhosphoSite® (Hornbeck et al., 2015). The data shown is for the human CNOT3 protein, which is 95% identical to mouse CNOT3. With the exception of one residue (E79), all of the residues that vary between human and mouse are located within the intrinsically disordered region between residues 329 and 604 of the mouse sequence. Blue circles indicate phosphorylated residues. The y- axis shows the number of records where the modification was assigned using proteomic discovery mass spectrometry. The numbered residues show a phosphorylation hotspot that extends from positions 291-299.

3 Figure S2

A

B WT

Clone 1

110bp ssODN

Clone 2

Cnot3 Exon 10

C

WT-DMSO WT-AZD1152 WT-AZD1152

24h 48h

Counts

Counts Counts

PI PI PI

D E G2/M CNOT3-DM S * ** NS G1 NS

** ** Counts * NS NS Percent of cells

PI

WT-DMSO F CNOT3-DM WT-AZD1152WT-AZD1152 24h 48h

V5-CNOT3-WTV5-CNOT3-T292AV5-CNOT3-S294AV5-CNOT3-DMV5-CNOT3-WTV5-CNOT3-T292AV5-CNOT3-S294AV5-CNOT3-DM

V5 (CNOT3) V5 (CNOT3)

-TUBULIN LAMIN B

Cytoplasm Nucleus HEK 293 cells

4

Supplementary Figure 2: Identification of Aurora B phosphorylation sites in the CNOT3 protein A. MS/MS fragmentation spectrum generated by ion trap collision-induced dissociation (CID) of doubly charged precursor ion with m/z 658.27 identifying CNOT3 tryptic peptide 291-302 (STDSEVSQSPAK). The spectrum presented two fragments which were as follows: Left panel: query number 2609 with assigned b- and y-ions; MASCOT score: -5 54; expectation value: 1.8x10 ; neutral loss of H3PO4 (-97.98 Da) from precursor ion detected at m/z 609.49; localisation probability for phosphorylation of Ser291: 3.0%, Thr292: 3.0%, Ser294 is 93.9%. Right panel: query number 2608 with assigned b- and y-ions; MASCOT score: 53; expectation value: 3.6x10-5; localisation probability for phosphorylation of Ser291: 25.6%, Thr292: 25.6%, Ser294: 48.4%. B. Left panel: Schematic representation showing use of CRISPR/Cas9 targeting of Cnot3 in exon 10 to mutate Threonine 292 and Serine 294 to alanine, using the guide RNA described in Methods and a 110 bp single-stranded donor oligonucleotides (ssODN) carrying the mutations. Right panel: Sequence images showing the wild-type (WT) sequence and the mutated sequences obtained from two CRISPR generated clones. C. FACS analysis of the cell cycle profiles of ESCs incubated for 24 and 48 hours with the Aurora B inhibitor AZD1152 or with vehicle (DMSO) and stained with propidium iodide (PI). D. FACS analysis of the cell cycle profile of Cnot3-DM ESCs. E. Histograms showing the percentage of cells measured at different phases of the cell cycle in (c) and (d). Mean ± SEM; P values (with respect to WT-DMSO) calculated by unpaired t-test, ns = non-significant, n=3. F. Representative immunoblot analysis of V5-tagged CNOT3 carried out on cytoplasmic and nuclear extracts from HEK 293 cells transfected with WT Cnot3, Cnot3-T292A, Cnot3-S294A, or Cnot3-DM expression constructs. Cells were harvested 48 hours post- transfection. a-Tubulin and Lamin B were used as loading controls for the cytoplasm and nucleus respectively.

5 Figure S3

A B

DAPI OTX2 WT CNOT3-DM WT

*** *** SURVIVAL RATIO SURVIVAL *

BMP4+FGF2 + - - CNOT3-DM ACTIVIN A - + - WNT3 - - +

C

Mesoderm Endoderm Late endoderm and Ectoderm Pluripotency trophectoderm *** ns *** *** * *** * *

**

** *** (RELATIVE TO UNDIFF) TO (RELATIVE FOLD CHANGE IN mRNA

(RELATIVE TO UNDIFF) TO (RELATIVE *** *** *** FOLD CHANGE IN mRNA (RELATIVE TO UNDIFF) TO (RELATIVE FOLD CHANGE IN mRNA

Cdx2 Fgf5 Mixl1Gata4Gata6Sox17 PdgfrCxcr4 Tead4 Sox1 Oct-4 Nanog Sox2 Brachury Lineage marker expression in differentiated wild-type (WT) mesendodermal cells 4 days mesendoderm differentiation

Supplementary Figure 3: Differentiation and survival of Cnot3-DM mesendoderm A. Representative confocal images of 8-day embryoid bodies (EBs) derived from wild- type (WT) and Cnot3-DM ESCs and subjected to immunostaining for the ectodermal marker OTX2 (green). Nuclei were counterstained with DAPI; scale bar: 100 µm. B. Survival of WT and Cnot3-DM ESCs over 4 days of differentiation induced in the presence of BMP4+FGF2 or Activin A or Wnt3 and the cell survival was determined using WST-1 reagent on day 4 of the differentiation. Histograms show survival ratios calculated relative to the values obtained for wild-type (WT) cells treated with BMP4+FGF2, which was assigned a value of 1. Mean ± SEM, significance calculated using unpaired t-test *P<0.05, ***P<0.001, n=3. C. Analysis of mRNA levels of lineage marker after BMP4+FGF2 induced differentiation of WT ESCs for 4 days. Values for the undifferentiated (undiff) ESCs were set at 1 as indicated by a horizontal line in each histogram. Cdx2 is a marker for late endoderm and trophectoderm. Tead4 is a marker for early trophectoderm. Mean ± SEM; n=3. P values are calculated by unpaired t-test, *P<0.05, **P<0.01, ***P<0.001, ns = non significant.

6 Figure S4

A DAPI BRACHURY GATA4 MERGED DAPI SMA WT WT CNOT3-DM CNOT3-DM

4 days mesendoderm differentiation B DAPI BRACHURY FOXA2 MERGED DAPI SMA WT CNOT3-DM

8 days mesendoderm differentiation

DAPI BRACHURY GATA4 MERGED WT CNOT3-DM

8 days mesendoderm differentiation

Supplementary Figure 4: CNOT3-T292/S294 phosphorylation promotes survival of differentiating mesendodermal cells. A. Staining for ME markers shows that differentiation of Cnot3-DM ESCs with BMP4 + FGF2 for 4 days generates reduced numbers of ME cells. Left panel: Representative IF analyses of staining of differentiated cells for the mesodermal marker Brachury (red) and endodermal marker GATA4 (green) carried out on ME cells differentiated from WT and Cnot3-DM ESCs by incubating them with BMP4 + FGF2. Right panel: staining for the additional mesodermal marker smooth muscle actin (SMA) (green). Nuclei were

7 stained with DAPI. Merged red and green images show the ME cells. Scale bars in (C) = 100 µm; scale bar in (D) = 60 µm. For IF analysis at 8 days of differentiation, see Figure S4A. B. Representative immunofluorescence analysis images of ME cells generated by inducing differentiation of WT and Cnot3-DM ESCs with BMP4 + FGF2 for 8 days (left panel: scale bar: 100 µm). Top left panel: staining for the mesodermal marker Brachury (red) and endodermal marker FOXA2 (green). Bottom panel: staining for Brachyury (red) and the endodermal marker GATA4. Merged images show dual staining of the mesendodermal cells. Top right panel: staining for the mesodermal marker a-SMA (scale bar: 60 µm). Nuclei were counterstained with DAPI.

8 Figure S5

A

DAPI NESTIN DAPI -SMA DAPI FOXA2 + WT WT WT CNOT3-DM CNOT3-DM CNOT3-DM

Ectoderm Mesoderm Endoderm

ns ns *** SURVIVAL RATIO SURVIVAL SURVIVAL RATIO SURVIVAL SURVIVAL RATIO SURVIVAL

WT WT WT CNOT3 CNOT3 CNOT3 -DM -DM -DM

B CNOT3-DM C WT CNOT3-DM G2/M WT Clone 2 Clone 3 ns S

G1 Counts Counts 8 days mesendoderm differentiation Percent of cells

WT PI PI CNOT3-DM

Supplementary Figure 5: Differentiation of Cnot3-DM ESCs into the three germ layers. A. Top panels: Representative confocal images of differentiation of wild-type (WT) and Cnot3-DM ESCs into ectoderm (3 days; immunostained for Nestin (green)), Activin A induced mesoderm (3 days; immunostained for a-SMA (green)) and FGF2+retinoic acid induced endoderm (3 days; immunostained for FOXA2 (green). Nuclei were counterstained with DAPI; scale bar: 12 µm. (See Methods). Bottom panels: Graphs represent the ratio of suriving Cnot3-DM cells versus WT cells after 3 days of differentiation into each lineage. Cell survival was measured using WST-1 reagent (see Methods). Mean ± SEM; (P values calculated by unpaired t-test, ***P<0.001. ns = non- significant, n=3). B. Crystal violet staining showing the efficiency of BMP4+FGF2 induced mesendoderm (ME) differentiation of wild-type (WT) ESCs and two additional clones of Cnot3-DM ESCs after 4 days of differentiation.

9 C. Cell cycle analysis of ME cells obtained by inducing differentiation of WT and Cnot3- DM ESCs with BMP4 + FGF2 for 4 days. Histograms represent the percentage of cells at different phases of the cell cycle. Mean ± SEM; P values calculated by unpaired t-test, ns = non-significant, n=3.

10 Figure S6

A B

P adj< 0.05 WT CNOT3-DM NON SIG

WT CNOT3-DM PC2 : 22.69 % VARIANCE P VALUE P 10 log

- PC1 : 40.13 % VARIANCE

-log2 FOLD CHANGE

C D ENRICHED ONTOLOGY TERMS (DOWNREGULATED) ENRICHED TERMS (UPREGULATED) Cell adhesion Epithelium development Developmental processes Signal transducer activity Mesenchymal cell development Molecular transducer activity Embryonic placenta development Signaling receptor activity Negative regulation of Wnt signaling pathway Cell surface receptor signaling pathway Apoptotic process involved in morphogenesis Growth factor binding Regulation of programmed cell death Regulation of cell proliferation Positive regulation of kinase activity -log P VALUE 10 Multicellular organism development Regulation of MAPK cascade

-log P VALUE E 10

11 Supplementary Figure 6: Transcriptome profiling by RNA sequencing in differentiating wild-type and Cnot3-DM mesendoderm cells. A. and B. RNA-seq analysis was used to compare patterns of ME cells obtained by differentiating wild-type (WT) and Cnot3-DM ESCs for 4 days with BMP4 + FGF2. A. Volcano plot of RNA-Seq transcriptome data showing differentially expressed genes after 4 days of BMP4 + FGF2 induced mesendoderm differentiation of WT and Cnot3- DM ESCs. Fold change ratio >1.5; adjusted P value <0.05. Three biological replicates for each were analyzed. B. Principal component plot of the RNA-seq data showing PC1 and PC2 and the percent variance obtained for each. This reveals separate clustering of the genes that are differentially expressed between the WT and Cnot3-DM-derived 4-day differentiated mesendodermal cells. C. and D. Selected gene ontology (GO) profiles obtained out using the GOSeq Bioconductor package and showing upregulated genes (C) and downregulated genes (D), which are involved in signalling pathways that have key roles in developmental processes. Adjusted P value <0.05. E. GSEA analysis showing Hallmark gene sets that were significantly altered in the differentiated cells. See also Supplementary Tables S3A and S3B for the differentially regulated genes. The GO and GSEA analysis both showed that expression of genes involved in cell death and apoptosis was upregulated (C and E). The negative regulator of Wnt signalling Axin2, and TGFb, which is part of the apoptosis program, were also elevated. Additionally, the GO analysis revealed downregulation of gene expression programs that are involved in cell adhesion and developmental processes (D). The changes included downregulation of signalling molecules such as Wnt6. Bcl2, a key cell survival factor downstream of ERK signalling was also downregulated.

12 Figure S7

A

WT CNOT3-DM

CNOT3

P-ERK T202/Y204

ERK1/2

-TUBULIN

Undifferentiated ES cells

B DAPI TUBULIN PLA:CNOT3+ERK MERGED

Cytoplasmic WT Nuclear PLA: CNOT3+ERK

*

* AVERAGE PLA DOTS PER CELL PLA AVERAGE

CNOT3-DM WT CNOT3-DM

C

AURKB CNOT3 ERK WT CNOT3-DM

13 Supplementary Figure 7: Phosphorylation of CNOT3-T292/S294 alters interaction of CNOT3 with ERK in mesendodermal cells. A.Immunoblotting analysis of the indicated from cell extracts prepared from undifferentiated wild-type (WT) and Cnot3-DM ESCs. a-Tubulin was used as a loading control. B. PLA detection of endogenous interaction between CNOT3 and ERK using a second pair of antibodies (results from the first antibody pair are shown in Figure 7b). WT and Cnot3-DM ESCs were induced to differentiate into mesendodermal cells by incubation for 4 days with BMP4 + FGF2. Positive PLA signals are visible as red dots. Nuclei were counterstained with DAPI and staining for Tubulin was used to mark the cell cytoplasm. Boxed areas represent an enlarged cell indicated by the arrows; scale bar: 10 µm. PLA dots were quantified from randomly chosen fields from at least 50 cells for each biological replicate. Histogram represents average interactions per cell (dots per cell) as well as the nuclear and cytoplasmic distribution. Error bars represent the Mean ± SEM; n=2. (P values calculated by unpaired t-test, *P<0.05). C. Negative control single PLAs for analysis shown in Figure 6. The control assays were conducted with conducted with anti-Aurora B only, anti-CNOT3 only or anti-ERK only.

14 Supplementary Table 1

List of primers used for qPCR

Oct-4 F-5′-CCAATCAGCTTGGGCTAGAG-3′ R-5′-CTGGGAAAGGTGTCCCTGTA-3′ Nanog F-5′-TACCTCAGCCTCCAGCAGAT-3′ R-5′-GCAATGGATGCTGGGATACT-3′ Sox2 F-5'-CACAACTCGGAGATCAGCAA-3' R-5'-CTCCGGGAAGCGTGTACTTA-3' Sox1 F-5'-CCTCGGATCTCTGGTCAAGT-3' R-5'-TACAGAGCCGGCAGTCATAC-3' Fgf5 F-5′-TCTGGATCTCCTTTGCGTTT-3′ R-5′-GGGCTTCGAAAGCACATTTA-3′ Cdx2 F-5′-TGGTGTACACAGACCATCAGC-3′ R-5′-CCTTGGCTCTGCGGTTCT-3′ Tead4 F-5’-ATCCTGACGGAGGAAGGCA-3’ R-5'-GCTTGATATGGCGTGCGAT-3' Brachury F-5′-CAGCTGTCGGGGAGCCTGG-3′ R-5′-TGCTGCCTGTGAGTCATAAC-3′ Mixl1 F-5′-AGTTGCTGGAGCTCGTCTTC-3′ R-5′-TTCTGGAACCACACCTGGAT-3′ Gata4 F-5′-TTCCTCTCCCAGGAACATCAA-3′ R-5′-GCTGCACAACTGGGCTCTACTT-3′ Gata6 F-5’-ACAGCCCACTTCTGTGTTCCC-3’ R-5’-GTGGGTTGGTCACGTGGTACAG-3’ Sox17 F-5’- TTCTGTACACTTTAATGAGGCTGTTC-3’ R-5’-TTGTGGGAAGTGGGATCAAG-3’ Pdgfr F-5’-CTGGTGCCTGCCTCCTATGAC-3’ R-5’-CACGATCGTTTCTCCTGCCTTAT-3’ Cxcr4 F-5’-GAGGCCAAGGAAACTGCTG-3’ R-5’-GCGGTCACAGATGTACCTGTC-3’ Mln51 b F-5’- ATGACGATGAGGATCGGAAAAAC-3’ R-5’-GTCCCCTTGGGTCGGACTTC-3’

15 Supplementary Table 2

List of antibodies

Antibody Source Catalouge number Rabbit polyclonal anti-H3 ABCAM Cat# ab1791 Rabbit polyclonal anti-CNOT3 Bethyl Laboratories,Inc Cat# A302-156A Mouse monoclonal anti-CNOT3 Abnova Cat# H00004849- M01 Mouse monoclonal anti-ERK1/2 Santa Cruz Cat# sc-514302 Biotechnology Rabbit polyclonal anti-ERK1/2 Cell Signaling Technology Cat# 9102 Rabbit polyclonal anti-phospho- Cell Signaling Technology Cat# 9101 ERK1/2(T202/Y204) Mouse monoclonal anti-Aurora B Santa Cruz Cat# sc-393357 Biotechnology Mouse monoclonal anti-Lamin B1 Santa Cruz Cat# sc-374015 Biotechnology Mouse monoclonal anti-a Tubulin Santa Cruz Cat# sc-5286 Biotechnology Mouse monoclonal anti-OTX2 Santa Cruz Cat# sc-514195 Biotechnology Rabbit polyclonal anti-GATA4 ABCAM Cat# ab84593 Rabbit polyclonal anti-a SMA ABCAM Cat# ab5694 Mouse monoclonal anti-V5 ABCAM Cat# ab27671 Goat polyclonal anti-Brachyury Santa Cruz Cat# sc-17743 Biotechnology Rabbit monoclonal anti-FoxA2 ABCAM Cat# ab108422 Mouse monoclonal anti-Nestin Santa Cruz Cat# sc-23927 Biotechnology Mouse monoclonal anti-a Tubulin-FITC Merck Cat# F2168 Aanti-rabbit IgG-HRP Santa Cruz Cat# sc-2357 Biotechnology Anti-mouse IgG-HRP Cell Signaling Technology Cat# 7076

Donkey anti-goat IgG-Alexa fluor 555 Thermo Fisher Scientific Cat# A-21432 Donkey anti-rabbit IgG-Alexa fluor 488 Thermo Fisher Scientific Cat# A-21206 Goat anti-mouse IgG-Alexa fluor 488 Thermo Fisher Scientific Cat# A-11029

16

Supplementary Table 3A

RNA sequencing: List of genes downregulated in 4-day differentiated mesendoderm cells from Cnot3-DM ESCs compared with differentiated wild-type ESCs

17

18 Supplementary Table 3B

RNA sequencing: List of genes upregulated in 4-day differentiated mesendoderm cells from Cnot3-DM ESCs compared with differentiated wild-type ESCs

19

20