PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools

Localization of NPFxD motif-containing proteins in Aspergillus nidulans

During growth, filamentous fungi produce polarized cells called hyphae. It is generally presumed that polarization of hyphae is dependent upon secretion through the Spitzenkörper, as well as a mechanism called apical recycling, which maintains a balance between the tightly coupled processes of endocytosis and exocytosis. Endocytosis predominates in an annular domain called the sub-apical endocytic collar, which is located in the region of plasma membrane 1-5 μm distal to the Spitzenkörper. It has previously been proposed that one function of the sub-apical endocytic collar is to maintain the apical localization of polarization proteins. These proteins mark areas of polarization at the apices of hyphae. However, as hyphae grow, these proteins are displaced along the membrane and some must then be removed at the sub-apical endocytic collar in order to maintain the hyphoid shape. While endocytosis is fairly well characterized in yeast, comparatively little is known about the process in filamentous fungi. Here, a bioinformatics approach was utilized to identify 39 Aspergillus nidulans proteins that are predicted to be cargo of endocytosis based on the presence of an NPFxD peptide motif. This motif is a necessary endocytic signal sequence first established in Saccharomyces cerevisiae, where it marks proteins for endocytosis through an interaction with the adapter protein Sla1p. It is hypothesized that some proteins that contain this NPFxD peptide sequence in A. nidulans will be potential targets for endocytosis, and therefore will localize either to the endocytic collar or to more proximal polarized regions of the cell, e.g. the apical dome or the Spitzenkörper. To test this, a subset of the motif-containing proteins in A. nidulans was tagged with GFP and the dynamic localization was evaluated. The documented localization patterns support the hypothesis that the motif marks proteins for localization to the polarized cell apex in growing hyphae.

Bioinformatics Analysis Identifies the Estrogen Receptor 1 (ESR1) Gene and hsa-miR-26a-5p as Potential Prognostic Biomarkers in Patients with Intrahepatic Cholangiocarcinoma

Background

Intrahepatic cholangiocarcinoma arises from the epithelial cells of the bile ducts and is associated with poor prognosis. This study aimed to use bioinformatics analysis to identify molecular biomarkers of intrahepatic cholangiocarcinoma and their potential mechanisms.

Material/Methods

MicroRNA (miRNA) and mRNA microarrays from GSE53870 and GSE32879 were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed miRNAs (DEMs) associated with prognosis were identified using limma software and Kaplan-Meier survival analysis. Predictive target genes of the DEMs were identified using miRWalk, miRTarBase, miRDB, and TargetScan databases of miRNA-binding sites and targets. Target genes underwent Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Hub genes were analyzed by constructing the protein-protein interaction (PPI) network using Cytoscape. DEMs validated the hub genes, followed by construction of the miRNA-gene regulatory network.

Results

Twenty-five DEMs were identified. Fifteen DEMs were upregulated, and ten were down-regulated. Kaplan-Meier survival analysis identified seven upregulated DEMs and nine down-regulated DEMs that were associated with the overall survival (OS), and 130 target genes were selected. GO analysis showed that target genes were mainly enriched for metabolism and development processes. KEGG analysis showed that target genes were mainly enriched for cancer processes and some signaling pathways. Fourteen hub genes identified from the PPI network were associated with the regulation of cell proliferation. The overlap between hub genes and DEMs identified the estrogen receptor 1 (ESR1) gene and hsa-miR-26a-5p.

Conclusions

Bioinformatics analysis identified ESR1 and hsa-miR-26a-5p as potential prognostic biomarkers for intrahepatic cholangiocarcinoma.

The FOXA1 transcriptional network coordinates key functions of primary human airway epithelial cells

The differentiated functions of the human airway epithelium are coordinated by a complex network of transcription factors. These include the pioneer factors Forkhead box A1 and A2 (FOXA1 and FOXA2), which are well studied in several tissues, but their role in airway epithelial cells is poorly characterized. Here, we define the cistrome of FOXA1 and FOXA2 in primary human bronchial epithelial (HBE) cells by chromatin immunoprecipitation with deep-sequencing (ChIP-seq). Next, siRNA-mediated depletion of each factor is used to investigate their transcriptome by RNA-seq. We found that, as predicted from their DNA-binding motifs, genome-wide occupancy of the two factors showed substantial overlap; however, their global impact on gene expression differed. FOXA1 is an abundant transcript in HBE cells, while FOXA2 is expressed at low levels, and both these factors likely exhibit autoregulation and cross-regulation. FOXA1 regulated loci are involved in cell adhesion and the maintenance of epithelial cell identity, particularly through repression of genes associated with epithelial to mesenchymal transition (EMT). FOXA1 also directly targets other transcription factors with a known role in the airway epithelium such as SAM-pointed domain-containing Ets-like factor (SPDEF). The intersection of the cistrome and transcriptome for FOXA1 revealed enrichment of genes involved in epithelial development and tissue morphogenesis. Moreover, depletion of FOXA1 was shown to reduce the transepithelial resistance of HBE cells, confirming the role of this factor in maintaining epithelial barrier integrity.

Translational control of one-carbon metabolism underpins ribosomal protein phenotypes in cell division and longevity

A long-standing problem is how cells that lack one of the highly similar ribosomal proteins (RPs) often display distinct phenotypes. Yeast and other organisms live longer when they lack specific ribosomal proteins, especially of the large 60S subunit of the ribosome. However, longevity is neither associated with the generation time of RP deletion mutants nor with bulk inhibition of protein synthesis. Here, we queried actively dividing RP mutants through the cell cycle. Our data link transcriptional, translational, and metabolic changes to phenotypes associated with the loss of paralogous RPs. We uncovered translational control of transcripts encoding enzymes of methionine and serine metabolism, which are part of one-carbon (1C) pathways. Cells lacking Rpl22Ap, which are long-lived, have lower levels of metabolites associated with 1C metabolism. Loss of 1C enzymes increased the longevity of wild type cells. 1C pathways exist in all organisms and targeting the relevant enzymes could represent longevity interventions.

The reactome pathway knowledgebase

