Complementing tissue characterization by integrating transcriptome profiling from the Human Protein Atlas and from the FANTOM5 consortium

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic β-cell identity and function. Elimination of Lkb1 from the β-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of β cell-selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from β-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors (TFs) important for β-cell identity, such as FOXA, MAFA or RFX6, and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates β-cell identity and function.

Anatomical and molecular insights into the antennal gland of the giant freshwater prawn Macrobrachium rosenbergii

In this study, the complex organization of the AnG in the giant freshwater prawn Macrobrachium rosenbergii was revealed using various techniques, including conventional histology, histochemistry, scanning electron microscopy, and X-ray tomography. The results showed the diversity of cells in the AnG and the detailed organization of the labyrinth's tubule into four radiated areas from the central to peripheral zones. The study also demonstrated the expression of some vertebrate kidney-associated homolog genes, aquaporin (AQP), solute carrier family 22 (SLC-22), nephrin, and uromodulin, in the AnG by qPCR. The result of in situ hybridization further showed the localization of SLC-22 and AQP transcript in the bladder and labyrinth's epithelium, specifically in regions 2, 3, and 4. Additionally, the study revealed neuropeptide expressions in the AnG by qPCR and in situ hybridization, i.e., crustacean hyperglycemic hormone (CHH) and molt inhibiting hormone (MIH), implying that the AnG may have a role in hormone production. Moreover, male and female prawns exhibited different levels of AQP, SLC-22, nephrin, and CHH expressions during the premolt and intermolt stages, suggesting a crucial role relevant to the molting stages. In conclusion, this study clarified the complex structure of the AnG in M. rosenbergii and demonstrated for the first time the expression of vertebrate kidney-associated genes and the possible endocrine role of the AnG. Further investigation is needed to clarify the role of these genes, particularly during ecdysis. The implications of these findings could significantly advance our understanding of the AnG in decapod crustaceans.

Evaluation of the impact of iPSC differentiation protocols on transcriptomic signatures

Human induced pluripotent stem cells (iPSC) have the potential to produce desired target cell types in vitro and allow for the high-throughput screening of drugs/chemicals at population level thereby minimising the cost of drug discovery and drug withdrawals after clinical trials. There is a substantial need for the characterisation of the iPSC derived models to better understand and utilise them for toxicological relevant applications. In our study, iPSC (SBAD2 or SBAD3 lines obtained from StemBANCC project) were differentiated towards toxicologically relevant cell types: alveolar macrophages, brain capillary endothelial cells, brain cells, endothelial cells, hepatocytes, lung airway epithelium, monocytes, podocytes and renal proximal tubular cells. A targeted transcriptomic approach was employed to understand the effects of differentiation protocols on these cell types. Pearson correlation and principal component analysis (PCA) separated most of the intended target cell types and undifferentiated iPSC models as distinct groups with a high correlation among replicates from the same model. Based on PCA, the intended target cell types could also be separated into the three germ layer groups (ectoderm, endoderm and mesoderm). Differential expression analysis (DESeq2) presented the upregulated genes in each intended target cell types that allowed the evaluation of the differentiation to certain degree and the selection of key differentiation markers. In conclusion, these data confirm the versatile use of iPSC differentiated cell types as standardizable and relevant model systems for in vitro toxicology.

Liver kinase B1 (LKB1) regulates the epigenetic landscape of mouse pancreatic beta cells

Liver kinase B1 (LKB1/STK11) is an important regulator of pancreatic β-cell identity and function. Elimination of Lkb1 from the β-cell results in improved glucose-stimulated insulin secretion and is accompanied by profound changes in gene expression, including the upregulation of several neuronal genes. The mechanisms through which LKB1 controls gene expression are, at present, poorly understood. Here, we explore the impact of β cell- selective deletion of Lkb1 on chromatin accessibility in mouse pancreatic islets. To characterize the role of LKB1 in the regulation of gene expression at the transcriptional level, we combine these data with a map of islet active transcription start sites and histone marks. We demonstrate that LKB1 elimination from β-cells results in widespread changes in chromatin accessibility, correlating with changes in transcript levels. Changes occurred in hundreds of promoter and enhancer regions, many of which were close to neuronal genes. We reveal that dysregulated enhancers are enriched in binding motifs for transcription factors important for β-cell identity, such as FOXA, MAFA or RFX6 and we identify microRNAs (miRNAs) that are regulated by LKB1 at the transcriptional level. Overall, our study provides important new insights into the epigenetic mechanisms by which LKB1 regulates β-cell identity and function.

Correlation between the RNA Expression and the DNA Methylation of Estrogen Receptor Genes in Normal and Malignant Human Tissues

Estrogen plays a multifaceted function in humans via interacting with the estrogen receptors ERα, ERβ, and G protein-coupled estrogen receptor 1 (GPER1). Previous research has predominantly concentrated on elucidating the signaling route of estrogen. However, the comprehensive understanding of the expression profile and control of these estrogen receptors in various human tissues is not well known. In the present study, the RNA levels of estrogen receptors in various normal and malignant human tissues were retrieved from the human protein atlas, the cancer genome atlas (TCGA), and the genotype-tissue expression (GTEx) databases for analyzing the expression profile of estrogen receptors through gene expression profiling interactive analysis (GEPIA). The status of DNA methylation of estrogen receptor genes from TCGA were analyzed through the software Wanderer and cBioPortal. The MethSurv tool was utilized to estimate the relevance between specific cytosine-guanine (CG) methylation and tumor survival. The expression profile analysis revealed that ERα, ERβ, and GPER1 have unique expression patterns in diverse tissues and malignancies. The interesting results were the higher expression of ERβ RNA in the male testis than in females and the positive association between the RNA level of ERα and the androgen receptor in different human normal tissues. Especially, the significant changes in GPER1 expression in multiple malignancies showed a consistent decrease with no exception, which indicates the role of GPER1 in common tumor inhibition. The finding on the expression profile provides clues for exploring novel potential physiological and pathophysiological functions of estrogen. The DNA methylation analysis manifested that the expression of GPER1 and ERα showed a substantial correlation with the methylation of specific CG sites in the cis-regulating region of the gene. However, no such association was observed for ERβ. When comparing tumor tissues to normal tissues, the DNA methylation of certain CG sites of estrogen receptors showed a correlation with tumor survival but did not always correlate with the expression of that gene or with the expression of DNA methyltransferases. We proposed that the variation in DNA methylation at different CG sites in estrogen receptor genes had other functions beyond its regulatory role in its gene expression, and this might be associated with the progression and therapy efficiency of the tumor based on the modulation of the chromatin configuration.

ICAM-1 nanoclusters regulate hepatic epithelial cell polarity by leukocyte adhesion-independent control of apical actomyosin

Epithelial intercellular adhesion molecule (ICAM)-1 is apically polarized, interacts with, and guides leukocytes across epithelial barriers. Polarized hepatic epithelia organize their apical membrane domain into bile canaliculi and ducts, which are not accessible to circulating immune cells but that nevertheless confine most of ICAM-1. Here, by analyzing ICAM-1_KO human hepatic cells, liver organoids from ICAM-1_KO mice and rescue-of-function experiments, we show that ICAM-1 regulates epithelial apicobasal polarity in a leukocyte adhesion-independent manner. ICAM-1 signals to an actomyosin network at the base of canalicular microvilli, thereby controlling the dynamics and size of bile canalicular-like structures. We identified the scaffolding protein EBP50/NHERF1/SLC9A3R1, which connects membrane proteins with the underlying actin cytoskeleton, in the proximity interactome of ICAM-1. EBP50 and ICAM-1 form nano-scale domains that overlap in microvilli, from which ICAM-1 regulates EBP50 nano-organization. Indeed, EBP50 expression is required for ICAM-1-mediated control of BC morphogenesis and actomyosin. Our findings indicate that ICAM-1 regulates the dynamics of epithelial apical membrane domains beyond its role as a heterotypic cell-cell adhesion molecule and reveal potential therapeutic strategies for preserving epithelial architecture during inflammatory stress.

