Search

Abstract Background Metastatic relapse of prostate cancer after radiotherapy is a significant cause of prostate cancer-related morbidity and mortality. PLOD2 is a mediator of invasion and metastasis that we identified as being upregulated in our highly aggressive radiorecurrent prostate cancer cell line. Methods Patient dataset analysis was conducted using a variety of prostate cancer cohorts. Prostate cancer cell lines were treated with siRNA, or the drug PX-478 prior to in vitro invasion, migration, or in vivo chick embryo (CAM) extravasation assay. Protein levels were detected by western blot. For RNA analysis, RNA sequencing was conducted on PLOD2 knockdown cells and validated by qRT-PCR. Results PLOD2 is a negative prognostic factor associated with biochemical relapse, driving invasion, migration, and extravasation in radiorecurrent prostate cancer. Mechanistically, HIF1α upregulation drives PLOD2 expression in our radiorecurrent prostate cancer cells, which is effectively inhibited by HIF1α inhibitor PX-478 to reduce invasion, migration, and extravasation. Finally, the long non-coding RNA LNCSRLR acts as a promoter of invasion downstream of PLOD2. Conclusions Together, our results demonstrate for the first time the role of PLOD2 in radiorecurrent prostate cancer invasiveness, and point towards its potential as a therapeutic target to reduce metastasis and improve survival outcomes in prostate cancer patients.

Abstract Urine is a complex biofluid that reflects both overall physiologic state and the state of the genitourinary tissues through which it passes. It contains both secreted proteins and proteins encapsulated in tissue-derived extracellular vesicles (EVs). To understand the population variability and clinical utility of urine, we quantified the secreted and EV proteomes from 190 men, including a subset with prostate cancer. We demonstrate that a simple protocol enriches prostatic proteins in urine. Secreted and EV proteins arise from different subcellular compartments. Urinary EVs are faithful surrogates of tissue proteomes, but secreted proteins in urine or cell line EVs are not. The urinary proteome is longitudinally stable over several years. It can accurately and non-invasively distinguish malignant from benign prostatic lesions and can risk-stratify prostate tumors. This resource quantifies the complexity of the urinary proteome and reveals the synergistic value of secreted and EV proteomes for translational and biomarker studies.

Abstract Chordomas are clinically aggressive tumors with a high rate of disease progression despite maximal therapy. Given the limited therapeutic options available, there remains an urgent need for the development of novel therapies to improve clinical outcomes. Cell surface proteins are attractive therapeutic targets yet are challenging to profile with common methods. Four chordoma cell lines were analyzed by quantitative proteomics using a differential ultracentrifugation organellar fractionation approach. A subtractive proteomics strategy was applied to select proteins that are plasma membrane enriched. Systematic data integration prioritized PLA2R1 (secretory phospholipase A2 receptor–PLA2R1) as a chordoma-enriched surface protein. The expression profile of PLA2R1 was validated across chordoma cell lines, patient surgical tissue samples, and normal tissue lysates via immunoblotting. PLA2R1 expression was further validated by immunohistochemical analysis in a richly annotated cohort of 25-patient tissues. Immunohistochemistry analysis revealed that elevated expression of PLA2R1 is correlated with poor prognosis. Using siRNA- and CRISPR/Cas9-mediated knockdown of PLA2R1, we demonstrated significant inhibition of 2D, 3D and in vivo chordoma growth. PLA2R1 depletion resulted in cell cycle defects and metabolic rewiring via the MAPK signaling pathway, suggesting that PLA2R1 plays an essential role in chordoma biology. We have characterized the proteome of four chordoma cell lines and uncovered PLA2R1 as a novel cell-surface protein required for chordoma cell survival and association with patient outcome.