The Reactome Knowledgebase (https://reactome.org) provides molecular details of signal transduction, transport, DNA replication, metabolism and other cellular processes as an ordered network of molecular transformations in a single consistent data model, an extended version of a classic metabolic map. Reactome functions both as an archive of biological processes and as a tool for discovering functional relationships in data such as gene expression profiles or somatic mutation catalogs from tumor cells. To extend our ability to annotate human disease processes, we have implemented a new drug class and have used it initially to annotate drugs relevant to cardiovascular disease. Our annotation model depends on external domain experts to identify new areas for annotation and to review new content. New web pages facilitate recruitment of community experts and allow those who have contributed to Reactome to identify their contributions and link them to their ORCID records. To improve visualization of our content, we have implemented a new tool to automatically lay out the components of individual reactions with multiple options for downloading the reaction diagrams and associated data, and a new display of our event hierarchy that will facilitate visual interpretation of pathway analysis results.

Mitochondrial tRNA methylation in Alzheimer's disease and progressive supranuclear palsy

BACKGROUND

Methylation of mitochondrial tRNAs (mt-tRNA) at the 9th position ("p9 site") is known to impact translational efficiency and downstream mitochondrial function; however, direct assessment of mt-RNA methylation is challenging. Recent RNA sequence-based methods have been developed to reliably identify post-transcriptional methylation. Though p9 methylation has been studied in healthy human populations and in the context of cancer, it has not yet been analyzed in neurodegenerative disease, where mitochondrial dysfunction is a prominent and early hallmark of disease progression.

METHODS

Mitochondrial p9 methylation was inferred from multi-allelic calls in RNA-seq data. Gene-based association studies were performed in FUMA. Correlations between nuclear gene expression and p9 methylation were tested using Spearman's rho. Fisher's Exact test was used in PANTHER and IPA to test for overrepresentation and enrichment of biological processes and pathways in the top nuclear genes correlated with p9 methylation.

RESULTS

Variable methylation was observed at 11 p9 sites in post-mortem cerebellar tissue of elderly subjects who were either healthy or diagnosed with Alzheimer's disease (AD), progressive supranuclear palsy (PSP) or pathological aging (PA). Similarities in degree of methylation were observed between AD and PSP. Certain nuclear encoded genes were identified as significantly associated with p9 methylation. Expression of 5300 nuclear encoded genes was significantly correlated with p9 methylation, with AD and PSP subjects exhibiting similar expression profiles. Overrepresentation and enrichment testing using the top transcripts revealed enrichment for a number of molecular processes, terms and pathways including many of which that were mitochondrial-related.

CONCLUSION

With mitochondrial dysfunction being an established hallmark of neurodegenerative disease pathophysiology, this work sheds light on the potential molecular underpinnings of this dysfunction. Here we show overlap in cerebellar pathophysiology between common tauopathies such as Alzheimer's disease and progressive supranuclear palsy. Whether p9 hypermethylation is a cause or consequence of pathology remains an area of focus.

miRDB: an online database for prediction of functional microRNA targets

MicroRNAs (miRNAs) are small noncoding RNAs that act as master regulators in many biological processes. miRNAs function mainly by downregulating the expression of their gene targets. Thus, accurate prediction of miRNA targets is critical for characterization of miRNA functions. To this end, we have developed an online database, miRDB, for miRNA target prediction and functional annotations. Recently, we have performed major updates for miRDB. Specifically, by employing an improved algorithm for miRNA target prediction, we now present updated transcriptome-wide target prediction data in miRDB, including 3.5 million predicted targets regulated by 7000 miRNAs in five species. Further, we have implemented the new prediction algorithm into a web server, allowing custom target prediction with user-provided sequences. Another new database feature is the prediction of cell-specific miRNA targets. miRDB now hosts the expression profiles of over 1000 cell lines and presents target prediction data that are tailored for specific cell models. At last, a new web query interface has been added to miRDB for prediction of miRNA functions by integrative analysis of target prediction and Gene Ontology data. All data in miRDB are freely accessible at http://mirdb.org.

Knockdown of the Ribosomal Protein eL29 in Mammalian Cells Leads to Significant Changes in Gene Expression at the Transcription Level

An imbalance in the synthesis of ribosomal proteins can lead to the disruption of various cellular processes. For mammalian cells, it has been shown that the level of the eukaryote-specific ribosomal protein eL29, also known as the one interacting with heparin/heparan sulfate, substantially affects their growth. Moreover, in animals lacking this protein, a number of anatomical abnormalities have been observed. Here, we applied next-generation RNA sequencing to HEK293 cells transfected with siRNAs specific for the mRNA of eL29 to determine what changes occur in the transcriptome profile with a decrease in the level of the target protein. We showed that an approximately 2.5-fold decrease in the content of eL29 leads to statistically significant changes in the expression of more than a thousand genes at the transcription level, without a noticeable effect on cell viability, rRNA level, and global translation. The set of eL29-dependent genes included both up-regulated and down-regulated ones, among which there are those previously identified as targets for proteins implicated in oncogenesis. Thus, our findings demonstrate that an insufficiency of eL29 in mammalian cells causes a significant reorganization of gene expression, thereby highlighting the relationship between the cellular balance of eL29 and the activities of certain genes.

Myocardial micro-biopsy procedure for molecular characterization with increased precision and reduced trauma

Endomyocardial biopsy is a valuable tool in cardiac diagnostics but is limited by low diagnostic yield and significant complication risks. Meanwhile, recent developments in transcriptomic and proteomic technologies promise a wealth of biological data from minimal tissue samples. To take advantage of the minimal tissue amount needed for molecular analyses, we have developed a sub-millimeter endovascular biopsy device, considerably smaller than current clinical equipment, and devised a low-input RNA-sequencing protocol for analyzing small tissue samples. In in vivo evaluation in swine, 81% of biopsy attempts (n = 157) were successful. High quality RNA-sequencing data was generated from 91% of the sequenced cardiac micro-biopsy samples (n = 32). Gene expression signatures of samples taken with the novel device were comparable with a conventional device. No major complications were detected either during procedures or during 7 days' follow-up, despite acquiring a relatively large number of biopsies (median 30) in each animal. In conclusion, the novel device coupled with RNA-sequencing provides a feasible method to obtain molecular data from the myocardium. The method is less traumatic and has a higher flexibility compared to conventional methods, enabling safer and more targeted sampling from different parts of the myocardium.

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

Haplotype-Based Genome-Wide Association Study and Identification of Candidate Genes Associated with Carcass Traits in Hanwoo Cattle

Hanwoo, is the most popular native beef cattle in South Korea. Due to its extensive popularity, research is ongoing to enhance its carcass quality and marbling traits. In this study we conducted a haplotype-based genome-wide association study (GWAS) by constructing haplotype blocks by three methods: number of single nucleotide polymorphisms (SNPs) in a haplotype block (nsnp), length of genomic region in kb (Len) and linkage disequilibrium (LD). Significant haplotype blocks and genes associated with them were identified for carcass traits such as BFT (back fat thickness), EMA (eye Muscle area), CWT (carcass weight) and MS (marbling score). Gene-set enrichment analysis and functional annotation of genes in the significantly-associated loci revealed candidate genes, including PLCB1 and PLCB4 present on BTA13, coding for phospholipases, which might be important candidates for increasing fat deposition due to their role in lipid metabolism and adipogenesis. CEL (carboxyl ester lipase), a bile-salt activated lipase, responsible for lipid catabolic process was also identified within the significantly-associated haplotype block on BTA11. The results were validated in a different Hanwoo population. The genes and pathways identified in this study may serve as good candidates for improving carcass traits in Hanwoo cattle.

Mitochondrial DNA alterations may influence the cisplatin responsiveness of oral squamous cell carcinoma

Cisplatin is the first-line chemotherapeutic agent for the treatment of oral squamous cell carcinoma (OSCC). However, the intrinsic or acquired resistance against cisplatin remains a major obstacle to treatment efficacy in OSCC. Recently, mitochondrial DNA (mtDNA) alterations have been reported in a variety of cancers. However, the role of mtDNA alterations in OSCC has not been comprehensively studied. In this study, we evaluated the correlation between mtDNA alterations (mtDNA content, point mutations, large-scale deletions, and methylation status) and cisplatin sensitivity using two OSCC cell lines, namely SAS and H103, and stem cell-like tumour spheres derived from SAS. By microarray analysis, we found that the tumour spheres profited from aberrant lipid and glucose metabolism and became resistant to cisplatin. By qPCR analysis, we found that the cells with less mtDNA were less responsive to cisplatin (H103 and the tumour spheres). Based on the findings, we theorised that the metabolic changes in the tumour spheres probably resulted in mtDNA depletion, as the cells suppressed mitochondrial respiration and switched to an alternative mode of energy production, i.e. glycolysis. Then, to ascertain the origin of the variation in mtDNA content, we used MinION, a nanopore sequencer, to sequence the mitochondrial genomes of H103, SAS, and the tumour spheres. We found that the lower cisplatin sensitivity of H103 could have been caused by a constellation of genetic and epigenetic changes in its mitochondrial genome. Future work may look into how changes in mtDNA translate into an impact on cell function and therefore cisplatin response.

Mitochondrial proteomics profile points oxidative phosphorylation as main target for beauvericin and enniatin B mixture

Beauvericin (BEA) and enniatin B (EN B) are non-legislated Fusarium mycotoxins usually found in cereal and cereal-based products all around the world. By the proteomic analysis of mitochondria enriched extracts from Jurkat cells exposed for 24 h to three concentrations of BEA:EN B (0.01-0.1-0.5 μM), a number of 1821 proteins (202 mitochondrial) were identified and relatively quantified. 340 proteins (59 mitochondrial) were statistically significant altered in our samples (Anova p-value ≤ 0.05 and fold change (FC) ≥1.5). The protein mitochondrial translational release factor 1 like (MTRF1L) was the most abundant protein in the three mycotoxin exposures studied. The mycotoxins mixture exposure induced concentration dependent changes at mitochondrial proteins levels that mainly involve inner and outer membrane complexes, Electron Transport Chain (ETC) and ribosomes. These results showed alteration of proteins levels related to oxidative phosphorylation, metabolic and neurodegenerative diseases related pathways.

Proteomic Profiling of Small Extracellular Vesicles Secreted by Human Pancreatic Cancer Cells Implicated in Cellular Transformation

Extracellular vesicles secreted from tumor cells are functional vehicles capable of contributing to intercellular communication and metastasis. A growing number of studies have focused on elucidating the role that tumor-derived extracellular vesicles play in spreading pancreatic cancer to other organs, due to the highly metastatic nature of the disease. We recently showed that small extracellular vesicles secreted from pancreatic cancer cells could initiate malignant transformation of healthy cells. Here, we analyzed the protein cargo contained within these vesicles using mass spectrometry-based proteomics to better understand their makeup and biological characteristics. Three different human pancreatic cancer cell lines were compared to normal pancreatic epithelial cells revealing distinct differences in protein cargo between cancer and normal vesicles. Vesicles from cancer cells contain an enrichment of proteins that function in the endosomal compartment of cells responsible for vesicle formation and secretion in addition to proteins that have been shown to contribute to oncogenic cell transformation. Conversely, vesicles from normal pancreatic cells were shown to be enriched for immune response proteins. Collectively, results contribute to what we know about the cargo contained within or excluded from cancer cell-derived extracellular vesicles, supporting their role in biological processes including metastasis and cancer progression.

Targeting Cellular Metabolism in Acute Myeloid Leukemia and The Role of Patient Heterogeneity

Acute myeloid leukemia (AML) is an aggressive blood cancer resulting in accumulation of immature, dysfunctional blood cells in the bone marrow. Changes in cell metabolism are features of many cancers, including AML and this may be exploited as a therapeutic target. In this study we investigated the in vitro antileukemic effects of seven metabolic inhibitors that target different metabolic pathways. The metabolic inhibitors were tested on AML cells derived from 81 patients using proliferation and viability assays; we also compared global gene expression and proteomic profiles for various patient subsets. Metformin, 2DG, 6AN, BPTES and ST1326 had strong antiproliferative and proapoptotic effects for most patients, whereas lonidamine and AZD3965 had an effect only for a minority. Antiproliferative effects on AML cells were additive when combined with the chemotherapeutic agent AraC. Using unsupervised hierarchical clustering, we identified a strong antiproliferative effect on AML cells after treatment with metabolic inhibitors for a subset of 29 patients. Gene expression and proteomic studies suggested that this subset was characterized by altered metabolic and transcriptional regulation. In addition, the Bcl-2 inhibitor venetoclax, in combination with 2DG or 6AN, increased the antiproliferative effects of these metabolic inhibitors on AML cells. Therapeutic targeting of cellular metabolism may have potential in AML, but the optimal strategy will likely differ between patients.

15q11.2 deletion is enriched in patients with total anomalous pulmonary venous connection

INTRODUCTION

CNV is a vital pathogenic factor of congenital heart disease (CHD). However, few CNVs have been reported for total anomalous pulmonary venous connection (TAPVC), which is a rare form of CHD. Using case-control study, we identified 15q11.2 deletion associated with TAPVC. We then used a TAPVC trio as model to reveal possible molecular basis of 15q11.2 microdeletion.

METHODS

CNVplex and Chromosomal Microarray were used to identify and validate CNVs in samples from 231 TAPVC cases and 200 healthy controls from Shanghai Children's Medical Center. In vitro cardiomyocyte differentiation of induced pluripotent stem cells from peripheral blood mononuclear cells for a TAPVC trio with paternal inherited 15q11.2 deletion was performed to characterise the effect of the deletion on cardiomyocyte differentiation and gene expression.

RESULTS

The 15q11.2 microdeletion was significantly enriched in patients with TAPVC compared with healthy control (13/231 in patients vs 0/200 in controls, p=5.872×10-2, Bonferroni adjusted) using Fisher's exact test. Induced pluripotent stem cells from the proband could not differentiate into normal cardiomyocyte. Transcriptomic analysis identified a number of differentially expressed genes in the 15q11.2 deletion carriers of the family. TAPVC disease-causing genes such as PITX2, NKX2-5 and ANKRD1 showed significantly higher expression in the proband compared with her healthy mother. Knockdown of TUBGCP5 could lead to abnormal cardiomyocyte differentiation.

CONCLUSION

We discovered that the 15q11.2 deletion is significantly associated with TAPVC. Gene expression profile that might arise from 15q11.2 deletion for a TAPVC family was characterised using cell experiments.

Applying knowledge-driven mechanistic inference to toxicogenomics

When considering toxic chemicals in the environment, a mechanistic, causal explanation of toxicity may be preferred over a statistical or machine learning-based prediction by itself. Elucidating a mechanism of toxicity is, however, a costly and time-consuming process that requires the participation of specialists from a variety of fields, often relying on animal models. We present an innovative mechanistic inference framework (MechSpy), which can be used as a hypothesis generation aid to narrow the scope of mechanistic toxicology analysis. MechSpy generates hypotheses of the most likely mechanisms of toxicity, by combining a semantically-interconnected knowledge representation of human biology, toxicology and biochemistry with gene expression time series on human tissue. Using vector representations of biological entities, MechSpy seeks enrichment in a manually curated list of high-level mechanisms of toxicity, represented as biochemically- and causally-linked ontology concepts. Besides predicting the canonical mechanism of toxicity for many well-studied compounds, we experimentally validated some of our predictions for other chemicals without an established mechanism of toxicity. This mechanistic inference framework is an advantageous tool for predictive toxicology, and the first of its kind to produce a mechanistic explanation for each prediction. MechSpy can be modified to include additional mechanisms of toxicity, and is generalizable to other types of mechanisms of human biology.

Therapeutic potential of lipids obtained from γ-irradiated PBMCs in dendritic cell-mediated skin inflammation

BACKGROUND

Since numerous pathological conditions are evoked by unwanted dendritic cell (DC) activity, therapeutic agents modulating DC functions are of great medical interest. In regenerative medicine, cellular secretomes have gained increasing attention and valuable immunomodulatory properties have been attributed to the secretome of γ-irradiated peripheral blood mononuclear cells (PBMCs). Potential effects of the PBMC secretome (PBMCsec) on key DC functions have not been elucidated so far.

METHODS

We used a hapten-mediated murine model of contact hypersensitivity (CH) to study the effects of PBMCsec on DCs in vivo. Effects of PBMCsec on human DCs were investigated in monocyte-derived DCs (MoDC) and ex vivo skin cultures. DCs were phenotypically characterised by transcriptomics analyses and flow cytometry. DC function was evaluated by cytokine secretion, antigen uptake, PBMC proliferation and T-cell priming.

FINDINGS

PBMCsec significantly alleviated tissue inflammation and cellular infiltration in hapten-sensitized mice. We found that PBMCsec abrogated differentiation of MoDCs, indicated by lower expression of classical DC markers CD1a, CD11c and MHC class II molecules. Furthermore, PBMCsec reduced DC maturation, antigen uptake, lipopolysaccharides-induced cytokine secretion, and DC-mediated immune cell proliferation. Moreover, MoDCs differentiated with PBMCsec displayed diminished ability to prime naïve CD4+T-cells into TH1 and TH2 cells. Furthermore, PBMCsec modulated the phenotype of DCs present in the skin in situ. Mechanistically, we identified lipids as the main biomolecule accountable for the observed immunomodulatory effects.

INTERPRETATION

Together, our data describe DC-modulatory actions of lipids secreted by stressed PBMCs and suggest PBMCsec as a therapeutic option for treatment of DC-mediated inflammatory skin conditions.

FUNDING

This research project was supported by the Austrian Research Promotion Agency (Vienna, Austria; grant "APOSEC" 862068; 2015-2019) and the Vienna Business Agency (Vienna, Austria; grant "APOSEC to clinic" 2343727).

Post-translational Modification Crosstalk and Hotspots in Sirtuin Interactors Implicated in Cardiovascular Diseases

Sirtuins are protein deacetylases that play a protective role in cardiovascular diseases (CVDs), as well as many other diseases. Absence of sirtuins can lead to hyperacetylation of both nuclear and mitochondrial proteins leading to metabolic dysregulation. The protein post-translational modifications (PTMs) are known to crosstalk among each other to bring about complex phenotypic outcomes. Various PTM types such as acetylation, ubiquitination, and phosphorylation, and so on, drive transcriptional regulation and metabolism, but such crosstalks are poorly understood. We integrated protein–protein interactions (PPI) and PTMs from several databases to integrate information on 1,251 sirtuin-interacting proteins, of which 544 are associated with cardiac diseases. Based on the ∼100,000 PTM sites obtained for sirtuin interactors, we observed that the frequency of PTM sites (83 per protein), as well as PTM types (five per protein), is higher than the global average for human proteome. We found that ∼60–70% PTM sites fall into ordered regions. Approximately 83% of the sirtuin interactors contained at least one competitive crosstalk (in situ) site, with half of the sites occurring in CVD-associated proteins. A large proportion of identified crosstalk sites were observed for acetylation and ubiquitination competition. We identified 614 proteins containing PTM hotspots (≥5 PTM sites) and 133 proteins containing crosstalk hotspots (≥3 crosstalk sites). We observed that a large proportion of disease-associated sequence variants were found in PTM motifs of CVD proteins. We identified seven proteins (TP53, LMNA, MAPT, ATP2A2, NCL, APEX1, and HIST1H3A) containing disease-associated variants in PTM and crosstalk hotspots. This is the first comprehensive bioinformatics analysis on sirtuin interactors with respect to PTMs and their crosstalks. This study forms a platform for generating interesting hypotheses that can be tested for a deeper mechanistic understanding gained or derived from big-data analytics.

Identification of Novel Targets of RBM5 in the Healthy and Injured Brain

The tumor suppressor RNA-binding motif 5 (RBM5) regulates the expression levels and cassette exon-definition (i.e. splicing) of a select set of mRNAs in a tissue-specific manner. Most RBM5-regulated targets were identified in oncological investigations and frequently involve genes which mediate apoptotic cell death. Little is known about the role of RBM5 in the brain. Also, it is unclear if a brain injury may be required to detect RBM5 mediated effects on pro-apoptotic genes due to their low expression levels in the healthy adult CNS at baseline. Conditional/floxed (brain-specific) gene deleter mice were generated to elucidate CNS-specific RBM5 mRNA targets. Male/female mice were subjected to a severe controlled cortical impact (CCI) traumatic brain injury (TBI) in order to increase the background expression of pro-death mRNAs and facilitate testing of the hypothesis that RBM5 inhibition decreases post-injury upregulation of caspases/FAS in the CNS. As expected, a CCI increased caspases/FAS mRNA in the injured cortex. RBM5 KO did not affect their levels or splicing. Surprisingly, KO increased the mRNA levels of novel targets including casein kinase 2 alpha prime interacting protein (Csnka2ip/CKT2) - a gene not thought to be expressed in the brain, contrary to findings here. Twenty-two unique splicing events were also detected in KOs including increased block-inclusion of cassette exons 20-22 in regulating synaptic membrane exocytosis 2 (Rims2). In conclusion, here we used genome-wide transcriptomic analysis on healthy and injured RBM5 KO mouse brain tissue to elucidate the first known gene targets of this enigmatic RBP in this CNS.

Reduction of protein phosphatase 2A (PP2A) complexity reveals cellular functions and dephosphorylation motifs of the PP2A/B'δ holoenzyme

Protein phosphatase 2A (PP2A) is a large enzyme family responsible for most cellular Ser/Thr dephosphorylation events. PP2A substrate specificity, localization, and regulation by second messengers rely on more than a dozen regulatory subunits (including B/R2, B'/R5, and B″/R3), which form the PP2A heterotrimeric holoenzyme by associating with a dimer comprising scaffolding (A) and catalytic (C) subunits. Because of partial redundancy and high endogenous expression of PP2A holoenzymes, traditional approaches of overexpressing, knocking down, or knocking out PP2A regulatory subunits have yielded only limited insights into their biological roles and substrates. To this end, here we sought to reduce the complexity of cellular PP2A holoenzymes. We used tetracycline-inducible expression of pairs of scaffolding and regulatory subunits with complementary charge-reversal substitutions in their interaction interfaces. For each of the three regulatory subunit families, we engineered A/B charge-swap variants that could bind to one another, but not to endogenous A and B subunits. Because endogenous Aα was targeted by a co-induced shRNA, endogenous B subunits were rapidly degraded, resulting in expression of predominantly a single PP2A heterotrimer composed of the A/B charge-swap pair and the endogenous catalytic subunit. Using B'δ/PPP2R5D, we show that PP2A complexity reduction, but not PP2A overexpression, reveals a role of this holoenzyme in suppression of extracellular signal-regulated kinase signaling and protein kinase A substrate dephosphorylation. When combined with global phosphoproteomics, the PP2A/B'δ reduction approach identified consensus dephosphorylation motifs in its substrates and suggested that residues surrounding the phosphorylation site play roles in PP2A substrate specificity.

Pilot study and bioinformatics analysis of differentially expressed genes in adipose tissues of rats with excess dietary intake

Excessive adipose tissue accumulation is an increasing health problem worldwide. The present study aimed to determine differentially expressed genes (DEGs) that are associated with the excessive accumulation of adipose tissues by PCR arrays in an excess dietary intake animal model. For this purpose, male Sprague Dawley rats were randomly assigned to 2 groups: Control (given an ordinary diet) and experimental (given twice the amount of the ordinary diet). After 2 months of feeding, the abdominal cavities of the rats from each group were opened, then subcutaneous and visceral adipose tissues were removed. The adipose tissues collected were then used for total RNA extraction and then reverse transcribed to cDNA, which was then used as a template to identify the DEGs of 84 transcripts for rat obesity by RT2 Profiler PCR Arrays. The results showed significant downregulation of bombesin‑like receptor 3 (BRS3) and uncoupling protein 1 (UCP1) in visceral adipose tissues of experimental rats compared with those of the control rats, and differential gene expression analysis showed an association with fat cell differentiation and regulation of triglyceride sequestration, as well as fatty acid binding. The gene expression patterns observed in the present study, which may be associated with peroxisome proliferator‑activated receptor‑γ (PPARG) on excessive visceral adipose tissue accumulation, may be useful in identifying a group of surrogate biomarkers for the early diet‑induced accumulation of visceral adipose tissue detection in humans. The biomarkers can also be the specific targets for drug development to reduce excessive visceral adipose tissue accumulation in the body and its associated diseases.

Decoy receptor 2 mediation of the senescent phenotype of tubular cells by interacting with peroxiredoxin 1 presents a novel mechanism of renal fibrosis in diabetic nephropathy

Premature senescence of renal tubular epithelial cell (RTEC), which is involved in kidney fibrosis, is a key event in the progression of diabetic nephropathy. However, the underlying mechanism remains unclear. Here we investigated the role and mechanism of decoy receptor 2 (DcR2) in kidney fibrosis and the senescent phenotype of RTEC. DcR2 was specifically expressed in senescent RTEC and associated with kidney fibrosis in patients with diabetic nephropathy and mice with streptozotocin-induced with diabetic nephropathy. Knockdown of DcR2 decreased the expression of α-smooth muscle actin, collagen I, fibronectin and serum creatinine levels in streptozotocin-induced mice. DcR2 knockdown also inhibited the expression of senescent markers p16, p21, senescence-associated beta-galactosidase and senescence-associated heterochromatic foci and promoted the secretion of a senescence-associated secretory phenotype including IL-6, TGF-β1, and matrix metalloproteinase 2 in vitro and in vivo. However, DcR2 overexpression showed the opposite effects. Quantitative proteomics and validation studies revealed that DcR2 interacted with peroxiredoxin 1 (PRDX1), which regulated the cell cycle and senescence. Knockdown of PRDX1 upregulated p16 and cyclin D1 while downregulating cyclin-dependent kinase 6 expression in vitro, resulting in RTEC senescence. Furthermore, PRDX1 knockdown promoted DcR2-induced p16, cyclin D1, IL-6, and TGF-β1 expression, whereas PRDX1 overexpression led to the opposite results. Subsequently, DcR2 regulated PRDX1 phosphorylation, which could be inhibited by the specific tyrosine kinase inhibitor genistein. Thus, DcR2 mediated the senescent phenotype of RTEC and kidney fibrosis by interacting with PRDX1. Hence, DcR2 may act as a potential therapeutic target for the amelioration of diabetic nephropathy progression.

AP-1 Contributes to Chromatin Accessibility to Promote Sarcomere Disassembly and Cardiomyocyte Protrusion During Zebrafish Heart Regeneration

RATIONALE

The adult human heart is an organ with low regenerative potential. Heart failure following acute myocardial infarction is a leading cause of death due to the inability of cardiomyocytes to proliferate and replenish lost cardiac muscle. While the zebrafish has emerged as a powerful model to study endogenous cardiac regeneration, the molecular mechanisms by which cardiomyocytes respond to damage by disassembling sarcomeres, proliferating, and repopulating the injured area remain unclear. Furthermore, we are far from understanding the regulation of the chromatin landscape and epigenetic barriers that must be overcome for cardiac regeneration to occur.

OBJECTIVE

To identify transcription factor regulators of the chromatin landscape, which promote cardiomyocyte regeneration in zebrafish, and investigate their function.

METHODS AND RESULTS

Using the Assay for Transposase-Accessible Chromatin coupled to high-throughput sequencing (ATAC-Seq), we first find that the regenerating cardiomyocyte chromatin accessibility landscape undergoes extensive changes following cryoinjury, and that activator protein-1 (AP-1) binding sites are the most highly enriched motifs in regions that gain accessibility during cardiac regeneration. Furthermore, using bioinformatic and gene expression analyses, we find that the AP-1 response in regenerating adult zebrafish cardiomyocytes is largely different from the response in adult mammalian cardiomyocytes. Using a cardiomyocyte-specific dominant negative approach, we show that blocking AP-1 function leads to defects in cardiomyocyte proliferation as well as decreased chromatin accessibility at the fbxl22 and ilk loci, which regulate sarcomere disassembly and cardiomyocyte protrusion into the injured area, respectively. We further show that overexpression of the AP-1 family members Junb and Fosl1 can promote changes in mammalian cardiomyocyte behavior in vitro.

CONCLUSIONS

AP-1 transcription factors play an essential role in the cardiomyocyte response to injury by regulating chromatin accessibility changes, thereby allowing the activation of gene expression programs that promote cardiomyocyte dedifferentiation, proliferation, and protrusion into the injured area.

Differential Response of the Chicken Trachea to Chronic Infection with Virulent Mycoplasma gallisepticum Strain Ap3AS and Vaxsafe MG (Strain ts-304): a Transcriptional Profile

Mycoplasma gallisepticum is the primary etiological agent of chronic respiratory disease in chickens. Live attenuated vaccines are most commonly used in the field to control the disease, but current vaccines have some limitations. Vaxsafe MG (strain ts-304) is a new vaccine candidate that is efficacious at a lower dose than the current commercial vaccine strain ts-11, from which it is derived. In this study, the transcriptional profiles of the trachea of unvaccinated chickens and chickens vaccinated with strain ts-304 were compared 2 weeks after challenge with M. gallisepticum strain Ap3AS during the chronic stage of infection. After challenge, genes, gene ontologies, pathways, and protein classes involved in inflammation, cytokine production and signaling, and cell proliferation were upregulated, while those involved in formation and motor movement of cilia, formation of intercellular junctional complexes, and formation of the cytoskeleton were downregulated in the unvaccinated birds compared to the vaccinated birds, reflecting immune dysregulation and the pathological changes induced in the trachea by infection with M. gallisepticum Vaccination appears to protect the structural and functional integrity of the tracheal mucosa 2 weeks after infection with M. gallisepticum.

Protein Synthesis by Day 16 Bovine Conceptuses during the Time of Maternal Recognition of Pregnancy

Interferon Tau (IFNT), the conceptus-derived pregnancy recognition signal in cattle, significantly modifies the transcriptome of the endometrium. However, the endometrium also responds to IFNT-independent conceptus-derived products. The aim of this study was to determine what proteins are produced by the bovine conceptus that may facilitate the pregnancy recognition process in cattle. We analysed by mass spectrometry the proteins present in conceptus-conditioned media (CCM) after 6 h culture of Day 16 bovine conceptuses (n = 8) in SILAC media (arginine- and lysine-depleted media supplemented with heavy isotopes) and the protein content of extracellular vesicles (EVs) isolated from uterine luminal fluid (ULF) of Day 16 pregnant (n = 7) and cyclic (n = 6) cross-bred heifers on day 16. In total, 11,122 proteins were identified in the CCM. Of these, 5.95% (662) had peptides with heavy labelled amino acids, i.e., de novo synthesised by the conceptuses. None of these proteins were detected in the EVs isolated from ULF. Pregnancy-associated glycoprotein 11, Trophoblast Kunitz domain protein 1 and DExD-Box Helicase 39A were de novo produced and present in the CCM from all conceptuses and in previously published CCM data following 6 and 24 h. A total of 463 proteins were present in the CCM from all the conceptuses in the present study, and after 6 and 24 h culture in a previous study, while expression of their transcripts was not detected in endometrium indicating that they are likely conceptus-derived. Of the proteins present in the EVs, 67 were uniquely identified in ULF from pregnant heifers; 35 of these had been previously reported in CCM from Day 16 conceptuses. This study has narrowed a set of conceptus-derived proteins that may be involved in EV-mediated IFNT-independent embryo–maternal communication during pregnancy recognition in cattle.

Transcriptome analysis of responses in Brachypodium distachyon overexpressing the BdbZIP26 transcription factor

BACKGROUND

Biotic and abiotic stresses are the major cause of reduced growth, persistence, and yield in agriculture. Over the past decade, RNA-Sequencing and the use of transgenics with altered expression of stress related genes have been utilized to gain a better understanding of the molecular mechanisms leading to salt tolerance in a variety of species. Identification of transcription factors that, when overexpressed in plants, improve multiple stress tolerance may be valuable for crop improvement, but sometimes overexpression leads to deleterious effects during normal plant growth.

RESULTS

Brachypodium constitutively expressing the BdbZIP26:GFP gene showed reduced stature compared to wild type plants (WT). RNA-Seq analysis comparing WT and bZIP26 transgenic plants revealed 7772 differentially expressed genes (DEGs). Of these DEGs, 987 of the DEGs were differentially expressed in all three transgenic lines. Many of these DEGs are similar to those often observed in response to abiotic and biotic stress, including signaling proteins such as kinases/phosphatases, calcium/calmodulin related proteins, oxidases/reductases, hormone production and signaling, transcription factors, as well as disease responsive proteins. Interestingly, there were many DEGs associated with protein turnover including ubiquitin-related proteins, F-Box and U-box related proteins, membrane proteins, and ribosomal synthesis proteins. Transgenic and control plants were exposed to salinity stress. Many of the DEGs between the WT and transgenic lines under control conditions were also found to be differentially expressed in WT in response to salinity stress. This suggests that the over-expression of the transcription factor is placing the plant in a state of stress, which may contribute to the plants diminished stature.

CONCLUSION

The constitutive expression of BdbZIP26:GFP had an overall negative effect on plant growth and resulted in stunted plants compared to WT plants under control conditions, and a similar response to WT plants under salt stress conditions. The results of gene expression analysis suggest that the transgenic plants are in a constant state of stress, and that they are trying to allocate resources to survive.

How tryptophan levels in plant-based aquafeeds affect fish physiology, metabolism and proteome

Fish meal replacement by plant-protein sources is a priority in aquaculture feeds. Within this framework, dietary supplementation with essential amino acids (EAA), as tryptophan (TRP), is strategic to ensure that the individual nutritional needs are met, besides promoting enhanced immunological status. The purpose of this study was to examine the beneficial effects of TRP incorporation in plant-protein source diets on fish growth performance and nutritional status. We tested diets with 20% lower (LTRP) and 27% higher (HTRP) of the putative requirements of TRP for seabream (Sparus aurata) and assessed its impact on fish physiology and liver metabolism and proteome. After 12 weeks, growth performance, body proximate, hepatic composition and liver metabolic profiling were similar between diets. Nevertheless, liver proteome analysis indicated a higher accumulation of proteins involved in acute-phase responses, typically triggered by infection, inflammation or trauma, in fish fed with HTRP diet as compared with those fed with LTRP. The overall results obtained suggest a potential beneficial effect of TRP supplementation in terms of immune stimulation, without compromising growth or feed intake. Moreover, proteomics and metabolic profiling demonstrate to be valuable tools in this endeavour. SIGNIFICANCE: Nutritional needs are hard to assess in aquaculture fisheries, and many times controversial depending on the methodology employed. The estimated amino acid requirements depend on both fish species and stage development, making it extremely hard to standardise. On the other hand, the substitution of fish-based to plant-based protein sources diets towards a sustainable aquaculture, may imbalance these requirements, being necessary further studies to assess the impact on fish growth and development. Finally, the incorporation of crystalized amino acids such as TRP into diets aims global better performance both at fish health/immune condition and growth development. This work focused on the potential beneficial effects of TRP supplementation into diets with a plant-based protein source, addressing the effects on the liver metabolism and proteome, and on growth performance of Gilthead seabream juveniles, a species with special relevance and economical importance in the Mediterranean region. The present study by employing proteomics together with metabolic profiling shows that TRP supplementation at the tested doses, does not compromise growth performance, and seems to stimulate the immune system. Our findings can contribute to the development of new feed formulations for Gilthead seabream species, therefore, reinforcing the resilience and competitiveness of the on-growing aquaculture industry and impact directly the sustainability of living resources with the decrease of the fisheries needs to fulfil the human search for quality proteins consume.

Loss of CHD1 Promotes Heterogeneous Mechanisms of Resistance to AR-Targeted Therapy via Chromatin Dysregulation

Metastatic prostate cancer is characterized by recurrent genomic copy number alterations that are presumed to contribute to resistance to hormone therapy. We identified CHD1 loss as a cause of antiandrogen resistance in an in vivo small hairpin RNA (shRNA) screen of 730 genes deleted in prostate cancer. ATAC-seq and RNA-seq analyses showed that CHD1 loss resulted in global changes in open and closed chromatin with associated transcriptomic changes. Integrative analysis of this data, together with CRISPR-based functional screening, identified four transcription factors (NR3C1, POU3F2, NR2F1, and TBX2) that contribute to antiandrogen resistance, with associated activation of non-luminal lineage programs. Thus, CHD1 loss results in chromatin dysregulation, thereby establishing a state of transcriptional plasticity that enables the emergence of antiandrogen resistance through heterogeneous mechanisms.

The secretome of skin cancer cells activates the mTOR/MYC pathway in healthy keratinocytes and induces tumorigenic properties

Cutaneous squamous cell carcinoma (cSCC) is the most prominent tumor of non-melanoma skin cancers and the most aggressive tumor among keratinocyte carcinoma of the skin, showing a high potential for local invasion and metastasis. The cSCC incidences increased dramatically in recent years and the disease occurs more commonly than any other malignancy. The secretome of cancer cells is currently the focus of many studies in order to identify new marker proteins for different types of cancer and to investigate its influence on the tumor microenvironment. In our study we evaluated whether the secretome of cSCC cells has an impact on keratinocytes, the surrounding tissue cells of cSCC. Therefore, we analyzed and compared the secretome of human A431 cancer cells and of HaCaT keratinocytes by mass spectrometry. In a second experiment, keratinocytes were exposed to the secretome of A431 cells and vice versa and the transcriptome was analyzed by next-generation sequencing. HaCaT cells incubated with A431 conditioned medium revealed a significantly activated mammalian target of rapamycin pathway with a concomitant increase in proliferation and migration. In conclusion, our data demonstrate the impact of the secretome of cancer cells on the transcription machinery of the cells surrounding the tumor, leading to a tumorigenic cell fate.

Quantitative SWATH-Based Proteomic Profiling for Identification of Mechanism-Driven Diagnostic Biomarkers Conferring in the Progression of Metastatic Prostate Cancer

Prostate cancer (PCa), the most frequently diagnosed malignancy in men is associated with significant mortality and morbidity. Therefore, demand exists for the identification of potential biomarkers for patient stratification according to prognostic risks and the mechanisms involved in cancer development and progression to avoid over/under treatment of patients and prevent relapse. Quantitative proteomic mass spectrometry profiling and gene enrichment analysis of TGF-β induced-EMT in human Prostate androgen-dependent (LNCaP) and androgen-independent (PC-3) adenocarcinoma cell lines was performed to investigate proteomics involved in Prostate carcinogenesis and their effect onto the survival of PCa patients. Amongst 1,795 proteins, which were analyzed, 474 proteins were significantly deregulated. These proteins contributed to apoptosis, gluconeogenesis, transcriptional regulation, RNA splicing, cell cycle, and MAPK cascade and hence indicating the crucial roles of these proteins in PCa initiation and progression. We have identified a panel of six proteins viz., GOT1, HNRNPA2B1, MAPK1, PAK2, UBE2N, and YWHAB, which contribute to cancer development, and the transition of PCa from androgen dependent to independent stages. The prognostic values of identified proteins were evaluated using UALCAN, GEPIA, and HPA datasets. The results demonstrate the utility of SWATH-LC-MS/MS for understanding the proteomics involved in EMT transition of PCa and identification of clinically relevant proteomic biomarkers.

Transcriptome Profile of Yeast Strain Used for Biological Wine Aging Revealed Dynamic Changes of Gene Expression in Course of Flor Development

Flor strains of Saccharomyces cerevisiae are principal microbial agents responsible for biological wine aging used for production of sherry-like wines. The flor yeast velum formed on the surface of fortified fermented must is a major adaptive and technological characteristic of flor yeasts that helps them to withstanding stressful winemaking conditions and ensures specific biochemical and sensory oxidative alterations typical for sherry wines. We have applied RNAseq technology for transcriptome analysis of an industrial flor yeast strain at different steps of velum development over 71 days under experimental winemaking conditions. Velum growth and maturation was accompanied by accumulation of aldehydes and acetales. We have identified 1490 differentially expressed genes including 816 genes upregulated and 674 downregulated more than 2-fold at mature biofilm stage as compared to the early biofilm. Distinct expression patterns of genes involved in carbon and nitrogen metabolism, respiration, cell cycle, DNA repair, cell adhesion, response to various stresses were observed. Many genes involved in response to different stresses, oxidative carbon metabolism, high affinity transport of sugars, glycerol utilization, sulfur metabolism, protein quality control and recycling, cell wall biogenesis, apoptosis were induced at the mature biofilm stage. Strong upregulation was observed for FLO11 flocculin while expression of other flocculins remained unaltered or moderately downregulated. Downregulated genes included those for proteins involved in glycolysis, transportation of ions, metals, aminoacids, sugars, indicating repression of some major transport and metabolic process at the mature biofilm stage. Presented results are important for in-depth understanding of cell response elicited by velum formation and sherry wine manufacturing conditions, and for the comprehension of relevant regulatory mechanisms. Such knowledge may help to better understand the molecular mechanisms that flor yeasts use to adapt to winemaking environments, establish the functions of previously uncharacterized genes, improve the technology of sherry- wine production, and find target genes for strain improvement.

Proteome dynamics during homeostatic scaling in cultured neurons

Protein turnover, the net result of protein synthesis and degradation, enables cells to remodel their proteomes in response to internal and external cues. Previously, we analyzed protein turnover rates in cultured brain cells under basal neuronal activity and found that protein turnover is influenced by subcellular localization, protein function, complex association, cell type of origin, and by the cellular environment (Dörrbaum et al., 2018). Here, we advanced our experimental approach to quantify changes in protein synthesis and degradation, as well as the resulting changes in protein turnover or abundance in rat primary hippocampal cultures during homeostatic scaling. Our data demonstrate that a large fraction of the neuronal proteome shows changes in protein synthesis and/or degradation during homeostatic up- and down-scaling. More than half of the quantified synaptic proteins were regulated, including pre- as well as postsynaptic proteins with diverse molecular functions.

A single-cell atlas of adult Drosophila ovary identifies transcriptional programs and somatic cell lineage regulating oogenesis

Oogenesis is a complex developmental process that involves spatiotemporally regulated coordination between the germline and supporting, somatic cell populations. This process has been modeled extensively using the Drosophila ovary. Although different ovarian cell types have been identified through traditional means, the large-scale expression profiles underlying each cell type remain unknown. Using single-cell RNA sequencing technology, we have built a transcriptomic data set for the adult Drosophila ovary and connected tissues. Using this data set, we identified the transcriptional trajectory of the entire follicle-cell population over the course of their development from stem cells to the oogenesis-to-ovulation transition. We further identify expression patterns during essential developmental events that take place in somatic and germline cell types such as differentiation, cell-cycle switching, migration, symmetry breaking, nurse-cell engulfment, egg-shell formation, and corpus luteum signaling. Extensive experimental validation of unique expression patterns in both ovarian and nearby, nonovarian cells also led to the identification of many new cell type-and stage-specific markers. The inclusion of several nearby tissue types in this data set also led to our identification of functional convergence in expression between distantly related cell types such as the immune-related genes that were similarly expressed in immune cells (hemocytes) and ovarian somatic cells (stretched cells) during their brief phagocytic role in nurse-cell engulfment. Taken together, these findings provide new insight into the temporal regulation of genes in a cell-type specific manner during oogenesis and begin to reveal the relatedness in expression between cell and tissues types.

No Dataset Left Behind: Mechanistic Insights into Thyroid Receptor Signaling Through Transcriptomic Consensome Meta-Analysis

Background: Discovery-scale omics datasets relevant to thyroid receptors (TRs) and their physiological and synthetic bioactive small-molecule ligands allow for genome-wide interrogation of TR-regulated genes. These datasets have considerable collective value as a reference resource to allow researchers to routinely generate hypotheses addressing the mechanisms underlying the cell biology and physiology of TR signaling in normal and disease states. Methods: Here, we searched the Gene Expression Omnibus database to identify a population of publicly archived transcriptomic datasets involving genetic or pharmacological manipulation of either TR isoform in a mouse tissue or cell line. After initial quality control, samples were organized into contrasts (experiments), and transcript differential expression values and associated measures of significance were generated and committed to a consensome (for consensus omics) meta-analysis pipeline. To gain insight into tissue-selective functions of TRs, we generated liver- and central nervous system (CNS)-specific consensomes and identified evidence for genes that were selectively responsive to TR signaling in each organ. Results: The TR transcriptomic consensome ranks genes based on the frequency of their significant differential expression over the entire group of experiments. The TR consensome assigns elevated rankings both to known TR-regulated genes and to genes previously uncharacterized as TR-regulated, which shed mechanistic light on known cellular and physiological roles of TR signaling in different organs. We identify evidence for unreported genomic targets of TR signaling for which it exhibits strikingly distinct regulatory preferences in the liver and CNS. Moreover, the intersection of the TR consensome with consensomes for other cellular receptors sheds light on transcripts potentially mediating crosstalk between TRs and these other signaling paradigms. Conclusions: The mouse TR datasets and consensomes are freely available in the Signaling Pathways Project website for hypothesis generation, data validation, and modeling of novel mechanisms of TR regulation of gene expression. Our results demonstrate the insights into the mechanistic basis of thyroid hormone action that can arise from an ongoing commitment on the part of the research community to the deposition of discovery-scale datasets.

Widespread Exon Junction Complex Footprints in the RNA Degradome Mark mRNA Degradation before Steady State Translation

Exon junction complexes (EJCs) are deposited on mRNAs during splicing and displaced by ribosomes during the pioneer round of translation. Nonsense-mediated mRNA decay (NMD) degrades EJC-bound mRNA, but the lack of suitable methodology has prevented the identification of other degradation pathways. Here, we show that the RNA degradomes of Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), worm (Caenorhabditis elegans), and human (Homo sapiens) cells exhibit an enrichment of 5' monophosphate (5'P) ends of degradation intermediates that map to the canonical EJC region. Inhibition of 5' to 3' exoribonuclease activity and overexpression of an EJC disassembly factor in Arabidopsis reduced the accumulation of these 5'P ends, supporting the notion that they are in vivo EJC footprints. Hundreds of Arabidopsis NMD targets possess evident EJC footprints, validating their degradation during the pioneer round of translation. In addition to premature termination codons, plant microRNAs can also direct the degradation of EJC-bound mRNAs. However, the production of EJC footprints from NMD but not microRNA targets requires the NMD factor SUPPRESSOR WITH MORPHOLOGICAL EFFECT ON GENITALIA PROTEIN7. Together, our results demonstrating in vivo EJC footprinting in Arabidopsis unravel the composition of the RNA degradome and provide a new avenue for studying NMD and other mechanisms targeting EJC-bound mRNAs for degradation before steady state translation.