High-throughput omics technologies in inflammatory bowel disease

Inflammatory bowel disease (IBD) is a chronic, relapsing intestinal disease. Elucidation of the pathogenic mechanisms of IBD requires high-throughput technologies (HTTs) to effectively obtain and analyze large amounts of data. Recently, HTTs have been widely used in IBD, including genomics, transcriptomics, proteomics, microbiomics, metabolomics and single-cell sequencing. When combined with endoscopy, the application of these technologies can provide an in-depth understanding on the alterations of intestinal microbe diversity and abundance, the abnormalities of signaling pathway-mediated immune responses and functionality, and the evaluation of therapeutic effects, improving the accuracy of early diagnosis and treatment of IBD. This review comprehensively summarizes the development and advancement of HTTs, and also highlights the challenges and future directions of these technologies in IBD research. Although HTTs have made striking breakthrough in IBD, more standardized methods and large-scale dataset processing are still needed to achieve the goal of personalized medicine.

A computational modeling of pri-miRNA expression

MicroRNAs (miRNAs) play crucial roles in gene regulation. Most studies focus on mature miRNAs, which leaves many unknowns about primary miRNAs (pri-miRNAs). To fill the gap, we attempted to model the expression of pri-miRNAs in 1829 primary cell types, cell lines, and tissues in this study. We demonstrated that the expression of pri-miRNAs can be modeled well by the expression of specific sets of mRNAs, which we termed their associated mRNAs. These associated mRNAs differ from their corresponding target mRNAs and are enriched with specific functions. Most associated mRNAs of a miRNA are shared across conditions, while on average, about one-fifth of the associated mRNAs are condition-specific. Our study shed new light on understanding miRNA biogenesis and general gene transcriptional regulation.

Omada: robust clustering of transcriptomes through multiple testing

BACKGROUND

Cohort studies increasingly collect biosamples for molecular profiling and are observing molecular heterogeneity. High-throughput RNA sequencing is providing large datasets capable of reflecting disease mechanisms. Clustering approaches have produced a number of tools to help dissect complex heterogeneous datasets, but selecting the appropriate method and parameters to perform exploratory clustering analysis of transcriptomic data requires deep understanding of machine learning and extensive computational experimentation. Tools that assist with such decisions without prior field knowledge are nonexistent. To address this, we have developed Omada, a suite of tools aiming to automate these processes and make robust unsupervised clustering of transcriptomic data more accessible through automated machine learning-based functions.

FINDINGS

The efficiency of each tool was tested with 7 datasets characterized by different expression signal strengths to capture a wide spectrum of RNA expression datasets. Our toolkit's decisions reflected the real number of stable partitions in datasets where the subgroups are discernible. Within datasets with less clear biological distinctions, our tools either formed stable subgroups with different expression profiles and robust clinical associations or revealed signs of problematic data such as biased measurements.

CONCLUSIONS

In conclusion, Omada successfully automates the robust unsupervised clustering of transcriptomic data, making advanced analysis accessible and reliable even for those without extensive machine learning expertise. Implementation of Omada is available at http://bioconductor.org/packages/omada/.

Disclosing the molecular profile of the human amniotic mesenchymal stromal cell secretome by filter-aided sample preparation proteomic characterization

BACKGROUND

The secretome of mesenchymal stromal cells isolated from the amniotic membrane (hAMSCs) has been extensively studied for its in vitro immunomodulatory activity as well as for the treatment of several preclinical models of immune-related disorders. The bioactive molecules within the hAMSCs secretome are capable of modulating the immune response and thus contribute to stimulating regenerative processes. At present, only a few studies have attempted to define the composition of the secretome, and several approaches, including multi-omics, are underway in an attempt to precisely define its composition and possibly identify key factors responsible for the therapeutic effect.

METHODS

In this study, we characterized the protein composition of the hAMSCs secretome by a filter-aided sample preparation (FASP) digestion and liquid chromatography-high resolution mass spectrometry (LC-MS) approach. Data were processed for gene ontology classification and functional protein interaction analysis by bioinformatics tools.

RESULTS

Proteomic analysis of the hAMSCs secretome resulted in the identification of 1521 total proteins, including 662 unique elements. A number of 157 elements, corresponding to 23.7%, were found as repeatedly characterizing the hAMSCs secretome, and those that resulted as significantly over-represented were involved in immunomodulation, hemostasis, development and remodeling of the extracellular matrix molecular pathways.

CONCLUSIONS

Overall, our characterization enriches the landscape of hAMSCs with new information that could enable a better understanding of the mechanisms of action underlying the therapeutic efficacy of the hAMSCs secretome while also providing a basis for its therapeutic translation.

Integrated analysis of the ubiquitination mechanism reveals the specific signatures of tissue and cancer

BACKGROUND

Ubiquitination controls almost all cellular processes. The dysregulation of ubiquitination signals is closely associated with the initiation and progression of multiple diseases. However, there is little comprehensive research on the interaction and potential function of ubiquitination regulators (UBRs) in spermatogenesis and cancer.

METHODS

We systematically characterized the mRNA and protein expression of UBRs across tissues and further evaluated their roles in testicular development and spermatogenesis. Subsequently, we explored the genetic alterations, expression perturbations, cancer hallmark-related pathways, and clinical relevance of UBRs in pan-cancer.

RESULTS

This work reveals heterogeneity in the expression patterns of UBRs across tissues, and the expression pattern in testis is the most distinct. UBRs are dynamically expressed during testis development, which are critical for normal spermatogenesis. Furthermore, UBRs have widespread genetic alterations and expression perturbations in pan-cancer. The expression of 79 UBRs was identified to be closely correlated with the activity of 32 cancer hallmark-related pathways, and ten hub genes were screened for further clinical relevance analysis by a network-based method. More than 90% of UBRs can affect the survival of cancer patients, and hub genes have an excellent prognostic classification for specific cancer types.

CONCLUSIONS

Our study provides a comprehensive analysis of UBRs in spermatogenesis and pan-cancer, which can build a foundation for understanding male infertility and developing cancer drugs in the aspect of ubiquitination.

Multivariate Evaluation of DNA Quality Differences in Different Preanalytical Procedures in Mouse Livers

Successful histogenetic research relies on proper handling of tissue samples to maximize DNA quality. As the largest gland in the body, the liver is particularly sensitive to sample mishandling owing to its enzymatic and transcriptional activity. However, the impact of preanalytical procedures on the quality of extracted liver DNA remains poorly understood. In this study, we assessed the impact of extraction methods, duration of ex vivo liver ischemia, liver storage time, and temperature on extracted DNA quality. Comprehensive parameters such as DNA yields, purity, DNA integrity number, the percentage of double-stranded DNA (%dsDNA), and PCR amplification of the GAPDH gene fragment were assessed to identify the quality of extracted DNA. Our results revealed that these preanalytical processes had little effect on DIN values and PCR efficiency of GAPDH gene fragments for each sample, whereas the DNA yields, purity, and %dsDNAs varied widely across different processes. For liver DNA extraction, RNase is necessary to isolate "pure" DNA, and the presence of RNase could significantly increase the %dsDNA. In addition, significant increases in the yields, purity, and %dsDNA of extracted DNA were observed in the TissueLyser-processed livers compared with the mortar and pestle or shear cell disruption methods. Further investigation revealed that livers experiencing longer periods of ex vivo ischemia resulted in significantly compromised DNA yields, and to obtain sufficient DNA, the ex vivo liver ischemia should be limited to within 30 minutes. Moreover, compared with storage of livers at -80°C, storage of livers in the vapor phase of liquid nitrogen yielded a higher quality of the extracted DNA. Our findings exhibited significant implications for liver-derived DNA quality assessment and management.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

Multi-scale systems genomics analysis predicts pathways, cell types, and drug targets involved in normative variation in peri-adolescent human cognition

