Ron Hause

Large datasets describing the quantitative effects of mutations on protein function are becoming increasingly available. Here, we leverage these datasets to develop Envision, which predicts the magnitude of a missense variant's molecular effect. Envision combines 21,026 variant effect measurements from nine large-scale experimental mutagenesis datasets, a hitherto untapped training resource, with a supervised, stochastic gradient boosting learning algorithm. Envision outperforms other missense… Show more

Large datasets describing the quantitative effects of mutations on protein function are becoming increasingly available. Here, we leverage these datasets to develop Envision, which predicts the magnitude of a missense variant's molecular effect. Envision combines 21,026 variant effect measurements from nine large-scale experimental mutagenesis datasets, a hitherto untapped training resource, with a supervised, stochastic gradient boosting learning algorithm. Envision outperforms other missense variant effect predictors both on large-scale mutagenesis data and on an independent test dataset comprising 2,312 TP53 variants whose effects were measured using a low-throughput approach. This dataset was never used for hyperparameter tuning or model training and thus serves as an independent validation set. Envision prediction accuracy is also more consistent across amino acids than other predictors. Finally, we demonstrate that Envision's performance improves as more large-scale mutagenesis data are incorporated. We precompute Envision predictions for every possible single amino acid variant in human, mouse, frog, zebrafish, fruit fly, worm, and yeast proteomes (https://envision.gs.washington.edu/). Show less

Determining the effect of chemotherapeutic treatment on changes in protein expression can provide important targets for overcoming resistance. Due to challenges in simultaneously measuring large numbers of proteins, a paucity of data exists on global changes. To overcome these challenges, we utilized microwestern arrays that allowed us to measure the abundance and modification state of hundreds of cell signaling and transcription factor proteins in cells following drug exposure. HapMap… Show more

Determining the effect of chemotherapeutic treatment on changes in protein expression can provide important targets for overcoming resistance. Due to challenges in simultaneously measuring large numbers of proteins, a paucity of data exists on global changes. To overcome these challenges, we utilized microwestern arrays that allowed us to measure the abundance and modification state of hundreds of cell signaling and transcription factor proteins in cells following drug exposure. HapMap lymphoblastoid cell lines (LCLs) were exposed to cisplatin, a chemotherapeutic agent commonly used to treat testicular, head and neck, non-small cell lung, and gynecological cancers. We evaluated the expression of 259 proteins following 2, 6, and 12 h of cisplatin treatment in two LCLs with discordant sensitivity to cisplatin. Of these 259 proteins, 66 displayed significantly different protein expression changes (p < 0.05). Fifteen of these proteins were evaluated in a second pair of LCLs with discordant sensitivities to cisplatin; six demonstrated significant differences in expression. We then evaluated a subset of 63 proteins in a second set of LCLs with discordant sensitivity, and 40% of those that were significant in the first pair were also significant in the second part with concordant directionality (p < 0.05). We functionally validated one of the top proteins identified, PDK1, and demonstrated a synergistic relationship between cisplatin and a PDK1 inhibitor in multiple lung cancer lines. This study highlights the potential for identifying novel targets through an understanding of cellular changes in protein expression and modification following drug treatments. Show less

Mutagenesis is a widely-used method for identifying protein positions that are important for function or ligand binding. Advances in high-throughput DNA sequencing and mutagenesis techniques have enabled measurement of the effects of nearly all possible amino acid substitutions in many proteins. The resulting large-scale mutagenesis data sets offer a unique opportunity to draw general conclusions about the effects of different amino acid substitutions. Thus, we analyzed 34,373 mutations in… Show more

Mutagenesis is a widely-used method for identifying protein positions that are important for function or ligand binding. Advances in high-throughput DNA sequencing and mutagenesis techniques have enabled measurement of the effects of nearly all possible amino acid substitutions in many proteins. The resulting large-scale mutagenesis data sets offer a unique opportunity to draw general conclusions about the effects of different amino acid substitutions. Thus, we analyzed 34,373 mutations in fourteen proteins whose effects were measured using large-scale mutagenesis approaches. Methionine was the most tolerated substitution while proline was the least tolerated. We found that several substitutions, including histidine and asparagine, best recapitulated the effects of other substitutions, even when the identity of the wild type amino acid was considered. The effects of histidine and asparagine substitutions also correlated best with the effects of other substitutions in different structural contexts. Furthermore, highly disruptive substitutions like aspartic and glutamic acid had the most discriminatory power for detecting ligand interface positions. Our work highlights the utility of large-scale mutagenesis data, and our conclusions can help guide future single substitution mutational scans. Show less

Colistin sulfate (polymixin E) is an antibiotic prescribed with increasing frequency for severe Gram-negative bacterial infections. As nephrotoxicity is a common side effect, the discovery of pharmacogenomic markers associated with toxicity would benefit the utility of this drug. Our objective was to identify genetic markers of colistin cytotoxicity that were also associated with expression of key proteins using an unbiased, whole genome approach and further evaluate the functional significance… Show more

Colistin sulfate (polymixin E) is an antibiotic prescribed with increasing frequency for severe Gram-negative bacterial infections. As nephrotoxicity is a common side effect, the discovery of pharmacogenomic markers associated with toxicity would benefit the utility of this drug. Our objective was to identify genetic markers of colistin cytotoxicity that were also associated with expression of key proteins using an unbiased, whole genome approach and further evaluate the functional significance in renal cell lines. To this end, we employed International HapMap lymphoblastoid cell lines (LCLs) of Yoruban ancestry with known genetic information to perform a genome-wide association study (GWAS) with cellular sensitivity to colistin. Further association studies revealed that single nucleotide polymorphisms (SNPs) associated with gene expression and protein expression were significantly enriched in SNPs associated with cytotoxicity (p ≤ 0.001 for gene and p = 0.015 for protein expression). The most highly associated SNP, chr18:3417240 (p = 6.49 × 10−8), was nominally a cis-expression quantitative trait locus (eQTL) of the gene TGIF1 (transforming growth factor β (TGFβ)-induced factor-1; p = 0.021) and was associated with expression of the protein HOXD10 (homeobox protein D10; p = 7.17 × 10−5). To demonstrate functional relevance in a murine colistin nephrotoxicity model, HOXD10 immunohistochemistry revealed upregulated protein expression independent of mRNA expression in response to colistin administration. Knockdown of TGIF1 resulted in decreased protein expression of HOXD10 and increased resistance to colistin cytotoxicity. Furthermore, knockdown of HOXD10 in renal cells also resulted in increased resistance to colistin cytotoxicity, supporting the physiological relevance of the initial genomic associations. Show less

Microsatellite instability (MSI), the spontaneous loss or gain of nucleotides from repetitive DNA tracts, is a diagnostic phenotype for gastrointestinal, endometrial, and colorectal tumors, yet the landscape of instability events across a wider variety of cancer types remains poorly understood. To explore MSI across malignancies, we examined 5,930 cancer exomes from 18 cancer types at more than 200,000 microsatellite loci and constructed a genomic classifier for MSI. We identified MSI-positive… Show more

Microsatellite instability (MSI), the spontaneous loss or gain of nucleotides from repetitive DNA tracts, is a diagnostic phenotype for gastrointestinal, endometrial, and colorectal tumors, yet the landscape of instability events across a wider variety of cancer types remains poorly understood. To explore MSI across malignancies, we examined 5,930 cancer exomes from 18 cancer types at more than 200,000 microsatellite loci and constructed a genomic classifier for MSI. We identified MSI-positive tumors in 14 of the 18 cancer types. We also identified loci that were more likely to be unstable in particular cancer types, resulting in specific instability signatures that involved cancer-associated genes, suggesting that instability patterns reflect selective pressures and can potentially identify novel cancer drivers. We also observed a correlation between survival outcomes and the overall burden of unstable microsatellites, suggesting that MSI may be a continuous, rather than discrete, phenotype that is informative across cancer types. These analyses offer insight into conserved and cancer-specific properties of MSI and reveal opportunities for improved methods of clinical MSI diagnosis and cancer gene discovery. Show less

With the advent of massively parallel sequencing, considerable work has gone into adapting chromosome conformation capture (3C) techniques to study chromosomal architecture at a genome-wide scale. We recently demonstrated that the inactive murine X chromosome adopts a bipartite structure using a novel 3C protocol, termed in situ DNase Hi-C. Like traditional Hi-C protocols, in situ DNase Hi-C requires that chromatin be chemically cross-linked, digested, end-repaired, and proximity-ligated with a… Show more

With the advent of massively parallel sequencing, considerable work has gone into adapting chromosome conformation capture (3C) techniques to study chromosomal architecture at a genome-wide scale. We recently demonstrated that the inactive murine X chromosome adopts a bipartite structure using a novel 3C protocol, termed in situ DNase Hi-C. Like traditional Hi-C protocols, in situ DNase Hi-C requires that chromatin be chemically cross-linked, digested, end-repaired, and proximity-ligated with a biotinylated bridge adaptor. The resulting ligation products are optionally sheared, affinity-purified via streptavidin bead immobilization, and subjected to traditional next-generation library preparation for Illumina paired-end sequencing. Importantly, in situ DNase Hi-C obviates the dependence on a restriction enzyme to digest chromatin, instead relying on the endonuclease DNase I. Libraries generated by in situ DNase Hi-C have a higher effective resolution than traditional Hi-C libraries, which makes them valuable in cases in which high sequencing depth is allowed for, or when hybrid capture technologies are expected to be used. The protocol described here, which involves ∼4 d of bench work, is optimized for the study of mammalian cells, but it can be broadly applicable to any cell or tissue of interest, given experimental parameter optimization. Show less

Interpreting variants of uncertain significance (VUS) is a central challenge in medical genetics. One approach is to experimentally measure the functional consequences of VUS, but to date this approach has been post hoc and low-throughput. Here we use massively parallel assays to measure the effects of nearly 2,000 missense substitutions in the RING domain of BRCA1 on its E3 ubiquitin ligase activity and its binding to the BARD1 RING domain. From the resulting scores, we generate a model to… Show more

Interpreting variants of uncertain significance (VUS) is a central challenge in medical genetics. One approach is to experimentally measure the functional consequences of VUS, but to date this approach has been post hoc and low-throughput. Here we use massively parallel assays to measure the effects of nearly 2,000 missense substitutions in the RING domain of BRCA1 on its E3 ubiquitin ligase activity and its binding to the BARD1 RING domain. From the resulting scores, we generate a model to predict the capacities of full-length BRCA1 variants to support homology-directed DNA repair, the essential role of BRCA1 in tumor suppression, and show that it outperforms widely used biological-effect prediction algorithms. We envision that massively parallel functional assays may facilitate the prospective interpretation of variants observed in clinical sequencing. Show less

Cell Preview for "Genetic Variation in Human DNA Replication Timing" (Koren et al. 2014)

Genome replication programs are highly orchestrated. In this issue, Koren and colleagues leverage whole-genome sequencing data to discover that DNA replication timing patterns differ between individuals and are associated with genetic variants. These so-called rtQTLs represent a new form of quantitative trait loci that may influence gene dosage and mutational frequency and provide mechanistic… Show more

Cell Preview for "Genetic Variation in Human DNA Replication Timing" (Koren et al. 2014)

Genome replication programs are highly orchestrated. In this issue, Koren and colleagues leverage whole-genome sequencing data to discover that DNA replication timing patterns differ between individuals and are associated with genetic variants. These so-called rtQTLs represent a new form of quantitative trait loci that may influence gene dosage and mutational frequency and provide mechanistic insights into the regulation of DNA replication. Show less

Saturation mutagenesis—coupled to an appropriate biological assay—represents a fundamental means of achieving a high-resolution understanding of regulatory and protein-coding nucleic acid sequences of interest. However, mutagenized sequences introduced in trans on episomes or via random or “safe-harbour” integration fail to capture the native context of the endogenous chromosomal locus. This shortcoming markedly limits the interpretability of the resulting measurements of mutational impact… Show more

Saturation mutagenesis—coupled to an appropriate biological assay—represents a fundamental means of achieving a high-resolution understanding of regulatory and protein-coding nucleic acid sequences of interest. However, mutagenized sequences introduced in trans on episomes or via random or “safe-harbour” integration fail to capture the native context of the endogenous chromosomal locus. This shortcoming markedly limits the interpretability of the resulting measurements of mutational impact. Here, we couple CRISPR/Cas9 RNA-guided cleavage with multiplex homology-directed repair using a complex library of donor templates to demonstrate saturation editing of genomic regions. In exon 18 of BRCA1, we replace a six-base-pair (bp) genomic region with all possible hexamers, or the full exon with all possible single nucleotide variants (SNVs), and measure strong effects on transcript abundance attributable to nonsense-mediated decay and exonic splicing elements. We similarly perform saturation genome editing of a well-conserved coding region of an essential gene, DBR1, and measure relative effects on growth that correlate with functional impact. Measurement of the functional consequences of large numbers of mutations with saturation genome editing will potentially facilitate high-resolution functional dissection of both cis-regulatory elements and trans-acting factors, as well as the interpretation of variants of uncertain significance observed in clinical sequencing. Show less

Many genetic variants associated with human disease have been found to be associated with alterations in mRNA expression. Although it is commonly assumed that mRNA expression changes will lead to consequent changes in protein levels, methodological challenges have limited our ability to test the degree to which this assumption holds true. Here, we further developed the micro-western array approach and globally examined relationships between human genetic variation and cellular protein levels… Show more

Many genetic variants associated with human disease have been found to be associated with alterations in mRNA expression. Although it is commonly assumed that mRNA expression changes will lead to consequent changes in protein levels, methodological challenges have limited our ability to test the degree to which this assumption holds true. Here, we further developed the micro-western array approach and globally examined relationships between human genetic variation and cellular protein levels. We collected more than 250,000 protein level measurements comprising 441 transcription factor and signaling protein isoforms across 68 Yoruba (YRI) HapMap lymphoblastoid cell lines (LCLs) and identified 12 cis and 160 trans protein level QTLs (pQTLs) at a false discovery rate (FDR) of 20%. Whereas up to two thirds of cis mRNA expression QTLs (eQTLs) were also pQTLs, many pQTLs were not associated with mRNA expression. Notably, we replicated and functionally validated a trans pQTL relationship between the KARS lysyl-tRNA synthetase locus and levels of the DIDO1 protein. This study demonstrates proof of concept in applying an antibody-based microarray approach to iteratively measure the levels of human proteins and relate these levels to human genome variation and other genomic data sets. Our results suggest that protein-based mechanisms might functionally buffer genetic alterations that influence mRNA expression levels and that pQTLs might contribute phenotypic diversity to a human population independently of influences on mRNA expression. Show less

Many human diseases are associated with aberrant regulation of phosphoprotein signaling networks. Src homology 2 (SH2) domains represent the major class of protein domains in metazoans that interact with proteins phosphorylated on the amino acid residue tyrosine. Although current SH2 domain prediction algorithms perform well at predicting the sequences of phosphorylated peptides that are likely to result in the highest possible interaction affinity in the context of random peptide library… Show more

Many human diseases are associated with aberrant regulation of phosphoprotein signaling networks. Src homology 2 (SH2) domains represent the major class of protein domains in metazoans that interact with proteins phosphorylated on the amino acid residue tyrosine. Although current SH2 domain prediction algorithms perform well at predicting the sequences of phosphorylated peptides that are likely to result in the highest possible interaction affinity in the context of random peptide library screens, these algorithms do poorly at predicting the interaction potential of SH2 domains with physiologically derived protein sequences. We employed a high-throughput interaction assay system to empirically determine the affinity between 93 human SH2 domains and phosphopeptides abstracted from several RTKs and signaling proteins. The resulting interaction experiments revealed over 1000 novel peptide-protein interactions and provided a glimpse into the common and specific interaction potentials of c-Met, c-Kit, GAB1, and the human androgen receptor. We used this data to build a permutation- based logistic (PEBL) regression classifier that performed considerably better than existing algorithms for predicting the interaction potential of several SH2 domains. Show less

Annotating and interpreting the results of genome-wide association studies (GWAS) remains challenging. Assigning function to genetic variants as expression quantitative trait loci is an expanding and useful approach, but focuses exclusively on mRNA rather than protein levels. Many variants remain without annotation. To address this problem, we measured the steady state abundance of 441 human signaling and transcription factor proteins from 68 Yoruba HapMap lymphoblastoid cell lines to identify… Show more

Annotating and interpreting the results of genome-wide association studies (GWAS) remains challenging. Assigning function to genetic variants as expression quantitative trait loci is an expanding and useful approach, but focuses exclusively on mRNA rather than protein levels. Many variants remain without annotation. To address this problem, we measured the steady state abundance of 441 human signaling and transcription factor proteins from 68 Yoruba HapMap lymphoblastoid cell lines to identify novel relationships between inter-individual protein levels, genetic variants, and sensitivity to chemotherapeutic agents. Proteins were measured using micro-western and reverse phase protein arrays from three independent cell line thaws to permit mixed effect modeling of protein biological replicates. We observed enrichment of protein quantitative trait loci (pQTLs) for cellular sensitivity to two commonly used chemotherapeutics: cisplatin and paclitaxel. We functionally validated the target protein of a genome-wide significant trans-pQTL for its relevance in paclitaxel-induced apoptosis. GWAS overlap results of drug-induced apoptosis and cytotoxicity for paclitaxel and cisplatin revealed unique SNPs associated with the pharmacologic traits (at p<0.001). Interestingly, GWAS SNPs from various regions of the genome implicated the same target protein (p<0.0001) that correlated with drug induced cytotoxicity or apoptosis (p≤0.05). Two genes were functionally validated for association with drug response using siRNA: SMC1A with cisplatin response and ZNF569 with paclitaxel response. This work allows pharmacogenomic discovery to progress from the transcriptome to the proteome and offers potential for identification of new therapeutic targets. This approach, linking targeted proteomic data to variation in pharmacologic response, can be generalized to other studies evaluating genotype-phenotype relationships and provide insight into chemotherapeutic mechanisms. Show less

To investigate the relationship between a protein's sequence and its biophysical properties, we studied the effects of more than 100 mutations in Avena sativa light-oxygen-voltage domain 2, a model protein of the Per-Arnt-Sim family. The A. sativa light-oxygen-voltage domain 2 undergoes a photocycle with a conformational change involving the unfolding of the terminal helices. Whereas selection studies typically search for winners in a large population and fail to characterize many sites, we… Show more

To investigate the relationship between a protein's sequence and its biophysical properties, we studied the effects of more than 100 mutations in Avena sativa light-oxygen-voltage domain 2, a model protein of the Per-Arnt-Sim family. The A. sativa light-oxygen-voltage domain 2 undergoes a photocycle with a conformational change involving the unfolding of the terminal helices. Whereas selection studies typically search for winners in a large population and fail to characterize many sites, we characterized the biophysical consequences of mutations throughout the protein using NMR, circular dichroism, and ultraviolet/visible spectroscopy. Despite our intention to introduce highly disruptive substitutions, most had modest or no effect on function, and many could even be considered to be more photoactive. Substitutions at evolutionarily conserved sites can have minimal effect, whereas those at nonconserved positions can have large effects, contrary to the view that the effects of mutations, especially at conserved positions, are predictable. Using predictive models, we found that the effects of mutations on biophysical function and allostery reflect a complex mixture of multiple characteristics including location, character, electrostatics, and chemistry. Show less

We developed a large-scale fluorescence polarization (FP) methodology that was able to characterize interactions between SH2 domains and ErbB receptor phosphosites with higher fidelity and sensitivity than was previously achieved with PMs. We used the FP assay to query the interaction of synthetic phosphopeptides corresponding to 89 ErbB receptor intracellular tyrosine sites against 93 human SH2 domains and 2 phosphotyrosine binding (PTB) domains. From 358,944 polarization measurements, the… Show more

We developed a large-scale fluorescence polarization (FP) methodology that was able to characterize interactions between SH2 domains and ErbB receptor phosphosites with higher fidelity and sensitivity than was previously achieved with PMs. We used the FP assay to query the interaction of synthetic phosphopeptides corresponding to 89 ErbB receptor intracellular tyrosine sites against 93 human SH2 domains and 2 phosphotyrosine binding (PTB) domains. From 358,944 polarization measurements, the affinities for 1,405 unique biological interactions were determined, 83% of which are novel. In contrast to data from previous reports, our analyses suggested that ErbB2 was not more promiscuous than the other ErbB receptors. Our results showed that each receptor displays unique preferences in the affinity and location of recruited SH2 domains that may contribute to differences in downstream signaling potential. ErbB1 was enriched versus the other receptors for recruitment of domains from RAS GEFs whereas ErbB2 was enriched for recruitment of domains from tyrosine and phosphatidyl inositol phosphatases. ErbB3, the kinase inactive ErbB receptor family member, was predictably enriched for recruitment of domains from phosphatidyl inositol kinases and surprisingly, was enriched for recruitment of domains from tyrosine kinases, cytoskeletal regulatory proteins, and RHO GEFs but depleted for recruitment of domains from phosphatidyl inositol phosphatases. Many novel interactions were also observed with phosphopeptides corresponding to ErbB receptor tyrosines not previously reported to be phosphorylated by mass spectrometry, suggesting the existence of many biologically relevant RTK sites that may be phosphorylated but below the detection threshold of standard mass spectrometry procedures. This dataset represents a rich source of testable hypotheses regarding the biological mechanisms of ErbB receptors. Show less

Caffeic acid phenethyl ester (CAPE) is a bioactive component derived from honeybee hive propolis. CAPE has been shown to have antimitogenic, anticarcinogenic, and other beneficial medicinal properties. Many of its effects have been shown to be mediated through its inhibition of NF-κB signaling pathways. We took a systematic approach to uncover the effects of CAPE from hours to days on the signaling networks in human prostate cancer cells. We observed that CAPE dosage dependently suppressed the… Show more

Caffeic acid phenethyl ester (CAPE) is a bioactive component derived from honeybee hive propolis. CAPE has been shown to have antimitogenic, anticarcinogenic, and other beneficial medicinal properties. Many of its effects have been shown to be mediated through its inhibition of NF-κB signaling pathways. We took a systematic approach to uncover the effects of CAPE from hours to days on the signaling networks in human prostate cancer cells. We observed that CAPE dosage dependently suppressed the proliferation of LNCaP, DU-145, and PC-3 human prostate cancer cells. Administration of CAPE by gavage significantly inhibited the tumor growth of LNCaP xenografts in nude mice. Using LNCaP cells as a model system, we examined the effect of CAPE on gene expression, protein signaling, and transcriptional regulatory networks using micro-Western arrays and PCR arrays. We built a model of the impact of CAPE on cell signaling which suggested that it acted through inhibition of Akt-related protein signaling networks. Overexpression of Akt1 or c-Myc, a downstream target of Akt signaling, significantly blocked the antiproliferative effects of CAPE. In summary, our results suggest that CAPE administration may be useful as an adjuvant therapy for prostate and potentially other types of cancers that are driven by the p70S6K and Akt signaling networks. Show less

Leaf growth and development determines a plant’s capacity for photosynthesis and carbon fixation. These morphological traits are the integration of genetic and environmental factors through time. Yet fine dissection of the developmental genetic basis of leaf expansion throughout a growing season is difficult, due to the complexity of the trait and the need for real time measurement. In this study, we developed a time-lapse image analysis approach, which traces leaf expansion under seasonal… Show more

Leaf growth and development determines a plant’s capacity for photosynthesis and carbon fixation. These morphological traits are the integration of genetic and environmental factors through time. Yet fine dissection of the developmental genetic basis of leaf expansion throughout a growing season is difficult, due to the complexity of the trait and the need for real time measurement. In this study, we developed a time-lapse image analysis approach, which traces leaf expansion under seasonal light variation. Three growth traits, rosette leaf area, circular area, and their ratio as compactness, were measured and normalized on a linear timescale to control for developmental heterogeneity. We found high heritability for all growth traits that changed over time. Our study highlights a cost-effective, high-throughput phenotyping approach that facilitates the dissection of genetic basis of plant shoot growth and development under dynamic environmental conditions. Show less

While proteomic methods have illuminated many areas of biological protein space, many fundamental questions remain with regard to systems-level relationships between mRNAs, proteins, and cell behaviors. While mass spectrometric methods offer a panoramic picture of the relative expression and modification of large numbers of proteins, they are neither optimal for the analysis of pre-defined targets across large numbers of samples nor for assessing differences in proteins between individual cells… Show more

While proteomic methods have illuminated many areas of biological protein space, many fundamental questions remain with regard to systems-level relationships between mRNAs, proteins, and cell behaviors. While mass spectrometric methods offer a panoramic picture of the relative expression and modification of large numbers of proteins, they are neither optimal for the analysis of pre-defined targets across large numbers of samples nor for assessing differences in proteins between individual cells or cell compartments. Conversely, traditional antibody-based methods are effective at sensitively analyzing small numbers of proteins across small numbers of conditions, and can be used to analyze relative differences in protein abundance and modification between cells and cell compartments. However, traditional antibody-based approaches are not optimal for analyzing large numbers of protein abundances and modifications across many samples. In this perspective article, we will review recent advances in methodologies and philosophies behind several microarray-based, intermediate-level, “protein-omic” methods including a focus on reverse phase lysate arrays and micro-western arrays that have been helpful for bridging gaps between large- and small-scale protein analysis approaches and have provided insight into the roles that protein systems play in several biological processes. Show less