Divergent northern and southern populations and demographic history of the pearl oyster in the western Pacific revealed with genomic SNPs

In the open ocean without terrain boundaries, marine invertebrates with pelagic larvae can migrate long distances using ocean currents, suggesting reduced genetic diversification. Contrary to this assumption, however, genetic differentiation is often observed in marine invertebrates. In the present study, we sought to explain how population structure is established in the western Pacific Ocean, where the strong Kuroshio Current maintains high levels of gene flow from south to north, presumably promoting genetic homogeneity. We determined the population structure of the pearl oyster, Pinctada fucata, in the Indo-Pacific Ocean using genome-wide genotyping data from multiple sampling localities. Cluster analysis showed that the western Pacific population is distinct from that of the Indian Ocean, and that it is divided into northern (Japanese mainland) and southern (Nansei Islands, China, and Cambodia) populations. Genetic differentiation of P. fucata can be explained by geographic barriers in the Indian Ocean and a local lagoon, and by environmental gradients of sea surface temperature (SST) and oxygen concentration in the western Pacific. A genome scan showed evidence of adaptive evolution in genomic loci, possibly associated with changes in environmental factors, including SST and oxygen concentration. Furthermore, Bayesian simulation demonstrated that the past population expansion and division are congruent with ocean warming after the last glacial period. It is highly likely that the environmental gradient forms a genetic barrier that diversifies P. fucata populations in the western Pacific. This hypothesis helps to explain genetic differentiation and possible speciation of marine invertebrates.