An open challenge in human genetics is to better understand the systems-level impact of genotype variation on developmental cognition. To characterize the genetic underpinnings of peri-adolescent cognition, we performed genotype-phenotype and systems analysis for binarized accuracy in nine cognitive tasks from the Philadelphia Neurodevelopmental Cohort (~2,200 individuals of European continental ancestry aged 8-21 years). We report a region of genome-wide significance within the 3' end of the Fibulin-1 gene (P = 4.6 × 10-8), associated with accuracy in nonverbal reasoning, a heritable form of complex reasoning ability. Diffusion tensor imaging data from a subset of these participants identified a significant association of white matter fractional anisotropy with FBLN1 genotypes (P < 0.025); poor performers show an increase in the C and A allele for rs77601382 and rs5765534, respectively, which is associated with increased fractional anisotropy. Integration of published human brain-specific 'omic maps, including single-cell transcriptomes of the developing human brain, shows that FBLN1 demonstrates greatest expression in the fetal brain, as a marker of intermediate progenitor cells, demonstrates negligible expression in the adolescent and adult human brain, and demonstrates increased expression in the brain in schizophrenia. Collectively these findings warrant further study of this gene and genetic locus in cognition, neurodevelopment, and disease. Separately, genotype-pathway analysis identified an enrichment of variants associated with working memory accuracy in pathways related to development and to autonomic nervous system dysfunction. Top-ranking pathway genes include those genetically associated with diseases with working memory deficits, such as schizophrenia and Parkinson's disease. This work advances the "molecules-to-behavior" view of cognition and provides a framework for using systems-level organization of data for other biomedical domains.

Altered DNA methylome profiles of blood leukocytes in Chinese patients with mild cognitive impairment and Alzheimer's disease

Objective: DNA methylation plays a potential role in the pathogenesis of Alzheimer's disease (AD). However, little is known about the global changes of blood leukocyte DNA methylome profiles from Chinese patients with mild cognitive impairment (MCI) and with AD, or the specific DNA methylation-based signatures associated with MCI and AD. In this study, we sought to dissect the characteristics of blood DNA methylome profiles in MCI- and AD-affected Chinese patients with the aim of identifying novel DNA methylation biomarkers for AD. Methods: In this study, we profiled the DNA methylome of peripheral blood leukocytes from 20 MCI- and 20 AD-affected Chinese patients and 20 cognitively healthy controls (CHCs) with the Infinium Methylation EPIC BeadChip array. Results: We identified significant alterations of the methylome profiles in MCI and AD blood leukocytes. A total of 2,582 and 20,829 CpG sites were significantly and differentially methylated in AD and MCI compared with CHCs (adjusted p < 0.05), respectively. Furthermore, 441 differentially methylated positions (DMPs), aligning to 213 unique genes, were overlapped by the three comparative groups of AD versus CHCs, MCI versus CHCs, and AD versus MCI, of which 6 and 5 DMPs were continuously hypermethylated and hypomethylated in MCI and AD relative to CHCs (adjusted p < 0.05), respectively, such as FLNC cg20186636 and AFAP1 cg06758191. The DMPs with an area under the curve >0.900, such as cg18771300, showed high potency for predicting MCI and AD. In addition, gene ontology and pathway enrichment results showed that these overlapping genes were mainly involved in neurotransmitter transport, GABAergic synaptic transmission, signal release from synapse, neurotransmitter secretion, and the regulation of neurotransmitter levels. Furthermore, tissue expression enrichment analysis revealed a subset of potentially cerebral cortex-enriched genes associated with MCI and AD, including SYT7, SYN3, and KCNT1. Conclusion: This study revealed a number of potential biomarkers for MCI and AD, also highlighted the presence of epigenetically dysregulated gene networks that may engage in the underlying pathological events resulting in the onset of cognitive impairment and AD progression. Collectively, this study provides prospective cues for developing therapeutic strategies to improve cognitive impairment and AD course.

Identification, expression profiling, and functional characterization of cystatin C from big-belly seahorse (Hippocampus abdominalis)

Cystatins are natural inhibitors of lysosomal cysteine proteases, including cathepsins B, L, H, and S. Cystatin C (CSTC) is a member of the type 2 cystatin family and is an essential biomarker in the prognosis of several diseases. Emerging evidence suggests the immune regulatory roles of CSTC in antigen presentation, the release of different inflammatory mediators, and apoptosis in various pathophysiologies. In this study, the 390-bp cystatin C (HaCSTC) cDNA from big-belly seahorse (Hippocampus abdominalis) was cloned and characterized by screening the pre-established cDNA library. Based on similarities in sequence, HaCSTC is a homolog of the teleost type 2 cystatin family with putative catalytic cystatin domains, signal peptides, and disulfide bonds. HaCSTC transcripts were ubiquitously expressed in all tested big-belly seahorse tissues, with the highest expression in ovaries. Immune challenge with lipopolysaccharides, polyinosinic:polycytidylic acid, Edwardsiella tarda, and Streptococcus iniae caused significant upregulation in HaCSTC transcript levels. Using a pMAL-c5X expression vector, the 14.29-kDa protein of recombinant HaCSTC (rHaCSTC) was expressed in Escherichia coli BL21 (DE3), and its protease inhibitory activity against papain cysteine protease was determined with the aid of a protease substrate. Papain was competitively blocked by rHaCSTC in a dose-dependent manner. In response to viral hemorrhagic septicemia virus (VHSV) infection, HaCSTC overexpression strongly decreased the expression of VHSV transcripts, pro-inflammatory cytokines, and pro-apoptotic genes; while increasing the expression of anti-apoptotic genes in fathead minnow (FHM) cells. Furthermore, HaCSTC overexpression protected VHSV-infected FHM cells against VHSV-induced apoptosis and increased cell viability. Our findings imply the profound role of HaCSTC against pathogen infections by modulating fish immune responses.

APPLICATION OF THE ACMG/AMP FRAMEWORK TO CAPTURE EVIDENCE RELEVANT TO PREDICTED AND OBSERVED IMPACT ON SPLICING: RECOMMENDATIONS FROM THE CLINGEN SVI SPLICING SUBGROUP

The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1 (null variant in a gene where loss-of-function is the mechanism of disease), PS3 (functional assays show damaging effect on splicing), PP3 (computational evidence supports a splicing effect), BS3 (functional assays show no damaging effect on splicing), BP4 (computational evidence suggests no splicing impact), and BP7 (silent change with no predicted impact on splicing). However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation (SVI) Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. Our study utilised empirically derived splicing evidence to: 1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, 2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and 3) exemplify methodology to calibrate bioinformatic splice prediction tools. We propose repurposing of the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely BP7 may be used to capture RNA results demonstrating no impact on splicing for both intronic and synonymous variants, and for missense variants if protein functional impact has been excluded. Furthermore, we propose that the PS3 and BS3 codes are applied only for well-established assays that measure functional impact that is not directly captured by RNA splicing assays. We recommend the application of PS1 based on similarity of predicted RNA splicing effects for a variant under assessment in comparison to a known Pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA assay evidence described aim to help standardise variant pathogenicity classification processes and result in greater consistency when interpreting splicing-based evidence.

The human protein atlas-Integrated omics for single cell mapping of the human proteome

Studying the spatial distribution of proteins provides the basis for understanding the biology, molecular repertoire, and architecture of every human cell. The Human Protein Atlas (HPA) has grown into one of the world's largest biological databases, and in the most recent version, a major update of the structure of the database was performed. The data has now been organized into 10 different comprehensive sections, each summarizing different aspects of the human proteome and the protein-coding genes. In particular, large datasets with information on the single cell type level have been integrated, refining the tissue and cell type specificity and detailing the expression in cell states with an increased resolution. The multi-modal data constitute an important resource for both basic and translational science, and hold promise for integration with novel emerging technologies at the protein and RNA level.

Spatial transcriptome analysis of long non-coding RNAs reveals tissue specificity and functional roles in cancer

