supporting information - proceedings of the national … relative expression dmrt1 ret e-cadherin...

12
Supporting Information Krentz et al. 10.1073/pnas.0905431106 SI Text Bisulfite Sequencing. GFP positive E13.5 germ cells from Dmrt1 / Oct4PE:GFP males were sorted into RNAprotect (Qiagen) by using a FACSDiva cell sorter. DNA was prepared by using AllPrep Micro kit (Qiagen) and bisulfate-treated by using EpiTect Bisulfite kit (Qiagen). H19, Lit1, and Xist loci were amplified and sequenced as described in ref. 1 and analyzed by using CpG Viewer (2). Array Comparative Genome Hybridization (CGH). DNA was phenol: chloroform:isoamyl alcohol-extracted from adult 129Sv wild-type and 129Sv Dmrt1 / testes and spleen tissue and from adult B6 wild-type and B6 Dmrt1 / testes. CGH was carried out on an Agilent mouse genome 4 44 K oligonucleotide microarray. Expression Array Data Analysis. Cel files were normalized by using the GC-RMA algorithm (3) as implemented in the Genedata Expressionist Software Package. Affymetrix probe set ID’s were converted to Entrez Gene Symbols by using Affymetrix map- pings for the array. Two group t tests were carried out to determine genes that were differentially expressed between genotypes. Genes found to have a P value 0.05 and an average fold change of greater than absolute value of 2 are included in the figure. Gene profiles in Dmrt1 / are shown relative to the expression pattern observed in wild-type. 1. Hajkova P, et al. (2002) Epigenetic reprogramming in mouse primordial germ cells. Mech Dev 117:15–23. 2. Carr IM, Valleley EM, Cordery SF, Markham AF, Bonthron DT (2007) Sequence analysis and editing for bisulfite genomic sequencing projects. Nucleic Acids Res 35:e79. 3. Wu Z, Irizarry R, Gentleman R, Martinez-Murillo F, Spencer F (2004) A model-based background adjustment for oligonucleotide expression arrays. J Am Stat Assoc 99:909 – 917. 4. Ancelin K, et al. (2006) Blimp1 associates with Prmt5 and directs histone arginine methylation in mouse germ cells. Nat Cell Biol 8:623– 630. 5. Looijenga LH, et al. (2006) Genomic and expression profiling of human spermatocytic seminomas: Primary spermatocyte as tumorigenic precursor and DMRT1 as candidate chromosome 9 gene. Cancer Res 66:290 –302. Krentz et al. www.pnas.org/cgi/content/short/0905431106 1 of 12

Upload: hoangthien

Post on 23-Mar-2018

222 views

Category:

Documents


1 download

TRANSCRIPT

Supporting InformationKrentz et al. 10.1073/pnas.0905431106SI TextBisulfite Sequencing. GFP positive E13.5 germ cells fromDmrt1�/� Oct4�PE:GFP� males were sorted into RNAprotect(Qiagen) by using a FACSDiva cell sorter. DNA was preparedby using AllPrep Micro kit (Qiagen) and bisulfate-treated byusing EpiTect Bisulfite kit (Qiagen). H19, Lit1, and Xist loci wereamplified and sequenced as described in ref. 1 and analyzed byusing CpG Viewer (2).

Array Comparative Genome Hybridization (CGH). DNA was phenol:chloroform:isoamyl alcohol-extracted from adult 129Sv wild-typeand 129Sv Dmrt1�/� testes and spleen tissue and from adult B6wild-type and B6 Dmrt1�/� testes. CGH was carried out on anAgilent mouse genome 4 � 44 K oligonucleotide microarray.

Expression Array Data Analysis. Cel files were normalized by usingthe GC-RMA algorithm (3) as implemented in the GenedataExpressionist Software Package. Affymetrix probe set ID’s wereconverted to Entrez Gene Symbols by using Affymetrix map-pings for the array. Two group t tests were carried out todetermine genes that were differentially expressed betweengenotypes. Genes found to have a P value �0.05 and an averagefold change of greater than absolute value of 2 are included inthe figure. Gene profiles in Dmrt1�/� are shown relative to theexpression pattern observed in wild-type.