Insights into the intracellular localization, protein associations and artemisinin resistance properties of Plasmodium falciparum K13

The emergence of artemisinin (ART) resistance in Plasmodium falciparum intra-erythrocytic parasites has led to increasing treatment failure rates with first-line ART-based combination therapies in Southeast Asia. Decreased parasite susceptibility is caused by K13 mutations, which are associated clinically with delayed parasite clearance in patients and in vitro with an enhanced ability of ring-stage parasites to survive brief exposure to the active ART metabolite dihydroartemisinin. Herein, we describe a panel of K13-specific monoclonal antibodies and gene-edited parasite lines co-expressing epitope-tagged versions of K13 in trans. By applying an analytical quantitative imaging pipeline, we localize K13 to the parasite endoplasmic reticulum, Rab-positive vesicles, and sites adjacent to cytostomes. These latter structures form at the parasite plasma membrane and traffic hemoglobin to the digestive vacuole wherein artemisinin-activating heme moieties are released. We also provide evidence of K13 partially localizing near the parasite mitochondria upon treatment with dihydroartemisinin. Immunoprecipitation data generated with K13-specific monoclonal antibodies identify multiple putative K13-associated proteins, including endoplasmic reticulum-resident molecules, mitochondrial proteins, and Rab GTPases, in both K13 mutant and wild-type isogenic lines. We also find that mutant K13-mediated resistance is reversed upon co-expression of wild-type or mutant K13. These data help define the biological properties of K13 and its role in mediating P. falciparum resistance to ART treatment.

