library_construction_protocol: Dneasy®  Plant  Maxi Kit (Qiagen) The experiment was carried out as described by Bartlett et al (2017). Briefly, leaf genomic DNA (gDNA) was extracted from tomato AC leaf , then the gDNA was fragmented into around 200bp fragments by Covaris S2 ultrasonicator; End-repair and A-tailing was performed sequentially, after that the truncated Illumina TruSeq adaptor was ligated and gDNA library was prepared; Meanwhile, full-length CDS of MBD5 was cloned into pIX-HALO tagged vector, and the fusion protein was expressed in vitro via TNT Coupled Wheat Germ Extract Systems; then the fusion protein was purified byMagne HaloTag Beads; Finally, beads coupled with protein were incubated with the gDNA library, then clean unbound gDNA fragments, MBD5 specific bound gDNA was purified and amplified for sequencing.