Summary: Deep mass spectrometry-based proteomic profiling of rare cell populations has been constrained by sample input requirements. Here, we present a protocol for droplet-based one-pot preparation for proteomic samples (DROPPS), an accessible low-input platform that generates high-fidelity proteomic profiles of 100–2,500 cells. We describe steps for depositing cellular material, cell lysis, and digesting proteins in the same microliter-droplet well. We anticipate DROPPS will accelerate biology-driven proteomic research for a multitude of rare cell populations.For complete details on the use and execution of this protocol, please refer to Waas et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Abstract Medulloblastoma (MB) is the most common type of malignant pediatric brain cancer. The current standard of care (SOC) involves maximal safe resection and chemoradiotherapy in individuals older than 3 years, often leading to devastating neurocognitive and developmental deficits. Out of the four distinct molecular subgroups, Group 3 and 4 have the poorest patient outcomes due to the aggressive nature of the tumor and propensity to metastasize and recur post therapy. The toxicity of the SOC and lack of response in specific subtypes to the SOC underscores the urgent need for developing and translating novel treatment options including immunotherapies. To identify differentially enriched surface proteins that could be evaluated for potential future immunotherapeutic interventions, we leveraged N-glycocapture surfaceome profiling on Group 3 MB cells from primary tumor, through therapy, to recurrence using our established therapy-adapted patient derived xenograft model. Integrin 𝛼5 (ITGA5) was one of the most differentially enriched targets found at recurrence when compared to engraftment and untreated timepoints. In addition to being enriched at recurrence, shRNA-mediated knockdown and small molecule inhibition of ITGA5 have resulted in marked decrease in proliferation and self-renewal in vitro and demonstrated a survival advantage in vivo. Together, our data highlights the value of dynamic profiling of cells as they evolve through therapy and the identification of ITGA5 as a promising therapeutic target for recurrent Group 3 MB.

The intercalated disc (ICD) is a membrane structure of the cardiac muscle involved in normal heart function. Here the authors report that knockdown of the ICD-bound transmembrane protein 65 results in impaired ICD structure, abnormal cardiac electrophysiology and cardiomyopathy in mice.

Abstract Multiparametric magnetic resonance imaging (mpMRI) is an emerging standard for diagnosing and prognosing prostate cancer, but ~ 20% of clinically significant tumors are invisible to mpMRI, as defined by the Prostate Imaging Reporting and Data System version 2 (PI-RADSv2) score of one or two. To understand the biological underpinnings of tumor visibility on mpMRI, we examined the proteomes of forty clinically significant tumors (i.e., International Society of Urological Pathology (ISUP) Grade Group 2)—twenty mpMRI-visible and twenty mpMRI-invisible, with matched histologically normal prostate. Normal prostate tissue was indistinguishable between patients with visible and invisible tumors, and invisible tumors closely resembled the normal prostate. These data indicate that mpMRI-visibility arises when tumor evolution leads to large-magnitude proteomic divergences from histologically normal prostate.

Abstract The development of functionally selective or biased ligands is a promising approach towards drugs with less side effects. Biased ligands for G protein-coupled receptors can selectively induce G protein activation or β-arrestin recruitment. The consequences of this selective action on cellular functions, however, are not fully understood. Here, we investigated the impact of five biased and balanced dopamine D2 receptor agonists and antagonists on the global protein expression in HEK293T cells by untargeted nanoscale liquid chromatography–tandem mass spectrometry. The proteome analysis detected 5290 protein groups. Hierarchical clustering and principal component analysis based on the expression levels of 1462 differential proteins led to a separation of antagonists and balanced agonist from the control treatment, while the biased ligands demonstrated larger similarities to the control. Functional analysis of affected proteins revealed that the antagonists haloperidol and sulpiride regulated exocytosis and peroxisome function. The balanced agonist quinpirole, but not the functionally selective agonists induced a downregulation of proteins involved in synaptic signaling. The β-arrestin-preferring agonist BM138, however, regulated several proteins related to neuron function and the dopamine receptor-mediated signaling pathway itself. The G protein-selective partial agonist MS308 influenced rather broad functional terms such as DNA processing and mitochondrial translation.