Shared genetic architecture and casual relationship between leptin levels and type 2 diabetes: large-scale cross-trait meta-analysis and Mendelian randomization analysis

OBJECTIVE

We aimed to estimate genetic correlation, identify shared loci and test causality between leptin levels and type 2 diabetes (T2D).

RESEARCH DESIGN AND METHODS

Our study consists of three parts. First, we calculated the genetic correlation of leptin levels and T2D or glycemic traits by using linkage disequilibrium score regression analysis. Second, we conducted a large-scale genome-wide cross-trait meta-analysis using cross-phenotype association to identify shared loci between trait pairs that showed significant genetic correlations in the first part. In the end, we carried out a bidirectional MR analysis to find out whether there is a causal relationship between leptin levels and T2D or glycemic traits.

RESULTS

We found positive genetic correlations between leptin levels and T2D (R_g=0.3165, p=0.0227), fasting insulin (FI) (R_g=0.517, p=0.0076), homeostasis model assessment-insulin resistance (HOMA-IR) (R_g=0.4785, p=0.0196), as well as surrogate estimates of β-cell function (HOMA-β) (R_g=0.4456, p=0.0214). We identified 12 shared loci between leptin levels and T2D, 1 locus between leptin levels and FI, 1 locus between leptin levels and HOMA-IR, and 1 locus between leptin levels and HOMA-β. We newly identified eight loci that did not achieve genome-wide significance in trait-specific genome-wide association studies. These shared genes were enriched in pancreas, thyroid gland, skeletal muscle, placenta, liver and cerebral cortex. In addition, we found that 1-SD increase in HOMA-IR was causally associated with a 0.329 ng/mL increase in leptin levels (β=0.329, p=0.001).

