Exploring Transcriptomic Landscapes in Red Blood Cells, in Their Extracellular Vesicles and on a Single-Cell Level.

Being enucleated, RBCs lack typical transcriptomes, but are known to contain small amounts of diverse long transcripts and microRNAs. However, the exact role and importance of these RNAs are lacking. Shedding of extracellular vesicles (EVs) from the plasma membrane constitutes an integral mechanism of RBC homeostasis, by which RBCs remove unnecessary cytoplasmic content and cell membrane. To study this further, we explored the transcriptomes of RBCs and extracellular vesicles (EVs) of RBCs using next-generation sequencing. Furthermore, we performed single-cell RNA sequencing on RBCs, which revealed that approximately 10% of the cells contained detectable levels of mRNA and cells formed three subpopulations based on their transcriptomes. A decrease in the mRNA quantity was observed across the populations. Qualitative changes included the differences in the globin transcripts and changes in the expression of ribosomal genes. A specific splice form of a long non-coding RNA, Metastasis Associated Lung Adenocarcinoma Transcript 1 (MALAT1), was the most enriched marker in one subpopulation of RBCs, co-expressing with ribosomal structural transcripts. MALAT1 expression was confirmed by qPCR in CD71-enriched reticulocytes, which were also characterized with imaging flow cytometry at the single cell level. Analysis of the RBC transcriptome shows enrichment of pathways and functional categories required for the maturation of reticulocytes and erythrocyte functions. The RBC transcriptome was detected in their EVs, making these transcripts available for intercellular communication in blood.

FinDonor 10 000 study: a cohort to identify iron depletion and factors affecting it in Finnish blood donors.

BACKGROUND AND OBJECTIVES: There is increasing evidence that frequent blood donation depletes the iron stores of some blood donors. The FinDonor 10 000 study was set up to study iron status and factors affecting iron stores in Finnish blood donors. In Finland, iron supplementation for at-risk groups has been in place since the 1980s. MATERIAL AND METHODS: A total of 2584 blood donors (N = 8003 samples) were recruited into the study alongside standard donation at three donation sites in the capital region of Finland between 5/2015 and 12/2017. All participants were asked to fill out a questionnaire about their health and lifestyle. Blood samples were collected from the sample pouch of whole blood collection set, kept in cool temperature and processed centrally. Whole blood count, CRP, ferritin and sTFR were measured from the samples, and DNA was isolated for GWAS studies. RESULTS: Participant demographics, albeit in general similar to the general blood donor population in Finland, indicated some bias towards older and more frequent donors. Participation in the study increased median donation frequency of the donors. Analysis of the effect of time lag from the sampling to the analysis and the time of day when sample was drawn revealed small but significant time-dependent changes. CONCLUSION: The FinDonor cohort now provides us with tools to identify potential donor groups at increased risk of iron deficiency and factors explaining this risk. The increase in donation frequency during the study suggests that scientific projects can be used to increase the commitment of blood donors.

The impact of RNA sequence library construction protocols on transcriptomic profiling of leukemia.

BACKGROUND: RNA sequencing (RNA-seq) has become an indispensable tool to identify disease associated transcriptional profiles and determine the molecular underpinnings of diseases. However, the broad adaptation of the methodology into the clinic is still hampered by inconsistent results from different RNA-seq protocols and involves further evaluation of its analytical reliability using patient samples. Here, we applied two commonly used RNA-seq library preparation protocols to samples from acute leukemia patients to understand how poly-A-tailed mRNA selection (PA) and ribo-depletion (RD) based RNA-seq library preparation protocols affect gene fusion detection, variant calling, and gene expression profiling. RESULTS: Overall, the protocols produced similar results with consistent outcomes. Nevertheless, the PA protocol was more efficient in quantifying expression of leukemia marker genes and showed better performance in the expression-based classification of leukemia. Independent qRT-PCR experiments verified that the PA protocol better represented total RNA compared to the RD protocol. In contrast, the RD protocol detected a higher number of non-coding RNA features and had better alignment efficiency. The RD protocol also recovered more known fusion-gene events, although variability was seen in fusion gene predictions. CONCLUSION: The overall findings provide a framework for the use of RNA-seq in a precision medicine setting with limited number of samples and suggest that selection of the library preparation protocol should be based on the objectives of the analysis.