1. Hajkova P, et al. (2002) Epigenetic reprogramming in mouse primordial germ cells.Mech Dev 117:15–23.

2. Carr IM, Valleley EM, Cordery SF, Markham AF, Bonthron DT (2007) Sequence analysisand editing for bisulfite genomic sequencing projects. Nucleic Acids Res 35:e79.

3. Wu Z, Irizarry R, Gentleman R, Martinez-Murillo F, Spencer F (2004) A model-basedbackground adjustment for oligonucleotide expression arrays. J Am Stat Assoc 99:909–917.

4. Ancelin K, et al. (2006) Blimp1 associates with Prmt5 and directs histone argininemethylation in mouse germ cells. Nat Cell Biol 8:623–630.

5. Looijenga LH, et al. (2006) Genomic and expression profiling of human spermatocyticseminomas: Primary spermatocyte as tumorigenic precursor and DMRT1 as candidatechromosome 9 gene. Cancer Res 66:290–302.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 1 of 12

TRA98 GM114 Merge

A B C

D E F

+/+

-/-

G H I

J K L

+/-

-/-

TRA98 BLIMP1 Merge

MVH

+/+

-/-

M

N

E13.5

E13.5

E13.5

E13.5

E13.5

E13.5

Fig. S1. Early gonocyte developmental markers are normal in Dmrt1 mutant germ cells. (A–F) GM114 is expressed at similar levels in germ cells of Dmrt1 mutantsand wild-type controls. (G–L) BLIMP1 translocation to the cytoplasm (4) is complete by E13.5. (M and N) Similar numbers of MVH positive germ cells are presentin Dmrt1 mutant testes compared to wild-type controls. In addition, the number of germ cells (MVH�) was counted for 50 tubules for Dmrt1�/� (avg � 5.0, SD �0.5) and wild-type (avg � 4.8, SD � 0.3) at P0; P � 0.69 (n � 2). (Scale bars, 50 �m.)

Krentz et al. www.pnas.org/cgi/content/short/0905431106 2 of 12

Lit1H19 Xist

-/- #1

-/- #2

-/- #3

+/+

3’ 5’Kcnq1

E11 E105’ 3’ 5’ 3’

E15’ 3’

G-rich

Fig. S2. Parental imprinting is mostly erased by E13.5 in Dmrt1 mutant germ cells. Methylation of the paternally imprinted locus H19 (Left) and the maternallyimprinted locus Lit1 within the Kcnq1 gene (Center) is erased in Dmrt1 mutant germ cells. The Xist locus (Right) underwent normal demethylation by E13.5 inDmrt1 mutant germ cells. Dmrt1�/� mutant #2 had not completed parental imprinting erasure by E13.5. Open circles, unmethylated; black filled circles,methylated; gray circles, sequence did not align to reference.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 3 of 12

TRA98 P27Kip1 MergeI J K

L M N

E16.5

E16.5

Dmrt1 +/-

Dmrt1 -/-

Dnd1Ter/+ Bax +/-

Dnd1Ter/Ter Bax -/-

A

B

E14.5

E14.5

E-cadherinDMRT1

;

;

TRA98 P27Kip1 Merge

E13.5

E13.5

Dmrt1 +/-

Dmrt1 -/-

C D E

F G H

-/-+/+

exp

ress

ion

Rel

ativ

e

E13.5

P27Kip1

E15.5

O

0

0.5

1

1.5

2

2.5

QP

P0.5

