Motif Tools

ETS class I (ETS1/ERG/FLI1)

ETSACCGGAAGT81 motifsaka ETS site, EBS, ETS-binding site
CTGATAGCGACGGATACAGAGCT
representative: ETS1.H14CORE.0.S.B · data consensus ACCGGAAGT
members tight · 0.82stated ACCGGAAGT · data ACCGGAAGT18 in vitro5 ChIP-only1 cross-species15 references, all PMIDs verified

ETS proteins are a family of 28 human transcription factors that each bind DNA through a single ~85-residue winged helix-turn-helix domain and read a nine-base site whose informative part is the four-base core GGA(A/T). The class I members - ETS1 and ETS2, ERG, FLI1, FEV, the PEA3 proteins ETV1/4/5, the ternary complex factors ELK1/3/4, the repressors ERF and ETV3, and GABPA - select essentially the same sequence, ACCGGAAGT, so a match names the class rather than a protein. These sites sit just upstream of a large fraction of human promoters and, in enhancers, next to RUNX, GATA, PAX and SRF sites. Class I also contains every ETS gene recurrently rearranged in cancer: ERG in prostate cancer, FLI1 and ERG in Ewing sarcoma, ETV1 and ETV4 in both.

The factors

The human genome contains 28 ETS genes. All encode one copy of the ETS domain: about 85 residues folded as a winged helix-turn-helix, three alpha-helices packed on a four-stranded antiparallel beta-sheet, that binds DNA as a monomer 8602247,21548782. Profiling the DNA-binding preference of every human and mouse ETS protein split the family into four classes, and class I is by far the largest: the ETS subfamily (ETS1, ETS2), the ERG subfamily (ERG, FLI1, FEV), the PEA3 proteins (ETV1, ETV4, ETV5), the ternary complex factors (ELK1, ELK3, ELK4), the repressors ERF and ETV3, and GABPA 20517297. GABPA is the structural exception: it functions only in complex with the ankyrin-repeat subunit GABPB, which it recruits using the ETS domain plus a carboxy-terminal extension 9461436. Several class I proteins carry autoinhibitory segments flanking the ETS domain that are relieved by partner binding 10075931,21548782. Expression ranges from near-ubiquitous (GABPA, ELK1, ERF, ETV3) through lymphoid-enriched (ETS1), endothelial and hematopoietic (ERG, FLI1), to narrowly restricted (FEV). SPDEF is grouped here in this atlas, but biochemically it is the sole class IV factor and prefers a GGAT core 20517297,21548782.

The motif

The site is nine base pairs, ACCGGAAGT, and nearly all of its information is in the GGA(A/T) core 21548782. Two invariant arginines in recognition helix H3 hydrogen-bond the two core guanines; a conserved tyrosine (Tyr395 in ETS1) supports the interface; a lysine in the wing contacts DNA 5' of the core and a lysine in the H2-H3 turn contacts the opposite strand 3' of it 8602247,21548782. The core is read base by base, the flanks indirectly, through sequence-dependent backbone geometry 21548782. The flanks separate subfamilies, and separate them weakly: class I prefers A three bases 5' of the core and disfavors C there, the ELF/EHF group prefers C at that position (CCCGGAAGT), SPI1/SPIB/SPIC use a purine-rich 5' flank (AGAGGAAGTG), and SPDEF alone prefers GGAT 20517297. Quantitatively, different ETS proteins vary only about ten-fold in affinity for the same consensus, whereas one protein's affinity across variant sites spans five- to eighty-fold 21548782; sequence therefore grades site quality much better than it identifies a paralogue, and no class I member is separable from another by motif alone. The CCGG contains a CpG, and methylating it blocks binding 20019798,21548782. Tandem arrangements are a distinct case: GGAA microsatellites of four or more consecutive repeats are the response elements of the EWS-FLI1 fusion 18626011,21548782, and some models grouped here (JASPAR MA0149.1, archetype AC0561) are GGAA repeats rather than single sites.

Where it acts

