Comprehensive analysis of amino acid and nucleotide composition in eukaryotic genomes, comparing genes and pseudogenes

Facile Access to Branched Multispecific Proteins

Approaches that leverage orthogonal chemical reactions to generate protein-protein conjugates have expanded access to bespoke chimeras. Although the literature is replete with examples of the semisynthesis of bispecific proteins, few methods exist for the semisynthesis of protein conjugates of higher complexity (i.e., greater than two-protein fusions). The recent emergence of trispecific cell engagers for immune cell redirection therapies necessitates the development of chemical methods for the construction of trispecific proteins that would otherwise be inaccessible via natural protein synthesis. Here, we demonstrate that 3-bromo-5-methylene pyrrolone (3Br-5MP) can be used to effect the facile chemical synthesis of trispecific peptides and proteins with exquisite control over the addition of each monomer. The multimeric complexes maintain epitope functionality both in human cells and upon immobilization. We anticipate that facile access to trispecific proteins using this 3Br-5MP will have broad utility in basic science research and will quicken the pace of research to establish novel, multimeric immune cell redirection therapies.

An Update on Protein Kinases as Therapeutic Targets-Part I: Protein Kinase C Activation and Its Role in Cancer and Cardiovascular Diseases

Protein kinases are one of the most significant drug targets in the human proteome, historically harnessed for the treatment of cancer, cardiovascular disease, and a growing number of other conditions, including autoimmune and inflammatory processes. Since the approval of the first kinase inhibitors in the late 1990s and early 2000s, the field has grown exponentially, comprising 98 approved therapeutics to date, 37 of which were approved between 2016 and 2021. While many of these small-molecule protein kinase inhibitors that interact orthosterically with the protein kinase ATP binding pocket have been massively successful for oncological indications, their poor selectively for protein kinase isozymes have limited them due to toxicities in their application to other disease spaces. Thus, recent attention has turned to the use of alternative allosteric binding mechanisms and improved drug platforms such as modified peptides to design protein kinase modulators with enhanced selectivity and other pharmacological properties. Herein we review the role of different protein kinase C (PKC) isoforms in cancer and cardiovascular disease, with particular attention to PKC-family inhibitors. We discuss translational examples and carefully consider the advantages and limitations of each compound (Part I). We also discuss the recent advances in the field of protein kinase modulators, leverage molecular docking to model inhibitor-kinase interactions, and propose mechanisms of action that will aid in the design of next-generation protein kinase modulators (Part II).

Use of pyridazinediones for tuneable and reversible covalent cysteine modification applied to peptides, proteins and hydrogels

Reversible cysteine modification has been found to be a useful tool for a plethora of applications such as selective enzymatic inhibition, activity-based protein profiling and/or cargo release from a protein or a material. However, only a limited number of reagents display reliable dynamic/reversible thiol modification and, in most cases, many of these reagents suffer from issues of stability, a lack of modularity and/or poor rate tunability. In this work, we demonstrate the potential of pyridazinediones as novel reversible and tuneable covalent cysteine modifiers. We show that the electrophilicity of pyridazinediones correlates to the rates of the Michael addition and retro-Michael deconjugation reactions, demonstrating that pyridazinediones provide an enticing platform for readily tuneable and reversible thiol addition/release. We explore the regioselectivity of the novel reaction and unveil the reason for the fundamental increased reactivity of aryl bearing pyridazinediones by using DFT calculations and corroborating findings with SCXRD. We also applied this fundamental discovery to making more rapid disulfide rebridging agents in related work. We finally provide the groundwork for potential applications in various areas with exemplification using readily functionalised "clickable" pyridazinediones on clinically relevant cysteine and disulfide conjugated proteins, as well as on a hydrogel material.

Methodology of stable peptide based on propargylated sulfonium

Peptides can be used as effective molecular tool for covalent modification of proteins and play important roles in ligand directed covalent modification. Tyr-selective protein modifications exert a profound impact on protein functionality. Here, we developed a general strategy that involves nucleophilic addition of alkyne for tyrosine modification. The terminal alkyne of propargyl sulfonium is motivated by the sulfonium center to react with phenolic hydroxyl. This approach provides a straightforward method for tyrosine modification due to its high yield in aqueous solution at physiological temperature. In addition, cyclic peptides could be obtained via adjusting pH to 8.0 from peptides consisting of tyrosine and methionine modified by propargyl bromide, and the resulting cyclic peptides are proved to have better stability, excellent 2-mercaptopyridine resistance and improved cellular uptakes. Furthermore, molecules made from the propargylated sulfonium have the potential to be used as warheads against tyrosine containing biomolecules. Collectively, we develop a direct and uncomplicated technique for modifying tyrosine residues, the strategy concerned can be widely utilized to construct stable peptides and biomolecules imaging.

Structural characterization of protein-material interfacial interactions using lysine reactivity profiling-mass spectrometry

Understanding how proteins and materials interact is useful for evaluating the safety of biomedical micro/nanomaterials, toxicity estimation and design of nano-drugs and catalytic activity improvement of bio-inorganic functional hybrids. However, characterizing the interfacial molecular details of protein-micro/nanomaterial hybrids remains a great challenge. This protocol describes the lysine reactivity profiling-mass spectrometry strategy for determining which parts of a protein are interacting with the micro/nanomaterials. Lysine residues occur frequently on hydrophilic protein surfaces, and their reactivity is dependent on the accessibility of their amine groups. The accessibility of a lysine residue is lower when it is in contact with another object; allosteric effects resulting from this interaction might reduce or increase the reactivity of remote lysine residues. Lysine reactivity is therefore a useful indicator of protein localization orientation, interaction sequence regions, binding sites and modulated protein structures in the protein-material hybrids. We describe the optimized two-step isotope dimethyl labeling strategy for protein-material hybrids under their native and denaturing conditions in sequence. The comparative quantification results of lysine reactivity are only dependent on the native microenvironments of lysine local structures. We also highlight other critical steps including protein digestion, elution from materials, data processing and interfacial structure analysis. The two-step isotope labeling steps need ~5 h, and the whole protocol including digestion, liquid chromatography-tandem mass spectrometry, data processing and structure analysis needs ~3-5 d.

Predicting small molecule binding pockets on diacylglycerol kinases using chemoproteomics and AlphaFold

