Search

Typical analyses of mass spectrometry data only identify amino acid sequences that exist in reference databases. This restricts the possibility of discovering new peptides such as those that contain uncharacterized mutations or originate from unexpected processing of RNAs and proteins. De novo peptide sequencing approaches address this limitation but often suffer from low accuracy and require extensive validation by experts. Here, we develop SMSNet, a deep learning-based de novo peptide sequencing framework that achieves >95% amino acid accuracy while retaining good identification coverage. Applications of SMSNet on landmark proteomics and peptidomics studies reveal over 10,000 previously uncharacterized HLA antigens and phosphopeptides, and in conjunction with database-search methods, expand the coverage of peptide identification by almost 30%. The power to accurately identify new peptides of SMSNet would make it an invaluable tool for any future proteomics and peptidomics studies, including tumor neoantigen discovery, antibody sequencing, and proteome characterization of non-model organisms.

Rabies is a deadly viral disease caused by the rabies virus (RABV), transmitted through a bite of an infected host, resulting in irreversible neurological symptoms and a 100% fatality rate in humans. Despite many aspects describing rabies neuropathogenesis, numerous hypotheses remain unanswered and concealed. Observations obtained from infected primary neurons or mouse brain samples are more relevant to human clinical rabies than permissive cell lines, however, limitations regarding the ethical issue and sample accessibility become a hurdle for discovering new insights into virus–host interplays. To better understand RABV pathogenesis in humans, we generated human-induced pluripotent stem cell (hiPSC)-derived neurons to offer the opportunity for an inimitable study of RABV infection at a molecular level in a pathologically relevant cell type. This study describes the characteristics and detailed proteomic changes of hiPSC-derived neurons in response to RABV infection using LC-MS/MS quantitative analysis. Gene ontology (GO) enrichment of differentially expressed proteins (DEPs) reveals temporal changes of proteins related to metabolic process, immune response, neurotransmitter transport/synaptic vesicle cycle, cytoskeleton organization, and cell stress response, demonstrating fundamental underlying mechanisms of neuropathogenesis in a time-course dependence. Lastly, we highlighted plausible functions of heat shock cognate protein 70 (HSC70 or HSPA8) that might play a pivotal role in regulating RABV replication and pathogenesis. Our findings acquired from this hiPSC-derived neuron platform help to define novel cellular mechanisms during RABV infection, which could be applicable to further studies to widen views of RABV-host interaction.

PIWI‐interacting RNAs (piRNAs) are germline‐specific small RNAs that form effector complexes with PIWI proteins (Piwi–piRNA complexes) and play critical roles for preserving genomic integrity by repressing transposable elements (TEs). Drosophila Piwi transcriptionally silences specific targets through heterochromatin formation and increases histone H3K9 methylation (H3K9me3) and histone H1 deposition at these loci, with nuclear RNA export factor variant Nxf2 serving as a co‐factor. Using ChEP and DamID‐seq, we now uncover a Piwi/Nxf2‐dependent target association with nuclear lamins. Hi‐C analysis of Piwi or Nxf2‐depleted cells reveals decreased intra‐TAD and increased inter‐TAD interactions in regions harboring Piwi–piRNA target TEs. Using a forced tethering system, we analyze the functional effects of Piwi–piRNA/Nxf2‐mediated recruitment of piRNA target regions to the nuclear periphery. Removal of active histone marks is followed by transcriptional silencing, chromatin conformational changes, and H3K9me3 and H1 association. Our data show that the Piwi–piRNA pathway can induce stepwise changes in nuclear architecture and chromatin state at target loci for transcriptional silencing., Formation of heterochromatin structure at piRNA target loci involves their Piwi/Nxf2‐dependent recruitment to the nuclear periphery.

Genome hybridization is an important evolutionary event that gives rise to species with novel capabilities. However, the merging of distinct genomes also brings together incompatible regulatory networks that must be resolved during the course of evolution. Understanding of the early stages of post-hybridization evolution is particularly important because changes in these stages have long-term evolutionary consequences. Here, via comparative transcriptomic analyses of two closely related, recently hybridized Trichosporon fungi, T. coremiiforme and T. ovoides, and three extant relatives, we show that early post-hybridization evolutionary processes occur separately at the gene sequence and gene expression levels but together contribute to the stabilization of hybrid genome and transcriptome. Our findings also highlight lineage-specific consequences of genome hybridization, revealing that the transcriptional regulatory dynamics in these hybrids responded completely differently to gene loss events: one involving both subgenomes and another that is strictly subgenome-specific. Sira Sriswasdi et al. use comparative transcriptomics to study the resolution of genomic shock, particularly in transcriptional regulation, in closely related Trichosporon hybrids and their extant relatives. They show that early post-hybridization evolution occurs separately at the gene sequence and gene expression levels.