The impact of variant calling algorithms on the analysis of intra-tumour heterogeneity has not been properly quantified. Here the authors measure the variability of 22 pipelines with different variant callers and clustering algorithms for subclonal reconstruction to inform future analyses.

Abstract The isolation and subsequent molecular analysis of extracellular vesicles (EVs) derived from patient samples is a widely used strategy to understand vesicle biology and to facilitate biomarker discovery. Expressed prostatic secretions in urine are a tumor proximal fluid that has received significant attention as a source of potential prostate cancer (PCa) biomarkers for use in liquid biopsy protocols. Standard EV isolation methods like differential ultracentrifugation (dUC) co‐isolate protein contaminants that mask lower‐abundance proteins in typical mass spectrometry (MS) protocols. Further complicating the analysis of expressed prostatic secretions, uromodulin, also known as Tamm‐Horsfall protein (THP), is present at high concentrations in urine. THP can form polymers that entrap EVs during purification, reducing yield. Disruption of THP polymer networks with dithiothreitol (DTT) can release trapped EVs, but smaller THP fibres co‐isolate with EVs during subsequent ultracentrifugation. To resolve these challenges, we describe here a dUC method that incorporates THP polymer reduction and alkaline washing to improve EV isolation and deplete both THP and other common protein contaminants. When applied to human expressed prostatic secretions in urine, we achieved relative enrichment of known prostate and prostate cancer‐associated EV‐resident proteins. Our approach provides a promising strategy for global proteomic analyses of urinary EVs.

Distributed multi-omic digitization of clinical specimen across multiple sites is a prerequisite for turning molecular precision medicine into reality. Here, the authors show that coordinated proteotype data acquisition is feasible using standardized MS data acquisition and analysis strategies.

Abstract High-grade serous carcinoma (HGSC) is the most prevalent and aggressive subtype of ovarian cancer. The large degree of clinical heterogeneity within HGSC has justified deviations from the traditional one-size-fits-all clinical management approach. However, the majority of HGSC patients still relapse with chemo-resistant cancer and eventually succumb to their disease, evidence that further work is needed to improve patient outcomes. Advancements in high-throughput technologies have enabled novel insights into biological complexity, offering a large potential for informing precision medicine efforts. Here, we review the current landscape of clinical management for HGSC and highlight applications of high-throughput biological approaches for molecular subtyping and the discovery of putative blood-based biomarkers and novel therapeutic targets. Additionally, we present recent improvements in model systems and discuss how their intersection with high-throughput platforms and technological advancements is positioned to accelerate the realization of precision medicine in HGSC.

The sarcoplasmic (SR) and endoplasmic reticulum (ER) are involved in heart development but how this arises is unclear. Here, the authors show that loss of a SR/ER protein REEP5 causes membrane destabilization and decreased cardiac myocyte contractility, with cardiac dysfunction in mutant mouse and zebrafish models.

The tumor microenvironment has recently emerged as a critical component of high-grade serous ovarian cancer (HGSC) disease progression. Specifically, cancer-associated fibroblasts (CAFs) have been recognized as key players in various pro-oncogenic processes. Here, we use mass-spectrometry (MS) to characterize the proteomes of HGSC patient-derived CAFs and compare them to those of the epithelial component of HGSC to gain a deeper understanding into their tumor-promoting phenotype. We integrate our data with primary tissue data to define a proteomic signature of HGSC CAFs and uncover multiple novel CAF proteins that are prognostic in an independent HGSC patient cohort. Our data represent the first MS-based global proteomic characterization of CAFs in HGSC and further highlights the clinical significance of HGSC CAFs.

Pancreatic progenitors (PPs) can be derived from human pluripotent stem cells in vitro but efficiency of differentiation varies, making it hard to sort for insulin-producing cells. Here, the authors use a proteomic approach to identify the secretory granule membrane glycoprotein 2 as a marker for PDX1+/NKX6-1+ PPs.