Fig. S3. Dmrt1 likely acts independently of Dnd1. (A and B) DMRT1 is expressed at similar levels in DndTer/Ter;Bax�/� and DndTer/�;Bax�/� germ cells at E14.5.Germ cells are marked with E-cadherin, which is expressed in all germ cells at this stage. (C–N) Double staining for the germ cell marker TRA98 and the cell cycleinhibitor P27Kip1 at E13.5 (C–H) and E16.5 (I–N). P27Kip1 is expressed in Dmrt1 mutant germ cells at a similar level compared to heterozygous control germ cellsat E13.5 and E16.5. (Scale bars, 50 �m.) (O) P27Kip1 is expressed in E13.5 and E15.5 Dmrt1�/� gonads at similar levels to wild-type controls. Expression wasnormalized to Hprt. Error bars: SD from three animals of each genotype. (P and Q) Serial sections of Dmrt1�/� gonad at P0.5. (P) PTEN and (Q) P-AKT Thr-308are expressed at similar levels in EC cluster, outlined in black, compared to adjacent tubules.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 4 of 12

Fo

ld e

nri

chm

ent

E13.5:DMRT1 ChIP

0

4

2

6B

-596 to 259P19INK4d Btg2

-564 to -357

A

P15INK4b0

0.5

1.0

1.5

exp

ress

ion

Rel

ativ

e

*

-/-+/+

***

E15.52.0

2.5

P18INK4c P19INK4d

*

Fig. S4. DMRT1 regulates P15INK4b, P18INK4c, and P19INK4d. (A) Dmrt1�/� testes have elevated expression of p15INK4b and have reduced expression of P18INK4c andP19INK4d relative to wild-type controls at E15.5 based on qRT-PCR. Expression is normalized to Hprt expression. Error bars: SD from three animals of each genotype.

*, P � 0.05; ***, P � 0.0005. (B) DMRT1 binds to P19INK4d promoter in E13.5 testes. ChIP-qPCR of P19INK4d and Btg2 regulatory regions, comparing enrichmentin chromatin immunoprecipitated with DMRT1 relative to input chromatin. Numbering below gene name indicates region covered by amplicons tested, relativeto start of transcription. Error bars: SD of duplicate qPCR of sample. The negative control promoter Btg2 was chosen because it is expressed at similar levels ingerm cells and Sertoli cells and does not change expression in Dmrt1 mutant testes.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 5 of 12

Mybl2 Klf2 Tdgf1 Zic3 Ss18 Prtg Jam2 Eomes FoxD3 Nr0b1

**

*

Rel

ativ

e ex

pre

ssio

n

E15.5 -/-+/+

0

5

10

15

20

25

** **

*

*

**

Fig. S5. Early direct response genes of OCT3/4 are overexpressed in Dmrt1�/� gonads. Seven out of 10 genes examined by qRT-PCR were significantlyup-regulated in Dmrt1�/� gonads compared to wild-type at E15.5. Expression was normalized to Hprt. Error bars: SD from three animals of each genotype.

*, P � 0.05; **, P � 0.005.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 6 of 12

Eras

Rel

ativ

e ex

pre

ssio

n

E15.5-/-+/+

0

2

4

6

**

B6 129Sv

8

10

Fig. S6. Eras is significantly up-regulated in 129Sv Dmrt1�/� gonads relative to 129Sv wild-type and B6 gonads of either genotype. Expression was normalizedto Hprt. Error bars: SD from three animals of each genotype. **, P � 0.005.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 7 of 12

NANOG

exp

ress

ion

Rel

ativ

eDMRT1 RET

E-cadherin (CDH1)

exp

ress

ion

Rel

ativ

eex

pre

ssio

nR

elat

ive

exp

ress

ion

Rel

ativ

eA B

C D

E

RetRetRet

SS1 SS2 SS3 SS4 SS5

SE1

SE2 SE3 SE4 DG1DG2 DG3

EC1 EC3EC2 EC4SS1 SS2 SS3 SS4 SS5 SE1 SE2 SE3

SE4 DG2DG3

EC1 EC3EC2 EC4

Dmrt1

SS1 SS2 SS3 SS4 SS5SE1

