New technologies enable single-cell transcriptome analysis, mapping genome-wide expression across the human body. Here, we present an extended analysis of protein-coding genes in all major human tissues and organs, combining single-cell and bulk transcriptomics. To enhance transcriptome depth, 31 tissues were analyzed using a pooling method, identifying 557 unique cell clusters, manually annotated by marker gene expression. Genes were classified by body-wide expression and validated through antibody-based profiling. All results are available in the updated open-access Single Cell Type section of the Human Protein Atlas for genome-wide exploration of genes, proteins, and their spatial distribution in cells.

Although most pituitary neuroendocrine tumors (PitNETs)/pituitary adenomas remain intrasellar, a significant proportion of tumors show parasellar invasive growth and 6% to 8% infiltrate the bone structures, thus affecting the prognosis. There is an unmet need to identify novel markers that can predict the parasellar growth of PitNETs. Furthermore, mechanisms that regulate bone invasiveness of PitNETs and factors related to tumor vascularization are largely unknown.We used genome-wide mRNA analysis in a cohort of 77 patients with PitNETs of different types to explore the differences in gene expression patterns between invasive and noninvasive tumors with respect to the parasellar growth and regarding the rare phenomenon of bone invasiveness. Additionally, we studied the genes correlated to the contrast enhancement quotient, a novel radiological parameter of tumor vascularization.Most of the genes differentially expressed related to the parasellar growth were genes involved in tumor invasiveness. Differentially expressed genes associated with bone invasiveness are involved in NF-kappa B pathway and antitumoral immune response. Lack of clear clustering regarding the parasellar and bone invasiveness may be explained by the influence of the cell lineage-related genes in this heterogeneous cohort of PitNETs.Our transcriptomics analysis revealed differences in the molecular fingerprints between invasive, including bone invasive, and noninvasive PitNETs, although without clear clustering. The contrast enhancement quotient emerged as a radiological parameter of tumor vascularization, correlating with several angiogenesis-related genes. Several of the top genes related to the PitNET invasiveness and vascularization have potential prognostic and therapeutic application requiring further research.