CONCLUSIONS

Our results have shown the shared genetic architecture between leptin levels and T2D and found causality of HOMA-IR on leptin levels, shedding light on the molecular mechanisms underlying the association between leptin levels and T2D.

Arabidopsis EDR1 Protein Kinase Regulates the Association of EDS1 and PAD4 to Inhibit Cell Death

ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) and PHYTOALEXIN DEFICIENT4 (PAD4) are sequence-related lipase-like proteins that function as a complex to regulate defense responses in Arabidopsis by both salicylic acid-dependent and independent pathways. Here, we describe a gain-of-function mutation in PAD4 (S135F) that enhances resistance and cell death in response to infection by the powdery mildew pathogen Golovinomyces cichoracearum. The mutant PAD4 protein accumulates to wild-type levels in Arabidopsis cells, thus these phenotypes are unlikely to be due to PAD4 over accumulation. The phenotypes are similar to loss-of-function mutations in the protein kinase EDR1 (Enhanced Disease Resistance1), and previous work has shown that loss of PAD4 or EDS1 suppresses edr1-mediated phenotypes, placing these proteins downstream of EDR1. Here, we show that EDR1 directly associates with EDS1 and PAD4 and inhibits their interaction in yeast and plant cells. We propose a model whereby EDR1 negatively regulates defense responses by interfering with the heteromeric association of EDS1 and PAD4. Our data indicate that the S135F mutation likely alters an EDS1-independent function of PAD4, potentially shedding light on a yet-unknown PAD4 signaling function.

