Christian D. Kelstrup

Recalibrating Equus evolution using the genome sequence of an early Middle Pleistocene horse

The rich fossil record of equids has made them a model for evolutionary processes. Here we present a 1.12-times coverage draft genome from a horse bone recovered from permafrost dated to approximately 560–780 thousand years before present (kyr bp). Our data represent the oldest full genome sequence determined so far by almost an order of magnitude. For comparison, we sequenced the genome of a Late Pleistocene horse (43 kyr bp), and modern genomes of five domestic horse breeds (Equus ferus caballus), a Przewalski’s horse (E. f. przewalskii) and a donkey (E. asinus). Our analyses suggest that the Equus lineage giving rise to all contemporary horses, zebras and donkeys originated 4.0–4.5 million years before present (Myr bp), twice the conventionally accepted time to the most recent common ancestor of the genus Equus. We also find that horse population size fluctuated multiple times over the past 2 Myr, particularly during periods of severe climatic changes. We estimate that the Przewalski’s and domestic horse populations diverged 38–72 kyr bp, and find no evidence of recent admixture between the domestic horse breeds and the Przewalski’s horse investigated. This supports the contention that Przewalski’s horses represent the last surviving wild horse population. We find similar levels of genetic variation among Przewalski’s and domestic populations, indicating that the former are genetically viable and worthy of conservation efforts. We also find evidence for continuous selection on the immune system and olfaction throughout horse evolution. Finally, we identify 29 genomic regions among horse breeds that deviate from neutrality and show low levels of genetic variation compared to the Przewalski’s horse. Such regions could correspond to loci selected early during domestication.

Proteomic Analysis of Lysine Acetylation Sites in Rat Tissues Reveals Organ Specificity and Subcellular Patterns

SUMMARY Lysine acetylation is a major posttranslational modification involved in a broad array of physiological functions. Here, we provide an organ-wide map of lysine acetylation sites from 16 rat tissues analyzed by high-resolution tandem mass spectrometry. We quantify 15,474 modification sites on 4,541 proteins and provide the data set as a web-based database. We demonstrate that lysine acetylation displays site-specific sequence motifs that diverge between cellular compartments, with a significant fraction of nuclear sites conforming to the consensus motifs G-AcK and AcK-P. Our data set reveals that the subcellular acetylation distribution is tissue-type dependent and that acetylation targets tissue-specific pathways involved in fundamental physiological processes. We compare lysine acetylation patterns for rat as well as human skeletal muscle biopsies and demonstrate its general involvement in muscle contraction. Furthermore, we illustrate that acetylation of fructose-bisphosphate aldolase and glycerol-3-phosphate dehydrogenase serves as a cellular mechanism to switch off enzymatic activity.

Proteomic Analyses Reveal Divergent Ubiquitylation Site Patterns in Murine Tissues

Posttranslational modifications of proteins increase the complexity of the cellular proteome and enable rapid regulation of protein functions in response to environmental changes. Protein ubiquitylation is a central regulatory posttranslational modification that controls numerous biological processes including proteasomal degradation of proteins, DNA damage repair and innate immune responses. Here we combine high-resolution mass spectrometry with single-step immunoenrichment of di-glycine modified peptides for mapping of endogenous putative ubiquitylation sites in murine tissues. We identify more than 20,000 unique ubiquitylation sites on proteins involved in diverse biological processes. Our data reveals that ubiquitylation regulates core signaling pathways common for each of the studied tissues. In addition, we discover that ubiquitylation regulates tissue-specific signaling networks. Many tissue-specific ubiquitylation sites were obtained from brain highlighting the complexity and unique physiology of this organ. We further demonstrate that different di-glycine-lysine-specific monoclonal antibodies exhibit sequence preferences, and that their complementary use increases the depth of ubiquitylation site analysis, thereby providing a more unbiased view of protein ubiquitylation.

Quantitative Phosphoproteomics Dissection of Seven-transmembrane Receptor Signaling Using Full and Biased Agonists