Long non-coding RNAs (lncRNAs) play a significant role in maintaining tissue morphology and functions, and their precise regulatory effectiveness is closely related to expression patterns. However, the spatial expression patterns of lncRNAs in humans are poorly characterized. Here, we constructed five comprehensive transcriptomic atlases of human lncRNAs covering thousands of major tissue samples in normal and disease states. The lncRNA transcriptomes exhibited high consistency within the same tissues across resources, and even higher complexity in specialized tissues. Tissue-elevated (TE) lncRNAs were identified in each resource and robust TE lncRNAs were refined by integrative analysis. We detected 1 to 4684 robust TE lncRNAs across tissues; the highest number was in testis tissue, followed by brain tissue. Functional analyses of TE lncRNAs indicated important roles in corresponding tissue-related pathways. Moreover, we found that the expression features of robust TE lncRNAs made them be effective biomarkers to distinguish tissues; TE lncRNAs also tended to be associated with cancer, and exhibited differential expression or were correlated with patient survival. In summary, spatial classification of lncRNAs is the starting point for elucidating the function of lncRNAs in both maintenance of tissue morphology and progress of tissue-constricted diseases.

Methylome-wide association study of different responses to risperidone in schizophrenia

Background: Accumulating evidence shows that DNA methylation plays a role in antipsychotic response. However, the mechanisms by which DNA methylation changes are associated with antipsychotic responses remain largely unknown. Methods: We performed a methylome-wide association study (MWAS) to evaluate the association between DNA methylation and the response to risperidone in schizophrenia. Genomic DNA methylation patterns were assessed using the Agilent Human DNA Methylation Microarray. Results: We identified numerous differentially methylated positions (DMPs) and regions (DMRs) associated with antipsychotic response. CYP46A1, SPATS2, and ATP6V1E1 had the most significant DMPs, with p values of 2.50 × 10^-6, 3.53 × 10^-6, and 5.71 × 10^-6, respectively. The top-ranked DMR was located on chromosome 7, corresponding to the PTPRN2 gene with a Šidák-corrected p-value of 9.04 × 10^-13. Additionally, a significant enrichment of synaptic function and neurotransmitters was found in the differentially methylated genes after gene ontology and pathway analysis. Conclusion: The identified DMP- and DMR-overlapping genes associated with antipsychotic response are related to synaptic function and neurotransmitters. These findings may improve understanding of the mechanisms underlying antipsychotic response and guide the choice of antipsychotic in schizophrenia.

Patient-derived parathyroid organoids as a tracer and drug-screening application model

Parathyroid diseases are characterized by dysregulation of calcium homeostasis and alterations in parathyroid hormone (PTH) excretion. The development of parathyroid-targeted treatment and imaging tracers could benefit from in vitro models. Therefore, we aim to establish a patient-derived parathyroid organoid model representing human parathyroid tissue. Hyperplastic parathyroid tissue was dispersed, and parathyroid organoids (PTOs) were cultured and characterized. PTO-derived cells exhibited self-renewal over several passages, indicative of the presence of putative stem cells. Immunofluorescence and RNA sequencing confirmed that PTOs phenocopy hyperplastic parathyroid tissue. Exposure of PTOs to increasing calcium concentrations and PTH-lowering drugs resulted in significantly reduced PTH excretion. PTOs showed specific binding of the imaging tracers ¹¹C-methionine and ^99mTc-sestamibi. These data show the functionality of PTOs resembling the parathyroid. This PTO model recapitulates the originating tissue on gene and protein expression and functionality, paving the way for future physiology studies and therapeutic target and tracer discovery.

Roles of dual specificity tyrosine-phosphorylation-regulated kinase 2 in nervous system development and disease

Dual specificity tyrosine-phosphorylation-regulated kinases (DYRKs) are a group of conserved eukaryotic kinases phosphorylating tyrosine, serine, and threonine residues. The human DYRK family comprises 5 members (DYRK1A, DYRK1B, DYRK2, DYRK3, and DYRK4). The different DYRKs have been implicated in neurological diseases, cancer, and virus infection. Specifically, DYRK2 has been mainly implicated in cancer progression. However, its role in healthy and pathological nervous system function has been overlooked. In this context, we review current available data on DYRK2 in the nervous system, where the available studies indicate that it has key roles in neuronal development and function. DYRK2 regulates neuronal morphogenesis (e.g., axon growth and branching) by phosphorylating cytoskeletal elements (e.g., doublecortin). Comparative data reveals that it is involved in the development of olfactory and visual systems, the spinal cord and possibly the cortex. DYRK2 also participates in processes such as olfaction, vision and, learning. However, DYRK2 could be involved in other brain functions since available expression data shows that it is expressed across the whole brain. High DYRK2 protein levels have been detected in basal ganglia and cerebellum. In adult nervous system, DYRK2 mRNA expression is highest in the cortex, hippocampus, and retina. Regarding nervous system disease, DYRK2 has been implicated in neuroblastoma, glioma, epilepsy, neuroinflammation, Alzheimer's disease, Parkinson's disease, spinal cord injury and virus infection. DYRK2 upregulation usually has a negative impact in cancer-related conditions and a positive impact in non-malignant conditions. Its role in axon growth makes DYRK2 as a promising target for spinal cord or brain injury and regeneration.

Protein Expression Profile of ACE2 in the Normal and COVID-19-Affected Human Brain

SARS-coronavirus 2 (SARS-CoV-2) that caused the coronavirus disease 2019 (COVID-19) pandemic has posed to be a global challenge. An increasing number of neurological symptoms have been linked to the COVID-19 disease, but the underlying mechanisms of such symptoms and which patients could be at risk are not yet established. The suggested key receptor for host cell entry is angiotensin I converting enzyme 2 (ACE2). Previous studies on limited tissue material have shown no or low protein expression of ACE2 in the normal brain. Here, we used stringently validated antibodies and immunohistochemistry to examine the protein expression of ACE2 in all major regions of the normal brain. The expression pattern was compared with the COVID-19-affected brain of patients with a varying degree of neurological symptoms. In the normal brain, the expression was restricted to the choroid plexus and ependymal cells with no expression in any other brain cell types. Interestingly, in the COVID-19-affected brain, an upregulation of ACE2 was observed in endothelial cells of certain patients, most prominently in the white matter and with the highest expression observed in the patient with the most severe neurological symptoms. The data shows differential expression of ACE2 in the diseased brain and highlights the need to further study the role of endothelial cells in COVID-19 disease in relation to neurological symptoms.

Advances in application of single-cell RNA sequencing in cardiovascular research

Single-cell RNA sequencing (scRNA-seq) provides high-resolution information on transcriptomic changes at the single-cell level, which is of great significance for distinguishing cell subtypes, identifying stem cell differentiation processes, and identifying targets for disease treatment. In recent years, emerging single-cell RNA sequencing technologies have been used to make breakthroughs regarding decoding developmental trajectories, phenotypic transitions, and cellular interactions in the cardiovascular system, providing new insights into cardiovascular disease. This paper reviews the technical processes of single-cell RNA sequencing and the latest progress based on single-cell RNA sequencing in the field of cardiovascular system research, compares single-cell RNA sequencing with other single-cell technologies, and summarizes the extended applications and advantages and disadvantages of single-cell RNA sequencing. Finally, the prospects for applying single-cell RNA sequencing in the field of cardiovascular research are discussed.

A systematic evaluation of the computational tools for ligand-receptor-based cell-cell interaction inference

Cell-cell interactions (CCIs) are essential for multicellular organisms to coordinate biological processes and functions. One classical type of CCI interaction is between secreted ligands and cell surface receptors, i.e. ligand-receptor (LR) interactions. With the recent development of single-cell technologies, a large amount of single-cell ribonucleic acid (RNA) sequencing (scRNA-Seq) data has become widely available. This data availability motivated the single-cell-resolution study of CCIs, particularly LR-based CCIs. Dozens of computational methods and tools have been developed to predict CCIs by identifying LR-based CCIs. Many of these tools have been theoretically reviewed. However, there is little study on current LR-based CCI prediction tools regarding their performance and running results on public scRNA-Seq datasets. In this work, to fill this gap, we tested and compared nine of the most recent computational tools for LR-based CCI prediction. We used 15 well-studied scRNA-Seq samples that correspond to approximately 100K single cells under different experimental conditions for testing and comparison. Besides briefing the methodology used in these nine tools, we summarized the similarities and differences of these tools in terms of both LR prediction and CCI inference between cell types. We provided insight into using these tools to make meaningful discoveries in understanding cell communications.