Diacylglycerol kinases (DGKs) are metabolic kinases involved in regulating cellular levels of diacylglycerol and phosphatidic lipid messengers. The development of selective inhibitors for individual DGKs would benefit from discovery of protein pockets available for inhibitor binding in cellular environments. Here we utilized a sulfonyl-triazole probe (TH211) bearing a DGK fragment ligand for covalent binding to tyrosine and lysine sites on DGKs in cells that map to predicted small molecule binding pockets in AlphaFold structures. We apply this chemoproteomics-AlphaFold approach to evaluate probe binding of DGK chimera proteins engineered to exchange regulatory C1 domains between DGK subtypes (DGKα and DGKζ). Specifically, we discovered loss of TH211 binding to a predicted pocket in the catalytic domain when C1 domains on DGKα were exchanged that correlated with impaired biochemical activity as measured by a DAG phosphorylation assay. Collectively, we provide a family-wide assessment of accessible sites for covalent targeting that combined with AlphaFold revealed predicted small molecule binding pockets for guiding future inhibitor development of the DGK superfamily.

Modular Diazo Compound for the Bioreversible Late-Stage Modification of Proteins

We introduce a versatile strategy for the bioreversible modification of proteins. Our strategy is based on a tricomponent molecule, synthesized in three steps, that incorporates a diazo moiety for chemoselective esterification of carboxyl groups, a pyridyl disulfide group for late-stage functionalization with thiolated ligands, and a self-immolative carbonate group for esterase-mediated cleavage. Using cytochrome c (Cyt c) and the green fluorescent protein (GFP) as models, we generated protein conjugates modified with diverse domains for cellular delivery that include a small molecule, targeting and cell-penetrating peptides (CPPs), and a large polysaccharide. As a proof of concept, we used our strategy to effect the delivery of proteins into the cytosol of live mammalian cells in the presence of serum. The cellular delivery of functional Cyt c, which induces apoptosis, highlighted the advantage of bioreversible conjugation on a carboxyl group versus irreversible conjugation on an amino group. The ease and utility of this traceless modification provide new opportunities for chemical biologists.

Ultrasensitive HPLC-MS Quantification of S-(2-Succino) Cysteine Based on Ethanol/Acetyl Chloride Derivatization in Fumarate Accumulation Cells

Succination is a nonenzymatic and irreversible post-translational modification (PTM) with important biological significance, yielding S-(2-succino) cysteine (2SC) residue. This PTM is low in abundance and often requires a large amount of protein samples for 2SC quantification. In this work, an efficient quantification method based on ethanol/acetyl chloride chemical derivatization was developed. The three carboxyl groups of 2SC were all esterified to increase hydrophobicity, greatly improving its ionization efficiency. The sensitivity was increased by 112 times; the limit of detection was reduced to 0.885 fmol, and the protein usage was reduced by at least 10 times. The established method was used to detect the overall concentration of 2SC in fumarate accumulation cells quantitatively.

Tyrosine bioconjugation with hypervalent iodine

Hypervalent iodine reagents have recently emerged as powerful tools for late-stage peptide and protein functionalization. Herein we report a tyrosine bioconjugation methodology for the introduction of hypervalent iodine onto biomolecules under physiological conditions. Tyrosine residues were engaged in a selective addition onto the alkynyl bond of ethynylbenziodoxolones (EBX), resulting in stable vinylbenziodoxolones (VBX) bioconjugates. The methodology was successfully applied to peptides and proteins and tolerated all other nucleophilic residues, with the exception of cysteine. The generated VBX were further functionalized by palladium-catalyzed cross-coupling and azide-alkyne cycloaddition reactions. The method could be successfully used to modify bioactive natural products and native streptavidin to enable thiol-mediated cellular uptake.

Context‐dependence of the Reactivity of Cysteine and Lysine Residues

The S‐alkylation of Cys residues with a maleimide and the N^ϵ‐acylation of Lys residues with an N‐hydroxysuccinimide (NHS) ester are common methods for bioconjugation. Using Cys and Lys derivatives as proxies, we assessed differences in reactivity depending on the position of Cys or Lys in a protein sequence. We find that Cys position is exploitable to improve site‐selectivity in maleimide‐based modifications. Reactivity decreases substantially in the order N‐terminal>in‐chain>C‐terminal Cys due to modulation of sulfhydryl pK_a by the α‐ammonium and carboxylate groups at the termini. A lower pK_a value yields a larger fraction thiolate, which promotes selectivity while somewhat decreasing thiolate nucleophilicity in accord with βnuc =0.41. Lowering pH and salt concentration enhances selectivity still further. In contrast, differences in the reactivity of Lys towards an NHS ester were modest due to an appreciable decrease in amino group nucleophilicity with a lower pK_a of its conjugate acid. Hence, site‐selective Lys modification protocols will require electrophiles other than NHS esters.

Reactive chemistry for covalent probe and therapeutic development

Bioactive small molecules that form covalent bonds with a target protein are important tools for basic research and can be highly effective drugs. This review highlights reactive groups found in a collection of thiophilic and oxophilic drugs that mediate pharmacological activity through a covalent mechanism of action (MOA). We describe the application of advanced proteomic and bioanalytical methodologies for assessing selectivity of these covalent agents to guide and inspire the search for additional electrophiles suitable for covalent probe and therapeutic development. While the emphasis is on chemistry for modifying catalytic serine, threonine or cysteine residues, we devote a substantial fraction of the review to a collection of exploratory reactive groups of understudied residues on proteins.

Protein Modification at Tyrosine with Iminoxyl Radicals

Post-translational modifications (PTMs) of proteins are a biological mechanism for reversibly controlling protein function. Synthetic protein modifications (SPMs) at specific canonical amino acids can mimic PTMs. However, reversible SPMs at hydrophobic amino acid residues in proteins are especially limited. Here, we report a tyrosine (Tyr)-selective SPM utilizing persistent iminoxyl radicals, which are readily generated from sterically hindered oximes via single-electron oxidation. The reactivity of iminoxyl radicals with Tyr was dependent on the steric and electronic demands of oximes; isopropyl methyl piperidinium oxime 1f formed stable adducts, whereas the reaction of tert-butyl methyl piperidinium oxime 1o was reversible. The difference in reversibility between 1f and 1o, differentiated only by one methyl group, is due to the stability of iminoxyl radicals, which is partly dictated by the bond dissociation energy of oxime O-H groups. The Tyr-selective modifications with 1f and 1o proceeded under physiologically relevant, mild conditions. Specifically, the stable Tyr-modification with 1f introduced functional small molecules, including an azobenzene photoswitch, to proteins. Moreover, masking critical Tyr residues by SPM with 1o, and subsequent deconjugation triggered by the treatment with a thiol, enabled on-demand control of protein functions. We applied this reversible Tyr modification with 1o to alter an enzymatic activity and the binding affinity of a monoclonal antibody with an antigen upon modification/deconjugation. The on-demand ON/OFF switch of protein functions through Tyr-selective and reversible covalent-bond formation will provide unique opportunities in biological research and therapeutics.

