Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited.

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues.

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Efficient protein production by yeast requires global tuning of metabolism.

The biotech industry relies on cell factories for production of pharmaceutical proteins, of which several are among the top-selling medicines. There is, therefore, considerable interest in improving the efficiency of protein production by cell factories. Protein secretion involves numerous intracellular processes with many underlying mechanisms still remaining unclear. Here, we use RNA-seq to study the genome-wide transcriptional response to protein secretion in mutant yeast strains. We find that many cellular processes have to be attuned to support efficient protein secretion. In particular, altered energy metabolism resulting in reduced respiration and increased fermentation, as well as balancing of amino-acid biosynthesis and reduced thiamine biosynthesis seem to be particularly important. We confirm our findings by inverse engineering and physiological characterization and show that by tuning metabolism cells are able to efficiently secrete recombinant proteins. Our findings provide increased understanding of which cellular regulations and pathways are associated with efficient protein secretion.

Improving the genotyping resolution of Cryptosporidium hominis subtype IbA10G2 using one step PCR-based amplicon sequencing.

Cryptosporidium hominis gp60 subtype IbA10G2 is a common cause of cryptosporidiosis. This subtype is responsible for many waterborne outbreaks as well as sporadic cases and is considered virulent and highly important in the epidemiology of cryptosporidiosis. Due to low heterogeneity within the genome of C. hominis it has been difficult to identify epidemiological markers with higher resolution than gp60. However, new markers are required in order to improve outbreak investigations and studies of the transmission dynamics of this clinically important subtype. Based on the whole genome sequences of 17 C. hominis isolates, we have identified several differential loci and developed a new sequence based typing panel with higher resolution than gp60. An amplicon sequencing method was also developed which is based on a one-step PCR which can be sequenced using a Next Generation Sequencing (NGS) platform. Such a system provides a rapid and high-throughput workflow. A panel of nine loci with 10 single nucleotide variants (SNV) was selected and evaluated using clinical IbA10G2 isolates from sporadic, cluster and outbreak associated cases. The specimens were separated into 10 different genetic profiles named sequence types (STs). All isolates within an outbreak or cluster belonged to the same ST, including several samples from the two large waterborne outbreaks which occurred in Sweden between 2010 and 2011 indicating that these outbreaks might be linked. The results demonstrate the methods suitability for improved genotyping of C. hominis IbA10G2.

Inflammatory infiltrates in parathyroid tumors.

CONTEXT: Inflammatory infiltrates are sometimes present in solid tumors and may be coupled to clinical behavior or etiology. Infectious viruses contribute to tumorigenesis in a significant fraction of human neoplasias. OBJECTIVE: Characterize inflammatory infiltrates and possible viral transcription in primary hyperparathyroidism. DESIGN: From the period 2007 to 2016, a total of 55 parathyroid tumors (51 adenomas and 4 hyperplasias) with prominent inflammatory infiltrates were identified from more than 2000 parathyroid tumors in the pathology archives, and investigated by immunohistochemistry for CD4, CD8, CD20 and CD45 and scored as +0, +1 or +2. Clinicopathological data were compared to 142 parathyroid adenomas without histological evidence of inflammation. Transcriptome sequencing was performed for 13 parathyroid tumors (four inflammatory, 9 non-inflammatory) to identify potential viral transcripts. RESULTS: Tumors had prominent germinal center-like nodular (+2) lymphocytic infiltrates consisting of T and B lymphocytes (31%) and/or diffuse (+1-2) infiltrates of predominantly CD8+ T lymphocytes (84%). In the majority of cases with adjacent normal parathyroid tissue, the normal rim was unaffected by the inflammatory infiltrates (96%). Presence of inflammatory infiltrates was associated with higher levels of serum-PTH (P = 0.007) and oxyphilic differentiation (P = 0.002). Co-existent autoimmune disease was observed in 27% of patients with inflammatory infiltrates, which in turn was associated with oxyphilic differentiation (P = 0.041). Additionally, prescription of anti-inflammatory drugs was associated with lower serum ionized calcium (P = 0.037). CONCLUSIONS: No evidence of virus-like sequences in the parathyroid tumors could be found by transcriptome sequencing, suggesting that other factors may contribute to attract the immune system to the parathyroid tumor tissue.