SE2 SE3SE4 DG1DG2DG3

EC1EC3EC2 EC4

SS1 SS2 SS3 SS4 SS5SE1

SE2SE3

SE4DG1DG2DG3

EC1EC3

EC2 EC4

exp

ress

ion

Rel

ativ

e

CDKN2C (P18INK4C)

CDKN2D (P19INK4D)

CDKN2CCDKN2C

CDKN2DCDKN2D

CDH1CDH1

SS1 SS2 SS3 SS4 SS5SE1

SE2 SE3SE4 DG1

DG2DG3EC1 EC3EC2 EC4

SS1 SS2 SS3 SS5SE2 EC1 EC3EC2 EC4

exp

ress

ion

Rel

ativ

e

F

NANOGNANOGNANOG

0

1

2

3

4

-1

-2

-3

-4

0

2

4

6

-2

-4

-6

-8

0

2

4

6

-2

-4

-6

-8

0

1

2

3

-5

-1

-2

-3

-4

0

1

2

3

4

-1

-2

-3

-4

-5

0

1

2

3

-1

-2

-3

-4

SS4SE1 SE3 SE4 DG1DG2DG3

DG1

Fig. S7. Gene expression in human germ cell tumors correlates with DMRT1 expression level. Expression profiling of spermatocytic seminomas (SS), seminomas(SE), dysgerminomas (DG), and embryonal carninomas (EC). (A) DMRT1, (B) RET, (C) E-CADHERIN/CDH1, (D) CDKN2C (P18INK4C), (E) NANOG, and (F) CDKN2D(P19INK4D). Tumor samples with high DMRT1 expression (SS1–5, SE2) also have high RET, CDKN2C, and CDKN2D expression and reduced E-CADHERIN expression.EC cells have low DMRT1 expression and have low RET and CDKN2C expression and elevated E-CADHERIN expression. NANOG expression was reduced in SS andelevated in SE, DG, and EC. Relative expression was determined in relation to the geometric mean of hybridization intensities for each patient sample. Differencesin relative expression represent deviations from the geometric mean (5).

Krentz et al. www.pnas.org/cgi/content/short/0905431106 8 of 12

L1 L1 L1 129Sv Dmrt1+ 129Sv Dmrt1-L1

Chromosome 19

Fig. S8. Diagram of complementation test with L1 and 129Sv Dmrt1. Blue bar represents region of MOLF-derived DNA substituted for 129Sv DNA in congenicL1 chromosome 19. Triangle indicates position of Dmrt1 null mutation.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 9 of 12

Table S1. Antibodies

Antibody Supplier Catalog no. Diluton

RET Neuromics GT15002 1:100TRA98 Bio Academia 73–003 1:500BLIMP1/PRDM1 Lifespan Biosciences LC-C2538 1:50Ki67 Thermo Scientific RM-9106-S1 1:250E-CADHERIN Zymed Laboratories 12–1900 1:500DMRT1 140A David Zarkower 1:500STRA8 Pierre Chambon 1:200SOX2 Chemicon AB5603 1:300OCT3/4 Santa Cruz SC-9081 1:150NANOG Cosmo Bio REC-RCAAB0002PF 1:250MVH Abcam AB13840 1:300P27Kip1 Santa Cruz Sc-528 1:100P-AKT Thr 308 Cell Signaling 244F9H2 1:50PTEN Cell Signaling 138G6 1:100GM114 Blanche Capel 1:500Goat � rat Alexa 594 Molecular Probes A11007 1:500Goat � rabbit Alexa 488 Molecular Probes A11034 1:500Goat � mouse Alexa Moleular Probes A11029 1:500Biotinylated � rabbit IgG Vector Labs BA-1000 1:200

Krentz et al. www.pnas.org/cgi/content/short/0905431106 10 of 12

Table S2. qRT-PCR primer sequences

qRT-PCR Forward Reverse