CBRPP: a new RNA-centric method to study RNA-protein interactions

RNA and protein are interconnected biomolecules that can influence each other's life cycles and functions through physical interactions. Abnormal RNA-protein interactions lead to cell dysfunctions and human diseases. Therefore, mapping networks of RNA-protein interactions is crucial for understanding cellular processes and pathogenesis of related diseases. Different practical protein-centric methods for studying RNA-protein interactions have been reported, but few robust RNA-centric methods exist. Here, we developed CRISPR-based RNA proximity proteomics (CBRPP), a new RNA-centric method to identify proteins associated with an endogenous RNA of interest in native cellular context without pre-editing of the target RNA, cross-linking or RNA-protein complexes manipulation in vitro. CBRPP is based on a fusion of dCas13 and proximity-based labelling (PBL) enzyme. dCas13 can deliver PBL enzyme to the target RNA with high specificity, while PBL enzyme labels the surrounding proteins of the target RNA, which are then identified by mass spectrometry.

Site-selective tyrosine bioconjugation via photoredox catalysis for native-to-bioorthogonal protein transformation

ABSTACT

The growing prevalence of synthetically modified proteins in pharmaceuticals and materials has exposed the need for efficient strategies to enable chemical modifications with high site-selectivity. While genetic engineering can incorporate non-natural amino acids into recombinant proteins, regioselective chemical modification of wild-type proteins remains a challenge. Herein, we use photoredox catalysis to develop a site-selective tyrosine bioconjugation pathway that incorporates bioorthogonal formyl groups, which subsequently allows for the synthesis of structurally defined fluorescent conjugates from native proteins. A water-soluble photocatalyst, lumiflavin, has been shown to induce oxidative coupling between a previously unreported phenoxazine dialdehyde tag and a single tyrosine site, even in the presence of multiple tyrosyl side chains, through the formation of a covalent C-N bond. A variety of native proteins, including those with multiple tyrosines, can successfully undergo both tyrosine-specific and single-site-selective labelling. This technology directly introduces aldehyde moieties onto native proteins, enabling rapid product diversification using an array of well-established bioorthogonal functionalization protocols including the alkyne-azide click reaction.

Two-Step Synthesis of α-Aryl-α-diazoamides as Modular Bioreversible Labels

α-Aryl-α-diazoamides were synthesized in two steps under mild conditions. This expeditious route employs Pd-catalyzed C-H arylation of N-succinimidyl 2-diazoacetate to obtain N-succinimidyl 2-aryl-2-diazoacetates, followed by aminolysis. The ensuing diazo compounds can esterify carboxyl groups in aqueous solution, and the ester products are substrates for an esterase. The broad scope of the synthetic route enables the continued development of diazo compounds in chemical biology.

Global Profiling of Lysine Accessibility to Evaluate Protein Structure Changes in Alzheimer's Disease

The linear sequence of amino acids in a protein folds into a 3D structure to execute protein activity and function, but it is still challenging to profile the 3D structure at the proteome scale. Here, we present a method of native protein tandem mass tag (TMT) profiling of Lys accessibility and its application to investigate structural alterations in human brain specimens of Alzheimer's disease (AD). In this method, proteins are extracted under a native condition, labeled by TMT reagents, followed by trypsin digestion and peptide analysis using two-dimensional liquid chromatography and tandem mass spectrometry (LC/LC-MS/MS). The method quantifies Lys labeling efficiency to evaluate its accessibility on the protein surface, which may be affected by protein conformations, protein modifications, and/or other molecular interactions. We systematically optimized the amount of TMT reagents, reaction time, and temperature and then analyzed protein samples under multiple conditions, including different labeling time (5 and 30 min), heat treatment, AD and normal human cases. The experiment profiled 15370 TMT-labeled peptides in 4475 proteins. As expected, the heat treatment led to extensive changes in protein conformations, with 17% of the detected proteome displaying differential labeling. Compared to the normal controls, AD brain showed different Lys accessibility of tau and RNA splicing complexes, which are the hallmarks of AD pathology, as well as proteins involved in transcription, mitochondrial, and synaptic functions. To eliminate the possibility that the observed differential Lys labeling was caused by protein level change, the whole proteome was quantified with standard TMT-LC/LC-MS/MS for normalization. Thus, this native protein TMT method enables the proteome-wide measurement of Lys accessibility, representing a straightforward strategy to explore protein structure in any biological system.

CBRPP: a new RNA-centric method to study RNA-protein interactions

RNA and protein are interconnected biomolecules that can influence each other's life cycles and functions through physical interactions. Abnormal RNA-protein interactions lead to cell dysfunctions and human diseases. Therefore, mapping networks of RNA-protein interactions is crucial for understanding cellular processes and pathogenesis of related diseases. Different practical protein-centric methods for studying RNA-protein interactions have been reported, but few robust RNA-centric methods exist. Here, we developed CRISPR-based RNA proximity proteomics (CBRPP), a new RNA-centric method to identify proteins associated with an endogenous RNA of interest in native cellular context without pre-editing of the target RNA, cross-linking or RNA-protein complexes manipulation in vitro. CBRPP is based on a fusion of dCas13 and proximity-based labelling (PBL) enzyme. dCas13 can deliver PBL enzyme to the target RNA with high specificity, while PBL enzyme labels the surrounding proteins of the target RNA, which are then identified by mass spectrometry.

Sequence Characterization and Molecular Modeling of Clinically Relevant Variants of the SARS-CoV-2 Main Protease

The SARS-CoV-2 main protease (Mpro) is essential to viral replication and cleaves highly specific substrate sequences, making it an obvious target for inhibitor design. However, as for any virus, SARS-CoV-2 is subject to constant neutral drift and selection pressure, with new Mpro mutations arising over time. Identification and structural characterization of Mpro variants is thus critical for robust inhibitor design. Here we report sequence analysis, structure predictions, and molecular modeling for seventy-nine Mpro variants, constituting all clinically observed mutations in this protein as of April 29, 2020. Residue substitution is widely distributed, with some tendency toward larger and more hydrophobic residues. Modeling and protein structure network analysis suggest differences in cohesion and active site flexibility, revealing patterns in viral evolution that have relevance for drug discovery.