AS3MT-mediated tolerance to arsenic evolved by multiple independent horizontal gene transfers from bacteria to eukaryotes.

Organisms have evolved the ability to tolerate toxic substances in their environments, often by producing metabolic enzymes that efficiently detoxify the toxicant. Inorganic arsenic is one of the most toxic and carcinogenic substances in the environment, but many organisms, including humans, metabolise inorganic arsenic to less toxic metabolites. This multistep process produces mono-, di-, and trimethylated arsenic metabolites, which the organism excretes. In humans, arsenite methyltransferase (AS3MT) appears to be the main metabolic enzyme that methylates arsenic. In this study, we examined the evolutionary origin of AS3MT and assessed the ability of different genotypes to produce methylated arsenic metabolites. Phylogenetic analysis suggests that multiple, independent horizontal gene transfers between different bacteria, and from bacteria to eukaryotes, increased tolerance to environmental arsenic during evolution. These findings are supported by the observation that genetic variation in AS3MT correlates with the capacity to methylate arsenic. Adaptation to arsenic thus serves as a model for how organisms evolve to survive under toxic conditions.

Personal model-assisted identification of NAD+ and glutathione metabolism as intervention target in NAFLD.

To elucidate the molecular mechanisms underlying non-alcoholic fatty liver disease (NAFLD), we recruited 86 subjects with varying degrees of hepatic steatosis (HS). We obtained experimental data on lipoprotein fluxes and used these individual measurements as personalized constraints of a hepatocyte genome-scale metabolic model to investigate metabolic differences in liver, taking into account its interactions with other tissues. Our systems level analysis predicted an altered demand for NAD+ and glutathione (GSH) in subjects with high HS Our analysis and metabolomic measurements showed that plasma levels of glycine, serine, and associated metabolites are negatively correlated with HS, suggesting that these GSH metabolism precursors might be limiting. Quantification of the hepatic expression levels of the associated enzymes further pointed to altered de novo GSH synthesis. To assess the effect of GSH and NAD+ repletion on the development of NAFLD, we added precursors for GSH and NAD+ biosynthesis to the Western diet and demonstrated that supplementation prevents HS in mice. In a proof-of-concept human study, we found improved liver function and decreased HS after supplementation with serine (a precursor to glycine) and hereby propose a strategy for NAFLD treatment.

Whole-Genome Sequencing of Invasion-Resistant Cells Identifies Laminin α2 as a Host Factor for Bacterial Invasion.

To understand the role of glycosaminoglycans in bacterial cellular invasion, xylosyltransferase-deficient mutants of Chinese hamster ovary (CHO) cells were created using clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated gene 9 (CRISPR-cas9) gene targeting. When these mutants were compared to the pgsA745 cell line, a CHO xylosyltransferase mutant generated previously using chemical mutagenesis, an unexpected result was obtained. Bacterial invasion of pgsA745 cells by group B Streptococcus (GBS), group A Streptococcus, and Staphylococcus aureus was markedly reduced compared to the invasion of wild-type cells, but newly generated CRISPR-cas9 mutants were only resistant to GBS. Invasion of pgsA745 cells was not restored by transfection with xylosyltransferase, suggesting that an additional mutation conferring panresistance to multiple bacteria was present in pgsA745 cells. Whole-genome sequencing and transcriptome sequencing (RNA-Seq) uncovered a deletion in the gene encoding the laminin subunit α2 (Lama2) that eliminated much of domain L4a. Silencing of the long Lama2 isoform in wild-type cells strongly reduced bacterial invasion, whereas transfection with human LAMA2 cDNA significantly enhanced invasion in pgsA745 cells. The addition of exogenous laminin-α2ß1γ1/laminin-α2ß2γ1 strongly increased bacterial invasion in CHO cells, as well as in human alveolar basal epithelial and human brain microvascular endothelial cells. Thus, the L4a domain in laminin α2 is important for cellular invasion by a number of bacterial pathogens. IMPORTANCE: Pathogenic bacteria penetrate host cellular barriers by attachment to extracellular matrix molecules, such as proteoglycans, laminins, and collagens, leading to invasion of epithelial and endothelial cells. Here, we show that cellular invasion by the human pathogens group B Streptococcus, group A Streptococcus, and Staphylococcus aureus depends on a specific domain of the laminin α2 subunit. This finding may provide new leads for the molecular pathogenesis of these bacteria and the development of novel antimicrobial drugs.