Somatic STAT3 mutations in large granular lymphocytic leukemia.

BACKGROUND: T-cell large granular lymphocytic leukemia is a rare lymphoproliferative disorder characterized by the expansion of clonal CD3+CD8+ cytotoxic T lymphocytes (CTLs) and often associated with autoimmune disorders and immune-mediated cytopenias. METHODS: We used next-generation exome sequencing to identify somatic mutations in CTLs from an index patient with large granular lymphocytic leukemia. Targeted resequencing was performed in a well-characterized cohort of 76 patients with this disorder, characterized by clonal T-cell-receptor rearrangements and increased numbers of large granular lymphocytes. RESULTS: Mutations in the signal transducer and activator of transcription 3 gene (STAT3) were found in 31 of 77 patients (40%) with large granular lymphocytic leukemia. Among these 31 patients, recurrent mutational hot spots included Y640F in 13 (17%), D661V in 7 (9%), D661Y in 7 (9%), and N647I in 3 (4%). All mutations were located in exon 21, encoding the Src homology 2 (SH2) domain, which mediates the dimerization and activation of STAT protein. The amino acid changes resulted in a more hydrophobic protein surface and were associated with phosphorylation of STAT3 and its localization in the nucleus. In vitro functional studies showed that the Y640F and D661V mutations increased the transcriptional activity of STAT3. In the affected patients, downstream target genes of the STAT3 pathway (IFNGR2, BCL2L1, and JAK2) were up-regulated. Patients with STAT3 mutations presented more often with neutropenia and rheumatoid arthritis than did patients without these mutations. CONCLUSIONS: The SH2 dimerization and activation domain of STAT3 is frequently mutated in patients with large granular lymphocytic leukemia; these findings suggest that aberrant STAT3 signaling underlies the pathogenesis of this disease. (Funded by the Academy of Finland and others.).

Molecular mechanisms of nasal epithelium in rhinitis and rhinosinusitis.

Allergic rhinitis, nonallergic rhinitis, and chronic rhinosinusitis are multifactorial upper airway diseases with high prevalence. Several genetic and environmental factors are proposed to predispose to the pathogenesis of the inflammatory upper airway diseases. Still, the molecular mechanisms leading toward the onset and progression of upper airway diseases are largely unknown. The upper airway epithelium has an important role in sensing the environment and regulating the inhaled air. As such, it links environmental insults to the host immunity. Human sinonasal epithelium serves as an excellent target for observing induced early-phase events, in vivo, and with a systems biological perspective. Actually, increasing number of investigations have provided evidence that altered homeostasis in the sinonasal epithelium might be important in the chronic upper airway inflammation.

Neurocan expression associates with better survival and viral positivity in Merkel cell carcinoma.

Merkel cell carcinoma (MCC) is a rare cutaneous neuroendocrine carcinoma that is frequently divided into Merkel cell polyomavirus negative and positive tumors due their distinct genomic and transcriptomic profiles, and disease outcomes. Although some prognostic factors in MCC are known, tumorigenic pathways, which that explain outcome differences in MCC are not fully understood. We investigated transcriptomes of 110 tissue samples of a formalin-fixed, paraffin-embedded MCC series by RNA sequencing to identify genes showing a bimodal expression pattern and predicting outcome in cancer and that potentially could play a role in tumorigenesis. We discovered 19 genes among which IGHM, IGKC, NCAN, OTOF, and USH2A were associated also with overall survival (all p-values < 0.05). From these genes, NCAN (neurocan) expression was detected in all 144 MCC samples by immunohistochemistry. Increased NCAN expression was associated with presence of Merkel cell polyomavirus DNA (p = 0.001) and viral large T antigen expression in tumor tissue (p = 0.004) and with improved MCC-specific survival (p = 0.027) and overall survival (p = 0.034). We conclude that NCAN expression is common in MCC, and further studies are warranted to investigate its role in MCC tumorigenesis.

