A Multiplexed Quantitative Proteomics Approach to the Human Plasma Protein Signature.

Despite the plasma proteome being able to provide a unique insight into the health and disease status of individuals, holding singular promise as a source of protein biomarkers that could be pivotal in the context of personalized medicine, only around 100 proteins covering a few human conditions have been approved as biomarkers by the US Food and Drug Administration (FDA) so far. Mass spectrometry (MS) currently has enormous potential for high-throughput analysis in clinical research; however, plasma proteomics remains challenging mainly due to the wide dynamic range of plasma protein abundances and the time-consuming procedures required. We applied a new MS-based multiplexed proteomics workflow to quantitate proteins, encompassing 67 FDA-approved biomarkers, in >1300 human plasma samples from a clinical cohort. Our results indicate that this workflow is suitable for large-scale clinical studies, showing good accuracy and reproducibility (coefficient of variation (CV) < 20 for 90% of the proteins). Furthermore, we identified plasma signature proteins (stable in time on an individual basis), stable proteins (exhibiting low biological variability and high temporal stability), and highly variable proteins (with low temporal stability) that can be used for personalized health monitoring and medicine.

In-gel protein digestion using acidic methanol produces a highly selective methylation of glutamic acid residues.

Mass-tolerant open search methods allow the high-throughput analysis of modified peptides by mass spectrometry. These techniques have paved the way to unbiased analysis of post-translational modifications in biological contexts, as well as of chemical modifications produced during the manipulation of protein samples. In this work, we have analyzed in-depth a wide variety of samples of different biological origin, including cells, extracellular vesicles, secretomes, centrosomes and tissue preparations, using Comet-ReCom, a recently improved version of the open search engine Comet-PTM. Our results demonstrate that glutamic acid residues undergo intensive methyl esterification when protein digestion is performed using in-gel techniques, but not using gel-free approaches. This effect was highly specific to Glu and was not found for other methylable residues such as Asp.

Metabolic changes contribute to maladaptive right ventricular hypertrophy in pulmonary hypertension beyond pressure overload: an integrative imaging and omics investigation.

Right ventricular (RV) failure remains the strongest determinant of survival in pulmonary hypertension (PH). We aimed to identify relevant mechanisms, beyond pressure overload, associated with maladaptive RV hypertrophy in PH. To separate the effect of pressure overload from other potential mechanisms, we developed in pigs two experimental models of PH (M1, by pulmonary vein banding and M2, by aorto-pulmonary shunting) and compared them with a model of pure pressure overload (M3, pulmonary artery banding) and a sham-operated group. Animals were assessed at 1 and 8 months by right heart catheterization, cardiac magnetic resonance and blood sampling, and myocardial tissue was analyzed. Plasma unbiased proteomic and metabolomic data were compared among groups and integrated by an interaction network analysis. A total of 33 pigs completed follow-up (M1, n = 8; M2, n = 6; M3, n = 10; and M0, n = 9). M1 and M2 animals developed PH and reduced RV systolic function, whereas animals in M3 showed increased RV systolic pressure but maintained normal function. Significant plasma arginine and histidine deficiency and complement system activation were observed in both PH models (M1&M2), with additional alterations to taurine and purine pathways in M2. Changes in lipid metabolism were very remarkable, particularly the elevation of free fatty acids in M2. In the integrative analysis, arginine-histidine-purines deficiency, complement activation, and fatty acid accumulation were significantly associated with maladaptive RV hypertrophy. Our study integrating imaging and omics in large-animal experimental models demonstrates that, beyond pressure overload, metabolic alterations play a relevant role in RV dysfunction in PH.

iSanXoT: A standalone application for the integrative analysis of mass spectrometry-based quantitative proteomics data.

Many bioinformatics tools are available for the quantitative analysis of proteomics experiments. Most of these tools use a dedicated statistical model to derive absolute quantitative protein values from mass spectrometry (MS) data. Here, we present iSanXoT, a standalone application that processes relative abundances between MS signals and then integrates them sequentially to upper levels using the previously published Generic Integration Algorithm (GIA). iSanXoT offers unique capabilities that complement conventional quantitative software applications, including statistical weighting and independent modeling of error distributions in each integration, aggregation of technical or biological replicates, quantification of posttranslational modifications, and analysis of coordinated protein behavior. iSanXoT is a standalone, user-friendly application that accepts output from popular proteomics pipelines and enables unrestricted creation of quantification workflows and fully customizable reports that can be reused across projects or shared among users. Numerous publications attest the successful application of diverse integrative workflows constructed using the GIA for the analysis of high-throughput quantitative proteomics experiments. iSanXoT has been tested with the main operating systems. Download links for the corresponding distributions are available at https://github.com/CNIC-Proteomics/iSanXoT/releases.