The Human Adrenal Gland Proteome Defined by Transcriptomics and Antibody-Based Profiling.

The adrenal gland is a composite endocrine organ with vital functions that include the synthesis and release of glucocorticoids and catecholamines. To define the molecular landscape that underlies the specific functions of the adrenal gland, we combined a genome-wide transcriptomics approach using messenger RNA sequencing of human tissues with immunohistochemistry-based protein profiling on tissue microarrays. Approximately two-thirds of all putative protein coding genes were expressed in the adrenal gland, and the analysis identified 253 genes with an elevated pattern of expression in the adrenal gland, with only 37 genes showing a markedly greater expression level (more than fivefold) in the adrenal gland compared with 31 other normal human tissue types analyzed. The analyses allowed for an assessment of the relative expression levels for well-known proteins involved in adrenal gland function but also identified previously poorly characterized proteins in the adrenal cortex, such as the FERM (4.1 protein, ezrin, radixin, moesin) domain containing 5 and the nephroblastoma overexpressed (NOV) protein homolog. We have provided a global analysis of the adrenal gland transcriptome and proteome, with a comprehensive list of genes with elevated expression in the adrenal gland and spatial information with examples of protein expression patterns for corresponding proteins. These genes and proteins constitute important starting points for an improved understanding of the normal function and pathophysiology of the adrenal glands.

Evolutionary engineering reveals divergent paths when yeast is adapted to different acidic environments.

Tolerance of yeast to acid stress is important for many industrial processes including organic acid production. Therefore, elucidating the molecular basis of long term adaptation to acidic environments will be beneficial for engineering production strains to thrive under such harsh conditions. Previous studies using gene expression analysis have suggested that both organic and inorganic acids display similar responses during short term exposure to acidic conditions. However, biological mechanisms that will lead to long term adaptation of yeast to acidic conditions remains unknown and whether these mechanisms will be similar for tolerance to both organic and inorganic acids is yet to be explored. We therefore evolved Saccharomyces cerevisiae to acquire tolerance to HCl (inorganic acid) and to 0.3M L-lactic acid (organic acid) at pH 2.8 and then isolated several low pH tolerant strains. Whole genome sequencing and RNA-seq analysis of the evolved strains revealed different sets of genome alterations suggesting a divergence in adaptation to these two acids. An altered sterol composition and impaired iron uptake contributed to HCl tolerance whereas the formation of a multicellular morphology and rapid lactate degradation was crucial for tolerance to high concentrations of lactic acid. Our findings highlight the contribution of both the selection pressure and nature of the acid as a driver for directing the evolutionary path towards tolerance to low pH. The choice of carbon source was also an important factor in the evolutionary process since cells evolved on two different carbon sources (raffinose and glucose) generated a different set of mutations in response to the presence of lactic acid. Therefore, different strategies are required for a rational design of low pH tolerant strains depending on the acid of interest.

The terpene synthase gene family in Tripterygium wilfordii harbors a labdane-type diterpene synthase among the monoterpene synthase TPS-b subfamily.