Quartet RNA reference materials improve the quality of transcriptomic data through ratio-based profiling.

Certified RNA reference materials are indispensable for assessing the reliability of RNA sequencing to detect intrinsically small biological differences in clinical settings, such as molecular subtyping of diseases. As part of the Quartet Project for quality control and data integration of multi-omics profiling, we established four RNA reference materials derived from immortalized B-lymphoblastoid cell lines from four members of a monozygotic twin family. Additionally, we constructed ratio-based transcriptome-wide reference datasets between two samples, providing cross-platform and cross-laboratory 'ground truth'. Investigation of the intrinsically subtle biological differences among the Quartet samples enables sensitive assessment of cross-batch integration of transcriptomic measurements at the ratio level. The Quartet RNA reference materials, combined with the ratio-based reference datasets, can serve as unique resources for assessing and improving the quality of transcriptomic data in clinical and biological settings.

Analysis of human brain tissue derived from DBS surgery.

BACKGROUND: Transcriptomic and proteomic profiling of human brain tissue is hindered by the availability of fresh samples from living patients. Postmortem samples usually represent the advanced disease stage of the patient. Furthermore, the postmortem interval can affect the transcriptomic and proteomic profiles. Therefore, fresh brain tissue samples from living patients represent a valuable resource of metabolically intact tissue. Implantation of deep brain stimulation (DBS) electrodes into the human brain is a neurosurgical treatment for, e.g., movement disorders. Here, we describe an improved approach to collecting brain tissues from surgical instruments used in implantation of DBS device for transcriptomics and proteomics analyses. METHODS: Samples were extracted from guide tubes and recording electrodes used in routine DBS implantation procedure to treat patients with Parkinson's disease, genetic dystonia and tremor. RNA sequencing was performed in tissues extracted from the recording microelectrodes and liquid chromatography-mass spectrometry (LC-MS) performed in tissues from guide tubes. To assess the performance of the current approach, the obtained datasets were compared with previously published datasets representing brain tissues. RESULTS: Altogether, 32,034 RNA transcripts representing the unique Ensembl gene identifiers were detected from eight samples representing both hemispheres of four patients. By using  LC-MS, we identified 734 unique proteins from 31 samples collected from 14 patients. The datasets are available in the BioStudies database (accession number S-BSST667). Our results indicate that surgical instruments used in DBS installation retain brain material sufficient for protein and gene expression studies. Comparison with previously published datasets obtained with similar approach proved the robustness and reproducibility of the protocol. CONCLUSIONS: The instruments used during routine DBS surgery are a useful source for obtaining fresh brain tissues from living patients. This approach overcomes the issues that arise from using postmortem tissues, such as the effect of postmortem interval on transcriptomic and proteomic landscape of the brain, and can be used for studying molecular aspects of DBS-treatable diseases.

A novel variant in SMG9 causes intellectual disability, confirming a role for nonsense-mediated decay components in neurocognitive development.