Oxidative Post-translational Protein Modifications upon Ischemia/Reperfusion Injury.

While reperfusion, or restoration of coronary blood flow in acute myocardial infarction, is a requisite for myocardial salvage, it can paradoxically induce a specific damage known as ischemia/reperfusion (I/R) injury. Our understanding of the precise pathophysiological molecular alterations leading to I/R remains limited. In this study, we conducted a comprehensive and unbiased time-course analysis of post-translational modifications (PTMs) in the post-reperfused myocardium of two different animal models (pig and mouse) and evaluated the effect of two different cardioprotective therapies (ischemic preconditioning and neutrophil depletion). In pigs, a first wave of irreversible oxidative damage was observed at the earliest reperfusion time (20 min), impacting proteins essential for cardiac contraction. A second wave, characterized by irreversible oxidation on different residues and reversible Cys oxidation, occurred at late stages (6-12 h), affecting mitochondrial, sarcomere, and inflammation-related proteins. Ischemic preconditioning mitigated the I/R damage caused by the late oxidative wave. In the mouse model, the two-phase pattern of oxidative damage was replicated, and neutrophil depletion mitigated the late wave of I/R-related damage by preventing both Cys reversible oxidation and irreversible oxidation. Altogether, these data identify protein PTMs occurring late after reperfusion as an actionable therapeutic target to reduce the impact of I/R injury.

MYH10 activation rescues contractile defects in arrhythmogenic cardiomyopathy (ACM).

The most prevalent genetic form of inherited arrhythmogenic cardiomyopathy (ACM) is caused by mutations in desmosomal plakophilin-2 (PKP2). By studying pathogenic deletion mutations in the desmosomal protein PKP2, here we identify a general mechanism by which PKP2 delocalization restricts actomyosin network organization and cardiac sarcomeric contraction in this untreatable disease. Computational modeling of PKP2 variants reveals that the carboxy-terminal (CT) domain is required for N-terminal domain stabilization, which determines PKP2 cortical localization and function. In mutant PKP2 cells the expression of the interacting protein MYH10 rescues actomyosin disorganization. Conversely, dominant-negative MYH10 mutant expression mimics the pathogenic CT-deletion PKP2 mutant causing actin network abnormalities and right ventricle systolic dysfunction. A chemical activator of non-muscle myosins, 4-hydroxyacetophenone (4-HAP), also restores normal contractility. Our findings demonstrate that activation of MYH10 corrects the deleterious effect of PKP2 mutant over systolic cardiac contraction, with potential implications for ACM therapy.

Incongruence between transcriptional and vascular pathophysiological cell states.

The Notch pathway is a major regulator of endothelial transcriptional specification. Targeting the Notch receptors or Delta-like ligand 4 (Dll4) dysregulates angiogenesis. Here, by analyzing single and compound genetic mutants for all Notch signaling members, we find significant differences in the way ligands and receptors regulate liver vascular homeostasis. Loss of Notch receptors caused endothelial hypermitogenic cell-cycle arrest and senescence. Conversely, Dll4 loss triggered a strong Myc-driven transcriptional switch inducing endothelial proliferation and the tip-cell state. Myc loss suppressed the induction of angiogenesis in the absence of Dll4, without preventing the vascular enlargement and organ pathology. Similarly, inhibition of other pro-angiogenic pathways, including MAPK/ERK and mTOR, had no effect on the vascular expansion induced by Dll4 loss; however, anti-VEGFA treatment prevented it without fully suppressing the transcriptional and metabolic programs. This study shows incongruence between single-cell transcriptional states, vascular phenotypes and related pathophysiology. Our findings also suggest that the vascular structure abnormalization, rather than neoplasms, causes the reported anti-Dll4 antibody toxicity.

ISG20L2: an RNA nuclease regulating T cell activation.

