This paper deals with the problem of the aggregation of implication functions. After characterizing those binary operations which merge two implication functions into a new implication function, we study some aspects of the aggregated implication functions in terms of the two given ones. In particular, we are interested in aggregation processes such that properties of the implication functions are preserved.