Pnas 1613199114 Sapp
Pnas 1613199114 Sapp
B
Gene ID Name Description
At2g18550 HB21 HD-ZIP Homeobox protein 21
At4g36740 HB40 HD-ZIP Homeobox protein 40
At5g66700 HB53 HD-ZIP Homeobox protein 53
At1g18330 RVE7 MYB protein
At4g11360 RHA1B RING-H2 finger A1B
At4g34000 ABF3/DPBF5 bZIP ABA responsive elements-binding factor 3
At2g46270 GBF3 bZIP G-box binding factor 3
At1g69490 NAP NAC domain containing protein activated by AP3/PI
At5g39610 ORE1/NAC6 NAC domain containing protein 6
At3g29575 AFP3 ABI five binding protein 3
At1g75490 ERF/AP2 DREB subfamily A-2 of ERF/AP2 protein
At2g36080 ABS2/NGAL1 AP2/B3-like transcriptional factor family protein
At2g18050 HIS1-3 Histone 1.3
At3g14440 NCED3 Nine-cis-epoxycarotenoid dioxygenase 3
At4g14930 Survival protein SurE-like phosphatase/nucleotidase
At2g47770 TSPO Outer membrane tryptophan-rich sensory protein
At1g79520 Cation efflux family protein
At4g17550 G3Pp4 Major facilitator superfamily protein
At1g73480 Alpha/beta-Hydrolases superfamily protein
At3g17520 LEA Late embryogenesis abundant protein
At2g35950 EDA12 Embryo sac development arrest 12
At1g79270 ECT8 Evolutionarily conserved C-terminal region 8
At1g11210 Protein of unknown function
At2g28400 Protein of unknown function
At2g36220 Protein of unknown function
At3g03870 Protein of unknown function
Fig. S1. A, ABA-related bud dormancy genes from (1) were compared with BRC1-dependent
genes induced in wild-type but not in brc1 W+FR-treated buds (FDR>0.06, (2)) using Venny (http://
bioinfogp.cnb.csic.es/tools/venny/) B, Gene list showing the 26 BRC1-dependent genes belonging
to the ABA-related GRN. HB21, HB40, HB53 are highlighted in blue. Additional genes analyzed in
this work are highlighted in green.
Fig. S2. Phylogenetic tree of Arabidopsis thaliana Class I HD-Zip genes. Class II HD-Zip HAT2
was included as an outgroup. HD-Zip domain aminoacid sequences were aligned with Clustal
Omega (3) and a phylogenetic tree was obtained with MrBayes (4). Numbers in branches indicate
the probability supporting the branch after generation of 100000 trees. The HB21, HB40, HB53
subclade is highlighted.
A
B
Fig. S3. HB21, HB40 and HB53 mRNA levels correlate with BRC1 levels and with the state of
bud activity. mRNA levels of BRC1, HB21, HB40 and HB53 analysed by qPCR in axillary buds. A,
buds 24 h after decapitation compared with buds of intact plants or plants in which 50 µM NAA was
applied to the decapitated stem. B, buds 8 h after application of 400 µM sucrose or 400 µM
mannitol (control). Error bars are SEM of three biological replicates. Asterisks are significant
differences (T-student, P<0.05) between control and treated plants. Letters denote significant
differences (one-way anova, P<0.05) among means.
Time after estradiol induction
A
75 KDa
α-HA
50 KDa
Ponceau
B
0
2 4
6
8
10 12 14 16 18
C
100
80
% GFP signal
60
40
20
0
0 2 4 6 8 10 12 14 16 18 20 22
Hours after induction
1.4
1.2
1
Relative mRNA levels
0.8
0.6
0.4
0.2
0
HB21 HB40 HB53
Fig. S5. HB21, HB40 and HB53 are not upregulated in estradiol-inducible GUS (GUSind)
control lines after an estradiol treatment. This is the negative control of Figure 1d. Seven-day-
old seedlings carrying an estradiol-inducible construct LEXA:GUS (GUSind). mRNA levels of HB21,
HB40 and HB53 analyzed by qPCR 8 h after a 10 µM estradiol treatment. Error bars are SEM of
three biological replicates.
Fig. S6. Conservation of BRC1 binding motifs in Brasicaceae HB genes. HB21, HB40 and HB53
orthologs in Arabidopsis lyrata, Thellungiella halophila, Caspela rubella and Brassica rapa were
obtained from phytozome (http://phytozome.jgi.doe.gov/) using the full protein sequence from
Arabidopsis thaliana as a query. AlHB21, Alyrata343486; ThHB21, Thhalv10022814; CrHB21,
Carubv10014603; BrHB21, Bra040941; AlHB40, Alyrata490960; ThHB40, Thhalv10026207;
CrHB40, Carubv10006436; BrHB40a, Bra017780; BrHB40b, Bra010573; AlHB53, Alyrata332898;
ThHB53, Thhalv10005612; CtHB53, Carubv10027040; BrHB53, Bra012093. Search for conserved
BRC1 binding sites was performed as in Figure 1f.
Fig. S7. Conservation of BRC1 binding motifs in Brasicaceae HB genes. AlHB21, Alyrata343486;
ThHB21, Thhalv10022814; CrHB21, Carubv10014603; BrHB21, Bra040941; AlHB40,
Alyrata490960; ThHB40, Thhalv10026207; CrHB40, Carubv10006436; BrHB40a, Bra017780;
BrHB40b, Bra010573; AlHB53, Alyrata332898; ThHB53, Thhalv10005612; CtHB53,
Carubv10027040; BrHB53, Bra012093. Alignment of the conserved motifs included in sites 1-6
analyzed by ChIP (Figs. 1f and 1g).
HA:BRC1ind BRC1:GFPind HA:BRC1ind
BRC1:GFPind
GFP:BRC1ind
BRC1ind
BRC1ind
GFP:BRC1ind
GUSind GUSind
Mock 10 µM estradiol
HB40p:GUS
D D E
HB53p:GUS
F G H
HB53p:GUS
Fig. S9. GUS staining of HB40p:GUS and HB53:GUS lines. A, HB40p:GUS flower showing GUS
staining. Close-up of stained HB40p:GUS sepal stomata (B) and pollen grains (C). D, GUS staining
of HB53p:GUS flower pedicel. E, Close-up of stained HB53p:GUS flower pedicel stomata. F to H,
GUS staining of HB53p:GUS lateral roots. GUS expression was noticeable from stages II to VI of
lateral root initiation.
A
B
HB21 1.6 HB40 1.2 HB53 1.4 BRC1
1.4
Relative mRNA levels
1.4 1.2
1
1.2
1.2
1
1 0.8
1
0.8
0.8
0.8 0.6
0.6 0.6
0.6
0.4
0.4 0.4
0.4
0.2 0.2
0.2 0.2
0 * 0
* 0 * * 0
* *
* * * *
* * *
Primer pair: 1 2 3 1 2 3 1 2 3 1 2 3
Fig. S10. HB21, HB40 and HB53 insertional mutants. A, Schematic representation of the
genomic structure of HB21, HB40 and HB53 and T-DNA insertions. Red triangles indicate the
insertion site for each mutation. Rectangles indicate exons; lines, introns. Primers used in qPCR
experiments shown in B and C are indicated. B, mRNA levels of BRC1, HB21, HB40 and HB53
analyzed by qPCR in axillary buds of mutants shown in A. Primer pair 2 was used for HB genes. C,
mRNA levels of BRC1, HB21, HB40 and HB53 analyzed by qPCR in axillary buds of mutants using
primer pairs indicated in A. Error bars are SEM of three biological replicates. Asterisks are
significant differences (T-student, P<0.05) between wild-type and mutant plants.
0.55 0.56 0.67
3 W
A a 0.54 a
0.63
a a
W+FR
Number of RI branches
2.5 ac
2 bc
b
b b
1.5 b
0.5
0
wildwt
type hb21 hb40 hb53-1 hb53-2
ac
2 c
0
wildwt
type hb21
hb214053-1 hb21
hb214053-2
hb40 hb40
hb53-1 hb53-2
C wild type hb21hb40hb53-1
10 10
Number
of
individuals
Reproductive
8 8 Vegetative 3
6 6 Vegetative 2
4 4 Vegetative 1
2 2 Meristem
Empty
0 0
c1 c2 L1 L2 L3 L4 L5 L6 L7 L8 L9 L10 c1 c2 L1 L2 L3 L4 L5 L6 L7 L8 L9 L10
Fig. S11. Shoot branching phenotype of hb21, hb40 and hb53 mutants. A, Number of primary
rosette branches (RI) of wild type and single hb mutants B, Number of secondary cauline branches
(CII) of wild type, hb21 hb40 hb53-1 and hb21hb40hb53-2 mutants. In A and B plants were grown
for 2 weeks after bolting in W or W+FR (N=28). Error bars are SEM. Letters denote significant
differences (one-way anova, P<0.05) among means. Rosette leaf number was undistinguishable in
the mutants and the wild type. C, Developmental stages of buds in the axils of cotyledons (c1 and
c2) and rosette leaves (L1 to L12) of wild-type (left) and hb21hb40hb53-1 (right) individuals
(N=10). Plants were grown in W light and long-day conditions after bolting. Developmental stages
are as defined in (6).
A 4000 HB21
*
300 HB40
*
500 HB53
*
Relative mRNA levels
3500
Control 250 400
3000 Estradiol
200
2500 300
2000 150
1500 200
100
1000
50 100
500
0 0 0
GUSind HB21indHB40indHB53ind GUSind HB21ind HB40ind HB53ind GUSind HB21ind HB40ind HB53ind
Ponceau
Fig. S12. Estradiol-Inducible HB lines. A, mRNA levels of HB21, HB40 and HB53 analyzed by
qPCR in 7-day-old HBind seedlings 8 h after a 10 µM estradiol treatment. Error bar are SEM of
three biological replicates. Asterisks are significant differences (T-student, P<0.05) between control
and treated plants. B, HA:HB protein accumulation analyzed by immunoblot using α-HA antibody
in estradiol-inducible HA:HBsind plants. Protein extracts (15 µg) from 7-day-old seedlings after an 8
h treatment with mock or 10 µM estradiol were separated by 10% SDS-PAGE. HA:HBs were
identified as a 30 KD band. Membranes were probed with α-HA-HRP antibody (1:1000; Roche).
A B
HB21
HB21
HB40
HB40
HB53
HB53
C
-843 TTCCAATTATTGTGT
At3g14440.1
-‐973
-‐843
-‐2747
-‐1849
-973 TCCGAATAATTGATA
-1849 ATTAAATTATTAGAC NCED3
-2747 TCTATATTATTGCTC
-100
-226
D 100 bp
NCED3
TSS
E
*
Fold enrichment
BRC1-no tag
BRC1-GFP
NCED3 ACT2
Fig. S13. HB21, HB40 and HB53 and BRC1 bind the NCED3 promoter. A. Sequence logos
representing the top motifs bound by HB21, HB40 and HB53 in the Arabidopsis genome, as
identified by DAP-seq (7). B. HB21, HB40 and HB53 binding in vitro to the 5´ genomic region of
NCED3 (At3g14440). Coordinates correspond to the distance from TSS to the middle position of
the peaks defined in (7). C. Sequence alignment of motifs found in the middle position of the peaks
shown in B. In red, the sequence corresponding to the logo consensus. D. NCED3 promoter region
analyzed. E. Relative enrichment of GFP:BRC1 binding to the NCED3 5´genomic region shown in
D. ACT2 was used as a negative control. Error bars are SEM. Asterisks are significant differences
(T-student, P<0.05) between untagged BRC1 ind and GFP:BRC1ind lines.
A
Relative mRNA levels
B
3
Control
ABA *
Relative mRNA levels
2
* *
1
0
BRC1 HB21 HB40 HB53
C
BRC1
Relative mRNA levels
0.5
0 wt hb21hb40hb53
Fig. S14. Response of BRC1, HB21, HB40 and HB53 in nced3 mutants and after ABA
application in axillary buds. A, mRNA levels of BRC1, HB21, HB40 and HB53 analyzed by qPCR
in axillary buds of wild-type or nced3 plants treated with W or W+FR for 8 h. B, mRNA levels of
same genes analyzed by qPCR in axillary buds, 8 h after application of mock or 50 µM ABA. C,
mRNA levels of BRC1 analyzed by qPCR in axillary buds of wild-type and hb21 hb40 hb53-1
plants. Error bars are SEM of three biological replicates. Asterisks are significant differences (T-
student, P<0.05) between control and treated plants.
A
wild type brc1
Fold Fold
Gene ID AGI logSignal FDR logSignal FDR
Change Change
ABF3 At4g34000 1.66 10.26 0.0139446 1.03 9.37 0.99997021
AFP3 At3g29575 2.81 10.915 0.02748892 1.755 9.96 0.19698109
GBF3 At2g46270 1.56 9.19 0.04363698 1.44 9.01 0.21433049
NAP At1g69490 3.165 9.345 0.00703062 1.89 8.425 0.33826567
B
C D
Fig. S15. Response of ABF3, AFP3, GBF3 and NAP. A, response to W+FR in wild-type and brc1
mutant buds. Data obtained from (2). Microarray was hybridized with mRNA from wild-type and
brc1 axillary buds treated for 8 h with W or W+FR. Data obtained from six biological replicates.
Fold change is W+FR vs W. logSignal, logarithm of signal intensity. FDR, false discovery rate. Fold
change of microarray data in this experiment was shown to underestimate the fold change
validated by q-RT-PCR (2). In brc1 mutants gene induction is not statistically significant,
FDR>0.05. B, response to W+FR in wild-type and nced3-2 mutant buds. Expression of GBF3 (C)
and AFP3 (D) in 7-day-old GUSind, HA:BRC1ind, HA:HB21ind, HA:HB40ind and HA:HB53ind seedlings
treated for 8 h with 10 µM estradiol. Error bars are SEM of three biological replicates. Asterisks are
significant differences (T-student, P<0.05) between control and treated plants. Letters denote
significant differences (one-way anova, P<0.05) among means.
A
W
W+FR
2 2 1.5 2
a a b
a a ab
1.5 1.5 1.5
1
a a b ab
1 a a 1 1 a
b b a
0.5
0.5 0.5 0.5
0 wt hb21hb40hb53
0 wt hb21hb40hb53
0 wt hb21hb40hb53
0 wt hb21hb40hb53
wild type hb21 wild type hb21 wild type hb21 wild type hb21
hb40 hb40 hb40 hb40
hb53-1 hb53-1 hb53-1 hb53-1
B
Mock
Estradiol
GRAS (At1g63100) KIN12B (At3g23670)
1.4 1.4
RelaGve
mRNA
levels
1.2 1.2
1 1
0.8 0.8
0.6 0.6
0.4 0.4
0.2 0.2
0 GUSind HB21ind HB40ind HB53ind
0 GUSind HB21ind HB40ind HB53ind
1.5
1
1
0.5
0.5
Fig. S16. HB21, HB40 and HB53 do not regulate other BRC1-dependent GRNs. A, mRNA
levels of bud activation genes (2) analyzed by qPCR in axillary buds of wild-type and hb21 hb40
hb53-1 plants treated with W or W+FR for 8 h. B, mRNA levels of same genes analyzed by qPCR
in 7-d-old seedlings of HBind lines treated for 8 hours with 10 µM estradiol. These genes were
downregulated in BRC1ind lines in the same conditions (2). Error bars are SEM of three biological
replicates. Asterisks are significant differences (T-student, P<0.05) between control and treated
plants. Letters denote significant differences (one-way anova, P<0.05) among means.
Fig. S17. Schematic representation of the proposed BRC1 function in the control of bud
dormancy. BRC1 represses (red) and activates (green) several GRNs. BRC1 effect on bud
dormancy would be the additive result of its efect on these GRNs. Those GRNs described in (2),
involved in cell cycle and protein synthesis, are included. In this paper we propose that BRC1
induces an ABA-related GRN at least in part via HB21, HB40 and HB53.
A
Fig. S18. Potential cross-regulations within the GRN. A. Summary of Dap-seq and ChIP-seq
data (7, 8) for the 26 BRC1-dependent co-regulated genes of Fig. S1B. A large proportion of these
genes are bound in vitro and/or in vivo (blue cells) by TFs of the list. No data is available for BRC1,
At1g75490 and At2g36080. B. Number and percentage of genes bound (7, 8) by BRC1-dependent
TFs of Fig. S1B: i) in the Arabidopsis genome, ii) in the 26 genes of Fig. S1B (Core) and iii) in their
co-regulated genes (GRN, Dataset S1). *,** Other binding motifs are found in (8). Most values are
significantly higher than those expected in a random gene list (Pearson X2 test).
Supplemental Materials and Methods
Plant material
The brc1-2 mutant was described (6). The hb21-1 (WiscDsLox468G4), hb40-1
wild-type Col-0 were studied. In all experiments wild-type, single, double and
triple mutant siblings segregating from the same families were compared. The
h dark photoperiod at 22°C in white light (PAR 120 µmol m−2 s−1). For short day
experiments, plants were grown in the same conditions, but with 8h light/16h dark.
Light sources
Light sources were like those described previously (2). W was provided by
was provided by lamp tubes carrying FR 735-nm LEDs (L735-03AU; Epitex). All
Plant treatments
to the light treatment for 8 h. In the decapitation assays, soil-grown plants with 1-2
cm-long main inflorescences were decapitated and the cut stump was inserted into
an inverted 0.5 mL tube with 500 µL of 0.6% (w/v) agarose and 50 µM NAA or
0.005% NaOH 1N (control), samples were collected 24 h after decapitation. For
sugar treatments, 500 µL of 400 mM sucrose or mannitol (control) and 0.1% (v/v)
Tween-20 were added to rosette axillary buds, and samples were collected 8 h
estradiol in DMSO or 0.0005% (v/v) DMSO (control) per plate). They were
collected 8 h after the treatment. For branching responses to ABA, after bolting,
400 µL of 50 µM ABA (SIGMA, A1049), 0.1% (v/v) Tween-20 were added directly
to the base of the stem, covering all the axillary buds, every day and branches were
Phenotypic analyses
transferred to W+FR after bolting. Two weeks after the main inflorescence became
visible, RI branches (shoots longer than 0.5 cm) were counted. For bud
each genotype were dissected and the developmental stage of each axillary bud
development analyses in short days, plants were grown for 8 weeks and rosette
leaves were removed. The number of plants with more than one axillary bud with
leaves longer than 1 mm was counted. The number of rosette leaves (that could
affect the number of buds and branches) was undistinguishable in all the mutants
analysed.
RNA preparation and qPCR analysis
Rosettes leaves with their petioles, main inflorescence, and roots were carefully
dissected out, leaving a tissue highly enriched in axillary buds and subtending
tissue. Material was harvested from at least eight individuals, and three biological
replicates per genotype/treatment and was stored in N2(l). For the estradiol
(2). Three technical replicates (reactions) were done for each biological replicate.
Pairs of primers described in S2 dataset were used. SAND was used as a reference
gene (10). Relative fold expression changes were calculated using qBASE software
(11).
Plasmids generated in this work are listed in S2 dataset. HB21, HB40 and
HB53 coding sequences (CDS) were obtained from axillary bud RNA using
polymerase (Roche) using ATG and STOP primers for each gene (S2 dataset) and
cloned in pGEM-T easy (Promega). For Gateway (Invitrogen) plasmids, PCRs were
(New England Biolabs). PCR fragments were cloned in pDONR207 using Gateway
BP Clonase (Invitrogen) and these plasmids were used as entry vectors to generate
binary vectors.
Generation of transgenic lines
floral dip method (12). At least 6 T3 homozygous lines were generated and
LB culture were induced with 1 mM IPTG, grown at 37°C for 4 h and then
centrifuged at 20000 g and frozen in N2(l). MBP:BRC1 protein extraction and PBM
binding was performed as described (13). The top 10 motifs with the best score
were aligned and a frequency matrix was generated. WebLogo was used to create
the logo for the BRC1 binding site. This frequency matrix was used in Rsat Pattern
Matching, Matrix-scan tool (14). Sites with a score >1 and P<0.05 were
represented.
levels were made relative to the first measurement after estradiol/DMSO addition.
GUS staining
described in (16).
Chromatin immunoprecipitation
Ten-day-old GFP:BRC1ind seedlings grown in 125 mm, MS, 0.7 (w/v) % agar
tissue were collected per sample. ChIP was performed as described (17), using
µMACS Anti-GFP (Miltenyi). Primers used for qPCR are in Dataset S2. Regions of
ACT2 and HSF1 genes without BRC1 binding sites were used as negative controls.
HSF1 was used as a normalizer for DNA quantity and ACT2 as a background
control for immunoprecipitated DNA. The same results were obtained when ACT2
ABA measurements
mixed by shaking during one hour at 4ºC. The extract was kept a -20ºC overnight
and then centrifuged and the supernatant dried in a vacuum evaporator. The dry
residue was dissolved in 1% acetic acid and passed through an Oasis HLB
curves and the Xcalibur 2.2 SP1 build 48 and TraceFinder programs.
DAP-Seq raw data in FASTQ format was obtained from GEO and mapped to
the TAIR10 genome sequence using Bowtie2 (18) with default parameters. Read
coverage across the genome was computed in 10 bp bins using deepTools (19)
Coordinates correspond to the distance from TSS to the middle position of the
enoLOGOS (20).
Supplemental Dataset 1
Name Sequence
attB1-BRC1 ggggacaagtttgtacaaaaaagcaggctTAATGAACAACAACATTTTCA
attB2-BRC1 ggggaccactttgtacaagaaagctgggtTCAATACATGTTTTGATAGTTGTGCATGAGGTC
Xba1-BRC1 GGGGGTCTAGAATGAACAACAACATTTTCAGTACTACTACCACC
BRC1-BamH1 GTGTGGGATCCTCAATACATGTTTTGATAGTTGTGCATGAGGTC
attB1-HA ggggacaagtttgtacaaaaaagcaggctATGGCATACCCATACGACGTTCCG
attB2-BRC1NOSTOP ggggaccactttgtacaagaaagctgggtAATACATGTTTTGATAGTTGTGCATG
attB1-HB21p1k ggggacaagtttgtacaaaaaagcaggctTCCTCTAATCCTCTCATCACA
attB2-HB21p ggggaccactttgtacaagaaagctgggtGTTTCTTCTTTTCGTTCTCT
attB1-HB40p1k ggggacaagtttgtacaaaaaagcaggctATGTATATATATCATTGAATT
attB2-HB40p ggggaccactttgtacaagaaagctgggtGTTCTTGGATCTCTCTCTGTC
attB1-HB53p2k ggggacaagtttgtacaaaaaagcaggctGTTCGTTGCCCACACATTACT
attB2-HB53p ggggaccactttgtacaagaaagctgggtTTTCTCTCTCTAGTTTTTCAG
attB1-HB53p1k ggggacaagtttgtacaaaaaagcaggctCTCGGGAAGATTCGTGGTTAT
attB2-HB21NOSTOP ggggaccactttgtacaagaaagctgggtTCATAAATTGGCTCATCCAGT
attB2-HB40NOSTOP ggggaccactttgtacaagaaagctgggtTTATGTATAGACTCATCCATT
attB2-HB53NOSTOP ggggaccactttgtacaagaaagctgggtTTACATACAAACCATCCCAAT
HB21-ATG ATGAATAACCAGAATGTAGATGA
HB21-STOP TTACATAAATTGGCTCATCC
HB40-ATG ATGAACTACACGGTGGATGAT
HB40-STOP TTATATGTATAGACTCATCCATT
HB53-ATG ATGGATCATGGTAGGTTAATG
HB53-STOP TTATACATACAAACCATCCCA
XhoI-GFP CCGCTCGAGatggtgagcaagggcgaggagctg
BRC1-PacI ccccTTAATTAAtcAATACATGTTTTGATAG
XhoI-BRC1 CCGCTCGAGatgaacaacaacattttcagtactactaccacc
GFP-PacI ccccTTAATTAAttacttgtacagctc
XbaI-HB21 TGCtctagaATGAATAACCAGAATGTAGATGA
HB21-BamHI GCGggatccTTACATAAATTGGCTCATCC
XbaI-HB40 TGCtctagaATGAACTACACGGTGGATGAT
HB40-BglII TCCagatctTTATATGTATAGACTCATCCATT
XbaI-HB53 TGCtctagaATGGATCATGGTAGGTTA
HB53-BamHI GCGggatccTTATACATACAAACCATCCCAAT
HB53-3' AATTTCACTTTGAGCTTCACTC
HB53-5' ACGTCGTACTCGGCCAATG
HB40-3' GGTTCTGAAACCAAACGGCA
HB21 genot5 TGGACACCAAGTGAAACCATCA
HB21 genot3 GGGTCCAAACTTGACCTCCAT
p745 AACGTCCGCAATGTGTTATTAAGTTGTC
LBa1 tggttcacgtagtgggccatcg
o8409 ATATTGACCATCATACTCATTGC
attB1 GUS ggggacaagtttgtacaaaaaagcaggcttaATGTTACGTCCTGTAGAAA
attB2 GUS ggggaccactttgtacaagaaagctgggtTCATTGTTTGCCTCCCTGCT
References