ISG20L2, a 3' to 5' exoribonuclease previously associated with ribosome biogenesis, is identified here in activated T cells as an enzyme with a preferential affinity for uridylated miRNA substrates. This enzyme is upregulated in T lymphocytes upon TCR and IFN type I stimulation and appears to be involved in regulating T cell function. ISG20L2 silencing leads to an increased basal expression of CD69 and induces greater IL2 secretion. However, ISG20L2 absence impairs CD25 upregulation, CD3 synaptic accumulation and MTOC translocation towards the antigen-presenting cell during immune synapsis. Remarkably, ISG20L2 controls the expression of immunoregulatory molecules, such as AHR, NKG2D, CTLA-4, CD137, TIM-3, PD-L1 or PD-1, which show increased levels in ISG20L2 knockout T cells. The dysregulation observed in these key molecules for T cell responses support a role for this exonuclease as a novel RNA-based regulator of T cell function.

γ-Linolenic acid in maternal milk drives cardiac metabolic maturation.

Birth presents a metabolic challenge to cardiomyocytes as they reshape fuel preference from glucose to fatty acids for postnatal energy production1,2. This adaptation is triggered in part by post-partum environmental changes3, but the molecules orchestrating cardiomyocyte maturation remain unknown. Here we show that this transition is coordinated by maternally supplied γ-linolenic acid (GLA), an 18:3 omega-6 fatty acid enriched in the maternal milk. GLA binds and activates retinoid X receptors4 (RXRs), ligand-regulated transcription factors that are expressed in cardiomyocytes from embryonic stages. Multifaceted genome-wide analysis revealed that the lack of RXR in embryonic cardiomyocytes caused an aberrant chromatin landscape that prevented the induction of an RXR-dependent gene expression signature controlling mitochondrial fatty acid homeostasis. The ensuing defective metabolic transition featured blunted mitochondrial lipid-derived energy production and enhanced glucose consumption, leading to perinatal cardiac dysfunction and death. Finally, GLA supplementation induced RXR-dependent expression of the mitochondrial fatty acid homeostasis signature in cardiomyocytes, both in vitro and in vivo. Thus, our study identifies the GLA-RXR axis as a key transcriptional regulatory mechanism underlying the maternal control of perinatal cardiac metabolism.

Quantification of Farnesylated Progerin in Hutchinson-Gilford Progeria Patient Cells by Mass Spectrometry.

Hutchinson-Gilford progeria syndrome (HGPS) is a rare fatal disorder characterized by premature aging and death at a median age of 14.5 years. The most common cause of HGPS (affecting circa 90% of patients) is a de novo heterozygous synonymous single-base substitution (c.1824C>T; p.G608G) in the LMNA gene that results in the accumulation of progerin, an aberrant form of lamin A that, unlike mature lamin A, remains permanently farnesylated. The ratio of progerin to mature lamin A correlates with disease severity in HGPS patients, and can be used to assess the effectiveness of therapies aimed at lessening aberrant splicing or progerin farnesylation. We recently showed that the endogenous content of lamin A and progerin can be measured by mass spectrometry (MS), providing an alternative to immunological methods, which lack the necessary specificity and quantitative accuracy. Here, we present the first non-immunological method that reliably quantifies the levels of wild-type lamin A and farnesylated progerin in cells from HGPS patients. This method, which is based on a targeted MS approach and the use of isotope-labeled internal standards, could be applied in ongoing clinical trials evaluating the efficacy of drugs that inhibit progerin farnesylation.

Extracellular vesicles from Listeria monocytogenes-infected dendritic cells alert the innate immune response.

Communication through cell-cell contacts and extracellular vesicles (EVs) enables immune cells to coordinate their responses against diverse types of pathogens. The function exerted by EVs in this context depends on the proteins and nucleic acids loaded into EVs, which elicit specific responses involved in the resolution of infection. Several mechanisms control protein and nucleic acid loading into EVs; in this regard, acetylation has been described as a mechanism of cellular retention during protein sorting to exosomes. HDAC6 is a deacetylase involved in the control of cytoskeleton trafficking, organelle polarity and cell migration, defense against Listeria monocytogenes (Lm) infection and other immune related functions. Here, we show that the protein content of dendritic cells (DCs) and their secreted EVs (DEVs) vary during Lm infection, is enriched in proteins related to antiviral functions compared to non-infected cells and depends on HDAC6 expression. Analyses of the post-translational modifications revealed an alteration of the acetylation and ubiquitination profiles upon Lm infection both in DC lysates and DEVs. Functionally, EVs derived from infected DCs upregulate anti-pathogenic genes (e.g. inflammatory cytokines) in recipient immature DCs, which translated into protection from subsequent infection with vaccinia virus. Interestingly, absence of Listeriolysin O in Lm prevents DEVs from inducing this anti-viral state. In summary, these data underscore a new mechanism of communication between bacteria-infected DC during infection as they alert neighboring, uninfected DCs to promote antiviral responses.