Non-invasive prenatal diagnosis (NIPD) of isolated cell-free DNA from maternal plasma has been applied to detect monogenic diseases in the fetus. Droplet digital PCR (ddPCR) is a sensitive and quantitative technique for NIPD. In the present study, the development and evaluation of ddPCR-based assays for common α and β-thalassemia variants amongst the Asian population was described; specifically, Southeast Asian (SEA) deletion, HbE, and 41/42 (-CTTT). SEA is caused by deletion of a 20 kb region surrounding the α-globin gene, whilst HbE and 41/42 (-CTTT) are caused by point mutations on the β-globin gene. Cell-free DNA samples from 46 singleton pregnant women who were carriers of these mutations were isolated and quantified using ddPCR with specially designed probes for each target allele. Allelic copy number calculation and likelihood ratio tests were used to classify fetal genotypes. Classification performances were evaluated against ground truth fetal genotypes obtained from conventional amniocentesis. Copy number variation analysis of SEA deletion accurately classified fetal genotypes in 20 out of 22 cases with an area under the receiver operating characteristic curve of 0.98 for detecting Hb Bart's hydrops fetalis. For HbE cases, 10 out of 16 samples were correctly classified, and three were inconclusive. For 41/42 (-CTTT) cases, 2 out of 8 were correctly classified, and four were inconclusive. The correct genotype was not rejected in any inconclusive case and may be resolved with additional ddPCR experiments. These results indicate that ddPCR-based analysis of maternal plasma can become an accurate and effective NIPD for SEA deletion α-(0) thalassemia. Although the performance of ddPCR on HbE and 41/42 (-CTTT) mutations were not sufficient for clinical application, these results may serve as a foundation for future works in this field.

Cholangiocarcinoma (CCA), a lethal malignancy of the biliary epithelium, is the second most common primary liver cancer. The poor prognosis of CCA is due to the high rate of tumour invasion and distant metastasis. We found that the RNA-binding protein LIN28B, a known regulator of microRNA biogenesis, stem cell maintenance, and oncogenesis, is expressed in a subpopulation of CCA patients. To further investigate the potential role of LIN28B in CCA pathogenesis, we studied the effect of LIN28B overexpression in the cholangiocyte cell line MMNK-1 and cholangiocarcinoma cell lines HuCCT-1 and KKU-214. Here, we show that enhanced LIN28B expression promoted cancer stem cell-like properties in CCA, including enhanced cell migration, epithelial-to-mesenchymal transition (EMT), increased cell proliferation and spheroid formation. Proteomic analysis revealed TGF-β-induced protein (TGFBI) as a novel LIN28B target gene, and further analysis showed upregulation of other components of the TGF-β signalling pathway, including TGF-β receptor type I (TGFBRI) expression and cytokine TGFB-I, II and III secretion. Importantly, the small molecule TGF-β inhibitor SB431542 negated the effects of LIN28B on both cell migration and clonogenic potential. Overexpression of TGFBI alone promoted cholangiocarcinoma cell migration and EMT changes, but not spheroid formation, suggesting that TGFBI partially contributes to LIN28B-mediated aggressive cell behaviour. These observations are consistent with a model in which TGF-β and LIN28B work together to form a positive feedback loop during cholangiocarcinoma metastasis and provide a therapeutic intervention opportunity.