ETS occupancy follows two logics. At promoters it is largely redundant: a strong CCGGAAGT sits 20-40 bp upstream of the transcription start site in CpG-island housekeeping promoters - roughly a quarter of human promoters carry an ETS site near the start site 17652178,21548782 - and different ETS proteins bind interchangeably, with most GABPA-occupied promoters in Jurkat cells also occupied by ETS1 17652178,20019798. These sites lie in hypomethylated CpG islands, so the CpG in the core stays unmethylated 21548782. At distal, DNase-sensitive enhancers, occupancy is factor-specific, the sequence departs from consensus 17652178,20019798, and specificity comes from composite elements: ETS1 with RUNX1 at AGGA(A/T)GTGG, which binds more than 30-fold cooperatively and marks over a thousand ETS1-specific regions in T cells 10075931,20019798,21548782; ETS with GATA in hemogenic endothelium, where FLI1, GATA2 and SCL cross-regulate through neighboring sites 17962413; ETS1 with PAX5 at the mb-1 promoter, where Gln22 of the PAX5 paired domain hydrogen-bonds Tyr395 of ETS1 11779502,21548782; the ternary complex factors with SRF at serum response elements 11715049,21548782; and ERG with the androgen receptor in prostate cancer, where the two share about 44% of sites 21548782. The ETS-IRF composite element (EICE) belongs to the SPI1 subfamily and is described there.

Biology and disease

Class I factors run housekeeping transcription, growth-factor responses and lineage programs. Gabpa-null mouse embryos die before implantation 15199140, consistent with GABPA's occupancy of housekeeping promoters. Ets1 loss impairs T-cell and natural killer cell development. ERG and FLI1 are required for endothelial and hematopoietic specification. The clinical weight of the class is oncogenic rearrangement: TMPRSS2-ERG and related fusions place ERG, ETV1, ETV4 or ETV5 under androgen control in roughly half of prostate cancers 16254181,21548782, and EWSR1-FLI1 or EWSR1-ERG defines Ewing sarcoma, where the fusion's prion-like EWSR1 domain converts GGAA microsatellites into de novo enhancers 18626011,28844694.

Reading a hit

Treat a single match as weak evidence. The four-base GGAA core is expected about 12 million times per strand in a 3.1-Gb genome and the full nine-mer ACCGGAAGT on the order of 10,000 times per strand by base-composition arithmetic; CCGG-containing words are rarer than that because the genome is CpG-depleted, but the site is still common enough that isolated hits mean little without context: conservation, accessibility, a nearby partner site, or a promoter position close to the start site. A hit identifies an ETS site, not an ETS protein: ETS1, ERG, FLI1, ELK1 and GABPA are not distinguishable from sequence. ChIP-derived models carry partner signal: HOCOMOCO's ChIP-based ETS1 model (ETS1.H14CORE.1.P.B) has consensus ACAGGAAGTGG, whose trailing GTGG is the RUNX1 half-site of the ETS1-RUNX1 composite 20019798, and several other ChIP-derived ETS models converge on CAGGAAG rather than the in vitro CCGGAAGT 20019798,21548782. Half of the JASPAR entries grouped here (16 of 32) are CAP-SELEX heterodimer models such as ERF::SREBF2 or ETV5::FIGLA, where much of the matrix is the partner's site, not an ETS site. Finally, consensus strings are stored on either strand: GABPA's JASPAR model reads CACTTCCTGT, the reverse complement of ACAGGAAGTG.