Plasmolipin regulates basolateral-to-apical transcytosis of ICAM-1 and leukocyte adhesion in polarized hepatic epithelial cells

Apical localization of Intercellular Adhesion Receptor (ICAM)-1 regulates the adhesion and guidance of leukocytes across polarized epithelial barriers. Here, we investigate the molecular mechanisms that determine ICAM-1 localization into apical membrane domains of polarized hepatic epithelial cells, and their effect on lymphocyte-hepatic epithelial cell interaction. We had previously shown that segregation of ICAM-1 into apical membrane domains, which form bile canaliculi and bile ducts in hepatic epithelial cells, requires basolateral-to-apical transcytosis. Searching for protein machinery potentially involved in ICAM-1 polarization we found that the SNARE-associated protein plasmolipin (PLLP) is expressed in the subapical compartment of hepatic epithelial cells in vitro and in vivo. BioID analysis of ICAM-1 revealed proximal interaction between this adhesion receptor and PLLP. ICAM-1 colocalized and interacted with PLLP during the transcytosis of the receptor. PLLP gene editing and silencing increased the basolateral localization and reduced the apical confinement of ICAM-1 without affecting apicobasal polarity of hepatic epithelial cells, indicating that ICAM-1 transcytosis is specifically impaired in the absence of PLLP. Importantly, PLLP depletion was sufficient to increase T-cell adhesion to hepatic epithelial cells. Such an increase depended on the epithelial cell polarity and ICAM-1 expression, showing that the epithelial transcytotic machinery regulates the adhesion of lymphocytes to polarized epithelial cells. Our findings strongly suggest that the polarized intracellular transport of adhesion receptors constitutes a new regulatory layer of the epithelial inflammatory response.

A network-based approach reveals the dysregulated transcriptional regulation in non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is a leading cause of chronic liver disease worldwide. We performed network analysis to investigate the dysregulated biological processes in the disease progression and revealed the molecular mechanism underlying NAFLD. Based on network analysis, we identified a highly conserved disease-associated gene module across three different NAFLD cohorts and highlighted the predominant role of key transcriptional regulators associated with lipid and cholesterol metabolism. In addition, we revealed the detailed metabolic differences between heterogeneous NAFLD patients through integrative systems analysis of transcriptomic data and liver-specific genome-scale metabolic model. Furthermore, we identified transcription factors (TFs), including SREBF2, HNF4A, SREBF1, YY1, and KLF13, showing regulation of hepatic expression of genes in the NAFLD-associated modules and validated the TFs using data generated from a mouse NAFLD model. In conclusion, our integrative analysis facilitates the understanding of the regulatory mechanism of these perturbed TFs and their associated biological processes.

•

Disease-associated gene modules are conserved across multiple NAFLD cohorts

•

The central genes in disease-associated modules are key enzymes in cholesterol synthesis

•

YY1 and KLF13 are potential key transcriptional regulators of NAFLD development

Elimination of acquired resistance to PD-1 blockade via the concurrent depletion of tumour cells and immunosuppressive cells

Antigen release resulting from the death of tumour cells induced by chemotherapies and targeted therapies can augment the antitumour responses induced by immune checkpoint blockade (ICB). However, tumours responding to ICB therapies often become resistant to them. Here we show that the specific targeting of tumour cells promotes the growth of tumour-cell variants that are resistant to ICB, and that the acquired resistance can be overcome via the concurrent depletion of tumour cells and of major types of immunosuppressive cell via a monoclonal antibody binding the enzyme CD73, which we identified as highly expressed on tumour cells and on regulatory T cells, myeloid-derived suppressor cells and tumour-associated macrophages, but not on cytolytic T lymphocytes, natural killer cells and dendritic cells. In mice with murine tumours, the systemic administration of anti-PD1 antibodies and anti-CD73 antibodies conjugated to a near-infrared dye prevented near-infrared-irradiated tumours from acquiring resistance to ICB and resulted in the eradication of advanced tumours. The elimination of immunosuppressive cells may overcome acquired resistance to ICB across a range of tumour types and combination therapies.

Genome-wide analysis of DNA methylation in 106 schizophrenia family trios in Han Chinese

Background

Schizophrenia (SCZ) is a severe psychiatric disorder that affects approximately 0.75% of the global population. Both genetic and environmental factors contribute to development of SCZ. SCZ tends to run in family while both genetic and environmental factor contribute to its etiology. Much evidence suggested that alterations in DNA methylations occurred in SCZ patients.

Methods

To investigate potential inheritable pattern of DNA methylation in SCZ family, we performed a genome-wide analysis of DNA methylation of peripheral blood samples from 106 Chinese SCZ family trios. Genome-wide DNA methylations were quantified by Agilent 1 × 244 k Human Methylation Microarray.

Findings

In this study, we proposed a loci inheritance frequency model that allows characterization of differential methylated regions as SCZ biomarkers. Based on this model, 112 hypermethylated and 125 hypomethylated regions were identified. Additionally, 121 hypermethylated and 139 hypomethylated genes were annotated. The results of functional enrichment analysis indicated that multiple differentially methylated genes (DMGs) involved in Notch/HH/Wnt signaling, MAPK signaling, GPCR signaling, immune response signaling. Notably, a number of hypomethylated genes were significantly enriched in cerebral cortex and functionally enriched in nervous system development.

Interpretation

Our findings not only validated previously discovered risk genes of SCZ but also identified novel candidate DMGs in SCZ. These results may further the understanding of altered DNA methylations in SCZ.

An information-theoretic approach for measuring the distance of organ tissue samples using their transcriptomic signatures

Motivation

Recapitulating aspects of human organ functions using in vitro (e.g. plates, transwells, etc.), in vivo (e.g. mouse, rat, etc.), or ex vivo (e.g. organ chips, 3D systems, etc.) organ models is of paramount importance for drug discovery and precision medicine. It will allow us to identify potential side effects and test the effectiveness of new therapeutic approaches early in their design phase, and will inform the development of better disease models. Developing mathematical methods to reliably compare the ‘distance/similarity’ of organ models from/to the real human organ they represent is an understudied problem with important applications in biomedicine and tissue engineering.

Results

We introduce the Transcriptomic Signature Distance (TSD), an information-theoretic distance for assessing the transcriptomic similarity of two tissue samples, or two groups of tissue samples. In developing TSD, we are leveraging next-generation sequencing data as well as information retrieved from well-curated databases providing signature gene sets characteristic for human organs. We present the justification and mathematical development of the new distance and demonstrate its effectiveness and advantages in different scenarios of practical importance using several publicly available RNA-seq datasets.

Availability and Implementation

The computation of both TSD versions (simple and weighted) has been implemented in R and can be downloaded from https://github.com/Cod3B3nd3R/Transcriptomic-Signature-Distance.

Contact

dimitris.manatakis@emulatebio.com

Supplementary information

Supplementary data are available at Bioinformatics online.

Circulating trophoblast cell clusters for early detection of placenta accreta spectrum disorders

Placenta accreta spectrum (PAS) is a high-risk obstetrical condition associated with significant morbidity and mortality. Current clinical screening modalities for PAS are not always conclusive. Here, we report a nanostructure-embedded microchip that efficiently enriches both single and clustered circulating trophoblasts (cTBs) from maternal blood for detecting PAS. We discover a uniquely high prevalence of cTB-clusters in PAS and subsequently optimize the device to preserve the intactness of these clusters. Our feasibility study on the enumeration of cTBs and cTB-clusters from 168 pregnant women demonstrates excellent diagnostic performance for distinguishing PAS from non-PAS. A logistic regression model is constructed using a training cohort and then cross-validated and tested using an independent cohort. The combined cTB assay achieves an Area Under ROC Curve of 0.942 (throughout gestation) and 0.924 (early gestation) for distinguishing PAS from non-PAS. Our assay holds the potential to improve current diagnostic modalities for the early detection of PAS.