TGF-ß-Upregulated Lnc-Nr6a1 Acts as a Reservoir of miR-181 and Mediates Assembly of a Glycolytic Complex.

Long non-coding RNAs (lncRNAs) have emerged as key regulators in a wide range of biological processes. Here, we identified a mouse miRNA-host gene lncRNA (lnc-Nr6a1) upregulated early during epithelial-to-mesenchymal transition (EMT). We show that when lncRNA is processed, it gives rise to two abundant polyadenylated isoforms, lnc-Nr6a1-1 and lnc-Nr6a1-2, and a longer non-polyadenylated microprocessor-driven lnc-pri-miRNA containing clustered pre-miR-181a2 and pre-miR-181b2 hairpins. Ectopic expression of the lnc-Nr6a1-1 or lnc-Nr6a1-2 isoform enhanced cell migration and the invasive capacity of the cells, whereas the expression of the isoforms and miR-181a2 and miR-181b2 conferred anoikis resistance. Lnc-Nr6a1 gene deletion resulted in cells with lower adhesion capacity and reduced glycolytic metabolism, which are restored by lnc-Nr6a1-1 isoform expression. We performed identification of direct RNA interacting proteins (iDRIP) to identify proteins interacting directly with the lnc-Nr6a1-1 isoform. We defined a network of interacting proteins, including glycolytic enzymes, desmosome proteins and chaperone proteins; and we demonstrated that the lnc-Nr6a1-1 isoform directly binds and acts as a scaffold molecule for the assembly of ENO1, ALDOA, GAPDH, and PKM glycolytic enzymes, along with LDHA, supporting substrate channeling for efficient glycolysis. Our results unveil a role of Lnc-Nr6a1 as a multifunctional lncRNA acting as a backbone for multiprotein complex formation and primary microRNAs.

Cardiovascular Progerin Suppression and Lamin A Restoration Rescue Hutchinson-Gilford Progeria Syndrome.

BACKGROUND: Hutchinson-Gilford progeria syndrome (HGPS) is a rare disorder characterized by premature aging and death mainly because of myocardial infarction, stroke, or heart failure. The disease is provoked by progerin, a variant of lamin A expressed in most differentiated cells. Patients look healthy at birth, and symptoms typically emerge in the first or second year of life. Assessing the reversibility of progerin-induced damage and the relative contribution of specific cell types is critical to determining the potential benefits of late treatment and to developing new therapies. METHODS: We used CRISPR-Cas9 technology to generate LmnaHGPSrev/HGPSrev (HGPSrev) mice engineered to ubiquitously express progerin while lacking lamin A and allowing progerin suppression and lamin A restoration in a time- and cell type-specific manner on Cre recombinase activation. We characterized the phenotype of HGPSrev mice and crossed them with Cre transgenic lines to assess the effects of suppressing progerin and restoring lamin A ubiquitously at different disease stages as well as specifically in vascular smooth muscle cells and cardiomyocytes. RESULTS: Like patients with HGPS, HGPSrev mice appear healthy at birth and progressively develop HGPS symptoms, including failure to thrive, lipodystrophy, vascular smooth muscle cell loss, vascular fibrosis, electrocardiographic anomalies, and precocious death (median lifespan of 15 months versus 26 months in wild-type controls, P<0.0001). Ubiquitous progerin suppression and lamin A restoration significantly extended lifespan when induced in 6-month-old mildly symptomatic mice and even in severely ill animals aged 13 months, although the benefit was much more pronounced on early intervention (84.5% lifespan extension in mildly symptomatic mice, P<0.0001, and 6.7% in severely ill mice, P<0.01). It is remarkable that major vascular alterations were prevented and lifespan normalized in HGPSrev mice when progerin suppression and lamin A restoration were restricted to vascular smooth muscle cells and cardiomyocytes. CONCLUSIONS: HGPSrev mice constitute a new experimental model for advancing knowledge of HGPS. Our findings suggest that it is never too late to treat HGPS, although benefit is much more pronounced when progerin is targeted in mice with mild symptoms. Despite the broad expression pattern of progerin and its deleterious effects in many organs, restricting its suppression to vascular smooth muscle cells and cardiomyocytes is sufficient to prevent vascular disease and normalize lifespan.