How this was worked out

  1. 1990mediumSequence comparison across ETS proteins, DNA-binding assays2253872

    The ETS domain is named as a new DNA-binding motif shared by the ETS proteins and shown to recognize a purine-rich core centred on GGA.

    Still open then: How one short domain read so small a core, and whether family members differed in what they bound, was unknown.

  2. 1996highX-ray crystallography8602247

    The first ETS domain-DNA crystal structure (PU.1) shows a winged helix-turn-helix contacting ten base pairs through a loop-helix-loop: two arginines of the recognition helix read the GGA core, a wing lysine binds 5' of it and a turn lysine 3' of it on the opposite strand.

    Still open then: No base contacts outside the core were seen, so the flanking preferences that separate subfamilies were left to indirect readout.

  3. 2005highCancer outlier profile analysis, FISH, RT-PCR, androgen stimulation of cell lines16254181

    Outlier-expression analysis finds ERG and ETV1 overexpressed in prostate cancer through fusion to the 5' end of the androgen-regulated TMPRSS2 gene; 23 of 29 tumors examined by FISH carry an ERG or ETV1 rearrangement.

    Still open then: Which genomic sites the overexpressed ETS proteins act through, and whether they differ from those of the endogenous ETS factors, was not addressed.

  4. 2007mediumChIP-chip on promoter microarrays, motif analysis17652178

    Promoter-array ChIP for three ETS proteins in a T-cell line finds that redundant occupancy is the rule: consensus sites near the start sites of housekeeping genes are shared, while factor-specific sites diverge from consensus, lie farther away, and include an ETS1-RUNX1 cooperative site.

    Still open then: Promoter arrays could not see enhancers, where specific occupancy turned out to concentrate.

  5. 2008highChIP-chip, in vitro binding, reporter assays, genome-wide sequence analysis18626011

    EWS-FLI1 regulates targets such as NR0B1 through GGAA microsatellites; the NR0B1 repeat binds the fusion in vitro and in vivo, is necessary and sufficient for regulation, and activity depends on the number of consecutive GGAA motifs.

    Still open then: Whether repeat binding was a property of the FLI1 ETS domain or required the EWS portion of the fusion was not settled.

  6. 2009mediumChIP-seq, DNase I sensitivity, in vitro binding20019798 21548782

    ETS1 ChIP-seq in Jurkat T cells separates redundant occupancy with GABPA at housekeeping promoters from ETS1-specific occupancy at T-cell enhancers, which carry either a CpG-free CAGGAAGT variant or an ETS-RUNX composite; more than a thousand specific regions show AGGA(A/T)GTGG.

    Still open then: FLI1 and ERG also co-occupy sites with RUNX1 in stem cells, so how specific the composite is for ETS1 was left open.

  7. 2010highHigh-throughput microwell binding assay, protein-binding microarrays, ChIP-seq20517297

    DNA-binding profiles of all human and mouse ETS proteins, from two in vitro methods, fall into four classes; every cancer-linked ETS factor (ERG, ETV1, ETV4, FLI1) is in class I, and ChIP-seq for one member of each class confirms that the small in vitro differences affect site choice in cells.

    Still open then: Members within class I remained indistinguishable, so paralogue specificity had to come from elsewhere.

  8. 2017mediumChIP-seq in mesenchymal stem cells and Ewing sarcoma cells, mutagenesis, proteomics, in vitro phase-transition assays28844694

    Wild-type FLI1 occupies canonical ETS sites but not GGAA repeats, whereas EWS-FLI1 binds repeats, recruits BAF complexes and opens de novo enhancers, and needs the tyrosines that drive phase transitions of the EWSR1 prion-like domain to bind them. This revises the earlier view that the EWS portion contributes only transcriptional activation.

Still unresolved

How much to trust a hit

A hit tells you the family, not the member: all class I ETS proteins select the same nine-mer in vitro and differ only in weak flanking preferences, so a "FLI1 site" and an "ERG site" are the same sequence. The core is short and frequent, so a match outside an accessible, conserved or composite context is expected by chance. Sixteen of the 32 JASPAR models here are CAP-SELEX heterodimers whose matrices are partly a bHLH, FOX, HOX or SREBF site; four models here are GGAA repeats rather than single sites and will hit any GGAA microsatellite; and the ChIP-derived HOCOMOCO models (evidence letter P) include RUNX-extended and CAGGAAG variants that reflect partners and indirect binding rather than intrinsic ETS specificity. SPDEF is placed in this family here but reads GGAT and is biochemically its own class.

References