Ret TTCCGAGGAAATTCCCACTT CAGGGCTTCCCAATCAGTTAE-cadherin TGAGCTGCCTCACAAAAACA AGCCTGAACCACCAGAGTGTSox2 GAGTGGAAACTTTTGTCCGAGA GAAGCGTGTACTTATCCTTCTTCATOct3/4 GAAGCAGAAGAGGATCACCTTG TTCTTAAGGCTGAGCTGCAAGNanog CCTCAGCCTCCAGCAGATGC CCGCTTGCACTTCATCCTTTGp15INK4b AATAACTTCCTACGCATTTTCTGC CCCTTGGCTTCAAGGTGAGP27Kip1 TGGAGGGCAGATACGAATG CGGGGGCCTGTAGTAGAACTP18INK4c CAGATTAACCATCCCAGTCCTT CCCCTTTCCTTTGCTCCTAAP19INK4d AATGTGACCCAAGGCCACT TTTCCTCTTTTGTTGAGAAGTAACCMybl2 AGGACAAGGAACAGCACCAG GCAGCTATGGCAATCTCCTCKlf2 ACCAAGAGCTCGCACCTAAA CTGTGACCTGTGTGCTTTCGTdgf1 GGAAGGCACAAACTGGAAAG GTTTGAATTTGGACCCGTTGZic3 TCCACAAGAGGACCCATACAG TATAGGGCTTGTCCGAGGTGss18 ACTGCTGACCTTGACCCTGA ACCTCCACAGCAAGGATACCEomes CAATGTTTTCGTGGAAGTGG GTGGGAGCCAGTGTTAGGAGPrtg GGGGAAGAACCTGGAGAGAG GACCGTAGCTATTGATGATCAGGJam2 CCTCCTGATGCTGCTGCT TGTGACTTCTTGACGGTGGTFoxD3 CAGCAACCGTTTTCCGTACT GGGTCCAGGGTCCAGTAGTTNr0b1 TATCTGAAAGGGACCGTGCT CAGCGGATCTGATCTGGTACTEras CCTCTGGAGATCTGGTGCAT GCCCCTCATCAGACTGCTACHPRT TCCTCCTCAGACCGCTTTT CCTGGTTCATCATCGCTAATC

Krentz et al. www.pnas.org/cgi/content/short/0905431106 11 of 12

Table S3. ChIP primer sequences

Distance from transcriptionstart

ChIP primers Forward Reverse Size F primer R primer

Sox2 4 CAGCAGCCACATCTCAGAAA TTTTCTCCTTAATCGGGGTGT 125 �2,104 �1,980Sox2 5 ATAAGTCCTTTCCGGGTTCC ATGCCCACTTTAGACCCAGA 148 �1,937 �1,790Sox2 2 AGCCAAGCCTGGGAGAAT GCTCCGCTCATTGTCCTTAC 108 �1,766 �1,659Sox2 6 GGGCAGAGAGGGAAGGATAA TTTTTCCCTGCAATACTCTCG 135 �1,710 �1,568Sox2 3 TGGCATCAGGACTTTCTTCC GCCCCGTCTAAGTTTCCTTC 100 �1,499 �1,400p19INK4d 2 GGGCCCATCTTGGATCTTTA GCGCTTTACCACTCAGAAGC 189 �408 �596p19INK4d 4 TTGTCTCTCCCGTAGCAGTG CGCTGGGGTCAGTTAAACC 141 �243 �383p19INK4d 1 GACTCACCCTCCCTCCTTTC GGCTCTCGCTACTCTGTTGC 160 �6 �165p19INK4d 3 CCTCCATTGTCACAAAAGCA GATCCAACTCCCCAAAATGA 119 259 157Btg2 GACACTGACAGAGCCGTTCA ACACTCCTCCCACCAAACAG 208 �357 �564

All ChIP qPCR reactions performed in the presence of 0.5 M Betain.

Krentz et al. www.pnas.org/cgi/content/short/0905431106 12 of 12