Biallelic loss-of-function variants in the SMG9 gene, encoding a regulatory subunit of the mRNA nonsense-mediated decay (NMD) machinery, are reported to cause heart and brain malformation syndrome. Here we report five patients from three unrelated families with intellectual disability (ID) and a novel pathogenic SMG9 c.551 T > C p.(Val184Ala) homozygous missense variant, identified using exome sequencing. Sanger sequencing confirmed recessive segregation in each family. SMG9 c.551T > C p.(Val184Ala) is most likely an autozygous variant identical by descent. Characteristic clinical findings in patients were mild to moderate ID, intention tremor, pyramidal signs, dyspraxia, and ocular manifestations. We used RNA sequencing of patients and age- and sex-matched healthy controls to assess the effect of the variant. RNA sequencing revealed that the SMG9 c.551T > C variant did not affect the splicing or expression level of SMG9 gene products, and allele-specific expression analysis did not provide evidence that the nonsense mRNA-induced NMD was affected. Differential gene expression analysis identified prevalent upregulation of genes in patients, including the genes SMOX, OSBP2, GPX3, and ZNF155. These findings suggest that normal SMG9 function may be involved in transcriptional regulation without affecting nonsense mRNA-induced NMD. In conclusion, we demonstrate that the SMG9 c.551T > C missense variant causes a neurodevelopmental disorder and impacts gene expression. NMD components have roles beyond aberrant mRNA degradation that are crucial for neurocognitive development.

Patient-tailored design for selective co-inhibition of leukemic cell subpopulations.

The extensive drug resistance requires rational approaches to design personalized combinatorial treatments that exploit patient-specific therapeutic vulnerabilities to selectively target disease-driving cell subpopulations. To solve the combinatorial explosion challenge, we implemented an effective machine learning approach that prioritizes patient-customized drug combinations with a desired synergy-efficacy-toxicity balance by combining single-cell RNA sequencing with ex vivo single-agent testing in scarce patient-derived primary cells. When applied to two diagnostic and two refractory acute myeloid leukemia (AML) patient cases, each with a different genetic background, we accurately predicted patient-specific combinations that not only resulted in synergistic cancer cell co-inhibition but also were capable of targeting specific AML cell subpopulations that emerge in differing stages of disease pathogenesis or treatment regimens. Our functional precision oncology approach provides an unbiased means for systematic identification of personalized combinatorial regimens that selectively co-inhibit leukemic cells while avoiding inhibition of nonmalignant cells, thereby increasing their likelihood for clinical translation.

Caveolar transport through nasal epithelium of birch pollen allergen Bet v 1 in allergic patients.

BACKGROUND: Previous work in type I pollen allergies has focused on aberrant immunoresponses. OBJECTIVE: Our systems-level analyses explore the role of epithelium in early pathogenesis of type I allergic reactions. METHODS: We began top-down analyses of differences in human nasal epithelial cells and biopsy specimens obtained from patients with birch allergy and healthy control subjects in the resting state and after intranasal in vivo birch pollen challenges. Immunohistochemistry, immunotransmission electron microscopy, mass spectrometry, transcriptomics, and integration of data to a pathway were conducted. RESULTS: Bet v 1 allergen bound to epithelium immediately after in vivo birch pollen challenge during winter only in allergic individuals. It also travelled through epithelium with caveolae to mast cells. Sixteen unique proteins were found to bind to the Bet v 1 column only in lysates from allergic epithelial cells; 6 of these were caveolar and 6 were cytoskeletal proteins. The nasal epithelial transcriptome analysis from allergic and healthy subjects differed during the winter season, and these subjects also responded differentially to birch pollen challenge. Within this pollen-induced response, the gene ontology categories of cytoskeleton and actin cytoskeleton were decreased in allergic patients, whereas the actin-binding category was enriched in healthy subjects. Integration of microscopic, mass spectrometric, and transcriptomic data to a common protein-protein binding network showed how these were connected to each other. CONCLUSION: We propose a hypothesis of caveolae-dependent uptake and transport of birch pollen allergen in the epithelium of allergic patients only. Application of discovery-driven methodologies can provide new hypotheses worth further analysis of complex multifactorial diseases, such as type I allergy.

Genomics of asthma, allergy and chronic rhinosinusitis: novel concepts and relevance in airway mucosa.