To construct a backbone tree consisting of basidiomycetous yeasts, draft genome sequences from 25 species of Trichosporonales (Tremellomycetes, Basidiomycota) were generated. In addition to the hybrid genomes of Trichosporon coremiiforme and Trichosporon ovoides that we described previously, we identified an interspecies hybrid genome in Cutaneotrichosporon mucoides (formerly Trichosporon mucoides). This hybrid genome had a gene retention rate of ~55%, and its closest haploid relative was Cutaneotrichosporon dermatis. After constructing the C. mucoides subgenomes, we generated a phylogenetic tree using genome data from the 27 haploid species and the subgenome data from the three hybrid genome species. It was a high-quality tree with 100% bootstrap support for all of the branches. The genome-based tree provided superior resolution compared with previous multi-gene analyses. Although our backbone tree does not include all Trichosporonales genera (e.g. Cryptotrichosporon), it will be valuable for future analyses of genome data. Interest in interspecies hybrid fungal genomes has recently increased because they may provide a basis for new technologies. The three Trichosporonales hybrid genomes described in this study are different from well-characterized hybrid genomes (e.g. those of Saccharomyces pastorianus and Saccharomyces bayanus) because these hybridization events probably occurred in the distant evolutionary past. Hence, they will be useful for studying genome stability following hybridization and speciation events. Copyright © 2017 John Wiley & Sons, Ltd.

Anion exchanger 1 (AE1) is a critical transporter and the primary structural scaffold for large macromolecular complexes responsible for erythrocyte membrane flexibility and integrity. We used zero-length crosslinking and mass spectrometry to probe AE1 structures and interactions in intact erythrocyte membranes. An experimentally verified full-length model of AE1 dimers was developed by combining crosslink-defined distance constraints with homology modeling. Previously unresolved cytoplasmic loops in the AE1 C-terminal domain are packed at the domain-domain interface on the cytoplasmic face of the membrane where they anchor the N-terminal domain’s location and prevent it from occluding the ion channel. Crosslinks between AE1 dimers and ankyrin-1 indicate the likely topology for AE1 tetramers and suggest that ankyrin-1 wraps around AE1 tetramers, which may stabilize this oligomer state. This interaction and interactions of AE1 with other major erythrocyte membrane proteins show that protein-protein contacts are often substantially more extensive than previously reported.

Typical analyses of mass spectrometry data only identify amino acid sequences that exist in reference databases. This restricts the possibility of discovering new peptides such as those that contain uncharacterized mutations or originate from unexpected processing of RNAs and proteins. De novo peptide sequencing approaches address this limitation but often suffer from low accuracy and require extensive validation by experts. Here, we develop SMSNet, a deep learning-based hybrid de novo peptide sequencing framework that achieves >95% amino acid accuracy while retaining good identification coverage. Applications of SMSNet on landmark proteomics and peptideomics studies reveal over 10,000 previously uncharacterized HLA antigens and phosphopeptides and in conjunction with database-search methods, expand the coverage of peptide identification by almost 30%. The power to accurately identify new peptides of SMSNet would make it an invaluable tool for any future proteomics and peptidomics studies – especially cancer neoantigen discovery and proteome characterization of non-model organisms.

Delineation and characterization of genera in Trichosporonales (Agaricomycotina, Basidiomycota) was performed using 24 haploid and 3 naturally occurring hybrid genomes, with 3 Tremellales genomes used as outgroups. Orthologous group analysis of those genomes showed presence-absence patterns of orthologs that were consistent with the genus classifications. Many shared unique orthologs were identified in the well-supported lineages (genera Apiotrichum and Trichosporon), supporting the definitions of the genera Apiotrichum and Trichosporon from a genomic perspective. Specifically, we obtained 24 and 285 genus-specific genes from eight Apiotrichum and five Trichosporon species, respectively, and propose that these genus-specific genes can be used for delineation of those genera. On the other hand, the genus Cutaneotrichosporon shared only one genus-specific gene among eight genomes, indicating that this genus definition might require re-examination based on genomic data. In addition, taxonomic revisions are presented in this study, including the proposal of two genera, Pascua and Prillingera. Because genomic data can be systematically obtained and analyzed to compare species from a comprehensive viewpoint, they can be used not only to reconstruct reliable phylogenetic trees, but also to re-examine the definitions of taxonomic classifications. To our knowledge, this is the first report to discuss the 'natural system' of genus level classification in fungi based on genomic data.