Tripterygium wilfordii (Celastraceae) is a medicinal plant with anti-inflammatory and immunosuppressive properties. Identification of a vast array of unusual sesquiterpenoids, diterpenoids and triterpenoids in T. wilfordii has spurred investigations of their pharmacological properties. The tri-epoxide lactone triptolide was the first of many diterpenoids identified, attracting interest due to the spectrum of bioactivities. To probe the genetic underpinning of diterpenoid diversity, an expansion of the class II diterpene synthase (diTPS) family was recently identified in a leaf transcriptome. Following detection of triptolide and simple diterpene scaffolds in the root, we sequenced and mined the root transcriptome. This allowed identification of the root-specific complement of TPSs and an expansion in the class I diTPS family. Functional characterization of the class II diTPSs established their activities in the formation of four C-20 diphosphate intermediates, precursors of both generalized and specialized metabolism and a novel scaffold for Celastraceae. Functional pairs of the class I and II enzymes resulted in formation of three scaffolds, accounting for some of the terpenoid diversity found in T. wilfordii. The absence of activity-forming abietane-type diterpenes encouraged further testing of TPSs outside the canonical class I diTPS family. TwTPS27, close relative of mono-TPSs, was found to couple with TwTPS9, converting normal-copalyl diphosphate to miltiradiene. The phylogenetic distance to established diTPSs indicates neo-functionalization of TwTPS27 into a diTPS, a function not previously observed in the TPS-b subfamily. This example of evolutionary convergence expands the functionality of TPSs in the TPS-b family and may contribute miltiradiene to the diterpenoids of T. wilfordii.

Gene-specific correlation of RNA and protein levels in human cells and tissues.

An important issue for molecular biology is to establish whether transcript levels of a given gene can be used as proxies for the corresponding protein levels. Here, we have developed a targeted proteomics approach for a set of human non-secreted proteins based on parallel reaction monitoring to measure, at steady-state conditions, absolute protein copy numbers across human tissues and cell lines and compared these levels with the corresponding mRNA levels using transcriptomics. The study shows that the transcript and protein levels do not correlate well unless a gene-specific RNA-to-protein (RTP) conversion factor independent of the tissue type is introduced, thus significantly enhancing the predictability of protein copy numbers from RNA levels. The results show that the RTP ratio varies significantly with a few hundred copies per mRNA molecule for some genes to several hundred thousands of protein copies per mRNA molecule for others. In conclusion, our data suggest that transcriptome analysis can be used as a tool to predict the protein copy numbers per cell, thus forming an attractive link between the field of genomics and proteomics.

Profiling cancer testis antigens in non-small-cell lung cancer.

Cancer testis antigens (CTAs) are of clinical interest as biomarkers and present valuable targets for immunotherapy. To comprehensively characterize the CTA landscape of non-small-cell lung cancer (NSCLC), we compared RNAseq data from 199 NSCLC tissues to the normal transcriptome of 142 samples from 32 different normal organs. Of 232 CTAs currently annotated in the Caner Testis Database (CTdatabase), 96 were confirmed in NSCLC. To obtain an unbiased CTA profile of NSCLC, we applied stringent criteria on our RNAseq data set and defined 90 genes as CTAs, of which 55 genes were not annotated in the CTdatabase, thus representing potential new CTAs. Cluster analysis revealed that CTA expression is histology dependent and concurrent expression is common. IHC confirmed tissue-specific protein expression of selected new CTAs (TKTL1, TGIF2LX, VCX, and CXORF67). Furthermore, methylation was identified as a regulatory mechanism of CTA expression based on independent data from The Cancer Genome Atlas. The proposed prognostic impact of CTAs in lung cancer was not confirmed, neither in our RNAseq cohort nor in an independent meta-analysis of 1,117 NSCLC cases. In summary, we defined a set of 90 reliable CTAs, including information on protein expression, methylation, and survival association. The detailed RNAseq catalog can guide biomarker studies and efforts to identify targets for immunotherapeutic strategies.

