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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

RNA sequencing of chorionic villi from recurrent pregnancy loss patients reveals impaired function of basic nuclear and cellular machinery.

Recurrent pregnancy loss (RPL) concerns ~3% of couples aiming at childbirth. In the current study, transcriptomes and miRNomes of 1st trimester placental chorionic villi were analysed for 2 RPL cases (≥6 miscarriages) and normal, but electively terminated pregnancies (ETP; n = 8). Sequencing was performed on Illumina HiSeq 2000 platform. Differential expression analyses detected 51 (27%) transcripts with increased and 138 (73%) with decreased expression in RPL compared to ETP (DESeq: FDR P < 0.1 and DESeq2: <0.05). RPL samples had substantially decreased transcript levels of histones, regulatory RNAs and genes involved in telomere, spliceosome, ribosomal, mitochondrial and intra-cellular signalling functions. Downregulated expression of HIST1H1B and HIST1H4A (Wilcoxon test, fc≤0.372, P≤9.37 × 10-4) was validated in an extended sample by quantitative PCR (RPL, n = 14; ETP, n = 24). Several upregulated genes are linked to placental function and pregnancy complications: ATF4, C3, PHLDA2, GPX4, ICAM1, SLC16A2. Analysis of the miRNA-Seq dataset identified no large disturbances in RPL samples. Notably, nearly 2/3 of differentially expressed genes have binding sites for E2F transcription factors, coordinating mammalian endocycle and placental development. For a conceptus destined to miscarriage, the E2F TF-family represents a potential key coordinator in reprogramming the placental genome towards gradually stopping the maintenance of basic nuclear and cellular functions.

Extensive shift in placental transcriptome profile in preeclampsia and placental origin of adverse pregnancy outcomes.

One in five pregnant women suffer from gestational complications, prevalently driven by placental malfunction. Using RNASeq, we analyzed differential placental gene expression in cases of normal gestation, late-onset preeclampsia (LO-PE), gestational diabetes (GD) and pregnancies ending with the birth of small-for-gestational-age (SGA) or large-for-gestational-age (LGA) newborns (n = 8/group). In all groups, the highest expression was detected for small noncoding RNAs and genes specifically implicated in placental function and hormonal regulation. The transcriptome of LO-PE placentas was clearly distinct, showing statistically significant (after FDR) expressional disturbances for hundreds of genes. Taqman RT-qPCR validation of 45 genes in an extended sample (n = 24/group) provided concordant results. A limited number of transcription factors including LRF, SP1 and AP2 were identified as possible drivers of these changes. Notable differences were detected in differential expression signatures of LO-PE subtypes defined by the presence or absence of intrauterine growth restriction (IUGR). LO-PE with IUGR showed higher correlation with SGA and LO-PE without IUGR with LGA placentas. Whereas changes in placental transcriptome in SGA, LGA and GD cases were less prominent, the overall profiles of expressional disturbances overlapped among pregnancy complications providing support to shared placental responses. The dataset represent a rich catalogue for potential biomarkers and therapeutic targets.

Expression of Toll-like receptors in nasal epithelium in allergic rhinitis.

Toll-like receptors (TLRs) are important in barrier homeostasis, but their role in airborne allergies is not fully understood. The aim was to evaluate baseline and allergen-induced expression of TLR proteins in nasal epithelium during allergic rhinitis. Nineteen otherwise healthy non-smoking volunteers both allergic to birch pollen and non-allergic controls were enrolled. We took nasal biopsies before and after off-seasonal intranasal birch pollen or diluent challenge. The expression of epithelial TLR1-7, TLR9-10, and MyD88 proteins was immunohistochemically evaluated from the nasal biopsies. The TLR1-3 and TLR5-10 mRNAs were observed by RNA-microarray. Baseline epithelial expression of TLR proteins was wide and identical in controls and atopics. After off-seasonal intranasal birch pollen challenge, a negative change in the expression score of TLR1 and TLR6 proteins was detected in the atopic group. TLR mRNA expression was not affected by birch pollen challenge. Nasal epithelium seems to express all known TLRs. The mechanisms by which TLR1, and TLR6 proteins could affect pollen allergen transport need further studies.