Missense mutations in small muscle protein X-linked (SMPX) cause distal myopathy with protein inclusions

Using deep phenotyping and high-throughput sequencing, we have identified a novel type of distal myopathy caused by mutations in the Small muscle protein X-linked (SMPX) gene. Four different missense mutations were identified in ten patients from nine families in five different countries, suggesting that this disease could be prevalent in other populations as well. Haplotype analysis of patients with similar ancestry revealed two different founder mutations in Southern Europe and France, indicating that the prevalence in these populations may be higher. In our study all patients presented with highly similar clinical features: adult-onset, usually distal more than proximal limb muscle weakness, slowly progressing over decades with preserved walking. Lower limb muscle imaging showed a characteristic pattern of muscle involvement and fatty degeneration. Histopathological and electron microscopic analysis of patient muscle biopsies revealed myopathic findings with rimmed vacuoles and the presence of sarcoplasmic inclusions, some with amyloid-like characteristics. In silico predictions and subsequent cell culture studies showed that the missense mutations increase aggregation propensity of the SMPX protein. In cell culture studies, overexpressed SMPX localized to stress granules and slowed down their clearance.

The CR9 element is a novel mechanical load-responsive enhancer that regulates natriuretic peptide genes expression

Enhancers regulate gene expressions in a tissue- and pathology-specific manner by altering its activities. Plasma levels of atrial and brain natriuretic peptides, encoded by the Nppa and Nppb, respectively, and synthesized predominantly in cardiomyocytes, vary depending on the severity of heart failure. We previously identified the noncoding conserved region 9 (CR9) element as a putative Nppb enhancer at 22-kb upstream from the Nppb gene. However, its regulatory mechanism remains unknown. Here, we therefore investigated the mechanism of CR9 activation in cardiomyocytes using different kinds of drugs that induce either cardiac hypertrophy or cardiac failure accompanied by natriuretic peptides upregulation. Chronic treatment of mice with either catecholamines or doxorubicin increased CR9 activity during the progression of cardiac hypertrophy to failure, which is accompanied by proportional increases in Nppb expression. Conversely, for cultured cardiomyocytes, doxorubicin decreased CR9 activity and Nppb expression, while catecholamines increased both. However, exposing cultured cardiomyocytes to mechanical loads, such as mechanical stretch or hydrostatic pressure, upregulate CR9 activity and Nppb expression even in the presence of doxorubicin. Furthermore, the enhancement of CR9 activity and Nppa and Nppb expressions by either catecholamines or mechanical loads can be blunted by suppressing mechanosensing and mechanotransduction pathways, such as muscle LIM protein (MLP) or myosin tension. Finally, the CR9 element showed a more robust and cell-specific response to mechanical loads than the -520-bp BNP promoter. We concluded that the CR9 element is a novel enhancer that responds to mechanical loads by upregulating natriuretic peptides expression in cardiomyocytes.

A Chinese hamster transcription start site atlas that enables targeted editing of CHO cells

Chinese hamster ovary (CHO) cells are widely used for producing biopharmaceuticals, and engineering gene expression in CHO is key to improving drug quality and affordability. However, engineering gene expression or activating silent genes requires accurate annotation of the underlying regulatory elements and transcription start sites (TSSs). Unfortunately, most TSSs in the published Chinese hamster genome sequence were computationally predicted and are frequently inaccurate. Here, we use nascent transcription start site sequencing methods to revise TSS annotations for 15 308 Chinese hamster genes and 3034 non-coding RNAs based on experimental data from CHO-K1 cells and 10 hamster tissues. We further capture tens of thousands of putative transcribed enhancer regions with this method. Our revised TSSs improves upon the RefSeq annotation by revealing core sequence features of gene regulation such as the TATA box and the Initiator and, as exemplified by targeting the glycosyltransferase gene Mgat3, facilitate activating silent genes by CRISPRa. Together, we envision our revised annotation and data will provide a rich resource for the CHO community, improve genome engineering efforts and aid comparative and evolutionary studies.

Conservation and diversity of the eukaryotic SAGA coactivator complex across kingdoms

The SAGA complex is an evolutionarily conserved transcriptional coactivator that regulates gene expression through its histone acetyltransferase and deubiquitylase activities, recognition of specific histone modifications, and interactions with transcription factors. Multiple lines of evidence indicate the existence of distinct variants of SAGA among organisms as well as within a species, permitting diverse functions to dynamically regulate cellular pathways. Our co-expression analysis of genes encoding human SAGA components showed enrichment in reproductive organs, brain tissues and the skeletal muscle, which corresponds to their established roles in developmental programs, emerging roles in neurodegenerative diseases, and understudied functions in specific cell types. SAGA subunits modulate growth, development and response to various stresses from yeast to plants and metazoans. In metazoans, SAGA further participates in the regulation of differentiation and maturation of both innate and adaptive immune cells, and is associated with initiation and progression of diseases including a broad range of cancers. The evolutionary conservation of SAGA highlights its indispensable role in eukaryotic life, thus deciphering the mechanisms of action of SAGA is key to understanding fundamental biological processes throughout evolution. To illuminate the diversity and conservation of this essential complex, here we discuss variations in composition, essentiality and co-expression of component genes, and its prominent functions across Fungi, Plantae and Animalia kingdoms.

Expression of Cancer Testis Antigens in Tumor-Adjacent Normal Liver Is Associated with Post-Resection Recurrence of Hepatocellular Carcinoma

Simple Summary

High recurrence rates after resection of liver cancer (hepatocellular carcinoma) with curative intent impair clinical outcomes of patients diagnosed with liver cancer. Cancer/testis antigens (CTAs) are expressed in cancer and can serve as therapeutic targets. We identified 12 CTAs expressed in 80% of liver cancer patients, and each one individually in at least 10%. Furthermore, we found that patients with expression of CTAs in macroscopically tumor-free liver tissue, experience more tumor recurrence and poor survival after surgical tumor removal. The increased risk of tumor recurrence in patients with CTA expression in tumor-free liver suggests that these patients already have micro-metastasis at the time of operation. These CTA-expressing (pre-)malignant cells may thus be a source of liver cancer recurrence, reflecting the relevance of targeting these to prevent liver cancer recurrence.

Abstract

High recurrence rates after resection of hepatocellular carcinoma (HCC) with curative intent impair clinical outcomes of HCC. Cancer/testis antigens (CTAs) are suitable targets for cancer immunotherapy if selectively expressed in tumor cells. The aims were to identify CTAs that are frequently and selectively expressed in HCC-tumors, and to investigate whether CTAs could serve as biomarkers for occult metastasis. Tumor and paired tumor-free liver (TFL) tissues of HCC-patients and healthy tissues were assessed for mRNA expression of 49 CTAs by RT-qPCR and protein expression of five CTAs by immunohistochemistry. Twelve CTA-mRNAs were expressed in ≥10% of HCC-tumors and not in healthy tissues except testis. In tumors, mRNA and protein of ≥ 1 CTA was expressed in 78% and 71% of HCC-patients, respectively. In TFL, CTA mRNA and protein was found in 45% and 30% of HCC-patients, respectively. Interestingly, CTA-expression in TFL was an independent negative prognostic factor for post-resection HCC-recurrence and survival. We established a panel of 12 testis-restricted CTAs expressed in tumors of most HCC-patients. The increased risk of HCC-recurrence in patients with CTA expression in TFL, suggests that CTA-expressing (pre-)malignant cells may be a source of HCC-recurrence, reflecting the relevance of targeting these to prevent HCC-recurrence.

Western Diet Induced Remodelling of the Tongue Proteome