Oxidation and cyclization of casbene in the biosynthesis of Euphorbia factors from mature seeds of Euphorbia lathyris L.

The seed oil of Euphorbia lathyris L. contains a series of macrocyclic diterpenoids known as Euphorbia factors. They are the current industrial source of ingenol mebutate, which is approved for the treatment of actinic keratosis, a precancerous skin condition. Here, we report an alcohol dehydrogenase-mediated cyclization step in the biosynthetic pathway of Euphorbia factors, illustrating the origin of the intramolecular carbon-carbon bonds present in lathyrane and ingenane diterpenoids. This unconventional cyclization describes the ring closure of the macrocyclic diterpene casbene. Through transcriptomic analysis of E. lathyris L. mature seeds and in planta functional characterization, we identified three enzymes involved in the cyclization route from casbene to jolkinol C, a lathyrane diterpene. These enzymes include two cytochromes P450 from the CYP71 clan and an alcohol dehydrogenase (ADH). CYP71D445 and CYP726A27 catalyze regio-specific 9-oxidation and 5-oxidation of casbene, respectively. When coupled with these P450-catalyzed monooxygenations, E. lathyris ADH1 catalyzes dehydrogenation of the hydroxyl groups, leading to the subsequent rearrangement and cyclization. The discovery of this nonconventional cyclization may provide the key link to complete elucidation of the biosynthetic pathways of ingenol mebutate and other bioactive macrocyclic diterpenoids.

Adaptation to different types of stress converge on mitochondrial metabolism.

Yeast cell factories encounter physical and chemical stresses when used for industrial production of fuels and chemicals. These stresses reduce productivity and increase bioprocess costs. Understanding the mechanisms of the stress response is essential for improving cellular robustness in platform strains. We investigated the three most commonly encountered industrial stresses for yeast (ethanol, salt, and temperature) to identify the mechanisms of general and stress-specific responses under chemostat conditions in which specific growth rate-dependent changes are eliminated. By applying systems-level analysis, we found that most stress responses converge on mitochondrial processes. Our analysis revealed that stress-specific factors differ between applied stresses; however, they are underpinned by an increased ATP demand. We found that when ATP demand increases to high levels, respiration cannot provide sufficient ATP, leading to onset of respirofermentative metabolism. Although stress-specific factors increase ATP demand for cellular growth under stressful conditions, increased ATP demand for cellular maintenance underpins a general stress response and is responsible for the onset of overflow metabolism.

Integrated Network Analysis Reveals an Association between Plasma Mannose Levels and Insulin Resistance.

To investigate the biological processes that are altered in obese subjects, we generated cell-specific integrated networks (INs) by merging genome-scale metabolic, transcriptional regulatory and protein-protein interaction networks. We performed genome-wide transcriptomics analysis to determine the global gene expression changes in the liver and three adipose tissues from obese subjects undergoing bariatric surgery and integrated these data into the cell-specific INs. We found dysregulations in mannose metabolism in obese subjects and validated our predictions by detecting mannose levels in the plasma of the lean and obese subjects. We observed significant correlations between plasma mannose levels, BMI, and insulin resistance (IR). We also measured plasma mannose levels of the subjects in two additional different cohorts and observed that an increased plasma mannose level was associated with IR and insulin secretion. We finally identified mannose as one of the best plasma metabolites in explaining the variance in obesity-independent IR.

Transcriptomics resources of human tissues and organs.

Quantifying the differential expression of genes in various human organs, tissues, and cell types is vital to understand human physiology and disease. Recently, several large-scale transcriptomics studies have analyzed the expression of protein-coding genes across tissues. These datasets provide a framework for defining the molecular constituents of the human body as well as for generating comprehensive lists of proteins expressed across tissues or in a tissue-restricted manner. Here, we review publicly available human transcriptome resources and discuss body-wide data from independent genome-wide transcriptome analyses of different tissues. Gene expression measurements from these independent datasets, generated using samples from fresh frozen surgical specimens and postmortem tissues, are consistent. Overall, the different genome-wide analyses support a distribution in which many proteins are found in all tissues and relatively few in a tissue-restricted manner. Moreover, we discuss the applications of publicly available omics data for building genome-scale metabolic models, used for analyzing cell and tissue functions both in physiological and in disease contexts.