Proximity-dependent labeling methods for proteomic profiling in living cells: An update

Characterizing the proteome composition of organelles and subcellular regions of living cells can facilitate the understanding of cellular organization as well as protein interactome networks. Proximity labeling-based methods coupled with mass spectrometry (MS) offer a high-throughput approach for systematic analysis of spatially restricted proteomes. Proximity labeling utilizes enzymes that generate reactive radicals to covalently tag neighboring proteins. The tagged endogenous proteins can then be isolated for further analysis by MS. To analyze protein-protein interactions or identify components that localize to discrete subcellular compartments, spatial expression is achieved by fusing the enzyme to specific proteins or signal peptides that target to particular subcellular regions. Although these technologies have only been introduced recently, they have already provided deep insights into a wide range of biological processes. Here, we provide an updated description and comparison of proximity labeling methods, as well as their applications and improvements. As each method has its own unique features, the goal of this review is to describe how different proximity labeling methods can be used to answer different biological questions. This article is categorized under: Technologies > Analysis of Proteins.

Transcriptional activity and strain-specific history of mouse pseudogenes

Pseudogenes are ideal markers of genome remodelling. In turn, the mouse is an ideal platform for studying them, particularly with the recent availability of strain-sequencing and transcriptional data. Here, combining both manual curation and automatic pipelines, we present a genome-wide annotation of the pseudogenes in the mouse reference genome and 18 inbred mouse strains (available via the mouse.pseudogene.org resource). We also annotate 165 unitary pseudogenes in mouse, and 303, in human. The overall pseudogene repertoire in mouse is similar to that in human in terms of size, biotype distribution, and family composition (e.g. with GAPDH and ribosomal proteins being the largest families). Notable differences arise in the pseudogene age distribution, with multiple retro-transpositional bursts in mouse evolutionary history and only one in human. Furthermore, in each strain about a fifth of all pseudogenes are unique, reflecting strain-specific evolution. Finally, we find that ~15% of the mouse pseudogenes are transcribed, and that highly transcribed parent genes tend to give rise to many processed pseudogenes.

Non-canonical functions of human cytoplasmic tyrosyl-, tryptophanyl- and other aminoacyl-tRNA synthetases

Aminoacyl-tRNA synthetases catalyze the aminoacylation of their cognate tRNAs. Here we review the accumulated knowledge of non-canonical functions of human cytoplasmic aminoacyl-tRNA synthetases, especially tyrosyl- (TyrRS) and tryptophanyl-tRNA synthetase (TrpRS). Human TyrRS and TrpRS have an extra domain. Two distinct cytokines, i.e., the core catalytic "mini TyrRS" and the extra C-domain, are generated from human TyrRS by proteolytic cleavage. Moreover, the core catalytic domains of human TyrRS and TrpRS function as angiogenic and angiostatic factors, respectively, whereas the full-length forms are inactive for this function. It is also known that many synthetases change their localization in response to a specific signal and subsequently exhibit alternative functions. Furthermore, some synthetases function as sensors for amino acids by changing their protein interactions in an amino acid-dependent manner. Further studies will be necessary to elucidate regulatory mechanisms of non-canonical functions of aminoacyl-tRNA synthetases in particular, by analyzing the effect of their post-translational modifications.

The Regulatory Role of Key Metabolites in the Control of Cell Signaling

Robust biological systems are able to adapt to internal and environmental perturbations. This is ensured by a thick crosstalk between metabolism and signal transduction pathways, through which cell cycle progression, cell metabolism and growth are coordinated. Although several reports describe the control of cell signaling on metabolism (mainly through transcriptional regulation and post-translational modifications), much fewer information is available on the role of metabolism in the regulation of signal transduction. Protein-metabolite interactions (PMIs) result in the modification of the protein activity due to a conformational change associated with the binding of a small molecule. An increasing amount of evidences highlight the role of metabolites of the central metabolism in the control of the activity of key signaling proteins in different eukaryotic systems. Here we review the known PMIs between primary metabolites and proteins, through which metabolism affects signal transduction pathways controlled by the conserved kinases Snf1/AMPK, Ras/PKA and TORC1. Interestingly, PMIs influence also the mitochondrial retrograde response (RTG) and calcium signaling, clearly demonstrating that the range of this phenomenon is not limited to signaling pathways related to metabolism.

Sequence characterization and molecular modeling of clinically relevant variants of the SARS-CoV-2 main protease

The SARS-CoV-2 main protease (M ^pro ) is essential to viral replication and cleaves highly specific substrate sequences, making it an obvious target for inhibitor design. However, as for any virus, SARS-CoV-2 is subject to constant selection pressure, with new M ^pro mutations arising over time. Identification and structural characterization of M ^pro variants is thus critical for robust inhibitor design. Here we report sequence analysis, structure predictions, and molecular modeling for seventy-nine M ^pro variants, constituting all clinically observed mutations in this protein as of April 29, 2020. Residue substitution is widely distributed, with some tendency toward larger and more hydrophobic residues. Modeling and protein structure network analysis suggest differences in cohesion and active site flexibility, revealing patterns in viral evolution that have relevance for drug discovery.

A suite of polymerase chain reaction-based peptide tagging plasmids for epitope-targeted enzymatic functionalization of yeast proteins