Local Pressure Drives Low-Density Lipoprotein Accumulation and Coronary Atherosclerosis in Hypertensive Minipigs.

BACKGROUND: The mechanisms by which hypertension accelerates coronary artery disease are poorly understood. Patients with hypertension often have confounding humoral changes, and to date, no experimental models have allowed analysis of the isolated effect of pressure on atherosclerosis in a setting that recapitulates the dimensions and biomechanics of human coronary arteries. OBJECTIVES: This study sought to analyze the effect of pressure on coronary atherosclerosis and explore the underlying mechanisms. METHODS: Using inflatable suprarenal aortic cuffs, we increased mean arterial pressure by >30 mm Hg in the cephalad body part of wild-type and hypercholesterolemic proprotein convertase subtilisin kexin type 9 (PCSK9)D374Y Yucatan minipigs for >1 year. Caudal pressures remained normal. RESULTS: Under hypercholesterolemic conditions in PCSK9D374Y transgenic minipigs, cephalad hypertension accelerated coronary atherosclerosis to almost 5-fold with consistent development of fibroatheromas that were sufficiently large to cause stenosis on computed tomography angiography. This was caused by local pressure forces, because vascular beds shielded from hypertension, but exposed to the same humoral factors, showed no changes in lesion formation. The same experiment was conducted under normocholesterolemic conditions in wild-type minipigs to examine the underlying mechanisms. Hypertension produced clear changes in the arterial proteome with increased abundance of mechanical strength proteins and reduced levels of infiltrating plasma macromolecules. This was paralleled by increased smooth muscle cells and increased intimal accumulation of low-density lipoproteins in the coronary arteries. CONCLUSIONS: Increased pressure per se facilitates coronary atherosclerosis. Our data indicate that restructuring of the artery to match increased tensile forces in hypertension alters the passage of macromolecules and leads to increased intimal accumulation of low-density lipoproteins.

Proteomic analysis of plasma proteins of high-flux haemodialysis and on-line haemodiafiltration patients reveals differences in transthyretin levels related with anaemia.

A large proportion of end-stage renal disease (ESRD) patients under long-term haemodialysis, have persistent anaemia and require high doses of recombinant human erythropoietin (rhEPO). However, the underlying mechanisms of renal anaemia have not been fully elucidated in these patients. In this study, we will be focusing on anaemia and plasma proteins in ESRD patients on high-flux haemodialysis (HF) and on-line haemodiafiltration (HDF), to investigate using two proteomic approaches if patients undergoing these treatments develop differences in their plasma protein composition and how this could be related to their anaemia. The demographic and biochemical data revealed that HDF patients had lower anaemia and much lower rhEPO requirements than HF patients. Regarding their plasma proteomes, HDF patients had increased levels of a protein highly similar to serotransferrin, trypsin-1 and immunoglobulin heavy constant chain alpha-1, and lower levels of alpha-1 antitrypsin, transthyretin, apolipoproteins E and C-III, and haptoglobin-related protein. Lower transthyretin levels in HDF patients were further confirmed by transthyretin-peptide quantification and western blot detection. Since ESRD patients have increased transthyretin, a protein that can aggregate and inhibit transferrin endocytosis and erythropoiesis, our finding that HDF patients have lower transthyretin and lower anaemia suggests that the decrease in transthyretin plasma levels would allow an increase in transferrin endocytosis, contributing to erythropoiesis. Thus, transthyretin could be a critical actor for anaemia in ESRD patients and a novel player for haemodialysis adequacy.

Identification of common cardiometabolic alterations and deregulated pathways in mouse and pig models of aging.

Aging is the main risk factor for cardiovascular and metabolic diseases, which have become a global concern as the world population ages. These diseases and the aging process are exacerbated in Hutchinson-Gilford progeria syndrome (HGPS or progeria). Here, we evaluated the cardiometabolic disease in animal models of premature and normal aging with the aim of identifying alterations that are shared or specific to each condition. Despite differences in body composition and metabolic markers, prematurely and normally aging mice developed heart failure and similar cardiac electrical abnormalities. High-throughput proteomics of the hearts of progeric and normally aged mice revealed altered protein oxidation and glycation, as well as dysregulated pathways regulating energy metabolism, proteostasis, gene expression, and cardiac muscle contraction. These results were corroborated in the hearts of progeric pigs, underscoring the translational potential of our findings, which could help in the design of strategies to prevent or slow age-related cardiometabolic disease.

