We have developed a Bayesian framework to infer phosphorylation networks from time series measurements of phosphosite concentrations upon ligand stimulation. To increase the prediction accuracy we integrated different types of data, e.g., amino acid sequence data, genomic context data (gene fusion, gene neighborhood, and phylogentic profiles)