Abstract Background Onco-proteogenomics aims to understand how changes in a cancer’s genome influences its proteome. One challenge in integrating these molecular data is the identification of aberrant protein products from mass-spectrometry (MS) datasets, as traditional proteomic analyses only identify proteins from a reference sequence database. Methods We established proteomic workflows to detect peptide variants within MS datasets. We used a combination of publicly available population variants (dbSNP and UniProt) and somatic variations in cancer (COSMIC) along with sample-specific genomic and transcriptomic data to examine proteome variation within and across 59 cancer cell-lines. Results We developed a set of recommendations for the detection of variants using three search algorithms, a split target-decoy approach for FDR estimation, and multiple post-search filters. We examined 7.3 million unique variant tryptic peptides not found within any reference proteome and identified 4771 mutations corresponding to somatic and germline deviations from reference proteomes in 2200 genes among the NCI60 cell-line proteomes. Conclusions We discuss in detail the technical and computational challenges in identifying variant peptides by MS and show that uncovering these variants allows the identification of druggable mutations within important cancer genes.

MHC class I and II molecules generally present endogenous and exogenous peptides, respectively, through distinct mechanisms. Here, the authors show that the class II molecule HLA-DP84Glyuses both class I and II mechanisms to constitutively present peptides.

The last decade has seen a sharp increase in the number of scientific publications describing physiological and pathological functions of extracellular vesicles (EVs), a collective term covering various subtypes of cell-released, membranous structures, called exosomes, microvesicles, microparticles, ectosomes, oncosomes, apoptotic bodies, and many other names. However, specific issues arise when working with these entities, whose size and amount often make them difficult to obtain as relatively pure preparations, and to characterize properly. The International Society for Extracellular Vesicles (ISEV) proposed Minimal Information for Studies of Extracellular Vesicles (“MISEV”) guidelines for the field in 2014. We now update these “MISEV2014” guidelines based on evolution of the collective knowledge in the last four years. An important point to consider is that ascribing a specific function to EVs in general, or to subtypes of EVs, requires reporting of specific information beyond mere description of function in a crude, potentially contaminated, and heterogeneous preparation. For example, claims that exosomes are endowed with exquisite and specific activities remain difficult to support experimentally, given our still limited knowledge of their specific molecular machineries of biogenesis and release, as compared with other biophysically similar EVs. The MISEV2018 guidelines include tables and outlines of suggested protocols and steps to follow to document specific EV-associated functional activities. Finally, a checklist is provided with summaries of key points.

Proteomic technologies are capable of identifying thousands of proteins in biological samples, but biomarker applications are lagging. Here the authors use Multiple Reaction Monitoring Mass Spectrometry to delineate peptide signatures that accurately distinguish between defined prostate cancer patient risk groups.

The contribution of basal and luminal cells to cancer progression and metastasis is poorly understood. We report generation of reporter systems driven by either keratin-14 (K14) or keratin-8 (K8) promoter that not only express a fluorescent protein but also an inducible suicide gene. Transgenic mice express the reporter genes in the right cell compartments of mammary gland epithelia and respond to treatment with toxins. In addition, we engineered the reporters into 4T1 metastatic mouse tumor cell line and demonstrate that K14+ cells, but not K14- or K8+, are both highly invasive in three-dimensional (3D) culture and metastatic in vivo. Treatment of cells in culture, or tumors in mice, with reporter-targeting toxin inhibited both invasive behavior and metastasis in vivo. RNA sequencing (RNA-seq), secretome, and epigenome analysis of K14+ and K14- cells led to the identification of amphoterin-induced protein 2 (Amigo2) as a new cell invasion driver whose expression correlated with decreased relapse-free survival in patients with TP53 wild-type (WT) breast cancer.