The tongue is a heavily innervated and vascularized striated muscle that plays an important role in vocalization, swallowing and digestion. The surface of the tongue is lined with papillae which contain gustatory cells expressing various taste receptors. There is growing evidence to suggest that our perceptions of taste and food preference are remodelled following chronic consumption of Western diets rich in carbohydrate and fats. Our sensitivity to taste and also to metabolising Western diets may be a key factor in the rising prevalence of obesity; however, a systems-wide analysis of the tongue is lacking. Here, we defined the proteomic landscape of the mouse tongue and quantified changes following chronic consumption of a chow or Western diet enriched in lipid, fructose and cholesterol for 7 months. We observed a dramatic remodelling of the tongue proteome including proteins that regulate fatty acid and mitochondrial metabolism. Furthermore, the expressions of several receptors, metabolic enzymes and hormones were differentially regulated, and are likely to provide novel therapeutic targets to alter taste perception and food preference to combat obesity.

Enhancers are activated by p300/CBP activity-dependent PIC assembly, RNAPII recruitment, and pause release

The metazoan-specific acetyltransferase p300/CBP is involved in activating signal-induced, enhancer-mediated transcription of cell-type-specific genes. However, the global kinetics and mechanisms of p300/CBP activity-dependent transcription activation remain poorly understood. We performed genome-wide, time-resolved analyses to show that enhancers and super-enhancers are dynamically activated through p300/CBP-catalyzed acetylation, deactivated by the opposing deacetylase activity, and kinetic acetylation directly contributes to maintaining cell identity at very rapid (minutes) timescales. The acetyltransferase activity is dispensable for the recruitment of p300/CBP and transcription factors but essential for promoting the recruitment of TFIID and RNAPII at virtually all enhancers and enhancer-regulated genes. This identifies pre-initiation complex assembly as a dynamically controlled step in the transcription cycle and reveals p300/CBP-catalyzed acetylation as the signal that specifically promotes transcription initiation at enhancer-regulated genes. We propose that p300/CBP activity uses a "recruit-and-release" mechanism to simultaneously promote RNAPII recruitment and pause release and thereby enables kinetic activation of enhancer-mediated transcription.

Interspecies chimeric conditions affect the developmental rate of human pluripotent stem cells

Human pluripotent stem cells hold significant promise for regenerative medicine. However, long differentiation protocols and immature characteristics of stem cell-derived cell types remain challenges to the development of many therapeutic applications. In contrast to the slow differentiation of human stem cells in vitro that mirrors a nine-month gestation period, mouse stem cells develop according to a much faster three-week gestation timeline. Here, we tested if co-differentiation with mouse pluripotent stem cells could accelerate the differentiation speed of human embryonic stem cells. Following a six-week RNA-sequencing time course of neural differentiation, we identified 929 human genes that were upregulated earlier and 535 genes that exhibited earlier peaked expression profiles in chimeric cell cultures than in human cell cultures alone. Genes with accelerated upregulation were significantly enriched in Gene Ontology terms associated with neurogenesis, neuron differentiation and maturation, and synapse signaling. Moreover, chimeric mixed samples correlated with in utero human embryonic samples earlier than human cells alone, and acceleration was dose-dependent on human-mouse co-culture ratios. The altered gene expression patterns and developmental rates described in this report have implications for accelerating human stem cell differentiation and the use of interspecies chimeric embryos in developing human organs for transplantation.

Longitudinal saliva omics responses to immune perturbation: a case study

Saliva omics has immense potential for non-invasive diagnostics, including monitoring very young or elderly populations, or individuals in remote locations. In this study, multiple saliva omics from an individual were monitored over three periods (100 timepoints) involving: (1) hourly sampling over 24 h without intervention, (2) hourly sampling over 24 h including immune system activation using the standard 23-valent pneumococcal polysaccharide vaccine, (3) daily sampling for 33 days profiling the post-vaccination response. At each timepoint total saliva transcriptome and proteome, and small RNA from salivary extracellular vesicles were profiled, including mRNA, miRNA, piRNA and bacterial RNA. The two 24-h periods were used in a paired analysis to remove daily variation and reveal vaccination responses. Over 18,000 omics longitudinal series had statistically significant temporal trends compared to a healthy baseline. Various immune response and regulation pathways were activated following vaccination, including interferon and cytokine signaling, and MHC antigen presentation. Immune response timeframes were concordant with innate and adaptive immunity development, and coincided with vaccination and reported fever. Overall, mRNA results appeared more specific and sensitive (timewise) to vaccination compared to other omics. The results suggest saliva omics can be consistently assessed for non-invasive personalized monitoring and immune response diagnostics.

Processing and Analysis of RNA-seq Data from Public Resources

Advances in next generation sequencing (NGS) technologies resulted in a broad array of large-scale gene expression studies and an unprecedented volume of whole messenger RNA (mRNA) sequencing data, or the transcriptome (also known as RNA sequencing, or RNA-seq). These include the Genotype Tissue Expression project (GTEx) and The Cancer Genome Atlas (TCGA), among others. Here we cover some of the commonly used datasets, provide an overview on how to begin the analysis pipeline, and how to explore and interpret the data provided by these publicly available resources.

Breast Cancer Risk with Progestin Subdermal Implants: A Challenge in Patients Counseling

There is a steady global rise in the use of progestin subdermal implants, where use has increased by more than 20 times in the past two decades. BC risk has been reported with the older progestin only methods such as oral pills, injectables, and intrauterine devices, however, little is known about the risk with subdermal implants. In this review, we aim to update clinicians and researchers on the current evidence to support patient counseling and to inform future research directions. The available evidence of the association between the use of progestin subdermal implants and BC risk is discussed. We provide an overview of the potential role of endogenous progesterone in BC development. The chemical structure and molecular targets of synthetic progestins of relevance are summarized together with the preclinical and clinical evidence on their association with BC risk. We review all studies that investigated the action of the specific progestins included in subdermal implants. As well, we discuss the potential effect of the use of subdermal implants in women at increased BC risk, including carriers of BC susceptibility genetic mutations.

The Human Protein Atlas-Spatial localization of the human proteome in health and disease

For a complete understanding of a system's processes and each protein's role in health and disease, it is essential to study protein expression with a spatial resolution, as the exact location of proteins at tissue, cellular, or subcellular levels is tightly linked to protein function. The Human Protein Atlas (HPA) project is a large-scale initiative aiming at mapping the entire human proteome using antibody-based proteomics and integration of various other omics technologies. The publicly available knowledge resource www.proteinatlas.org is one of the world's most visited biological databases and has been extensively updated during the last few years. The current version is divided into six main sections, each focusing on particular aspects of the human proteome: (a) the Tissue Atlas showing the distribution of proteins across all major tissues and organs in the human body; (b) the Cell Atlas showing the subcellular localization of proteins in single cells; (c) the Pathology Atlas showing the impact of protein levels on survival of patients with cancer; (d) the Blood Atlas showing the expression profiles of blood cells and actively secreted proteins; (e) the Brain Atlas showing the distribution of proteins in human, mouse, and pig brain; and (f) the Metabolic Atlas showing the involvement of proteins in human metabolism. The HPA constitutes an important resource for further understanding of human biology, and the publicly available datasets hold much promise for integration with other emerging efforts focusing on single cell analyses, both at transcriptomic and proteomic level.

Pathophysiology and pharmacological management of pulmonary and cardiovascular features of COVID-19

The first confirmed case of novel Coronavirus Disease 2019 (COVID-19) in the United States was reported on January 20, 2020. As of November 24, 2020, close to 12.2 million cases of COVID-19 was confirmed in the US, with over 255,958 deaths. The rapid transmission of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), its unusual and divergent presentation has strengthened the status of COVID-19 as a major public health threat. In this review, we aim to 1- discuss the epidemiological data from various COVID-19 patient cohorts around the world and the USA as well the associated risk factors; 2- summarize the pathophysiology of SARS-CoV-2 infection and the underlying molecular mechanisms for the respiratory and cardiovascular manifestations; 3- highlight the potential treatments and vaccines as well as current clinical trials for COVID-19.

Proteomic signatures of 16 major types of human cancer reveal universal and cancer-type-specific proteins for the identification of potential therapeutic targets

Background

Proteomic characterization of cancers is essential for a comprehensive understanding of key molecular aberrations. However, proteomic profiling of a large cohort of cancer tissues is often limited by the conventional approaches.

Methods