Immunotherapy is an emerging approach in cancer treatment that activates the host immune system to destroy cancer cells that express unique peptide signatures (neoepitopes). Administrations of cancer-specific neoepitopes in the form of synthetic peptide vaccine have been proven effective in both mouse models and human patients. Because only a tiny fraction of cancer-specific neoepitopes actually elicit immune response, selection of potent, immunogenic neoepitopes remains a challenging step in cancer vaccine development. A basic approach for immunogenicity prediction is based on the premise that effective neoepitope should bind with the Major Histocompatibility Complex (MHC) with high affinity. In this study, we developed MHCSeqNet, an open-source deep learning model, which not only competes favorably with state-of-the-art predictors but also exhibits promising generalization to new MHC class I and class II alleles. MHCSeqNet employed neural network architectures developed for natural language processing to model amino acid sequence representations of MHC allele and epitope peptide as sentences with amino acids as individual words. This consideration allows MHCSeqNet to accept new MHC alleles as well as peptides of any length. The flexibility offered by MHCSeqNet should make it a valuable tool for screening effective neoepitopes in cancer vaccine development.Availability and implementationhttps://github.com/cmbcu/MHCSeqNet

Red blood cells have been extensively studied but many questions regarding membrane properties and pathophysiology remain unanswered. Proteome analysis of red cell membranes is complicated by a very wide dynamic range of protein concentrations as well as the presence of proteins that are very large, very hydrophobic, or heterogeneously glycosylated. This study investigated the removal of other blood cell types, red cell membrane extraction, differing degrees of fractionation using 1-D SDS gels, and label-free quantitative methods to determine optimized conditions for proteomic comparisons of clinical blood samples. The results showed that fractionation of red cell membranes on 1-D SDS gels was more efficient than low-ionic-strength extractions followed by 1-D gel fractionation. When gel lanes were sliced into 30 uniform slices, a good depth of analysis that included the identification of most well-characterized, low-abundance red cell membrane proteins including those present at 500 to 10,000 copies per cell was obtained. Furthermore, the size separation enabled detection of changes due to proteolysis or in vivo protein crosslinking. A combination of Rosetta Elucidator quantitation and subsequent statistical analysis enabled the robust detection of protein differences that could be used to address unresolved questions in red cell disorders. This article is part of a Special Issue entitled: Integrated omics.

Polyploidization events such as whole-genome duplication and inter-species hybridization are major evolutionary forces that shape genomes. Although long-term effects of polyploidization have been well-characterized, early molecular evolutionary consequences of polyploidization remain largely unexplored. Here, we report the discovery of two recent and independent genome hybridizations within a single clade of a fungal genus, Trichosporon. Comparative genomic analyses revealed that redundant genes are experiencing decelerations, not accelerations, of evolutionary rates. We identified a relationship between gene conversion and decelerated evolution suggesting that gene conversion may improve the genome stability of young hybrids by restricting gene functional divergences. Furthermore, we detected large-scale gene losses from transcriptional and translational machineries that indicate a global compensatory mechanism against increased gene dosages. Overall, our findings illustrate counteracting mechanisms during an early phase of post-genome hybridization and fill a critical gap in existing theories on genome evolution.

Ectopic pregnancy (EP) and normal intrauterine pregnancy (IUP) serum proteomes were quantitatively compared to systematically identify candidate biomarkers. A 3-D biomarker discovery strategy consisting of abundant protein immunodepletion, SDS gels, LC-MS/MS, and label-free quantitation of MS signal intensities identified 70 candidate biomarkers with differences between groups greater than 2.5-fold. Further statistical analyses of peptide quantities were used to select the most promising 12 biomarkers for further study, which included known EP biomarkers, novel EP biomarkers (ADAM12 and ISM2), and five specific isoforms of the pregnancy specific beta-1-glycoprotein family. Technical replicates showed good reproducibility and protein intensities from the label-free discovery analysis compared favorably with reported abundance levels of several known reference serum proteins over at least three orders of magnitude. Similarly, relative abundances of candidate biomarkers from the label-free discovery analysis were consistent with relative abundances from pilot validation assays performed for five of the 12 most promising biomarkers using label-free multiple reaction monitoring of both the patient serum pools used for discovery and the individual samples that constituted these pools. These results demonstrate robust, reproducible, in-depth 3-D serum proteome discovery, and subsequent pilot-scale validation studies can be achieved readily using label-free quantitation strategies.

