We leverage existing state-of-the-art transcriptomics and proteomics datasets from the GTEx project and the Human Protein Atlas to compute the protein-to-mRNA ratios of 36 human tissues. Using this as a proxy of translational efficiency, we build a machine learning model that identifies codons enriched or depleted in specific tissues. In particular, we detect two clusters of tissues with an opposite pattern of codon preferences. We then use the identified patterns for the development of CUSTOM (https://custom.crg.eu), a codon optimizer algorithm which suggests a synonymous codon design in order to optimize protein production in a tissue-specific manner. In a human cell model, we provide evidence that codon optimization should indeed take into account particularities of the translational machinery of the tissues in which the target proteins are expressed and that our approach can design genes with tissue-optimized expression profiles.