Genome wide association studies (GWASs) have revealed several airway disease-associated risk loci. Their role in the onset of asthma, allergic rhinitis (AR) or chronic rhinosinusitis (CRS), however, is not yet fully understood. The aim of this review is to evaluate the airway relevance of loci and genes identified in GWAS studies. GWASs were searched from databases, and a list of loci associating significantly (p < 10-8) with asthma, AR and CRS was created. This yielded a total of 267 significantly asthma/AR-associated loci from 31 GWASs. No significant CRS -associated loci were found in this search. A total of 170 protein coding genes were connected to these loci. Of these, 76/170 (44%) showed bronchial epithelial protein expression in stained microscopic figures of Human Protein Atlas (HPA), and 61/170 (36%) had a literature report of having airway epithelial function. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analyses were performed, and 19 functional protein categories were found as significantly (p < 0.05) enriched among these genes. These were related to cytokine production, cell activation and adaptive immune response, and all were strongly connected in network analysis. We also identified 15 protein pathways that were significantly (p < 0.05) enriched in these genes, related to T-helper cell differentiation, virus infection, JAK-STAT signaling pathway, and asthma. A third of GWAS-level risk loci genes of asthma or AR seemed to have airway epithelial functions according to our database and literature searches. In addition, many of the risk loci genes were immunity related. Some risk loci genes also related to metabolism, neuro-musculoskeletal or other functions. Functions overlapped and formed a strong network in our pathway analyses and are worth future studies of biomarker and therapeutics.

Immunogenomic Landscape of Hematological Malignancies.

Understanding factors that shape the immune landscape across hematological malignancies is essential for immunotherapy development. We integrated over 8,000 transcriptomes and 2,000 samples with multilevel genomics of hematological cancers to investigate how immunological features are linked to cancer subtypes, genetic and epigenetic alterations, and patient survival, and validated key findings experimentally. Infiltration of cytotoxic lymphocytes was associated with TP53 and myelodysplasia-related changes in acute myeloid leukemia, and activated B cell-like phenotype and interferon-γ response in lymphoma. CIITA methylation regulating antigen presentation, cancer type-specific immune checkpoints, such as VISTA in myeloid malignancies, and variation in cancer antigen expression further contributed to immune heterogeneity and predicted survival. Our study provides a resource linking immunology with cancer subtypes and genomics in hematological malignancies.

Concordant gene regulation related to perturbations of three GDP-mannose-related genes.

Glycosylation of proteins is one of the most crucial post-translational modifications. In order to access system-level and state-dependent data related to the regulation of glycosylation events, we cultivated yeast cell strains each harboring a selected conditional knockdown construct for a gene (either SEC53, VRG4 or DPM1) related to GDP-mannose synthesis or its utilization in glycan biosynthesis. In order to carry this out efficiently, we developed automated sampling from bioreactor cultivations, a collection of in silico workflows for data analysis as well as their integration into a large data warehouse. Using the above-mentioned approaches, we could show that conditional knocking down of transcripts related to GDP-mannose synthesis or transportation led to altered levels of over 300 transcripts. These transcripts and their corresponding proteins were characterized by their gene ontology (GO) annotations, and their putative transcriptional regulation was analyzed. Furthermore, novel pathways were generated indicating interactions between GO categories with common proteins, putative transcriptional regulators of such induced GO categories, and the large protein-protein interaction network among the proteins whose transcripts indicated altered expression levels. When these results are always added to an ever-expanding data warehouse as annotations, they will incrementally increase the knowledge of biological systems.

Absence of NEFL in patient-specific neurons in early-onset Charcot-Marie-Tooth neuropathy.