Seven-transmembrane receptors (7TMRs) signal through the well described heterotrimeric G proteins but can also activate G protein-independent signaling pathways of which the impact and complexity are less understood. The angiotensin II type 1 receptor (AT1R) is a prototypical 7TMR and an important drug target in cardiovascular diseases. “Biased agonists” with intrinsic “functional selectivity” that simultaneously blocks Gαq protein activity and activates G protein-independent pathways of the AT1R confer important perspectives in treatment of cardiovascular diseases. In this study, we performed a global quantitative phosphoproteomics analysis of the AT1R signaling network. We analyzed ligand-stimulated SILAC (stable isotope labeling by amino acids in cell culture) cells by high resolution (LTQ-Orbitrap) MS and compared the phosphoproteomes of the AT1R agonist angiotensin II and the biased agonist [Sar1,Ile4,Ile8]angiotensin II (SII angiotensin II), which only activates the Gαq protein-independent signaling. We quantified more than 10,000 phosphorylation sites of which 1183 were regulated by angiotensin II or its analogue SII angiotensin II. 36% of the AT1R-regulated phosphorylations were regulated by SII angiotensin II. Analysis of phosphorylation site patterns showed a striking distinction between protein kinases activated by Gαq protein-dependent and -independent mechanisms, and we now place protein kinase D as a key protein involved in both Gαq-dependent and -independent AT1R signaling. This study provides substantial novel insight into angiotensin II signal transduction and is the first study dissecting the differences between a full agonist and a biased agonist from a 7TMR on a systems-wide scale. Importantly, it reveals a previously unappreciated diversity and quantity of Gαq protein-independent signaling and uncovers novel signaling pathways. We foresee that the amount and diversity of G protein-independent signaling may be more pronounced than previously recognized for other 7TMRs as well. Quantitative mass spectrometry is a promising tool for evaluation of the signaling properties of biased agonists to other receptors in the future.

Rapid and deep proteomes by faster sequencing on a benchtop quadrupole ultra-high-field Orbitrap mass spectrometer.

Shotgun proteomics is a powerful technology for global analysis of proteins and their post-translational modifications. Here, we investigate the faster sequencing speed of the latest Q Exactive HF mass spectrometer, which features an ultra-high-field Orbitrap mass analyzer. Proteome coverage is evaluated by four different acquisition methods and benchmarked across three generations of Q Exactive instruments (ProteomeXchange data set PXD001305). We find the ultra-high-field Orbitrap mass analyzer to be capable of attaining a sequencing speed above 20 Hz, and it routinely exceeds 10 peptide spectrum matches per second or up to 600 new peptides sequenced per gradient minute. We identify 4400 proteins from 1 μg of HeLa digest using a 1 h gradient, which is an approximately 30% improvement compared to that with previous instrumentation. In addition, we show that very deep proteome coverage can be achieved in less than 24 h of analysis time by offline high-pH reversed-phase peptide fractionation, from which we identify more than 140,000 unique peptide sequences. This is comparable to state-of-the-art multiday, multienzyme efforts. Finally, the acquisition methods are evaluated for single-shot phosphoproteomics, where we identify 7600 unique HeLa phosphopeptides in one gradient hour and find the quality of fragmentation spectra to be more important than quantity for accurate site assignment.

In Vivo Phosphoproteomics Analysis Reveals the Cardiac Targets of β-Adrenergic Receptor Signaling

Analysis of phosphorylated proteins from the hearts of mice given drugs targeting β-adrenergic receptors may aid in treating heart disease. Getting to the Heart of Signaling Patients with high blood pressure and other heart-related conditions routinely take inhibitors of β-adrenergic receptors (βARs) to prevent cardiac dysfunction. βAR signaling leads to the increased contractility of cardiomyocytes, among other effects; however, the number of downstream targets of βARs is unclear. Lundby et al. treated mice with combinations of specific β1AR and β2AR agonists and antagonists to activate each receptor isoform individually before harvesting the hearts and subjecting them to phosphoproteomics analysis. The authors identified previously uncharacterized peptides and sites phosphorylated in response to β1AR signaling, as well as characterized the activation of a potassium channel important for increasing heart rate. This in vivo approach provides insights into βAR signaling pathways that may help in understanding how heart diseases develop and how they may be treated. β-Blockers are widely used to prevent cardiac arrhythmias and to treat hypertension by inhibiting β-adrenergic receptors (βARs) and thus decreasing contractility and heart rate. βARs initiate phosphorylation-dependent signaling cascades, but only a small number of the target proteins are known. We used quantitative in vivo phosphoproteomics to identify 670 site-specific phosphorylation changes in murine hearts in response to acute treatment with specific βAR agonists. The residues adjacent to the regulated phosphorylation sites exhibited a sequence-specific preference (R-X-X-pS/T), and integrative analysis of sequence motifs and interaction networks suggested that the kinases AMPK (adenosine 5′-monophosphate–activated protein kinase), Akt, and mTOR (mammalian target of rapamycin) mediate βAR signaling, in addition to the well-established pathways mediated by PKA (cyclic adenosine monophosphate–dependent protein kinase) and CaMKII (calcium/calmodulin-dependent protein kinase type II). We found specific regulation of phosphorylation sites on six ion channels and transporters that mediate increased ion fluxes at higher heart rates, and we showed that phosphorylation of one of these, Ser92 of the potassium channel KV7.1, increased current amplitude. Our data set represents a quantitative analysis of phosphorylated proteins regulated in vivo upon stimulation of seven-transmembrane receptors, and our findings reveal previously unknown phosphorylation sites that regulate myocardial contractility, suggesting new potential targets for the treatment of heart disease and hypertension.