The Urinary Bladder Transcriptome and Proteome Defined by Transcriptomics and Antibody-Based Profiling.

To understand functions and diseases of urinary bladder, it is important to define its molecular constituents and their roles in urinary bladder biology. Here, we performed genome-wide deep RNA sequencing analysis of human urinary bladder samples and identified genes up-regulated in the urinary bladder by comparing the transcriptome data to those of all other major human tissue types. 90 protein-coding genes were elevated in the urinary bladder, either with enhanced expression uniquely in the urinary bladder or elevated expression together with at least one other tissue (group enriched). We further examined the localization of these proteins by immunohistochemistry and tissue microarrays and 20 of these 90 proteins were localized to the whole urothelium with a majority not yet described in the context of the urinary bladder. Four additional proteins were found specifically in the umbrella cells (Uroplakin 1a, 2, 3a, and 3b), and three in the intermediate/basal cells (KRT17, PCP4L1 and ATP1A4). 61 of the 90 elevated genes have not been previously described in the context of urinary bladder and the corresponding proteins are interesting targets for more in-depth studies. In summary, an integrated omics approach using transcriptomics and antibody-based profiling has been used to define a comprehensive list of proteins elevated in the urinary bladder.

The Human Endometrium-Specific Proteome Defined by Transcriptomics and Antibody-Based Profiling.

The human uterus includes the complex endometrial mucosa, the endometrium that undergoes dynamic, hormone-dependent alterations throughout the life of fertile females. Here we have combined a genome-wide transcriptomics analysis with immunohistochemistry-based protein profiling to analyze gene expression patterns in the normal endometrium. Human endometrial tissues from five women were used for deep sequencing (RNA-Seq). The mRNA and protein expression data from the endometrium were compared to 31 (RNA) and 44 (protein) other normal tissue types, to identify genes with elevated expression in the endometrium and to localize the expression of corresponding proteins at a cellular resolution. Based on the expression levels of transcripts, we could classify all putative human protein coding genes into categories defined by expression patterns and found altogether 101 genes that showed an elevated pattern of expression in the endometrium, with only four genes showing more than five-fold higher expression levels in the endometrium compared to other tissues. In conclusion, our analysis based on transcriptomics and antibody-based protein profiling reports here comprehensive lists of genes with elevated expression levels in the endometrium, providing important starting points for a better molecular understanding of human reproductive biology and disease.

Analysis of the Human Prostate-Specific Proteome Defined by Transcriptomics and Antibody-Based Profiling Identifies TMEM79 and ACOXL as Two Putative, Diagnostic Markers in Prostate Cancer.

To better understand prostate function and disease, it is important to define and explore the molecular constituents that signify the prostate gland. The aim of this study was to define the prostate specific transcriptome and proteome, in comparison to 26 other human tissues. Deep sequencing of mRNA (RNA-seq) and immunohistochemistry-based protein profiling were combined to identify prostate specific gene expression patterns and to explore tissue biomarkers for potential clinical use in prostate cancer diagnostics. We identified 203 genes with elevated expression in the prostate, 22 of which showed more than five-fold higher expression levels compared to all other tissue types. In addition to previously well-known proteins we identified two poorly characterized proteins, TMEM79 and ACOXL, with potential to differentiate between benign and cancerous prostatic glands in tissue biopsies. In conclusion, we have applied a genome-wide analysis to identify the prostate specific proteome using transcriptomics and antibody-based protein profiling to identify genes with elevated expression in the prostate. Our data provides a starting point for further functional studies to explore the molecular repertoire of normal and diseased prostate including potential prostate cancer markers such as TMEM79 and ACOXL.