Bioinformatics tools are imperative for the in depth analysis of heterogeneous high-throughput data. Most of the software tools are developed by specific laboratories or groups or companies wherein they are designed to perform the required analysis for the group. However, such software tools may fail to capture “what the community needs in a tool”. Here, we describe a novel community-driven approach to build a comprehensive functional enrichment analysis tool. Using the existing FunRich tool as a template, we invited researchers to request additional features and/or changes. Remarkably, with the enthusiastic participation of the community, we were able to implement 90% of the requested features. FunRich enables plugin for extracellular vesicles wherein users can download and analyse data from Vesiclepedia database. By involving researchers early through community needs software development, we believe that comprehensive analysis tools can be developed in various scientific disciplines.

Alternative splicing acts on transcripts from almost all human multi-exon genes. Notwithstanding its ubiquity, fundamental ramifications of splicing on protein expression remain unresolved. The number and identity of spliced transcripts that form stably folded proteins remain the sources of considerable debate, due largely to low coverage of experimental methods and the resulting absence of negative data. We circumvent this issue by developing a semi-supervised learning algorithm, positive unlabeled learning for splicing elucidation (PULSE; http://www.kimlab.org/software/pulse), which uses 48 features spanning various categories. We validated its accuracy on sets of bona fide protein isoforms and directly on mass spectrometry (MS) spectra for an overall AU-ROC of 0.85. We predict that around 32% of “exon skipping” alternative splicing events produce stable proteins, suggesting that the process engenders a significant number of previously uncharacterized proteins. We also provide insights into the distribution of positive isoforms in various functional classes and into the structural effects of alternative splicing.

Condensin is a central regulator of mitotic genome structure with mutants showing poorly condensed chromosomes and profound segregation defects. Here, we identify NCT, a complex comprising the Nrc1 BET-family tandem bromodomain protein (SPAC631.02), casein kinase II (CKII), and several TAFs, as a regulator of condensin function. We show that NCT and condensin bind similar genomic regions but only briefly colocalize during the periods of chromosome condensation and decondensation. This pattern of NCT binding at the core centromere, the region of maximal condensin enrichment, tracks the abundance of acetylated histone H4, as regulated by the Hat1-Mis16 acetyltransferase complex and recognized by the first Nrc1 bromodomain. Strikingly, mutants in NCT or Hat1-Mis16 restore the formation of segregation-competent chromosomes in cells containing defective condensin. These results are consistent with a model where NCT targets CKII to chromatin in a cell-cycle-directed manner in order to modulate the activity of condensin during chromosome condensation and decondensation.

Renal cell carcinoma comprises 2 to 3% of malignancies in adults with the most prevalent subtype being clear-cell RCC (ccRCC). This type of cancer is well characterized at the genomic and transcriptomic level and is associated with a loss of VHL that results in stabilization of HIF1. The current study focused on evaluating ccRCC stage dependent changes at the proteome level to provide insight into the molecular pathogenesis of ccRCC progression. To accomplish this, label-free proteomics was used to characterize matched tumor and normal-adjacent tissues from 84 patients with stage I to IV ccRCC. Using pooled samples 1551 proteins were identified, of which 290 were differentially abundant, while 783 proteins were identified using individual samples, with 344 being differentially abundant. These 344 differentially abundant proteins were enriched in metabolic pathways and further examination revealed metabolic dysfunction consistent with the Warburg effect. Additionally, the protein data indicated activation of ESRRA and ESRRG, and HIF1A, as well as inhibition of FOXA1, MAPK1 and WISP2. A subset analysis of complementary gene expression array data on 47 pairs of these same tissues indicated similar upstream changes, such as increased HIF1A activation with stage, though ESRRA and ESRRG activation and FOXA1 inhibition were not predicted from the transcriptomic data. The activation of ESRRA and ESRRG implied that HIF2A may also be activated during later stages of ccRCC, which was confirmed in the transcriptional analysis. This combined analysis highlights the importance of HIF1A and HIF2A in developing the ccRCC molecular phenotype as well as the potential involvement of ESRRA and ESRRG in driving these changes. In addition, cofilin-1, profilin-1, nicotinamide N-methyltransferase, and fructose-bisphosphate aldolase A were identified as candidate markers of late stage ccRCC. Utilization of data collected from heterogeneous biological domains strengthened the findings from each domain, demonstrating the complementary nature of such an analysis. Together these results highlight the importance of the VHL/HIF1A/HIF2A axis and provide a foundation and therapeutic targets for future studies. (Data are available via ProteomeXchange with identifier PXD003271 and MassIVE with identifier MSV000079511.).

Proliferating eukaryotic cells undergo a finite number of cell divisions before irreversibly exiting mitosis. Yet pathways that normally limit the number of cell divisions remain poorly characterized. Here we describe a screen of a collection of 3762 single gene mutants in the yeast Saccharomyces cerevisiae, accounting for 2/3 of annotated yeast ORFs, to search for mutants that undergo an atypically high number of cell divisions. Many of the potential longevity genes map to cellular processes not previously implicated in mitotic senescence, suggesting that regulatory mechanisms governing mitotic exit may be broader than currently anticipated. We focused on an ER-Golgi gene cluster isolated in this screen to determine how these ubiquitous organelles integrate into mitotic longevity. We report that a chronic increase in ER protein load signals an expansion in the assembly of autophagosomes in an Ire1-independent manner, accelerates trafficking of high molecular weight protein aggregates from the cytoplasm to the vacuoles, and leads to a profound enhancement of daughter cell production. We demonstrate that this catabolic network is evolutionarily conserved, as it also extends reproductive lifespan in the nematode Caenorhabditis elegans. Our data provide evidence that catabolism of protein aggregates, a natural byproduct of high protein synthesis and turn over in dividing cells, is among the drivers of mitotic longevity in eukaryotes.

Advances in mass spectrometry-based proteomics have allowed researchers to quantify the abundances of the different forms of three closely related proteins in the neurons of mice.

H2A.Z is an essential histone variant that has been implicated to have multiple chromosomal functions. To understand how H2A.Z participates in such diverse activities, we sought to identify downstream effector proteins that are recruited to chromatin via H2A.Z. For this purpose, we developed a nucleosome purification method to isolate H2A.Z-containing nucleosomes from human cells and used mass spectrometry to identify the co-purified nuclear proteins. Through stringent filtering, we identified the top 21 candidates, many of which have conserved structural motifs that bind post-translationally modified histones. We further validated the biological significance of one such candidate, Brd2, which is a double-bromodomain-containing protein known to function in transcriptional activation. We found that Brd2's preference for H2A.Z nucleosomes is mediated through a combination of hyperacetylated H4 on these nucleosomes, as well as additional features on H2A.Z itself. In addition, comparison of nucleosomes containing either H2A.Z-1 or H2A.Z-2 isoforms showed that significantly more Brd2 co-purifies with the former, suggesting these two isoforms engage different downstream effector proteins. Consistent with these biochemical analyses, we found that Brd2 is recruited to AR-regulated genes in an H2A.Z-dependent manner and that chemical inhibition of Brd2 recruitment greatly inhibits AR-regulated gene expression. Taken together, we propose that Brd2 is a key downstream mediator that links H2A.Z and transcriptional activation of AR-regulated genes. Moreover, this study validates the approach of using proteomics to identify nucleosome-interacting proteins in order to elucidate downstream mechanistic functions associated with the histone variant H2A.Z.

Phospholamban (PLN) is an effective inhibitor of the sarco(endo)plasmic reticulum Ca(2+)-ATPase, which transports Ca(2+) into the SR lumen, leading to muscle relaxation. A mutation of PLN in which one of the di-arginine residues at positions 13 and 14 was deleted led to a severe, early onset dilated cardiomyopathy. Here we were interested in determining the cellular mechanisms involved in this disease-causing mutation.Mutations deleting codons for either or both Arg13 or Arg14 resulted in the mislocalization of PLN from the ER. Our data show that PLN is recycled via the retrograde Golgi to ER membrane traffic pathway involving COP-I vesicles, since co-immunoprecipitation assays determined that COP I interactions are dependent on an intact di-arginine motif as PLN RDelta14 did not co-precipitate with COP I containing vesicles. Bioinformatic analysis determined that the di-arginine motif is present in the first 25 residues in a large number of all ER/SR Gene Ontology (GO) annotated proteins. Mutations in the di-arginine motif of the Sigma 1-type opioid receptor, the beta-subunit of the signal recognition particle receptor, and Sterol-O-acyltransferase, three proteins identified in our bioinformatic screen also caused mislocalization of these known ER-resident proteins.We conclude that PLN is enriched in the ER due to COP I-mediated transport that is dependent on its intact di-arginine motif and that the N-terminal di-arginine motif may act as a general ER retrieval sequence.

Nuclear receptors (NRs) belong to a superfamily of transcription factors that regulate numerous homeostatic, metabolic and reproductive processes. Taken together with their modulation by small lipophilic molecules, they also represent an important and successful class of drug targets. Although many NRs have been targeted successfully, the majority have not, and one third are still orphans. Here we report the development of an in vivo GFP-based reporter system suitable for monitoring NR activities in all cells and tissues using live zebrafish (Danio rerio). The human NR fusion proteins used also contain a new affinity tag cassette allowing the purification of receptors with bound molecules from responsive tissues. We show that these constructs 1) respond as expected to endogenous zebrafish hormones and cofactors, 2) facilitate efficient receptor and cofactor purification, 3) respond robustly to NR hormones and drugs and 4) yield readily quantifiable signals. Transgenic lines representing the majority of human NRs have been established and are available for the investigation of tissue- and isoform-specific ligands and cofactors.