We present a proteomic landscape of 16 major types of human cancer, based on the analysis of 126 treatment-naïve primary tumor tissues, 94 tumor-matched normal adjacent tissues, and 12 normal tissues, using mass spectrometry-based data-independent acquisition approach.

Results

In our study, a total of 8527 proteins were mapped to brain, head and neck, breast, lung (both small cell and non-small cell lung cancers), esophagus, stomach, pancreas, liver, colon, kidney, bladder, prostate, uterus and ovary cancers, including 2458 tissue-enriched proteins. Our DIA-based proteomic approach has characterized major human cancers and identified universally expressed proteins as well as tissue-type-specific and cancer-type-specific proteins. In addition, 1139 therapeutic targetable proteins and 21 cancer/testis (CT) antigens were observed.

Conclusions

Our discoveries not only advance our understanding of human cancers, but also have implications for the design of future large-scale cancer proteomic studies to assist the development of diagnostic and/or therapeutic targets in multiple cancers.

Intestinal region-specific Wnt signalling profiles reveal interrelation between cell identity and oncogenic pathway activity in cancer development

BACKGROUND

Cancer results from the accumulation of mutations leading to the acquisition of cancer promoting characteristics such as increased proliferation and resistance to cell death. In colorectal cancer, an early mutation leading to such features usually occurs in the APC or CTNNB1 genes, thereby activating Wnt signalling. However, substantial phenotypic differences between cancers originating within the same organ, such as molecular subtypes, are not fully reflected by differences in mutations. Indeed, the phenotype seems to result from a complex interplay between the cell-intrinsic features and the acquired mutations, which is difficult to disentangle when established tumours are studied.

METHODS

We use a 3D in vitro organoid model to study the early phase of colorectal cancer development. From three different murine intestinal locations we grow organoids. These are transformed to resemble adenomas after Wnt activation through lentiviral transduction with a stable form of β-Catenin. The gene expression before and after Wnt activation is compared within each intestinal origin and across the three locations using RNA sequencing. To validate and generalize our findings, we use gene expression data from patients.

RESULTS

In reaction to Wnt activation we observe downregulation of location specific genes and differentiation markers. A similar effect is seen in patient data, where genes with significant differential expression between the normal left and right colon are downregulated in the cancer samples. Furthermore, the signature of Wnt target genes differs between the three intestinal locations in the organoids. The location specific Wnt signatures are dominated by genes which have been lowly expressed in the tissue of origin, and are the targets of transcription factors that are activated following enhanced Wnt signalling.

CONCLUSION

We observed that the region-specific cell identity has a substantial effect on the reaction to Wnt activation in a simple intestinal adenoma model. These findings provide a way forward in resolving the distinct biology between left- and right-sided human colon cancers with potential clinical relevance.

ACE2: Evidence of role as entry receptor for SARS-CoV-2 and implications in comorbidities

Pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus 19 disease (COVID-19) which presents a large spectrum of manifestations with fatal outcomes in vulnerable people over 70-years-old and with hypertension, diabetes, obesity, cardiovascular disease, COPD, and smoking status. Knowledge of the entry receptor is key to understand SARS-CoV-2 tropism, transmission and pathogenesis. Early evidence pointed to angiotensin-converting enzyme 2 (ACE2) as SARS-CoV-2 entry receptor. Here, we provide a critical summary of the current knowledge highlighting the limitations and remaining gaps that need to be addressed to fully characterize ACE2 function in SARS-CoV-2 infection and associated pathogenesis. We also discuss ACE2 expression and potential role in the context of comorbidities associated with poor COVID-19 outcomes. Finally, we discuss the potential co-receptors/attachment factors such as neuropilins, heparan sulfate and sialic acids and the putative alternative receptors, such as CD147 and GRP78.

Molecular Portraits of Early Rheumatoid Arthritis Identify Clinical and Treatment Response Phenotypes

There is a current imperative to unravel the hierarchy of molecular pathways that drive the transition of early to established disease in rheumatoid arthritis (RA). Herein, we report a comprehensive RNA sequencing analysis of the molecular pathways that drive early RA progression in the disease tissue (synovium), comparing matched peripheral blood RNA-seq in a large cohort of early treatment-naive patients, namely, the Pathobiology of Early Arthritis Cohort (PEAC). We developed a data exploration website (https://peac.hpc.qmul.ac.uk/) to dissect gene signatures across synovial and blood compartments, integrated with deep phenotypic profiling. We identified transcriptional subgroups in synovium linked to three distinct pathotypes: fibroblastic pauci-immune pathotype, macrophage-rich diffuse-myeloid pathotype, and a lympho-myeloid pathotype characterized by infiltration of lymphocytes and myeloid cells. This is suggestive of divergent pathogenic pathways or activation disease states. Pro-myeloid inflammatory synovial gene signatures correlated with clinical response to initial drug therapy, whereas plasma cell genes identified a poor prognosis subgroup with progressive structural damage.

•

Deep phenotyping and RNA-seq of early rheumatoid arthritis individuals pre-treatment

•

Synovial plasma cell gene expression predicts future progressive joint damage on X-ray

•

Blood interferon gene signature associates with synovial B and plasma cell infiltration

•

Interactive website enables RNA-seq and clinical data to be fully explored

Elevated mRNA Expression Levels of NCAPG are Associated with Poor Prognosis in Ovarian Cancer

BACKGROUND

Ovarian cancer is a major gynecologic malignancy that is often detected at a late stage due to the lack of detailed studies on its pathogenesis and reliable biomarkers for predicting its prognosis.

MATERIALS AND METHODS

Four ovarian cancer data sets GSE18520, GSE27651, GSE40595, and GSE52037 were downloaded from the Gene Expression Omnibus (GEO) database and the robust rank aggregation approach was used to find common differentially expressed genes (DEGs). Cytoscape software was used to construct and detect key models of protein-protein interaction (PPI) network. While the expression, prognostic value and potential mechanism of the hub gene non-SMC condensin I complex subunit G (NCAPG) was carried out through Gene Expression Profiling Interactive Analysis, Kaplan-Meier plotter online dataset and gene set enrichment analysis. To further investigate the role of NCAPG in ovarian cancer, in vitro experiments were carried out.

RESULTS

A total of 232 DEGs were identified in the four GEO datasets; and we detected 32 hub genes from the PPI network and 21 of these genes were associated with ovarian cancer prognosis, one of which was NCAPG. NCAPG was significantly upregulated in most of the ovarian cancer samples. High NCAPG expression was mainly involved in homologous recombination, DNA replication, proteasome, and more correlated pathways. NCAPG knockdown arrested the cell cycle, inhibited the proliferation, and attenuated the migration ability of A2780 cells. Meanwhile, silencing of NCAPG significantly promoted cisplatin-induced apoptosis thus increased the sensitivity to cisplatin.

CONCLUSION

NCAPG together with the other 31 hub genes play a vital role in the tumorigenesis of ovarian, meanwhile, the cell cycle pathway may be a potential pathway contributing to progression in OC; and NCAPG expression can be used as a promising target for the treatment of OC.

The protein expression profile of ACE2 in human tissues

The novel SARS-coronavirus 2 (SARS-CoV-2) poses a global challenge on healthcare and society. For understanding the susceptibility for SARS-CoV-2 infection, the cell type-specific expression of the host cell surface receptor is necessary. The key protein suggested to be involved in host cell entry is angiotensin I converting enzyme 2 (ACE2). Here, we report the expression pattern of ACE2 across > 150 different cell types corresponding to all major human tissues and organs based on stringent immunohistochemical analysis. The results were compared with several datasets both on the mRNA and protein level. ACE2 expression was mainly observed in enterocytes, renal tubules, gallbladder, cardiomyocytes, male reproductive cells, placental trophoblasts, ductal cells, eye, and vasculature. In the respiratory system, the expression was limited, with no or only low expression in a subset of cells in a few individuals, observed by one antibody only. Our data constitute an important resource for further studies on SARS-CoV-2 host cell entry, in order to understand the biology of the disease and to aid in the development of effective treatments to the viral infection.