OBJECTIVE: We used patient-specific neuronal cultures to characterize the molecular genetic mechanism of recessive nonsense mutations in neurofilament light (NEFL) underlying early-onset Charcot-Marie-Tooth (CMT) disease. METHODS: Motor neurons were differentiated from induced pluripotent stem cells of a patient with early-onset CMT carrying a novel homozygous nonsense mutation in NEFL. Quantitative PCR, protein analytics, immunocytochemistry, electron microscopy, and single-cell transcriptomics were used to investigate patient and control neurons. RESULTS: We show that the recessive nonsense mutation causes a nearly total loss of NEFL messenger RNA (mRNA), leading to the complete absence of NEFL protein in patient's cultured neurons. Yet the cultured neurons were able to differentiate and form neuronal networks and neurofilaments. Single-neuron gene expression fingerprinting pinpointed NEFL as the most downregulated gene in the patient neurons and provided data of intermediate filament transcript abundancy and dynamics in cultured neurons. Blocking of nonsense-mediated decay partially rescued the loss of NEFL mRNA. CONCLUSIONS: The strict neuronal specificity of neurofilament has hindered the mechanistic studies of recessive NEFL nonsense mutations. Here, we show that such mutation leads to the absence of NEFL, causing childhood-onset neuropathy through a loss-of-function mechanism. We propose that the neurofilament accumulation, a common feature of many neurodegenerative diseases, mimics the absence of NEFL seen in recessive CMT if aggregation prevents the proper localization of wild-type NEFL in neurons. Our results suggest that the removal of NEFL as a proposed treatment option is harmful in humans.

Genetic Control of Myelin Plasticity after Chronic Psychosocial Stress.

Anxiety disorders often manifest in genetically susceptible individuals after psychosocial stress, but the mechanisms underlying these gene-environment interactions are largely unknown. We used the chronic social defeat stress (CSDS) mouse model to study resilience and susceptibility to chronic psychosocial stress. We identified a strong genetic background effect in CSDS-induced social avoidance (SA) using four inbred mouse strains: 69% of C57BL/6NCrl (B6), 23% of BALB/cAnNCrl, 19% of 129S2/SvPasCrl, and 5% of DBA/2NCrl (D2) mice were stress resilient. Furthermore, different inbred mouse strains responded differently to stress, suggesting they use distinct coping strategies. To identify biological pathways affected by CSDS, we used RNA-sequencing (RNA-seq) of three brain regions of two strains, B6 and D2: medial prefrontal cortex (mPFC), ventral hippocampus (vHPC), and bed nucleus of the stria terminalis (BNST). We discovered overrepresentation of oligodendrocyte (OLG)-related genes in the differentially expressed gene population. Because OLGs myelinate axons, we measured myelin thickness and found significant region and strain-specific differences. For example, in resilient D2 mice, mPFC axons had thinner myelin than controls, whereas susceptible B6 mice had thinner myelin than controls in the vHPC. Neither myelin-related gene expression in several other regions nor corpus callosum thickness differed between stressed and control animals. Our unbiased gene expression experiment suggests that myelin plasticity is a substantial response to chronic psychosocial stress, varies across brain regions, and is genetically controlled. Identification of genetic regulators of the myelin response will provide mechanistic insight into the molecular basis of stress-related diseases, such as anxiety disorders, a critical step in developing targeted therapy.

Retrospective analysis of capillary hemoglobin recovery in nearly 1 200 000 blood donor returns.

Measuring the concentration of capillary hemoglobin (cHb) is a standard procedure before blood donation. To further assess the time period needed for cHb recovery after blood donation and to have a more in-depth understanding of features of recovery, we used data-mining tools in a large, retrospective data pool containing all 1 163 524 donor returns that took place in Finland in 2010 to 2015. The results show that the average recovery times for cHb to return back to the level preceding donation were substantially longer, over 200 days in all age groups, than were the minimum allowed donation intervals. cHb recovery was especially poor in women under the age of 30 who returned to donate soon after the minimum allowed donation interval. It was of interest that frequent donors recovered substantially faster, with the average recovery times of ∼100 days in men and ∼200 days in women, than did infrequent donors, suggesting that there is a subpopulation of donors who can donate frequently without fear of iron deficiency. Return interval in fact explained only 1% of the variation in cHb recovery, which points to unknown, individual features, such as genetic or lifestyle factors, warranting further studies and suggesting that simply extending the allowed donation intervals may not suffice to improve cHb recovery. The study demonstrates that data mining of blood bank records is a powerful tool for depicting features of blood donor population.