The budding yeast and model eukaryote Saccharomyces cerevisiae has been invaluable for purification and analysis of numerous evolutionarily conserved proteins and multisubunit complexes that cannot be readily reconstituted in Escherichia coli. For many studies, it is desirable to functionalize a particular protein or subunit of a complex with a ligand, fluorophore or other small molecule. Enzyme-catalysed site-specific modification of proteins bearing short peptide tags is a powerful strategy to overcome the limitations associated with traditional nonselective labelling chemistries. Towards this end, we developed a suite of template plasmids for C-terminal tagging with short peptide sequences that can be site-specifically functionalized with high efficiency and selectivity. We have also combined these sequences with the FLAG tag as a handle for purification or immunological detection of the modified protein. We demonstrate the utility of these plasmids by site-specifically labelling the 28-subunit core particle subcomplex of the 26S proteasome with the small molecule fluorophore Cy5. The full set of plasmids has been deposited in the non-profit plasmid repository Addgene (http://www.addgene.org).

PseudoFuN: Deriving functional potentials of pseudogenes from integrative relationships with genes and microRNAs across 32 cancers

Background

Long thought “relics” of evolution, not until recently have pseudogenes been of medical interest regarding regulation in cancer. Often, these regulatory roles are a direct by-product of their close sequence homology to protein-coding genes. Novel pseudogene-gene (PGG) functional associations can be identified through the integration of biomedical data, such as sequence homology, functional pathways, gene expression, pseudogene expression, and microRNA expression. However, not all of the information has been integrated, and almost all previous pseudogene studies relied on 1:1 pseudogene–parent gene relationships without leveraging other homologous genes/pseudogenes.

Results

We produce PGG families that expand beyond the current 1:1 paradigm. First, we construct expansive PGG databases by (i) CUDAlign graphics processing unit (GPU) accelerated local alignment of all pseudogenes to gene families (totaling 1.6 billion individual local alignments and >40,000 GPU hours) and (ii) BLAST-based assignment of pseudogenes to gene families. Second, we create an open-source web application (PseudoFuN [Pseudogene Functional Networks]) to search for integrative functional relationships of sequence homology, microRNA expression, gene expression, pseudogene expression, and gene ontology. We produce four “flavors” of CUDAlign-based databases (>462,000,000 PGG pairwise alignments and 133,770 PGG families) that can be queried and downloaded using PseudoFuN. These databases are consistent with previous 1:1 PGG annotation and also are much more powerful including millions of de novo PGG associations. For example, we find multiple known (e.g., miR-20a-PTEN-PTENP1) and novel (e.g., miR-375-SOX15-PPP4R1L) microRNA-gene-pseudogene associations in prostate cancer. PseudoFuN provides a “one stop shop” for identifying and visualizing thousands of potential regulatory relationships related to pseudogenes in The Cancer Genome Atlas cancers.

Conclusions

Thousands of new PGG associations can be explored in the context of microRNA-gene-pseudogene co-expression and differential expression with a simple-to-use online tool by bioinformaticians and oncologists alike.

Expression fusion immunogen by live attenuated Escherichia coli against enterotoxins infection in mice

Previous epidemiological studies have shown that enterotoxins from enterotoxigenic Escherichia coli (ETEC) appear to be the most important causes of neonatal piglet and porcine post-weaning diarrhoea (PWD). Thus, it is necessary to develop an effective vaccine against ETEC infection. In the present study, the Kil cassette was inserted into the pseudogene yaiT by homologous recombination to create an attenuated E. coli double selection platform O142(yaiT-Kil). After that, PRPL-Kil was replaced with a fusion gene (LTA1-STa13 -STb-LTA2-LTB-STa13 -STb) to establish oral vaccines O142(yaiT::LTA1-STa13 -STb-LTA2-LTB-STa13 -STb) (ER-T). Subsequently, BALB/c mice were orally immunized with ER-T. Results showed that serum IgG and faecal sIgA responded against all ETEC enterotoxins and induced F41 antibody in BALB/c mice by orogastrically inoculation with recombinant E. coli ER-T. Moreover, the determination of cellular immune response demonstrated that the stimulation index (SI) was significantly higher in immunized mice than in control mice, and a clear trend in the helper T-cell (Th) response was Th2-cell (IL-4) exceed Th1-cell (IFN-γ).Our results indicated that recombinant E. coli ER-T provides effective protection against ETEC infection.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.

dreamBase: DNA modification, RNA regulation and protein binding of expressed pseudogenes in human health and disease

Although thousands of pseudogenes have been annotated in the human genome, their transcriptional regulation, expression profiles and functional mechanisms are largely unknown. In this study, we developed dreamBase (http://rna.sysu.edu.cn/dreamBase) to facilitate the investigation of DNA modification, RNA regulation and protein binding of potential expressed pseudogenes from multidimensional high-throughput sequencing data. Based on ∼5500 ChIP-seq and DNase-seq datasets, we identified genome-wide binding profiles of various transcription-associated factors around pseudogene loci. By integrating ∼18 000 RNA-seq data, we analysed the expression profiles of pseudogenes and explored their co-expression patterns with their parent genes in 32 cancers and 31 normal tissues. By combining microRNA binding sites, we demonstrated complex post-transcriptional regulation networks involving 275 microRNAs and 1201 pseudogenes. We generated ceRNA networks to illustrate the crosstalk between pseudogenes and their parent genes through competitive binding of microRNAs. In addition, we studied transcriptome-wide interactions between RNA binding proteins (RBPs) and pseudogenes based on 458 CLIP-seq datasets. In conjunction with epitranscriptome sequencing data, we also mapped 1039 RNA modification sites onto 635 pseudogenes. This database will provide insights into the transcriptional regulation, expression, functions and mechanisms of pseudogenes as well as their roles in biological processes and diseases.

Glycolytic flux in Saccharomyces cerevisiae is dependent on RNA polymerase III and its negative regulator Maf1

Protein biosynthesis is energetically costly, is tightly regulated and is coupled to stress conditions including glucose deprivation. RNA polymerase III (RNAP III)-driven transcription of tDNA genes for production of tRNAs is a key element in efficient protein biosynthesis. Here we present an analysis of the effects of altered RNAP III activity on the Saccharomyces cerevisiae proteome and metabolism under glucose-rich conditions. We show for the first time that RNAP III is tightly coupled to the glycolytic system at the molecular systems level. Decreased RNAP III activity or the absence of the RNAP III negative regulator, Maf1 elicit broad changes in the abundance profiles of enzymes engaged in fundamental metabolism in S. cerevisiae In a mutant compromised in RNAP III activity, there is a repartitioning towards amino acids synthesis de novo at the expense of glycolytic throughput. Conversely, cells lacking Maf1 protein have greater potential for glycolytic flux.

De Novo Gene Evolution of Antifreeze Glycoproteins in Codfishes Revealed by Whole Genome Sequence Data

New genes can arise through duplication of a pre-existing gene or de novo from non-coding DNA, providing raw material for evolution of new functions in response to a changing environment. A prime example is the independent evolution of antifreeze glycoprotein genes (afgps) in the Arctic codfishes and Antarctic notothenioids to prevent freezing. However, the highly repetitive nature of these genes complicates studies of their organization. In notothenioids, afgps evolved from an extant gene, yet the evolutionary origin of afgps in codfishes is unknown. Here, we demonstrate that afgps in codfishes have evolved de novo from non-coding DNA 13-18 Ma, coinciding with the cooling of the Northern Hemisphere. Using whole-genome sequence data from several codfishes and notothenioids, we find higher copy number of afgp in species exposed to more severe freezing suggesting a gene dosage effect. Notably, antifreeze function is lost in one lineage of codfishes analogous to the afgp losses in non-Antarctic notothenioids. This indicates that selection can eliminate the antifreeze function when freezing is no longer imminent. In addition, we show that evolution of afgp-assisting antifreeze potentiating protein genes (afpps) in notothenioids coincides with origin and lineage-specific losses of afgp. The origin of afgps in codfishes is one of the first examples of an essential gene born from non-coding DNA in a non-model species. Our study underlines the power of comparative genomics to uncover past molecular signatures of genome evolution, and further highlights the impact of de novo gene origin in response to a changing selection regime.

Identification of Conflicting Selective Effects on Highly Expressed Genes

Many different selective effects on DNA and proteins influence the frequency of codons and amino acids in coding sequences. Selection is often stronger on highly expressed genes. Hence, by comparing high- and low-expression genes it is possible to distinguish the factors that are selected by evolution. It has been proposed that highly expressed genes should (i) preferentially use codons matching abundant tRNAs (translational efficiency), (ii) preferentially use amino acids with low cost of synthesis, (iii) be under stronger selection to maintain the required amino acid content, and (iv) be selected for translational robustness. These effects act simultaneously and can be contradictory. We develop a model that combines these factors, and use Akaike's Information Criterion for model selection. We consider pairs of paralogues that arose by whole-genome duplication in Saccharmyces cerevisiae. A codon-based model is used that includes asymmetric effects due to selection on highly expressed genes. The largest effect is translational efficiency, which is found to strongly influence synonymous, but not non-synonymous rates. Minimization of the cost of amino acid synthesis is implicated. However, when a more general measure of selection for amino acid usage is used, the cost minimization effect becomes redundant. Small effects that we attribute to selection for translational robustness can be identified as an improvement in the model fit on top of the effects of translational efficiency and amino acid usage.

The Architecture of the Rag GTPase Signaling Network

The evolutionarily conserved target of rapamycin complex 1 (TORC1) couples an array of intra- and extracellular stimuli to cell growth, proliferation and metabolism, and its deregulation is associated with various human pathologies such as immunodeficiency, epilepsy, and cancer. Among the diverse stimuli impinging on TORC1, amino acids represent essential input signals, but how they control TORC1 has long remained a mystery. The recent discovery of the Rag GTPases, which assemble as heterodimeric complexes on vacuolar/lysosomal membranes, as central elements of an amino acid signaling network upstream of TORC1 in yeast, flies, and mammalian cells represented a breakthrough in this field. Here, we review the architecture of the Rag GTPase signaling network with a special focus on structural aspects of the Rag GTPases and their regulators in yeast and highlight both the evolutionary conservation and divergence of the mechanisms that control Rag GTPases.

A Computational Methodology to Overcome the Challenges Associated With the Search for Specific Enzyme Targets to Develop Drugs Against Leishmania major

We present an approach for detecting enzymes that are specific of Leishmania major compared with Homo sapiens and provide targets that may assist research in drug development. This approach is based on traditional techniques of sequence homology comparison by similarity search and Markov modeling; it integrates the characterization of enzymatic functionality, secondary and tertiary protein structures, protein domain architecture, and metabolic environment. From 67 enzymes represented by 42 enzymatic activities classified by AnEnPi (Analogous Enzymes Pipeline) as specific for L major compared with H sapiens, only 40 (23 Enzyme Commission [EC] numbers) could actually be considered as strictly specific of L major and 27 enzymes (19 EC numbers) were disregarded for having ambiguous homologies or analogies with H sapiens. Among the 40 strictly specific enzymes, we identified sterol 24-C-methyltransferase, pyruvate phosphate dikinase, trypanothione synthetase, and RNA-editing ligase as 4 essential enzymes for L major that may serve as targets for drug development.

Sequence statistics of tertiary structural motifs reflect protein stability

The Protein Data Bank (PDB) has been a key resource for learning general rules of sequence-structure relationships in proteins. Quantitative insights have been gained by defining geometric descriptors of structure (e.g., distances, dihedral angles, solvent exposure, etc.) and observing their distributions and sequence preferences. Here we argue that as the PDB continues to grow, it may become unnecessary to reduce structure into a set of elementary descriptors. Instead, it could be possible to deduce quantitative sequence-structure relationships in the context of precisely-defined complex structural motifs by mining the PDB for closely matching backbone geometries. To validate this idea, we turned to the the task of predicting changes in protein stability upon amino-acid substitution—a difficult problem of broad significance. We defined non-contiguous tertiary motifs (TERMs) around a protein site of interest and extracted sequence preferences from ensembles of closely-matching substructures in the PDB to predict mutational stability changes at the site, ΔΔG_m. We demonstrate that these ensemble statistics predict ΔΔG_m on par with state-of-the-art statistical and machine-learning methods on large thermodynamic datasets, and outperform these, along with a leading structure-based modeling approach, when tested in the context of unbiased diverse mutations. Further, we show that the performance of the TERM-based method is directly related to the amount of available relevant structural data, automatically improving with the growing PDB. This enables a means of estimating prediction accuracy. Our results clearly demonstrate that: 1) statistics of non-contiguous structural motifs in the PDB encode fundamental sequence-structure relationships related to protein thermodynamic stability, and 2) the PDB is now large enough that such statistics are already useful in practice, with their accuracy expected to continue increasing as the database grows. These observations suggest new ways of using structural data towards addressing problems of computational structural biology.