Integrative -omics and HLA-ligandomics analysis to identify novel drug targets for ccRCC immunotherapy

Background

Clear cell renal cell carcinoma (ccRCC) is the dominant subtype of renal cancer. With currently available therapies, cure of advanced and metastatic ccRCC is achieved only in rare cases. Here, we developed a workflow integrating different -omics technologies to identify ccRCC-specific HLA-presented peptides as potential drug targets for ccRCC immunotherapy.

Methods

We analyzed HLA-presented peptides by MS-based ligandomics of 55 ccRCC tumors (cohort 1), paired non-tumor renal tissues, and 158 benign tissues from other organs. Pathways enriched in ccRCC compared to its cell type of origin were identified by transcriptome and gene set enrichment analyses in 51 tumor tissues of the same cohort. To retrieve a list of candidate targets with involvement in ccRCC pathogenesis, ccRCC-specific pathway genes were intersected with the source genes of tumor-exclusive peptides. The candidates were validated in an independent cohort from The Cancer Genome Atlas (TCGA KIRC, n = 452). DNA methylation (TCGA KIRC, n = 273), somatic mutations (TCGA KIRC, n = 392), and gene ontology (GO) and correlations with tumor metabolites (cohort 1, n = 30) and immune-oncological markers (cohort 1, n = 37) were analyzed to characterize regulatory and functional involvements. CD8⁺ T cell priming assays were used to identify immunogenic peptides. The candidate gene EGLN3 was functionally investigated in cell culture.

Results

A total of 34,226 HLA class I- and 19,325 class II-presented peptides were identified in ccRCC tissue, of which 443 class I and 203 class II peptides were ccRCC-specific and presented in ≥ 3 tumors. One hundred eighty-five of the 499 corresponding source genes were involved in pathways activated by ccRCC tumors. After validation in the independent cohort from TCGA, 113 final candidate genes remained. Candidates were involved in extracellular matrix organization, hypoxic signaling, immune processes, and others. Nine of the 12 peptides assessed by immunogenicity analysis were able to activate naïve CD8⁺ T cells, including peptides derived from EGLN3. Functional analysis of EGLN3 revealed possible tumor-promoting functions.

Conclusions

Integration of HLA ligandomics, transcriptomics, genetic, and epigenetic data leads to the identification of novel functionally relevant therapeutic targets for ccRCC immunotherapy. Validation of the identified targets is recommended to expand the treatment landscape of ccRCC.

Demographic assessment of the Dalmatian dog - effective population size, linkage disequilibrium and inbreeding coefficients

Background

The calculation of demographic measures is a useful tool for evaluating the genomic architecture of dog breeds and enables ranking dog breeds in terms of genetic diversity. To achieve this for the German Dalmatian dog population, 307 purebred animals of this breed were genotyped on the Illumina Canine high density BeadChip. The analysis of pedigree-based inbreeding was performed based on a pedigree with 25,761 dogs including the genotyped dogs.

Results

The effective population size derived from squared correlation coefficients between SNP alleles (r²) was 69. The maximum value of r² was 0.56, resulting in a 50% decay value of 0.28 at a marker distance of 37.5 kb. The effective population size calculated from pedigree data using individual increase in inbreeding over equivalent generations was 116. The pedigree inbreeding coefficient was 0.026. The genomic inbreeding coefficient based on the length of runs of homozygosity (ROH) was calculated for seven length categories of ROHs, and ranged from 0.08 to 0.28. The fixation coefficients F_{IS_PED} and F_{IS_GENO} were at 0.017 and 0.004. PANTHER statistical overrepresentation analysis of genes located in consensus ROHs revealed highly underrepresented biological processes in 50% of the investigated dogs. One of those is the 0.28 fold enriched “immune response”, which might be associated to the high prevalence of allergic dermatitis in the breed. Candidate genes for congenital sensorineural deafness (CCSD, a highly prevalent disease in the Dalmatian) were discovered in consensus ROHs.

Conclusions

The fast decay of r² and the moderate inbreeding coefficients indicate that the German Dalmatian dog population is rather diverse. Pedigree- and genomic-based inbreeding measures were highly correlated and therefore prove good reliability for the given population. Analyses of consensus ROHs with genes coding for deafness and other breed-defining traits, such as hyperuricosuria, indicate that those ROH became fixed in the Dalmatian population about 500 years ago. In case of the Dalmatian dog, a ROH of 40 SNPs length is enough to investigate signatures of selection (e.g. the ROH with the fixed hyperuricosuria mutation) as far back as the breed formation point approximately 500 years ago.

Protein Extraction, Enrichment and MALDI MS and MS/MS Analysis from Bitter Orange Leaves (Citrus aurantium)

Citrus aurantium is a widespread tree in the Mediterranean area, and it is mainly used as rootstock for other citrus. In the present study, a vacuum infiltration centrifugation procedure, followed by solid phase extraction matrix-assisted laser desorption ionization tandem mass spectrometry (SPE MALDI MS/MS) analysis, was adopted to isolate proteins from leaves. The results of mass spectrometry (MS) profiling, combined with the top-down proteomics approach, allowed the identification of 78 proteins. The bioinformatic databases TargetP, SignalP, ChloroP, WallProtDB, and mGOASVM-Loc were used to predict the subcellular localization of the identified proteins. Among 78 identified proteins, 20 were targeted as secretory pathway proteins and 36 were predicted to be in cellular compartments including cytoplasm, nucleus, and cell membrane. The largest subcellular fraction was the secretory pathway, accounting for 25% of total proteins. Gene Ontology (GO) of Citrus sinensis was used to simplify the functional annotation of the proteins that were identified in the leaves. The Kyoto Encyclopedia of Genes and Genomes (KEGG) showed the enrichment of metabolic pathways including glutathione metabolism and biosynthesis of secondary metabolites, suggesting that the response to a range of environmental factors is the key processes in citrus leaves. Finally, the Lipase GDSL domain-containing protein GDSL esterase/lipase, which is involved in plant development and defense response, was for the first time identified and characterized in Citrus aurantium.

Identification of marker proteins of muscular dystrophy in the urine proteome from the mdx-4cv model of dystrophinopathy