Complement C5 Protein as a Marker of Subclinical Atherosclerosis.

BACKGROUND: The mechanisms underlying early atherosclerotic plaque formation are not completely understood. Moreover, plasma biomarkers of subclinical atherosclerosis are lacking. OBJECTIVES: The purpose of this study was to analyze the temporal and topologically resolved protein changes taking place in human aortas with early atherosclerosis to find new potential diagnostic and/or therapeutic targets. METHODS: The protein composition of healthy aortas (media layer) or with early atheroma (fatty streak and fibrolipidic, media and intima layers) was analyzed by deep quantitative multiplexed proteomics. Further analysis was performed by Western blot, immunohistochemistry, real-time polymerase chain reaction, and enzyme-linked immunosorbent assay. Plasma levels of complement C5 were analyzed in relation to the presence of generalized (>2 plaques) or incipient (0 to 2 plaques) subclinical atherosclerosis in 2 independent clinical cohorts (PESA [Progression of Early Subclinical Atherosclerosis] [n = 360] and NEFRONA [National Observatory of Atherosclerosis in Nephrology] [n = 394]). RESULTS: Proteins involved in lipid transport, complement system, immunoglobulin superfamily, and hemostasis are increased in early plaques. Components from the complement activation pathway were predominantly increased in the intima of fibrolipidic plaques. Among them, increased C5 protein levels were further confirmed by Western blot, enzyme-linked immunosorbent assay and immunohistochemistry, and associated with in situ complement activation. Plasma C5 was significantly increased in individuals with generalized subclinical atherosclerosis in both PESA and NEFRONA cohorts, independently of risk factors. Moreover, in the PESA study, C5 plasma levels positively correlated with global plaque volume and coronary calcification. CONCLUSIONS: Activation of the complement system is a major alteration in early atherosclerotic plaques and is reflected by increased C5 plasma levels, which have promising value as a novel circulating biomarker of subclinical atherosclerosis.

Improved integrative analysis of the thiol redox proteome using filter-aided sample preparation.

Changes in the oxidation state of protein Cys residues are involved in cell signalling and play a key role in a variety of pathophysiological states. We had previously developed GELSILOX, an in-gel method that enables the large-scale, parallel analysis of dynamic alterations to the redox state of Cys sites and protein abundance changes. Here we present FASILOX, a further development of the GELSILOX approach featuring: i) significantly increased peptide recovery, ii) enhanced sensitivity for the detection of Cys oxidative alterations, and iii) streamlined workflow that results in shortened assay duration. In mitochondria isolated from the adipose tissue of obese, diabetic patients, FASILOX revealed a sexually dimorphic trait of Cys oxidation involving mainly mitochondrial oxidative phosphorylation complexes. These results provide the first evidence for a decreased efficiency in the antioxidant response of men as compared to women.

A Proteomics Signature of Mild Hypospadias: A Pilot Study.

Background and Objective: Mild hypospadias is a birth congenital condition characterized by the relocation of the male urethral meatus from its typical anatomical position near the tip of the glans penis, to a lower ventral position up to the brim of the glans corona, which can also be accompanied by foreskin ventral deficiency. For the most part, a limited number of cases have known etiology. We have followed a high-throughput proteomics approach to study the proteome in mild hypospadias patients. Methods: Foreskin samples from patients with mild hypospadias were collected during urethroplasty, while control samples were collected during elective circumcision (n = 5/group). A high-throughput, quantitative proteomics approach based on multiplexed peptide stable isotope labeling (SIL) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis was used to ascertain protein abundance changes in hypospadias patients when compared to control samples. Results: A total of 4,815 proteins were quantitated (2,522 with at least two unique peptides). One hundred and thirty-three proteins from patients with mild hypospadias showed significant abundance changes with respect to control samples, where 38 proteins were increased, and 95 proteins were decreased. Unbiased functional biological analysis revealed that both mitochondrial energy production and apoptotic signaling pathways were enriched in mild hypospadias. Conclusions: This first comprehensive proteomics characterization of mild hypospadias shows molecular changes associated with essential cellular processes related to energy production and apoptosis. Further evaluation of the proteome may expand the search of novel candidates in the etiology of mild hypospadias and could also lead to the identification of biomarkers for this congenital urogenital condition.