Proximity-dependent labeling methods for proteomic profiling in living cells

Characterizing the proteome composition of organelles and subcellular regions of living cells can facilitate the understanding of cellular organization as well as protein interactome networks. Proximity labeling-based methods coupled with mass spectrometry (MS) offer a high-throughput approach for systematic analysis of spatially restricted proteomes. Proximity labeling utilizes enzymes that generate reactive radicals to covalently tag neighboring proteins with biotin. The biotinylated endogenous proteins can then be isolated for further analysis by MS. To analyze protein-protein interactions or identify components that localize to discrete subcellular compartments, spatial expression is achieved by fusing the enzyme to specific proteins or signal peptides that target to particular subcellular regions. Although these technologies have only been introduced recently, they have already provided deep insights into a wide range of biological processes. Here, we describe and compare current methods of proximity labeling as well as their applications. As each method has its own unique features, the goal of this review is to describe how different proximity labeling methods can be used to answer different biological questions. WIREs Dev Biol 2017, 6:e272. doi: 10.1002/wdev.272 For further resources related to this article, please visit the WIREs website.

Filling the Void: Proximity-Based Labeling of Proteins in Living Cells

There are inherent limitations with traditional methods to study protein behavior or to determine the constituency of proteins in discrete subcellular compartments. In response to these limitations, several methods have recently been developed that use proximity-dependent labeling. By fusing proteins to enzymes that generate reactive molecules, most commonly biotin, proximate proteins are covalently labeled to enable their isolation and identification. In this review we describe current methods for proximity-dependent labeling in living cells and discuss their applications and future use in the study of protein behavior.