BACKGROUND: MicroRNAs (miRNAs) are known to be important regulators of both organ development and tumorigenesis. MiRNA networks and their regulation of messenger RNA (mRNA) translation and protein expression in specific biological processes are poorly understood. METHODS: We explored the dynamic regulation of miRNAs in mouse lung organogenesis. Comprehensive miRNA and mRNA profiling was performed encompassing all recognized stages of lung development beginning at embryonic day 12 and continuing to adulthood. We analyzed the expression patterns of dynamically regulated miRNAs and mRNAs using a number of statistical and computational approaches, and in an integrated manner with protein levels from an existing mass-spectrometry derived protein database for lung development. RESULTS: In total, 117 statistically significant miRNAs were dynamically regulated during mouse lung organogenesis and clustered into distinct temporal expression patterns. 11,220 mRNA probes were also shown to be dynamically regulated and clustered into distinct temporal expression patterns, with 3 major patterns accounting for 75% of all probes. 3,067 direct miRNA-mRNA correlation pairs were identified involving 37 miRNAs. Two defined correlation patterns were observed upon integration with protein data: 1) increased levels of specific miRNAs directly correlating with downregulation of predicted mRNA targets; and 2) increased levels of specific miRNAs directly correlating with downregulation of translated target proteins without detectable changes in mRNA levels. Of 1345 proteins analyzed, 55% appeared to be regulated in this manner with a direct correlation between miRNA and protein level, but without detectable change in mRNA levels. CONCLUSION: Systematic analysis of microRNA, mRNA, and protein levels over the time course of lung organogenesis demonstrates dynamic regulation and reveals 2 distinct patterns of miRNA-mRNA interaction. The translation of target proteins affected by miRNAs independent of changes in mRNA level appears to be a prominent mechanism of developmental regulation in lung organogenesis.