Peroxiredoxin-6 (PRDX6) is an unusual member of the peroxiredoxin family of antioxidant enzymes that has only one evolutionarily conserved cysteine. It reduces oxidized lipids and reactive oxygen species (ROS) by oxidation of the active-site cysteine (Cys47) to a sulfenic acid, but the mechanism for conversion back to a thiol is not completely understood. Moreover, it has phospholipase A2 (PLA2) activity in addition to its peroxidase activity. Interestingly, some biochemical data are inconsistent with a known high-resolution crystal structure of the catalytic intermediate of the protein, and biophysical data indicate that the protein undergoes conformational changes that affect enzyme activity. In order to further elucidate the solution structure of this important enzyme, we used chemical cross-linking coupled with high-resolution MS (CX–MS), with an emphasis on zero-length cross-links. Distance constraints from high confidence cross-links were used in homology modelling experiments to determine a solution structure of the reduced form of the protein. This structure was further evaluated using chemical cross-links produced by several homo-bifunctional amine-reactive cross-linking reagents, which helped to confirm the solution structure. The results show that several regions of the reduced version of human PRDX6 are in a substantially different conformation from that shown for the crystal structure of the peroxidase catalytic intermediate. The differences between these two structures are likely to reflect catalysis-related conformational changes. These studies also demonstrate that CX–MS using zero-length cross-linking is a powerful strategy for probing protein conformational changes that is complementary to alternative methods such as crystallographic, NMR and biophysical studies.

Pleiotropy refers to the phenomenon in which a single gene controls several distinct, and seemingly unrelated, phenotypic effects. We use C. elegans early embryogenesis as a model to conduct systematic studies of pleiotropy. We analyze high-throughput RNA interference (RNAi) data from C. elegans and identify “phenotypic signatures”, which are sets of cellular defects indicative of certain biological functions. By matching phenotypic profiles to our identified signatures, we assign genes with complex phenotypic profiles to multiple functional classes. Overall, we observe that pleiotropy occurs extensively among genes involved in early embryogenesis, and a small proportion of these genes are highly pleiotropic. We hypothesize that genes involved in early embryogenesis are organized into partially overlapping functional modules, and that pleiotropic genes represent “connectors” between these modules. In support of this hypothesis, we find that highly pleiotropic genes tend to reside in central positions in protein-protein interaction networks, suggesting that pleiotropic genes act as connecting points between different protein complexes or pathways., Author Summary In a biological system, some genes play single roles while others perform multiple functions. How can we determine which genes are multi-functional? An informative way for probing gene functions is to eliminate the expression of a given gene and observe the phenotypic consequences. RNAi techniques have enabled the generation of genome-wide phenotypic data. Conventionally, genes are clustered into mutually exclusive categories according to the observed defects following RNAi. However, assigning genes that may play multiple roles exclusively into a single category is arbitrary. This paper works out a computational approach that categorizes genes while allowing assignment of genes with complex phenotypes into multiple categories. We apply this approach to genes involved in cell divisions of C. elegans early embryos, and find that about half of these genes can be assigned to more than one functional category. This approach has allowed the identification of previously undiscovered gene functions. We also find that genes playing many roles in early embryos tend to reside in central positions in protein networks. Our approach can be used to perform functional annotations based on phenotypic data in other systems and to identify genes that coordinate multiple biological functions.

Diamond Blackfan Anemia (DBA) is a congenital anemia usually caused by diverse mutations in ribosomal proteins. Although the genetics of DBA are well characterized, the mechanisms that lead to macrocytic anemia remain unclear. We systematically analyzed the proteomes of red blood cell membranes from multiple DBA patients to determine whether abnormalities in protein translation or erythropoiesis contribute to the observed macrocytosis or alterations in the mature red blood cell membrane. In depth proteome analysis of red cell membranes enabled highly reproducible identification and quantitative comparisons of 1100 or more proteins. These comparisons revealed clear differences between red cell membrane proteomes in DBA patients and healthy controls that were consistent across DBA patients with different ribosomal gene mutations. Proteins exhibiting changes in abundance included those known to be increased in DBA such as fetal hemoglobin and a number of proteins not normally found in mature red cell membranes, including proteins involved in the major histocompatibility complex class I pathway. Most striking was the presence of dysferlin in the red blood cell membranes of DBA patients but absent in healthy controls. Immunoblot validation using red cell membranes isolated from additional DBA patients and healthy controls confirmed a distinct membrane protein signature specific to patients with DBA.