The emerging role of long non-coding RNAs in endometrial cancer

The human genome is pervasively transcribed and approximately 98% of the genome is non-coding. Long non-coding RNAs (lncRNAs) are a heterogeneous group of RNA transcripts that are >200 nucleotides in length with minimal to no protein-coding potential. Similar to proteins, lncRNAs have important biological functions in both normal cells and disease states including many types of cancer. This review summarizes recent advances in our understanding of lncRNAs in cancer biology and highlights the potential for lncRNA as diagnostic biomarkers and therapeutics. Herein we focus on the poorly understood role of lncRNAs in endometrial cancer, the most common gynecologic malignancy in the developed world.

Chemoselective modifications for the traceless ligation of thioamide-containing peptides and proteins

Thioamides are single-atom substitutions of canonical amide bonds, and have been proven to be versatile and minimally perturbing probes in protein folding studies. Previously, our group showed that thioamides can be incorporated into proteins by native chemical ligation (NCL) with Cys as a ligation handle. In this study, we report the expansion of this strategy into non-Cys ligation sites, utilizing radical initiated desulfurization to "erase" the side chain thiol after ligation. The reaction exhibited high chemoselectivity against thioamides, which can be further enhanced with thioacetamide as a sacrificial scavenger. As a proof-of-concept example, we demonstrated the incorporation of a thioamide probe into a 56 amino acid protein, the B1 domain of Protein G (GB1). Finally, we showed that the method can be extended to β-thiol amino acid analogs and selenocysteine.

H2A/K pseudogene mutation may promote cell proliferation

Little attention has been paid to the histone H2A/K pseudogene. Results from our laboratory showed that 7 of 10 kidney cancer patients carried a mutant H2A/K pseudogene; therefore, we were interested in determining the relationship between mutant H2A/K and cell proliferation. We used shotgun and label-free proteomics methods to study whether mutant H2A/K lncRNAs affected cell proliferation. Quantitative proteomic analysis indicated that the expression of mutant H2A/K lncRNAs resulted in the upregulation of many oncogenes, which promoted cell proliferation. Further interaction analyses revealed that a proliferating cell nuclear antigen (PCNA)-protein interaction network, with PCNA in the center, contributes to cell proliferation in cells expressing the mutant H2A/K lncRNAs. Western blotting confirmed the critical upregulation of PCNA by mutant H2A/K lncRNA expression. Finally, the promotion of cell proliferation by mutant H2A/K lncRNAs (C290T, C228A and A45G) was confirmed using cell proliferation assays. Although we did not determine the exact mechanism by which the oncogenes were upregulated by the mutant H2A/K lncRNAs, we confirmed that the mutant H2A/K lncRNAs promoted cell proliferation by upregulating PCNA and other oncogenes. The hypothesis that cell proliferation is promoted by the mutant H2A/K lncRNAs was supported by the protein expression and cell proliferation assay results. Therefore, mutant H2A/K lncRNAs may be a new factor in renal carcinogenesis.

Exploring G protein-coupled receptor signaling networks using SILAC-based phosphoproteomics

The type 1 parathyroid hormone receptor (PTH1R) is a key regulator of calcium homeostasis and bone turnover. Here, we employed SILAC-based quantitative mass spectrometry and bioinformatic pathways analysis to examine global changes in protein phosphorylation following short-term stimulation of endogenously expressed PTH1R in osteoblastic cells in vitro. Following 5min exposure to the conventional agonist, PTH(1-34), we detected significant changes in the phosphorylation of 224 distinct proteins. Kinase substrate motif enrichment demonstrated that consensus motifs for PKA and CAMK2 were the most heavily upregulated within the phosphoproteome, while consensus motifs for mitogen-activated protein kinases were strongly downregulated. Signaling pathways analysis identified ERK1/2 and AKT as important nodal kinases in the downstream network and revealed strong regulation of small GTPases involved in cytoskeletal rearrangement, cell motility, and focal adhesion complex signaling. Our data illustrate the utility of quantitative mass spectrometry in measuring dynamic changes in protein phosphorylation following GPCR activation.

Comparative genomic and phenomic analysis of Clostridium difficile and Clostridium sordellii, two related pathogens with differing host tissue preference

Background

Clostridium difficile and C. sordellii are two anaerobic, spore forming, gram positive pathogens with a broad host range and the ability to cause lethal infections. Despite strong similarities between the two Clostridial strains, differences in their host tissue preference place C. difficile infections in the gastrointestinal tract and C. sordellii infections in soft tissues.

Results

In this study, to improve our understanding of C. sordellii and C. difficile virulence and pathogenesis, we have performed a comparative genomic and phenomic analysis of the two. The global phenomes of C. difficile and C. sordellii were compared using Biolog Phenotype microarrays. When compared to C. difficile, C. sordellii was found to better utilize more complex sources of carbon and nitrogen, including peptides. Phenotype microarray comparison also revealed that C. sordellii was better able to grow in acidic pH conditions. Using next generation sequencing technology, we determined the draft genome of C. sordellii strain 8483 and performed comparative genome analysis with C. difficile and other Clostridial genomes. Comparative genome analysis revealed the presence of several enzymes, including the urease gene cluster, specific to the C. sordellii genome that confer the ability of expanded peptide utilization and survival in acidic pH.

Conclusions

The identified phenotypes of C. sordellii might be important in causing wound and vaginal infections respectively. Proteins involved in the metabolic differences between C. sordellii and C. difficile should be targets for further studies aimed at understanding C. difficile and C. sordellii infection site specificity and pathogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1663-5) contains supplementary material, which is available to authorized users.

Construction and evaluation of a fluorescence-based live attenuated Escherichia coli delivery system for generating oral vaccine candidate