81 motifs in this family — JASPAR 32 · HOCOMOCO 23 · CIS-BP 16 · Vierstra 10
motifdatabaseconsensusbits
EWSR1-FLI1 (MA0149.1)JASPARGGAAGGAAGGAAGGAAGG32.9
M08111_2.00CIS-BPGGAAGGAAGGAAGGAAGG32.9
M04778_2.00CIS-BPACCGGAAGtaCTTCCGGT26.7
M07950_2.00CIS-BPGCcTkCTGGGAatTGTAGTtC25
ERF::SREBF2 (MA1939.2)dimerJASPARCGGAAATCACGTGAT20.8
ERF::NHLH1 (MA1938.2)dimerJASPARaGCAGCTGCCGGAAgT20.6
AC0561:FLI/ZNF:EtsVierstraTTCCTTCCtTCCTTCC19.9
AC0079:ERF::NHLH:['Tryptophan_cluster_factors',_'Basic_helix-loop-helix_factors_(bHLH)']VierstraGCAGCTGCCGGAAg19.4
AC0365:SPDEF:EtsVierstraTaATCCGGGACCA19
ELK1::HOXB13 (MA1932.2)dimerJASPARACCGGAAGThGTAAA18.9
AC0643:ETS/ETV:EtsVierstracGGAAGTACTTCCg18.6
AC0646:ERG/FLI:EtsVierstraCCGGAtATCCGG18.4
ELK1::SREBF2 (MA1933.2)dimerJASPARCCGGAAGTggCGTGA18.3
AC0105:SPDEF:EtsVierstraCAGdAAGAAGTAaa18.1
ERF::FIGLA (MA1934.2)dimerJASPARCCGGAAaCAcCTG18
ELK1::HOXA1 (MA1931.1)dimerJASPARACCGGAAGTAATT17.4
ETV5::FIGLA (MA1945.2)dimerJASPARcGGAAgCAGGTG17.4
ETV4.H14CORE.1.SM.BSMBHOCOMOCOACCGGAAGT15.7
M04752_2.00CIS-BPACCGGAAGT15.2
AC0076:ERF::FIGLA/ETV2::FIGLA:['Tryptophan_cluster_factors',_'Basic_helix-loop-helix_factors_(bHLH)']VierstraCGGAAgCAGgTG15.2
AC0541:ETS/SIX:EtsVierstraACTACAacTCCCAG14.8
SPDEF.H14CORE.1.SM.BSMBHOCOMOCOAaCCGGATGTdvnnhccg14.7
M02973_2.00CIS-BPAcCCGGATGT14.5
FEV (MA0156.4)JASPARACCGGAAGT14.4
SPDEF (MA0686.2)JASPARAcCCGGATGT14.3
GABPA.H14CORE.0.PSM.APSMAHOCOMOCOACCGGAAGT14.3
ELK1 (MA0028.3)JASPARACCGGAAGT14.2
ELK1.H14CORE.0.PSM.APSMAHOCOMOCOACCGGAAGT14.2
ERG.H14CORE.1.SM.BSMBHOCOMOCOACCGGAAGT14.2
ETV5 (MA0765.4)JASPARACCGGAAGT14.1
ETV5::HOXA2 (MA1948.2)dimerJASPARCGGAaGTaATTA14.1
FLI1::DRGX (MA1949.2)dimerJASPARaCCGGAAGTAATTA14.1
ELK3.H14CORE.0.PSM.APSMAHOCOMOCOaCCGGAAGT14.1
ELK3 (MA0759.3)JASPARACCGGAAGT14
ETV5::FOXI1 (MA1946.2)dimerJASPARGTAAACAGGAaG14
M02734_2.00CIS-BPCGGAtaTCCG14
AC0054:ELK/ETV:EtsVierstraTCGTAAATGC13.9
FEV.H14CORE.0.S.BSBHOCOMOCOaCCGGAAGT13.8
ETS2.H14CORE.0.S.CSCHOCOMOCOaCCGGAAGT13.5
ETS2 (MA1484.2)JASPARACCGGAAGt13.3
FLI1 (MA0475.3)JASPARACCGGAAgT13.2
ETV5::DRGX (MA1944.2)dimerJASPARCGGAAGtAATTA13.2
ETV1.H14CORE.1.PM.APMAHOCOMOCOgCAGGAAGTg13.2
ETS1.H14CORE.1.P.BPBHOCOMOCOaCaGGAAGTGG13.1
FLI1::FOXI1 (MA1950.2)dimerJASPARtAAACAGGAAg13
FLI1.H14CORE.1.P.BPBHOCOMOCOGgcAGGAAGbgg13
M07952_2.00CIS-BPCCGGAAGT13
ERF::FOXO1 (MA1936.2)dimerJASPARgtaAACAGGAAG12.9
ERG.H14CORE.0.P.BPBHOCOMOCOaCCGGAAGTg12.9
M02968_2.00CIS-BPACCGGAAGT12.9
ELK4 (MA0076.3)JASPARaCTTCCGG12.7
ERF (MA0760.2)JASPARACCGGAAGT12.7
FLI1.H14CORE.0.PSM.APSMAHOCOMOCOaCCGGAAGTg12.7
M04730_2.00CIS-BPaCCGGAAGT12.6
ETV1.H14CORE.0.PSM.APSMAHOCOMOCOACCGGAAGT12.5
SPDEF.H14CORE.0.PSM.APSMAHOCOMOCOamCCGGATgt12.5
M04791_2.00CIS-BPaCCGGAAGT12.5
ETS1.H14CORE.0.S.BSBHOCOMOCOACCGGAAGT12.4
M07957_2.00CIS-BPCCGGAAGt12.2
Erg (MA0474.4)JASPARACAGGAAGTg12.1
ETV3.H14CORE.0.SM.BSMBHOCOMOCOACCGGAAGT12.1
ETS1 (MA0098.4)JASPARACCGGAAgT12
ERF::HOXB13 (MA1937.2)dimerJASPARCGGAAgtngTaAA11.9
GABPA (MA0062.4)JASPARCACTTCCTGT11.7
ETV1 (MA0761.3)JASPARACAGGAAGT11.6
ETV5::FOXO1 (MA1947.2)dimerJASPARGtAAACAGGA11.6
M02970_2.00CIS-BPACCGGAAGT11.6
ETV3 (MA0763.2)JASPARACCGGAAGT11.5
ETV4 (MA0764.4)JASPARaCCGGAAGT11.5
ELK4.H14CORE.0.PSM.APSMAHOCOMOCOCCGGAAGt11.4
ETV5.H14CORE.0.PS.APSAHOCOMOCOCCGGAAGt11.4
ETS2.H14CORE.1.P.BPBHOCOMOCOccGGAAGtG11.2
ETV4.H14CORE.0.P.BPBHOCOMOCOCaGGAAG11
ETV5.H14CORE.1.PM.APMAHOCOMOCOCCGGAAGt10.9
M03006_2.00CIS-BPaCCGGAAGT10.6
ERF.H14CORE.0.PS.APSAHOCOMOCOAGGAAgT9.8
M04800_2.00CIS-BPcCGGAaG9.6
ERF::FOXI1 (MA1935.2)dimerJASPARACaGGAAg9.5
M04708_2.00CIS-BPCGGAaG8.8
AC0641:ETV/ETS:EtsVierstraGGaAGGAAG7.8
M01477_2.00Mus musculus 1CIS-BPGGAa5.9

HOCOMOCO tags read as assay + quality grade: P ChIP-seq, S HT-SELEX, M methyl-HT-SELEX, G GHT-SELEX, I SMiLE-seq, B PBM; grade A means two or more assay types agreed. A P-only matrix is a record of where the protein was crosslinked, which need not be a sequence it binds itself.

The factors, elsewhere

how this page was made Drafted by a language model reading the primary literature, then checked: every cited PMID is resolved against PubMed and its citation text regenerated from PubMed's own record, and every claim about the motifs themselves is checked against the matrices in the four source databases. Methods and limits →

Other ETS families