Since the protein constituents of urine present a dynamic proteome that can reflect a variety of disease-related alterations in the body, the mass spectrometric survey of proteome-wide changes in urine promises new insights into pathogenic mechanisms. Urine can be investigated in a completely non-invasive way and provides valuable biomedical information on body-wide changes. In this report, we have focused on the urine proteome in X-linked muscular dystrophy using the established mdx-4cv mouse model of dystrophinopathy. In order to avoid potential artefacts due to the manipulation of the biofluid proteome prior to mass spectrometry, crude urine specimens were analyzed without the prior usage of centrifugation steps or concentration procedures. Comparative proteomics revealed 21 increased and 8 decreased proteins out of 870 identified urinary proteoforms using 50 μl of biofluid per investigated sample, i.e. 14 wild type versus 14 mdx-4cv specimens. Promising marker proteins that were almost exclusively found in mdx-4cv urine included nidogen, parvalbumin and titin. Interestingly, the mass spectrometric identification of urine-associated titin revealed a wide spread of peptides over the sequence of this giant muscle protein. The newly established urinomic signature of dystrophinopathy might be helpful for the design of non-invasive assays to improve diagnosis, prognosis, therapy-monitoring and evaluation of potential harmful side effects of novel treatments in the field of muscular dystrophy research.

Temporal patterns of macrophage- and neutrophil-related markers are associated with clinical outcome in heart failure patients

AIMS

Evidence on the association of macrophage- and neutrophil-related blood biomarkers with clinical outcome in heart failure patients is limited, and, with the exception of C-reactive protein, no data exist on their temporal evolution. We aimed to investigate whether temporal patterns of these biomarkers are related to clinical outcome in patients with stable chronic heart failure (CHF).

METHODS AND RESULTS

In 263 patients with CHF, we performed serial plasma measurements of scavenger receptor cysteine-rich type 1 protein M130 (CD163), tartrate-resistant acid phosphatase type 5 (TRAP), granulins (GRN), spondin-1 (SPON1), peptidoglycan recognition protein 1 (PGLYRP1), and tissue factor pathway inhibitor (TFPI). The Cardiovascular Panel III (Olink Proteomics AB, Uppsala, Sweden) was used. During 2.2 years of follow-up, we collected 1984 samples before the occurrence of the composite primary endpoint (PE) or censoring. For efficiency, we selected 567 samples for the measurements (all baseline samples, the last two samples preceding the PE, and the last sample before censoring in event-free patients). The relationship between repeatedly measured biomarker levels and the PE was evaluated by joint models. Mean (±standard deviation) age was 67 ± 13 years; 189 (72%) were men; left ventricular ejection fraction (%) was 32 ± 11. During follow-up, 70 (27%) patients experienced the PE. Serially measured biomarkers predicted the PE in a multivariable model adjusted for baseline clinical characteristics [hazard ratio (95% confidence interval) per 1-standard deviation change in biomarker]: CD163 [2.07(1.47-2.98), P < 0.001], TRAP [0.62 (0.43-0.90), P = 0.009], GRN [2.46 (1.64-3.84), P < 0.001], SPON1 [3.94 (2.50-6.50), P < 0.001], and PGLYRP1 [1.62 (1.14-2.31), P = 0.006].

CONCLUSIONS

Changes in plasma levels of CD163, TRAP, GRN, SPON1, and PGLYRP1 precede adverse cardiovascular events in patients with CHF.

Arginine-Enriched Mixed-Charge Domains Provide Cohesion for Nuclear Speckle Condensation

Low-complexity protein domains promote the formation of various biomolecular condensates. However, in many cases, the precise sequence features governing condensate formation and identity remain unclear. Here, we investigate the role of intrinsically disordered mixed-charge domains (MCDs) in nuclear speckle condensation. Proteins composed exclusively of arginine-aspartic acid dipeptide repeats undergo length-dependent condensation and speckle incorporation. Substituting arginine with lysine in synthetic and natural speckle-associated MCDs abolishes these activities, identifying a key role for multivalent contacts through arginine's guanidinium ion. MCDs can synergize with a speckle-associated RNA recognition motif to promote speckle specificity and residence. MCD behavior is tunable through net-charge: increasing negative charge abolishes condensation and speckle incorporation. Contrastingly, increasing positive charge through arginine leads to enhanced condensation, speckle enlargement, decreased splicing factor mobility, and defective mRNA export. Together, these results identify key sequence determinants of MCD-promoted speckle condensation and link the dynamic material properties of speckles with function in mRNA processing.

Clinical and molecular characterization of virus-positive and virus-negative Merkel cell carcinoma

BACKGROUND

Merkel cell carcinoma (MCC) is a highly aggressive neuroendocrine carcinoma of the skin caused by either the integration of Merkel cell polyomavirus (MCPyV) and expression of viral T antigens or by ultraviolet-induced damage to the tumor genome from excessive sunlight exposure. An increasing number of deep sequencing studies of MCC have identified significant differences between the number and types of point mutations, copy number alterations, and structural variants between virus-positive and virus-negative tumors. However, it has been challenging to reliably distinguish between virus positive and UV damaged MCC.

METHODS

In this study, we assembled a cohort of 71 MCC patients and performed deep sequencing with OncoPanel, a clinically implemented, next-generation sequencing assay targeting over 400 cancer-associated genes. To improve the accuracy and sensitivity for virus detection compared to traditional PCR and IHC methods, we developed a hybrid capture baitset against the entire MCPyV genome and software to detect integration sites and structure.

RESULTS

Sequencing from this approach revealed distinct integration junctions in the tumor genome and generated assemblies that strongly support a model of microhomology-initiated hybrid, virus-host, circular DNA intermediate that promotes focal amplification of host and viral DNA. Using the clear delineation between virus-positive and virus-negative tumors from this method, we identified recurrent somatic alterations common across MCC and alterations specific to each class of tumor, associated with differences in overall survival. Finally, comparing the molecular and clinical data from these patients revealed a surprising association of immunosuppression with virus-negative MCC and significantly shortened overall survival.

CONCLUSIONS

These results demonstrate the value of high-confidence virus detection for identifying molecular mechanisms of UV and viral oncogenesis in MCC. Furthermore, integrating these data with clinical data revealed features that could impact patient outcome and improve our understanding of MCC risk factors.

Quantitative global studies reveal differential translational control by start codon context across the fungal kingdom

Eukaryotic protein synthesis generally initiates at a start codon defined by an AUG and its surrounding Kozak sequence context, but the quantitative importance of this context in different species is unclear. We tested this concept in two pathogenic Cryptococcus yeast species by genome-wide mapping of translation and of mRNA 5' and 3' ends. We observed thousands of AUG-initiated upstream open reading frames (uORFs) that are a major contributor to translation repression. uORF use depends on the Kozak sequence context of its start codon, and uORFs with strong contexts promote nonsense-mediated mRNA decay. Transcript leaders in Cryptococcus and other fungi are substantially longer and more AUG-dense than in Saccharomyces. Numerous Cryptococcus mRNAs encode predicted dual-localized proteins, including many aminoacyl-tRNA synthetases, in which a leaky AUG start codon is followed by a strong Kozak context in-frame AUG, separated by mitochondrial-targeting sequence. Analysis of other fungal species shows that such dual-localization is also predicted to be common in the ascomycete mould, Neurospora crassa. Kozak-controlled regulation is correlated with insertions in translational initiation factors in fidelity-determining regions that contact the initiator tRNA. Thus, start codon context is a signal that quantitatively programs both the expression and the structures of proteins in diverse fungi.

Global profile changes in transcripts induced with a phosphate analogue: implications for gene regulation

Previous work has shown that thiophosphate, a phosphate analogue, leads to a global shift in the distribution of cellular proteins in a variety of organisms. Thiophosphate, when added to culture media, gets incorporated into the nucleic acids of cells resulting in nuclease-resistant phosphorothioate linkages. Using Escherichia coli, as a model organism, it was found that the global changes in protein expression induced with thiophosphate could be accounted for by significant changes in the absolute transcription levels of more than 1500 genes detected via RNA seq analysis. In fact, 58% of transcripts detected in RNA seq studies using total RNA were increased an average of 44 × fold while the remaining 42% were decreased an average of 20 × fold in thiophosphate-treated cells. Furthermore, microarray results showed no correlation between the transcriptional changes observed and the known stability of the corresponding mRNAs measured. Overall, the total amount of non-ribosomal RNA accumulated in TP-treated cells was increased relative to rRNA ~ 4 × fold (1.5-6 ×). The results further indicated that metabolic changes may play a role in inducing the transcriptional profiles observed with thiophosphate. Indeed, pathway analysis of transcripts showed an increase in routes for phosphoribosyl pyrophosphate (PRPP) synthesis and related derivatives, presumably due to a reduction in RNA turnover. These results raise the possibility that the energy savings with reduced RNA turnover could lead to an increased energy charge in the cell that induces transcriptional changes leading to an increase in biosynthetic processes.

Frameshifting preserves key physicochemical properties of proteins

Frameshifts in protein coding sequences are widely perceived as resulting in either nonfunctional or even deleterious protein products. Indeed, frameshifts typically lead to markedly altered protein sequences and premature stop codons. By analyzing complete proteomes from all three domains of life, we demonstrate that, in contrast, several key physicochemical properties of protein sequences exhibit significant robustness against +1 and -1 frameshifts. In particular, we show that hydrophobicity profiles of many protein sequences remain largely invariant upon frameshifting. For example, over 2,900 human proteins exhibit a Pearson's correlation coefficient R between the hydrophobicity profiles of the original and the +1-frameshifted variants greater than 0.7, despite an average sequence identity between the two of only 6.5% in this group. We observe a similar effect for protein sequence profiles of affinity for certain nucleobases as well as protein sequence profiles of intrinsic disorder. Finally, analysis of significance and optimality demonstrates that frameshift stability is embedded in the structure of the universal genetic code and may have contributed to shaping it. Our results suggest that frameshifting may be a powerful evolutionary mechanism for creating new proteins with vastly different sequences, yet similar physicochemical properties to the proteins from which they originate.