The Ndc80 internal loop is required for recruitment of the Ska complex to establish end-on microtubule attachment to kinetochores

Summary The Ndc80 complex establishes end-on attachment of kinetochores to microtubules, which is essential for chromosome segregation. The Ndc80 subunit is characterized by an N-terminal region that binds directly to microtubules, and a long coiled-coil region that interacts with Nuf2. A loop region in Ndc80 that generates a kink in the structure disrupts the long coiled-coil region but the exact function of this loop, has until now, not been clear. Here we show that this loop region is essential for end-on attachment of kinetochores to microtubules in human cells. Cells expressing loop mutants of Ndc80 are unable to align the chromosomes, and stable kinetochore fibers are absent. Through quantitative mass spectrometry and immunofluorescence we found that the binding of the spindle and kinetochore associated (Ska) complex depends on the loop region, explaining why end-on attachment is defective. This underscores the importance of the Ndc80 loop region in coordinating chromosome segregation through the recruitment of specific proteins to the kinetochore.

Global phosphoproteome profiling reveals unanticipated networks responsive to cisplatin treatment of embryonic stem cells

ABSTRACT Cellular responses to DNA-damaging agents involve the activation of various DNA damage signaling and transduction pathways. Using quantitative and high-resolution tandem mass spectrometry, we determined global changes in protein level and phosphorylation site profiles following treatment of SILAC (stable isotope labeling by amino acids in cell culture)-labeled murine embryonic stem cells with the anticancer drug cisplatin. Network and pathway analyses indicated that processes related to the DNA damage response and cytoskeleton organization were significantly affected. Although the ATM (ataxia telangiectasia mutated) and ATR (ATM and Rad3-related) consensus sequence (S/T-Q motif) was significantly overrepresented among hyperphosphorylated peptides, about half of the >2-fold-upregulated phosphorylation sites based on the consensus sequence were not direct substrates of ATM and ATR. Eleven protein kinases mainly belonging to the mitogen-activated protein kinase (MAPK) family were identified as being regulated in their kinase domain activation loop. The biological importance of three of these kinases (cyclin-dependent kinase 7 [CDK7], Plk1, and KPCD1) in the protection against cisplatin-induced cytotoxicity was demonstrated by small interfering RNA (siRNA)-mediated knockdown. Our results indicate that the cellular response to cisplatin involves a variety of kinases and phosphatases not only acting in the nucleus but also regulating cytoplasmic targets, resulting in extensive cytoskeletal rearrangements. Integration of transcriptomic and proteomic data revealed a poor correlation between changes in the relative levels of transcripts and their corresponding proteins, but a large overlap in affected pathways at the levels of mRNA, protein, and phosphoprotein. This study provides an integrated view of pathways activated by genotoxic stress and deciphers kinases that play a pivotal role in regulating cellular processes other than the DNA damage response.

qcML: An Exchange Format for Quality Control Metrics from Mass Spectrometry Experiments

Quality control is increasingly recognized as a crucial aspect of mass spectrometry based proteomics. Several recent papers discuss relevant parameters for quality control and present applications to extract these from the instrumental raw data. What has been missing, however, is a standard data exchange format for reporting these performance metrics. We therefore developed the qcML format, an XML-based standard that follows the design principles of the related mzML, mzIdentML, mzQuantML, and TraML standards from the HUPO-PSI (Proteomics Standards Initiative). In addition to the XML format, we also provide tools for the calculation of a wide range of quality metrics as well as a database format and interconversion tools, so that existing LIMS systems can easily add relational storage of the quality control data to their existing schema. We here describe the qcML specification, along with possible use cases and an illustrative example of the subsequent analysis possibilities. All information about qcML is available at http://code.google.com/p/qcml.

An Optimized Shotgun Strategy for the Rapid Generation of Comprehensive Human Proteomes

Summary This study investigates the challenge of comprehensively cataloging the complete human proteome from a single-cell type using mass spectrometry (MS)-based shotgun proteomics. We modify a classical two-dimensional high-resolution reversed-phase peptide fractionation scheme and optimize a protocol that provides sufficient peak capacity to saturate the sequencing speed of modern MS instruments. This strategy enables the deepest proteome of a human single-cell type to date, with the HeLa proteome sequenced to a depth of ∼584,000 unique peptide sequences and ∼14,200 protein isoforms (∼12,200 protein-coding genes). This depth is comparable with next-generation RNA sequencing and enables the identification of post-translational modifications, including ∼7,000 N-acetylation sites and ∼10,000 phosphorylation sites, without the need for enrichment. We further demonstrate the general applicability and clinical potential of this proteomics strategy by comprehensively quantifying global proteome expression in several different human cancer cell lines and patient tissue samples.