Enter toxigenic Escherichia coli (ETEC) is a major pathogen of swine industry that can have a substantial impact on morbidity and mortality. Therefore, it is necessary to develop effective vaccines for the prevention of ETEC infection. Live attenuated bacteria delivery system are effective tools for mucosal immunization. The purpose of this study was to construct a novel delivery system that can present the LTR192G-STb fusion protein as oral vaccine candidate. Firstly, the PRPL-mKate2 fluorescent cassette was inserted into the genome (yaiT pseudogene) of an attenuated E. coli by homologous recombination methods to construct the delivery system O142(yaiT::PRPL-mKate2). Secondly, the oral vaccine O142(yaiT:: LT192-STb) (ER-B) was derived for replacing the PRPL-mKate2 by LT192-STb fusion gene, and then it was tested for its feasibility as oral vaccine candidate. Subsequently, BALB/c mice were orogastrically immunized with ER-B. Results showed that mice orally immunized with ER-B produced high levels of specific IgA and IgG antibodies. The induced antibodies demonstrated neutralizing effects to enter toxins LT and STb. In addition, results of cellular immune responses showed that stimulation index values of immunized mice were significantly higher than the control group (P < 0.05) and with a marked shift towards Th 2 immunity. These data indicated that the recombinant E. coli ER-B could be a valuable candidate of future vaccines against ETEC infection.

Comparative analysis of pseudogenes across three phyla

Pseudogenes are degraded fossil copies of genes. Here, we report a comparison of pseudogenes spanning three phyla, leveraging the completed annotations of the human, worm, and fly genomes, which we make available as an online resource. We find that pseudogenes are lineage specific, much more so than protein-coding genes, reflecting the different remodeling processes marking each organism's genome evolution. The majority of human pseudogenes are processed, resulting from a retrotranspositional burst at the dawn of the primate lineage. This burst can be seen in the largely uniform distribution of pseudogenes across the genome, their preservation in areas with low recombination rates, and their preponderance in highly expressed gene families. In contrast, worm and fly pseudogenes tell a story of numerous duplication events. In worm, these duplications have been preserved through selective sweeps, so we see a large number of pseudogenes associated with highly duplicated families such as chemoreceptors. However, in fly, the large effective population size and high deletion rate resulted in a depletion of the pseudogene complement. Despite large variations between these species, we also find notable similarities. Overall, we identify a broad spectrum of biochemical activity for pseudogenes, with the majority in each organism exhibiting varying degrees of partial activity. In particular, we identify a consistent amount of transcription (∼15%) across all species, suggesting a uniform degradation process. Also, we see a uniform decay of pseudogene promoter activity relative to their coding counterparts and identify a number of pseudogenes with conserved upstream sequences and activity, hinting at potential regulatory roles.

Ultradeep human phosphoproteome reveals a distinct regulatory nature of Tyr and Ser/Thr-based signaling

Regulatory protein phosphorylation controls normal and pathophysiological signaling in eukaryotic cells. Despite great advances in mass-spectrometry-based proteomics, the extent, localization, and site-specific stoichiometry of this posttranslational modification (PTM) are unknown. Here, we develop a stringent experimental and computational workflow, capable of mapping more than 50,000 distinct phosphorylated peptides in a single human cancer cell line. We detected more than three-quarters of cellular proteins as phosphoproteins and determined very high stoichiometries in mitosis or growth factor signaling by label-free quantitation. The proportion of phospho-Tyr drastically decreases as coverage of the phosphoproteome increases, whereas Ser/Thr sites saturate only for technical reasons. Tyrosine phosphorylation is maintained at especially low stoichiometric levels in the absence of specific signaling events. Unexpectedly, it is enriched on higher-abundance proteins, and this correlates with the substrate KM values of tyrosine kinases. Our data suggest that P-Tyr should be considered a functionally separate PTM of eukaryotic proteomes.

Comparative genomics suggests that the human pathogenic fungus Pneumocystis jirovecii acquired obligate biotrophy through gene loss

Pneumocystis jirovecii is a fungal parasite that colonizes specifically humans and turns into an opportunistic pathogen in immunodeficient individuals. The fungus is able to reproduce extracellularly in host lungs without eliciting massive cellular death. The molecular mechanisms that govern this process are poorly understood, in part because of the lack of an in vitro culture system for Pneumocystis spp. In this study, we explored the origin and evolution of the putative biotrophy of P. jirovecii through comparative genomics and reconstruction of ancestral gene repertoires. We used the maximum parsimony method and genomes of related fungi of the Taphrinomycotina subphylum. Our results suggest that the last common ancestor of Pneumocystis spp. lost 2,324 genes in relation to the acquisition of obligate biotrophy. These losses may result from neutral drift and affect the biosyntheses of amino acids and thiamine, the assimilation of inorganic nitrogen and sulfur, and the catabolism of purines. In addition, P. jirovecii shows a reduced panel of lytic proteases and has lost the RNA interference machinery, which might contribute to its genome plasticity. Together with other characteristics, that is, a sex life cycle within the host, the absence of massive destruction of host cells, difficult culturing, and the lack of virulence factors, these gene losses constitute a unique combination of characteristics which are hallmarks of both obligate biotrophs and animal parasites. These findings suggest that Pneumocystis spp. should be considered as the first described obligate biotrophs of animals, whose evolution has been marked by gene losses.

A genetic variant in pseudogene E2F3P1 contributes to prognosis of hepatocellular carcinoma

Certain pseudogenes may regulate their protein-coding cousins by competing for miRNAs and play an active biological role in cancer. However, few studies have focused on the association of genetic variations in pseudogenes with cancer prognosis. We selected six potentially functional single nucleotide polymorphisms (SNPs) in cancer-related pseudogenes, and performed a case-only study to assess the association between those SNPs and the prognosis of hepatocellular carcinoma (HCC) in 331 HBV-positive HCC patients without surgical treatment. Log-rank test and Cox proportional hazard models were used for survival analysis. We found that the A allele of rs9909601 in E2F3P1 was significantly associated with a better prognosis compared with the G allele [adjusted hazard ratio (HR)  =  0.69, 95% confidence interval (CI)  =  0.56–0.86, P  =  0.001]. Additionally, this protective effect was more predominant for patients without chemotherapy and transcatheter hepatic arterial chemoembolization (TACE) treatment. Interestingly, we also detected a statistically significant multiplicative interaction between genotypes of rs9909601 and chemotherapy or TACE status on HCC survival (P for multiplicative interaction < 0.001). These findings indicate that rs9909601 in the pseudogene E2F3P1 may be a genetic marker for HCC prognosis in Chinese.