Background: There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs. Results: An open-access pig expression map (www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented. Conclusions: Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource (www.rnaatlas.org), including a comparison to the expression of human orthologs.

Pituitary neuroendocrine tumors (PitNETs) are common, generally benign tumors with complex clinical characteristics related to hormone hypersecretion and/or growing sellar tumor mass. PitNETs can be classified based on the expression pattern of anterior pituitary hormones and three main transcriptions factors (TF), SF1, PIT1 and TPIT that regulate differentiation of adenohypophysial cells. Here, we have extended this classification based on the global transcriptomics landscape using tumor tissue from a well-defined cohort comprising 51 PitNETs of different clinical and histological types. The molecular profiles were compared with current classification schemes based on immunohistochemistry. Our results identified three main clusters of PitNETs that were aligned with the main pituitary TFs expression patterns. Our analyses enabled further identification of specific genes and expression patterns, including both known and unknown genes, that could distinguish the three different classes of PitNETs. We conclude that the current classification of PitNETs based on the expression of SF1, PIT1 and TPIT reflects three distinct subtypes of PitNETs with different underlying biology and partly independent from the expression of corresponding hormones. The transcriptomic analysis reveals several potentially targetable tumor-driving genes with previously unknown role in pituitary tumorigenesis.

Simple Summary Pyruvate kinase muscle type (PKM) is a key enzyme in glycolysis and is a mediator of the Warburg effect in tumors. The association of PKM with survival of cancer patients is controversial. In this study, we investigated the associations of the alternatively spliced transcripts of PKM with cancer patients' survival outcomes and explained the conflicts in previous studies. We discovered three poorly studied alternatively spliced PKM transcripts that exhibited opposite prognostic indications in different human cancers based on integrative systems analysis. We also detected their protein products and explored their potential biological functions based on in-vitro experiments. Our analysis demonstrated that alternatively spliced transcripts of not only PKM but also other genes should be considered in cancer studies, since it may enable the discovery and targeting of the right protein product for development of the efficient treatment strategies. Pyruvate kinase muscle type (PKM) is a key enzyme in glycolysis and plays an important oncological role in cancer. However, the association of PKM expression and the survival outcome of patients with different cancers is controversial. We employed systems biology methods to reveal prognostic value and potential biological functions of PKM transcripts in different human cancers. Protein products of transcripts were shown and detected by western blot and mass spectrometry analysis. We focused on different transcripts of PKM and investigated the associations between their mRNA expression and the clinical survival of the patients in 25 different cancers. We find that the transcripts encoding PKM2 and three previously unstudied transcripts, namely ENST00000389093, ENST00000568883, and ENST00000561609, exhibited opposite prognostic indications in different cancers. Moreover, we validated the prognostic effect of these transcripts in an independent kidney cancer cohort. Finally, we revealed that ENST00000389093 and ENST00000568883 possess pyruvate kinase enzymatic activity and may have functional roles in metabolism, cell invasion, and hypoxia response in cancer cells. Our study provided a potential explanation to the controversial prognostic indication of PKM, and could invoke future studies focusing on revealing the biological and oncological roles of these alternative spliced variants of PKM.

Using the Cap Analysis of Gene Expression (CAGE) technology, the FANTOM5 consortium provided one of the most comprehensive maps of transcription start sites (TSSs) in several species. Strikingly, ~72% of them could not be assigned to a specific gene and initiate at unconventional regions, outside promoters or enhancers. Here, we probe these unassigned TSSs and show that, in all species studied, a significant fraction of CAGE peaks initiate at microsatellites, also called short tandem repeats (STRs). To confirm this transcription, we develop Cap Trap RNA-seq, a technology which combines cap trapping and long read MinION sequencing. We train sequence-based deep learning models able to predict CAGE signal at STRs with high accuracy. These models unveil the importance of STR surrounding sequences not only to distinguish STR classes, but also to predict the level of transcription initiation. Importantly, genetic variants linked to human diseases are preferentially found at STRs with high transcription initiation level, supporting the biological and clinical relevance of transcription initiation at STRs. Together, our results extend the repertoire of non-coding transcription associated with DNA tandem repeats and complexify STR polymorphism.

Objective. The COVID-19 pandemic poses an immense need for accurate, sensitive and high-throughput clinical tests, and serological assays are needed for both overarching epidemiological studies and evaluating vaccines. Here, we present the development and validation of a high-throughput multiplex bead-based serological assay. Methods. More than 100 representations of SARS-CoV-2 proteins were included for initial evaluation, including antigens produced in bacterial and mammalian hosts as well as synthetic peptides. The five best-performing antigens, three representing the spike glycoprotein and two representing the nucleocapsid protein, were further evaluated for detection of IgG antibodies in samples from 331 COVID-19 patients and convalescents, and in 2090 negative controls sampled before 2020. Results. Three antigens were finally selected, represented by a soluble trimeric form and the S1-domain of the spike glycoprotein as well as by the C-terminal domain of the nucleocapsid. The sensitivity for these three antigens individually was found to be 99.7%, 99.1% and 99.7%, and the specificity was found to be 98.1%, 98.7% and 95.7%. The best assay performance was although achieved when utilising two antigens in combination, enabling a sensitivity of up to 99.7% combined with a specificity of 100%. Requiring any two of the three antigens resulted in a sensitivity of 99.7% and a specificity of 99.4%. Conclusion. These observations demonstrate that a serological test based on a combination of several SARS-CoV-2 antigens enables a highly specific and sensitive multiplex serological COVID-19 assay.

Simple Summary: Pituitary neuroendocrine tumours originate from the endocrine cells of the anterior pituitary gland and may develop from any of the cell lineages responsible for producing the different pituitary hormones. The details related to tumour differentiation and hormone production in these tumours are not fully understood. The aim of our study was to investigate an uncharacterised pituitary enriched protein, transforming growth factor beta-receptor 3 like (TGFBR3L). The TGFBR3L protein is highly expressed in the pituitary compared to other organs. We found the protein to be gonadotroph-specific, i.e., detected in the cells that express follicle-stimulating and luteinizing hormones (FSH/LH). The gonadotroph-specific nature of TGFBR3L, a correlation to both FSH and LH as well as an inverse correlation to membranous E-cadherin and oestrogen receptor beta suggests a role in gonadotroph cell development and function and, possibly, tumour progression. Here, we report the investigation of transforming growth factor beta-receptor 3 like (TGFBR3L), an uncharacterised pituitary specific membrane protein, in non-neoplastic anterior pituitary gland and pituitary neuroendocrine tumours. A polyclonal antibody produced within the Human Protein Atlas project (HPA074356) was used for TGFBR3L staining and combined with SF1 and FSH for a 3-plex fluorescent protocol, providing more details about the cell lineage specificity of TGFBR3L expression. A cohort of 230 pituitary neuroendocrine tumours were analysed. In a subgroup of previously characterised gonadotroph tumours, correlation with expression of FSH/LH, E-cadherin, oestrogen (ER) and somatostatin receptors (SSTR) was explored. TGFBR3L showed membranous immunolabeling and was found to be gonadotroph cell lineage-specific, verified by co-expression with SF1 and FSH/LH staining in both tumour and non-neoplastic anterior pituitary tissues. TGFBR3L immunoreactivity was observed in gonadotroph tumours only and demonstrated intra-tumour heterogeneity with a perivascular location. TGFBR3L immunostaining correlated positively to both FSH (R = 0.290) and LH (R = 0.390) immunostaining, and SSTR3 (R = 0.315). TGFBR3L correlated inversely to membranous E-cadherin staining (R = -0.351) and oestrogen receptor beta mRNA (R = -0.274). In conclusion, TGFBR3L is a novel pituitary gland specific protein, located in the membrane of gonadotroph cells in non-neoplastic anterior pituitary gland and in a subset of gonadotroph pituitary tumours.

The localization of proteins at a tissue- or cell-type-specific level is tightly linked to the protein function. To better understand each protein's role in cellular systems, spatial information constitutes an important complement to quantitative data. The standard methods for determining the spatial distribution of proteins in single cells of complex tissue samples make use of antibodies. For a stringent analysis of the human proteome, we used orthogonal methods and independent antibodies to validate 5981 antibodies that show the expression of 3775 human proteins across all major human tissues. This enhanced validation uncovered 56 proteins corresponding to the group of "missing proteins" and 171 proteins of unknown function. The presented strategy will facilitate further discussions around criteria for evidence of protein existence based on immunohistochemistry and serves as a useful guide to identify candidate proteins for integrative studies with quantitative proteomics methods.

The proteins secreted by human tissues and blood cells, the secretome, are important both for the basic understanding of human biology and for identification of potential targets for future diagnosis and therapy. Here, a high-throughput mammalian cell factory is presented that was established to create a resource of recombinant full-length proteins covering the majority of those annotated as 'secreted' in humans. The full-length DNA sequences of each of the predicted secreted proteins were generated by gene synthesis, the constructs were transfected into Chinese hamster ovary (CHO) cells and the recombinant proteins were produced, purified and analyzed. Almost 1,300 proteins were successfully generated and proteins predicted to be secreted into the blood were produced with a success rate of 65%, while the success rates for the other categories of secreted proteins were somewhat lower giving an overall one-pass success rate of ca. 58%. The proteins were used to generate targeted proteomics assays and several of the proteins were shown to be active in a phenotypic assay involving pancreatic beta-cell dedifferentiation. Many of the proteins that failed during production in CHO cells could be rescued in human embryonic kidney (HEK 293) cells suggesting that a cell factory of human origin can be an attractive alternative for production in mammalian cells. In conclusion, a high-throughput protein production and purification system has been successfully established to create a unique resource of the human secretome.