Next Article in Journal
Neuroinflammation in Dementia—Therapeutic Directions in a COVID-19 Pandemic Setting
Next Article in Special Issue
Peripheral Immunophenotype in IgG4-Related Disease and Its Association with Clinical Phenotypes and Disease Activity
Previous Article in Journal
CORK1, A LRR-Malectin Receptor Kinase, Is Required for Cellooligomer-Induced Responses in Arabidopsis thaliana
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Timing of Blood Sample Processing Affects the Transcriptomic and Epigenomic Profiles in CD4+ T-cells of Atopic Subjects

1
Translational Inflammation Research Division & Core Facility for Single Cell Multiomics, Medical Faculty, Philipps University of Marburg, Member of the German Center for Lung Research (DZL) and the Universities of Giessen and Marburg Lung Center, D-35043 Marburg, Germany
2
Department of Medicine, Pulmonary and Critical Care Medicine, Brigham and Women’s Hospital, Harvard Medical School, Boston, MA 02115, USA
3
Pediatric Allergology, Department of Pediatrics, Dr. von Hauner Children’s Hospital, University Hospital, Ludwig Maximilian University (LMU) of Munich, Comprehensive Pneumology Center Munich (CPC-M), Member of the German Center for Lung Research (DZL), D-80337 Munich, Germany
4
Department of Pediatric Pneumology, Allergology and Neonatology, Hannover Medical School, Biomedical Research in Endstage and Obstructive Lung Disease Hannover (BREATH), Member of the German Center for Lung Research (DZL) and the Cluster of Excellence RESIST (EXC 2155), D-30625 Hannover, Germany
5
Institute of Medical Bioinformatics and Biostatistics, Medical Faculty, Philipps University of Marburg, D-35032 Marburg, Germany
6
Bioscientia MVZ Labor Mittelhessen GmbH, D-35394 Gießen, Germany
*
Author to whom correspondence should be addressed.
Cells 2022, 11(19), 2958; https://doi.org/10.3390/cells11192958
Submission received: 4 July 2022 / Revised: 10 August 2022 / Accepted: 9 September 2022 / Published: 22 September 2022
(This article belongs to the Special Issue T Cell Responses in Human Health and Disease - Second Edition)

Abstract

:
Optimal pre-analytical conditions for blood sample processing and isolation of selected cell populations for subsequent transcriptomic and epigenomic studies are required to obtain robust and reproducible results. This pilot study was conducted to investigate the potential effects of timing of CD4+ T-cell processing from peripheral blood of atopic and non-atopic adults on their transcriptomic and epigenetic profiles. Two heparinized blood samples were drawn from each of three atopic and three healthy individuals. For each individual, CD4+ T-cells were isolated from the first blood sample within 2 h (immediate) or from the second blood sample after 24 h storage (delayed). RNA sequencing (RNA-Seq) and histone H3K27 acetylation chromatin immunoprecipitation sequencing (ChIP-Seq) analyses were performed. A multiplicity of genes was shown to be differentially expressed in immediately processed CD4+ T-cells from atopic versus healthy subjects. These differences disappeared when comparing delayed processed cells due to a drastic change in expression levels of atopy-related genes in delayed processed CD4+ T-cells from atopic donors. This finding was further validated on the epigenomic level by examining H3K27 acetylation profiles. In contrast, transcriptomic and epigenomic profiles of blood CD4+ T-cells of healthy donors remained rather unaffected. Taken together, for successful transcriptomics and epigenomics studies, detailed standard operation procedures developed on the basis of samples from both healthy and disease conditions are implicitly recommended.

1. Introduction

In recent years there has been a growing interest in diverse “omics” analyses utilizing different high-throughput techniques such as RNA sequencing (RNA-Seq) and chromatin immunoprecipitation sequencing (ChIP-Seq) [1,2]. These techniques have been employed to develop a better understanding of underlying molecular mechanisms and the pathobiology of a wealth of disease conditions including atopy [3,4,5]. The latter involves an increased risk of sensitization associated with the generation of inflammatory immune responses to harmless substances, based on a genetic predisposition to develop allergic diseases such as allergic rhinitis and asthma and/or mostly modern lifestyle-associated environmental influences [6,7,8].
Peripheral blood provides a minimally invasive liquid biopsy that comprises a vital source of molecular information and thus serves as an easily accessible diagnostic tool [9]. CD4+ T-cells in the blood mediate different adaptive immune responses to environmental structures including pathogens and allergens. Central to their function is the activation of other immune and non-immune cells by secreting many types of cytokines [10].
Various ways of handling peripheral blood samples during collection, shipping and processing prior to downstream analyses have a significant impact on the robustness and reproducibility of the results [11,12]. Processing of preselected cells at a given time point following blood collection in biobank sampling for big cohort studies represents logistical challenges that might overwhelm the subsequent biological associations [13].
In this pilot study, we intended to better define optimal experimental conditions for subsequent larger cohort studies by investigating potential impacts of timing of CD4+ T-cell processing from peripheral blood of atopic and non-atopic (healthy) adults on their transcriptomic (RNA-Seq) profiles. Moreover, we validated the RNA-Seq results on the epigenomic level by examining histone H3K27 acetylation (H3K27ac) profiles, known to be associated with active gene transcription, by ChIP-Seq [14]. 

2. Materials and Methods

2.1. Study Population and CD4+ T-Cell Sorting

In this pilot study, we recruited three atopic and three non-atopic (healthy) adult donors from two research centers of the German Center for Lung Research (DZL): two atopic and one healthy subject from the Department of Pediatric Pneumology, Allergology and Neonatology, Hannover and two healthy and one atopic subject from the Pediatric Allergology, Department of Pediatrics, Munich. All atopic subjects were diagnosed with sensitization to house dust mite, with either a history of hayfever (2 patients) or allergic rhinitis (1 patient). The patients were symptom-free and did not receive anti-inflammatory treatment at the time of analysis. Two heparinized blood samples were drawn from each of the six study subjects. The first blood sample was processed and CD4+ T-cells were FACS-sorted within two hours upon blood withdrawal (immediate processing), while the second blood sample was kept at room temperature (according to common protocols in biobank sampling) for 24 h prior to CD4+ T-cell isolation (delayed processing), followed by performing RNA-Seq and histone H3K27 acetylation ChIP-Seq applications and conducting bioinformatic analysis (Figure 1). In detail, heparinized venous blood was collected from each of the three atopic and healthy volunteers. Of each sample, 10 mL were immediately processed and 10 mL were left at room temperature for 24 h for delayed processing. At the two time points, peripheral blood mononuclear cells (PBMCs) were isolated by Biocoll (Biochrom, Berlin, Germany) density-gradient centrifugation and stained with fluorescent-labeled antibodies: CD8-FITC (RPA-T8), CD4-APC-H7 (RPA-T4) (Becton Dickinson, Heidelberg, Germany), CD19-PerCP-Cy5.5 (HIB19) (Biolegend, San Diego, CA, USA) and CD3-PE-Cy7 (UCHT1) (eBioscience, San Diego, CA, USA). CD3+CD4+CD8CD19 T-cells were sorted by fluorescent-activated cell sorting (FACS) using a FACSAria cell sorter (Becton Dickinson). The underlying gating strategy is shown in Figure S1. CD4+ T-cells were not cryopreserved at any stage of the analytical process. About 400,000 CD4+ T-cells were pelleted and re-suspended in RLT-Buffer (Qiagen, Hilden, Germany) including 1% β-Mercaptoethanol (Carl Roth, Karlsruhe, Germany) for transcriptome analysis and immediately transferred to −80 °C for storage. For ChIP-Seq analysis, about 600,000 CD4+ T-cells were washed and fixed in 1% formaldehyde solution (Merck, Darmstadt, Germany) followed by an addition of 1.3 M Glycine (Carl Roth) to crosslink DNA-associated proteins. Cells were washed once with ice-cold Dulbecco’s phosphate-buffered saline (pH 7.2) containing 1x protease inhibitor mix (Diagenode, Seraing, Belgium) and pelleted. After supernatant aspiration, pellets were immediately snap-frozen in liquid nitrogen and stored at −80 °C.

2.2. RNA Extraction, Library Preparation, and RNA Sequencing (RNA-Seq) Analysis

Total RNA was extracted from CD4+ T-cells using miRNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol. RNA quality was assessed by Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). Library preparation was conducted using QuantSeq 3’ mRNA-Seq Library Prep Kit REV for Illumina (Lexogen, Greenland, NH, USA) according to the manufacturer’s protocol. Sequencing was performed on an Illumina NextSeq550 platform, High Output Kit v2.5, 50 bases single-reads (Illumina, San Diego, CA, USA) according to the manufacturer’s instructions.

2.3. Chromatin Immunoprecipitation Sequencing (ChIP-Seq) Analysis

ChIP-Seq analysis was performed using the True MicroChIP kit (Diagenode). Briefly, 100,000 fixed CD4+ T-cells were lysed and sheared using the Bioruptor Pico sonication device (Diagenode) for 15 cycles using a 30” [ON] 30” [OFF] setting. A total of 10% of the sheared chromatin was kept aside as an input control. For immunoprecipitation, sheared chromatin was incubated with 1 µg of an anti-H3K27ac polyclonal rabbit antibody. For quality check purposes, sheared chromatin was also incubated with 1 µg of unspecific rabbit IgG as isotype control (both antibodies from Diagenode). Samples were incubated with Protein A-coated magnetic beads and washed 4 times. Afterwards, the immunoprecipitated chromatin was de-crosslinked and DNA was purified. Quality check was performed for all samples by quantitative polymerase chain reaction (qPCR) using primer pairs (Diagenode) identifying a positive region, i.e., house-keeping gene GAPDH transcription starting site (TSS) control region, and myoglobin exon 2 as a negative control region using the Rotor-Gene SYBR Green PCR Kit (Qiagen). Afterwards, library preparation for H3K27ac and input control samples was performed using the MicroPlex Library Preparation Kit v3 with dual indexes (Diagenode) according to the manufacturer’s protocol. Sequencing was performed as described above.

2.4. Bioinformatic Analysis

RNA-Seq and ChIP-Seq bioinformatics workflows were performed using the Galaxy platform [15]. For RNA-Seq, quality of reads was assessed by FastQC v0.72 and reads were aligned to the human reference genome (Homo sapiens hg38 full) with RNA STAR v2.7.8a. Mapped reads were counted with featureCounts v2.0.1, followed by exon DESeq2 v2.11.40.6 analysis to determine differentially expressed genes. Enrichr BioPlanet and Gene ontologies databases were utilized to perform functional pathway enrichment analyses and investigate biological processes, respectively [16,17,18]. Gene expression values in Heatmaps were calculated as z-scores. ChIP-Seq reads were aligned to the same human reference genome using Bowtie2 v2.4.2. deepTools was utilized to generate enriched score genomic heatmaps and overlapped genes profiles. Differential whole gene analysis was conducted by DESeq2. MACS2 callpeak v2.2.7.1 and DiffBind v2.10.0 were used to perform peak calling and differential binding analyses, respectively. 

3. Results

3.1. Comparison of Transcriptome Profiles of Immediate Versus Delayed Processed CD4+ T-Cells from Atopic and Healthy Subjects

When comparing transcriptomic profiles of immediately processed CD4+ T-cells of atopic versus healthy subjects, a total of 1303 genes were identified to be differentially expressed at a false discovery rate (FDR) < 0.1 between these two groups (Figure 2A). We observed a multiplicity of genes involved in biological pathways known to play a role in atopy-associated T-cell activation processes (Figure S2A–C). Surprisingly, when performing the same analysis in delayed processed CD4+ T-cells, only three genes (RPS27, RPL41, and RPL13), all coding for ribosomal proteins, remained differentially expressed (Figure 2B). We then intended to explore the changes in the transcriptomes of the samples from the two time points of the same physiological condition. The above-described difference observed between delayed and immediate cell processing was predominantly caused by transcriptomic changes in CD4+ T-cells from atopic subjects since 320 genes were significantly differentially expressed when comparing immediately and delayed processed blood samples from atopic donors (Figure 2C). In contrast, the same comparison yielded only 17 differentially expressed genes in cells from healthy subjects (Figure 2D). Of note, expression of genes involved in deadenylation and degradation of mRNAs such as CCR4-NOT Transcription Complex Subunit 6 Like (CNOT6L) and CCR4-NOT Transcription Complex Subunit 1 (CNOT1) were found to be significantly up-regulated in delayed cell preparations from atopic subjects. In accordance, only pathways involved in the regulation of RNA stability were found to be upregulated in cells after delayed processing while a variety of atopy-related T-cell activation pathways (e.g., T helper cell surface molecules, T-cell receptor and IL-2 signaling pathways) were significantly down-regulated upon delayed processing (Figure 3A–C).

3.2. Epigenetic Profiles at the Level of Histone H3K27 Acetylation of Immediate Versus Delayed Processed CD4+ T-Cells from Atopic and Healthy Subjects

To verify the previous findings on the epigenetic level, we carried out ChIP-Seq analysis of H3K27 acetylation, an epigenetic mark strongly correlating with active gene transcription [14]. When comparing delayed versus immediate cell processing, acetylation levels of H3K27 flanking around the TSSs were remarkably confined in the delayed processed CD4+ T-cells of atopic subjects (Figure 4A). In contrast, acetylation of H3K27 remained at comparable levels between the two processing time points in healthy donors (Figure 4B). To comprehensively demonstrate the different changes in the H3K27ac profiles, we combined the coverage of the H3K27ac histone marks at both time points and for both conditions (Figure 4C, D). Generally, these results confirmed the observations of the transcriptome analyses. In addition, we performed differential analysis on the genes encompassed by H3K27ac histone signals and obtained concordant results (Figure S3A,B). Biological pathway analysis of the differentially acetylated genes in the delayed versus immediate cell processing comparison of atopics revealed a variety of pathways affected by downregulated genes (Table 1) of which T-cell receptor and IL 2 signaling pathways were shared between the transcriptomic and epigenomic profiles. Additionally, ribosomal gene transcription and cell cycle-related pathways such as RNA polymerase I transcription and chromosome maintenance, respectively, were diminished during the 24 h ex vivo blood storage (Figure 4E,F). To verify the previous results at the level of H3K27ac peaks, we conducted a differential analysis of the called peaks using the DiffBind analysis tool. As a result, samples from atopic subjects segregated into two clusters demonstrating a discrete separation between the delayed and immediate time points of CD4+ T-cells processing (Figure 5A). The related volcano plot showed a greater gain of peaks in immediately processed samples compared to those prepared at the later time point (Figure 5B). Accordingly, the binding affinity analysis showed a higher mean read concentration to be related with the immediate time point (Figure 5C). On the other side, a similar analysis was applied to the healthy samples of both time points of CD4+ T-cells processing, yet no significant differences were observed in all parts of the analysis as shown in Figure 5D–F. 

4. Discussion

Pre-analytical processes represent a critical step in clinical biomaterial sampling [19]. Therefore, controlling and standardizing as many variables as possible in the experimental design of translational clinical and cohort studies is a critical prerequisite to achieve most reliable outcomes. Experimental pitfalls including timing of blood cell processing in biobank sampling are often not sufficiently considered [20,21]. In the pilot study presented here, we focused on atopy as a target disease condition that involves various activated immune mechanisms and is associated with prominent inflammatory characteristics [22,23]. CD4+ T-cells were selected for this study as a major cell type playing a key role in the underlying adaptive immune response [24].
The most striking result emerging from our data is not only that significant differences in gene expression were observed as a consequence of the ex vivo 24 h rest of blood samples, but also that such changes occurred quite specifically in CD4+ T-cells from atopic donors. Histone acetylation profiles and downstream gene expression associated with atopy was changed in CD4+ T cells following delayed processing. This resulted in the disappearance of significant differences observed in immediately processed CD4+ T-cells derived from atopic versus healthy individuals. This remark was supported by the downregulation of a variety of well-known atopy-related genes and their associated pathways such as IL-2, IL-17 and T cell receptor signaling pathways after 24 h storage of blood samples of atopic individuals [25]. Although too few studies have been conducted to investigate such pre-analytical effects, our results corroborated very well with observations by Baechler et al. [26]. In their study they utilized PBMCs of rheumatoid arthritis and systemic lupus erythematosus patients and observed time-dependent expression alterations in genes involved in inflammatory response, immune signaling and others. The emergence of RNA degradation genes and pathways after delayed CD4+ T-cells processing could propose evidence for the previous results [27,28]. However, the overall RNA quality did not vary much among all samples. This implies that key posttranscriptional control mechanisms/pathways involved in cell homeostasis are specifically triggered in CD4+ T-cells from atopic subjects upon delayed processing, presumably due to discontinuation of signals required for the maintenance of cellular activation during ex vivo blood storage [29]. However, it remains unclear why and how products of RNA degradation genes such as CNOT1 and CNOT6L became specifically upregulated and how such mechanisms were selectively directed against atopy-related mRNA transcripts.
In addition, we investigated the epigenetic landscape as an upstream regulation mechanism of gene expression known to be involved in different allergic and atopic diseases [30]. Concordant results between the transcriptomic and epigenomic profiles were obtained, including a strong decrease in H3K27 acetylation signals spanning 3 kilobases (kb) around the TSSs in delayed processed CD4+ T-cells of atopics. These reduced histone acetylation signals were associated with genes involved in several molecular pathways (see Table 1), among them a variety of pathways involved in the regulation of general molecular and cellular processes such as control of gene expression, e.g., RNA polymerase I transcription or chromosome maintenance pathways [31].
Certainly, our work has some limitations such as small sample size and few time points of cell processing. Nevertheless, it outperformed previous studies in which similar timing effects were interrogated on other cell types derived from healthy donors only [32].
Taken together, our work demonstrates that delayed processing of blood-derived CD4+ T-cells has drastically changed their transcriptomic and epigenomic profiles specifically in atopic patients, which remained widely unaffected in healthy donors. RNA degradation mechanisms were activated in cells from atopic donors during prolonged 24 h ex vivo blood storage. To circumvent these time and condition-dependent pre-analytical effects, it is highly recommended to develop and apply standard operation procedures (SOPs) and protocols established on the basis of biomaterial samples derived from both patients and healthy individuals.

Supplementary Materials

The following supporting information can be downloaded online at: https://www.mdpi.com/article/10.3390/cells11192958/s1, Figure S1: Gating strategy for sorting of CD3+CD4+ T-cells. Stepwise gating is exemplarily indicated from left to right showing plots of related samples from one atopic subject with immediate and delayed processing. The purity of the sorted CD3+CD4+ T-cells used for the analyses is shown in the last plot to the right; Figure S2: Biological processes and pathways affected by the differentially expressed genes in immediately processed CD4+ T-cells in atopic versus healthy. (A) Heatmap depicting the top 10 significant up- and down-regulated genes of immediately processed CD4+ T-cells in atopic versus healthy subjects alongside with their associated biological processes (GO Terms). (B) Biological pathways affected by up- and down-regulated genes of the previous comparison using a significance cut-off of p < 0.05, and (C) heatmap showing three selected biological pathways (indicated by colored lines) and their associated genes; Figure S3: Differential gene analysis of H3K27ac ChIP-Seq analysis in delayed versus immediately processed blood-derived CD4+ T-cells. Volcano plots depicting genes associated with differential H3K27 acetylation identified by ChIP-Seq analyses (significance cut-off at FDR < 0.1) of delayed versus immediately processed blood-derived CD4+ T-cells from (A) atopic and (B) healthy subjects.

Author Contributions

Conceptualization, F.A., C.H., D.P.P., B.S., R.G. and H.G.; Methodology, F.A., K.L., A.H., A.C.J. and R.C.; Software, F.A., C.T. and H.-R.C.; Validation, F.A., C.T., H.-R.C., B.S., R.G. and H.G.; Formal Analysis, F.A., K.L., C.T., H.-R.C., B.S., R.G. and H.G.; Investigation, F.A., K.L., A.H., A.C.J. and R.C.; Resources, G.H., B.S., R.G. and H.G.; Data Curation, F.A., K.L., B.S., R.G. and H.G.; Writing—Original Draft Preparation, F.A.; Writing—Review & Editing, F.A., K.L., C.H., A.H., A.C.J., C.T., R.C., H.-R.C., G.H., D.P.P., B.S., R.G. and H.G.; Visualization, F.A.; Supervision, B.S., R.G. and H.G.; Project Administration, F.A., B.S., R.G. and H.G.; Funding Acquisition, F.A., C.H., G.H., D.P.P., B.S., R.G. and H.G. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Deutsches Zentrum für Lungenforschung (DZL, German Center for Lung Research)—Marburg (UGLMC), Munich (CPC-M) and Hannover (BREATH) sites. F.A. received a personal grant (reference number 91726294) from the Deutscher Akademischer Austauschdienst (DAAD, German Academic Exchange Service). Open Access funding provided by the Open Acess Publication Fund of Philipps-Universität Marburg with support of the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation).

Institutional Review Board Statement

The study was part of the DZL ALLIANCE cohort activities and was approved by the Ethics Committee of the University of Lübeck (Approval No. 12–215).

Informed Consent Statement

Informed consent was obtained from all individual participants involved in the study.

Data Availability Statement

The data presented in this study are openly available in Dryad at doi:10.5061/dryad.zw3r22890.

Acknowledgments

The authors would like to thank Andrea Nist and Thorsten Stiewe from the Genomics Core Facility of the Philipps University of Marburg for sample sequencing.

Conflicts of Interest

Authors declared no conflict of interest.

References

  1. Pei, J.; Schuldt, M.; Nagyova, E.; Gu, Z.; Bouhaddani, S.; Yiangou, L.; Jansen, M.; Calis, J.J.A.; Dorsch, L.M.; Blok, C.S.; et al. Multi-omics integration identifies key upstream regulators of pathomechanisms in hypertrophic cardiomyopathy due to truncating MYBPC3 mutations. Clin. Epigenetics 2021, 13, 1–20. [Google Scholar] [CrossRef] [PubMed]
  2. Andrysik, Z.; Bender, H.; Galbraith, M.D.; Espinosa, J.M. Multi-omics analysis reveals contextual tumor suppressive and oncogenic gene modules within the acute hypoxic response. Nat. Commun. 2021, 12, 1–18. [Google Scholar] [CrossRef]
  3. Jiang, Y.; Gruzieva, O.; Wang, T.; Forno, E.; Boutaoui, N.; Sun, T.; Merid, S.K.; Acosta-Pérez, E.; Kull, I.; Canino, G.; et al. Transcriptomics of atopy and atopic asthma in white blood cells from children and adolescents. Eur. Respir. J. 2019, 53, 1900102. [Google Scholar] [CrossRef]
  4. Park, J.; Lutz, S.M.; Choi, S.; Lee, S.; Park, S.; Kim, K.; Choi, H.; Park, H.; Lee, S.Y.; Weiss, S.T.; et al. Multi-omics analyses implicate EARS2 in the pathogenesis of atopic dermatitis. Allergy 2021, 76, 2602–2604. [Google Scholar] [CrossRef] [PubMed]
  5. Alhamdan, F.; Marsh, L.M.; Pedersen, F.; Alhamwe, B.A.; Thölken, C.; Pfefferle, P.I.; Bahmer, T.; Greulich, T.; Potaczek, D.P.; Garn, H. Differential Regulation of Interferon Signaling Pathways in CD4+ T Cells of the Low Type-2 Obesity-Associated Asthma Phenotype. Int. J. Mol. Sci. 2021, 22, 10144. [Google Scholar] [CrossRef] [PubMed]
  6. Garn, H.; Potaczek, D.P.; Pfefferle, P.I. The Hygiene Hypothesis and New Perspectives—Current Challenges Meeting an Old Postulate. Front. Immunol. 2021, 12, 637087. [Google Scholar] [CrossRef] [PubMed]
  7. Brown, P.; Jarreau, P. Atopy: Marching with allergies. Nat. 2011, 479, S14–S15. [Google Scholar] [CrossRef]
  8. Karakioulaki, M.; Papakonstantinou, E.; Goulas, A.; Stolz, D. The Role of Atopy in COPD and Asthma. Front. Med. 2021, 8, 674742. [Google Scholar] [CrossRef]
  9. Lianidou, E.; Pantel, K. Liquid biopsies. Genes Chromosom. Cancer 2019, 58, 219–232. [Google Scholar] [CrossRef]
  10. Luckheeram, R.V.; Zhou, R.; Verma, A.D.; Xia, B. CD4+T Cells: Differentiation and Functions. Clin. Dev. Immunol. 2012, 2012, 925135. [Google Scholar] [CrossRef] [Green Version]
  11. Narayanan, S. The Preanalytic Phase: An Important Component of Laboratory Medicine. Am. J. Clin. Pathol. 2000, 113, 429–452. [Google Scholar] [CrossRef] [PubMed]
  12. Yin, P.; Lehmann, R.; Xu, G. Effects of pre-analytical processes on blood samples used in metabolomics studies. Anal. Bioanal. Chem. 2015, 407, 4879–4892. [Google Scholar] [CrossRef] [PubMed]
  13. Savage, A.K.; Gutschow, M.V.; Chiang, T.; Henderson, K.; Green, R.; Chaudhari, M.; Swanson, E.; Heubeck, A.T.; Kondza, N.; Burley, K.C.; et al. Multimodal analysis for human ex vivo studies shows extensive molecular changes from delays in blood processing. iScience 2021, 24, 102404. [Google Scholar] [CrossRef]
  14. Paauw, N.D.; Lely, T.; Joles, J.A.; Franx, A.; Nikkels, P.G.; Mokry, M.; Van Rijn, B.B. H3K27 acetylation and gene expression analysis reveals differences in placental chromatin activity in fetal growth restriction. Clin. Epigenetics 2018, 10, 1–11. [Google Scholar] [CrossRef]
  15. Afgan, E.; Baker, D.; Batut, B.; van den Beek, M.; Bouvier, D.; Čech, M.; Chilton, J.; Clements, D.; Coraor, N.; Grüning, B.A.; et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018, 46, W537–W544. [Google Scholar] [CrossRef]
  16. Xie, Z.; Bailey, A.; Kuleshov, M.V.; Clarke, D.J.B.; Evangelista, J.E.; Jenkins, S.L.; Lachmann, A.; Wojciechowicz, M.L.; Kropiwnicki, E.; Jagodnik, K.M.; et al. Gene Set Knowledge Discovery with Enrichr. Curr. Protoc. 2021, 1, e90. [Google Scholar] [CrossRef] [PubMed]
  17. Kuleshov, M.V.; Jones, M.R.; Rouillard, A.D.; Fernandez, N.F.; Duan, Q.; Wang, Z.; Koplev, S.; Jenkins, S.L.; Jagodnik, K.M.; Lachmann, A.; et al. Enrichr: A comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016, 44, W90–W97. [Google Scholar] [CrossRef]
  18. Chen, E.Y.; Tan, C.M.; Kou, Y.; Duan, Q.; Wang, Z.; Meirelles, G.V.; Clark, N.R.; Ma’Ayan, A. Enrichr: Interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinform. 2013, 14, 128. [Google Scholar] [CrossRef] [PubMed]
  19. Dvinge, H.; Ries, R.E.; Ilagan, J.O.; Stirewalt, D.L.; Meshinchi, S.; Bradley, R.K. Sample processing obscures cancer-specific alterations in leukemic transcriptomes. Proc. Natl. Acad. Sci. 2014, 111, 16802–16807. [Google Scholar] [CrossRef]
  20. Lommen, K.; Odeh, S.; de Theije, C.C.; Smits, K.M. Biobanking in Molecular Biomarker Research for the Early Detection of Cancer. Cancers 2020, 12, 776. [Google Scholar] [CrossRef] [Green Version]
  21. Coppola, L.; Cianflone, A.; Grimaldi, A.M.; Incoronato, M.; Bevilacqua, P.; Messina, F.; Baselice, S.; Soricelli, A.; Mirabelli, P.; Salvatore, M. Biobanking in health care: Evolution and future directions. J. Transl. Med. 2019, 17, 1–18. [Google Scholar] [CrossRef] [PubMed]
  22. Lyons, J.J.; Milner, J.D. Primary Atopic Disorders. J. Exp. Med. 2018, 215, 1009–1022. [Google Scholar] [CrossRef] [PubMed]
  23. Lauzon-Joset, J.-F.; Jones, A.C.; Mincham, K.T.; Thomas, J.A.; Rosenthal, L.A.; Bosco, A.; Holt, P.G.; Strickland, D.H. Atopy-Dependent and Independent Immune Responses in the Heightened Severity of Atopics to Respiratory Viral Infections: Rat Model Studies. Front. Immunol. 2018, 9, 1805. [Google Scholar] [CrossRef] [PubMed]
  24. Geginat, J.; Paroni, M.; Maglie, S.; Alfen, J.S.; Kastirr, I.; Gruarin, P.; De Simone, M.; Pagani, M.; Abrignani, S. Plasticity of Human CD4 T Cell Subsets. Front. Immunol. 2014, 5, 630. [Google Scholar] [CrossRef]
  25. Howrylak, J.A.; Moll, M.; Weiss, S.T.; Raby, B.A.; Wu, W.; Xing, E.P. Gene expression profiling of asthma phenotypes demonstrates molecular signatures of atopy and asthma control. J. Allergy Clin. Immunol. 2016, 137, 1390–1397.e6. [Google Scholar] [CrossRef]
  26. Baechler, E.C.; Batliwalla, F.M.; Karypis, G.; Gaffney, P.; Moser, K.; A Ortmann, W.; Espe, K.J.; Balasubramanian, S.; Hughes, K.M.; Chan, J.P.; et al. Expression levels for many genes in human peripheral blood cells are highly sensitive to ex vivo incubation. Genes Immun. 2004, 5, 347–353. [Google Scholar] [CrossRef]
  27. Romero, I.G.; A Pai, A.; Tung, J.; Gilad, Y. RNA-seq: Impact of RNA degradation on transcript quantification. BMC Biol. 2014, 12, 42. [Google Scholar] [CrossRef]
  28. Shirai, Y.-T.; Suzuki, T.; Morita, M.; Takahashi, A.; Yamamoto, T. Multifunctional roles of the mammalian CCR4-NOT complex in physiological phenomena. Front. Genet. 2014, 5, 286. [Google Scholar] [CrossRef]
  29. Mino, T.; Takeuchi, O. Post-Transcriptional Regulation of Immune Responses by RNA Binding Proteins. Proc. Jpn. Acad. Ser. B Phys. Biol. Sci. 2018, 94, 248–258. [Google Scholar] [CrossRef]
  30. Alhamwe, B.A.; Alhamdan, F.; Ruhl, A.; Potaczek, D.P.; Renz, H. The role of epigenetics in allergy and asthma development. Curr. Opin. Allergy Clin. Immunol. 2020, 20, 48–55. [Google Scholar] [CrossRef]
  31. McStay, B.; Grummt, I. The Epigenetics of rRNA Genes: From Molecular to Chromosome Biology. Annu. Rev. Cell Dev. Biol. 2008, 24, 131–157. [Google Scholar] [CrossRef] [PubMed]
  32. Navas, A.; Giraldo-Parra, L.; Prieto, M.D.; Cabrera, J.; Gómez, M.A. Phenotypic and functional stability of leukocytes from human peripheral blood samples: Considerations for the design of immunological studies. BMC Immunol. 2019, 20, 5. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of the general study design.
Figure 1. Schematic representation of the general study design.
Cells 11 02958 g001
Figure 2. Differential transcriptome profiles in immediately versus delayed processed CD4+ T-cells from atopic and healthy subjects. Volcano plots presenting differentially regulated genes of (A) immediately and (B) delayed processed blood-derived CD4+ T-cells of atopic versus healthy subjects, and delayed versus immediately processed cells of (C) atopic subjects and (D) healthy subjects, all with significance threshold at FDR < 0.1.
Figure 2. Differential transcriptome profiles in immediately versus delayed processed CD4+ T-cells from atopic and healthy subjects. Volcano plots presenting differentially regulated genes of (A) immediately and (B) delayed processed blood-derived CD4+ T-cells of atopic versus healthy subjects, and delayed versus immediately processed cells of (C) atopic subjects and (D) healthy subjects, all with significance threshold at FDR < 0.1.
Cells 11 02958 g002
Figure 3. Biological processes and pathways affected by the differentially expressed genes in immediately processed CD4+ T-cells from atopic versus healthy subjects. (A) Heatmap depicting the top 10 significant up- and down-regulated genes of delayed versus immediately processed CD4+ T-cells from atopic subjects alongside their associated biological processes (GO Terms). (B) Biological pathways affected by up- and down-regulated genes of the previous comparison using a significance cut-off of p < 0.05, and (C) heatmap showing three selected biological pathways (indicated by colored lines) and their associated genes.
Figure 3. Biological processes and pathways affected by the differentially expressed genes in immediately processed CD4+ T-cells from atopic versus healthy subjects. (A) Heatmap depicting the top 10 significant up- and down-regulated genes of delayed versus immediately processed CD4+ T-cells from atopic subjects alongside their associated biological processes (GO Terms). (B) Biological pathways affected by up- and down-regulated genes of the previous comparison using a significance cut-off of p < 0.05, and (C) heatmap showing three selected biological pathways (indicated by colored lines) and their associated genes.
Cells 11 02958 g003
Figure 4. Selective histone modification changes within immediately versus delayed processed CD4+ T-cells from atopic subjects. (A,B) Genomic heatmaps visualizing H3K27 acetylation ChIP-Seq analyses scores associated with genomic regions of 3 kilobases (kb) around the gene transcription start site (TSS) and (C,D) average enrichments over all genomic regions of immediately and delayed processed blood-derived CD4+ T-cells from atopic (A,C) and healthy (B,D) subjects. (E,F) Gene set enrichment analysis (GSEA) of downregulated genes in delayed versus immediately processed blood-derived CD4+ T-cells from atopic subjects.
Figure 4. Selective histone modification changes within immediately versus delayed processed CD4+ T-cells from atopic subjects. (A,B) Genomic heatmaps visualizing H3K27 acetylation ChIP-Seq analyses scores associated with genomic regions of 3 kilobases (kb) around the gene transcription start site (TSS) and (C,D) average enrichments over all genomic regions of immediately and delayed processed blood-derived CD4+ T-cells from atopic (A,C) and healthy (B,D) subjects. (E,F) Gene set enrichment analysis (GSEA) of downregulated genes in delayed versus immediately processed blood-derived CD4+ T-cells from atopic subjects.
Cells 11 02958 g004
Figure 5. DiffBind analysis of called peaks of H3K27 acetylation ChIP-Seq analysis in delayed versus immediately processed blood-derived CD4+ T-cells of atopic and healthy subjects. (A) PCA plot showing the distance between samples based on the called peaks of delayed vs. immediate (atopic) comparison. (B) Volcano plot showing the differentially called peaks of the previous comparison (significance cut-off at FDR < 0.1). (C) Box plots of read distributions for significantly differentially called peaks of the previous comparison (binding affinity) for all sites (all), sites with increased affinity in the atopic delayed (+), and sites with increased affinity the atopic immediate (−) condition. (DF) This is the same as the previous description for delayed vs. immediate (healthy) comparison.
Figure 5. DiffBind analysis of called peaks of H3K27 acetylation ChIP-Seq analysis in delayed versus immediately processed blood-derived CD4+ T-cells of atopic and healthy subjects. (A) PCA plot showing the distance between samples based on the called peaks of delayed vs. immediate (atopic) comparison. (B) Volcano plot showing the differentially called peaks of the previous comparison (significance cut-off at FDR < 0.1). (C) Box plots of read distributions for significantly differentially called peaks of the previous comparison (binding affinity) for all sites (all), sites with increased affinity in the atopic delayed (+), and sites with increased affinity the atopic immediate (−) condition. (DF) This is the same as the previous description for delayed vs. immediate (healthy) comparison.
Cells 11 02958 g005
Table 1. Top significant 30 biological pathways affected by down-regulated genes in H3K27ac ChIP-Seq analysis of delayed versus immediately processed blood-derived CD4+ T-cells from atopic subjects.
Table 1. Top significant 30 biological pathways affected by down-regulated genes in H3K27ac ChIP-Seq analysis of delayed versus immediately processed blood-derived CD4+ T-cells from atopic subjects.
Termp-ValueFDR
Translation3.36 × 10−203.04 × 10−17
Influenza viral RNA transcription and replication6.97 × 10−203.16 × 10−17
Cytoplasmic ribosomal proteins1.98 × 10−196.00 × 10−17
Systemic lupus erythematosus5.34 × 10−191.21 × 10−16
Influenza infection6.30 × 10−171.14 × 10−14
Packaging of telomere ends3.26 × 10−134.93 × 10−11
RNA polymerase I promoter opening4.53 × 10−135.87 × 10−11
Cap-dependent translation initiation4.55 × 10−125.16 × 10−10
T cell receptor regulation of apoptosis4.31 × 10−114.34 × 10−09
Activation of mRNA upon binding of the cap-binding complex and eIFs, and subsequent binding to 43S6.91 × 10−116.27 × 10−09
Meiotic recombination1.91 × 10−101.57 × 10−08
Protein metabolism6.45 × 10−104.87 × 10−08
Amyloids1.79 × 10−091.25 × 10−07
Telomere maintenance2.18 × 10−091.41 × 10−07
Gene expression2.99 × 10−091.81 × 10−07
Meiotic synapsis6.98 × 10−093.76 × 10−07
RNA polymerase I transcription7.46 × 10−093.76 × 10−07
Meiosis7.46 × 10−093.76 × 10−07
Deposition of new CENP-A-containing nucleosomes at the centromere2.34 × 10−081.12 × 10−06
Transcription2.49 × 10−081.13 × 10−06
RNA polymerase I, RNA polymerase III, and mitochondrial transcription4.82 × 10−062.04 × 10−04
Type II interferon signaling (interferon-gamma)4.95 × 10−062.04 × 10−04
Chromosome maintenance5.32 × 10−062.10 × 10−04
Interleukin-2 signaling pathway1.52 × 10−055.73 × 10−04
Messenger RNA splicing: major pathway1.95 × 10−057.09 × 10−04
Diurnally regulated genes with circadian orthologs3.06 × 10−050.001069
Respiratory electron transport, ATP biosynthesis by chemiosmotic coupling, and heat production by uncoupling proteins1.18 × 10−040.003957
p75 neurotrophin receptor signaling via NF-kB7.34 × 10−040.023787
Clathrin derived vesicle budding8.94 × 10−040.027101
mRNA stability regulation by proteins that bind AU-rich elements8.96 × 10−040.027101
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alhamdan, F.; Laubhahn, K.; Happle, C.; Habener, A.; Jirmo, A.C.; Thölken, C.; Conca, R.; Chung, H.-R.; Hansen, G.; Potaczek, D.P.; et al. Timing of Blood Sample Processing Affects the Transcriptomic and Epigenomic Profiles in CD4+ T-cells of Atopic Subjects. Cells 2022, 11, 2958. https://doi.org/10.3390/cells11192958

AMA Style

Alhamdan F, Laubhahn K, Happle C, Habener A, Jirmo AC, Thölken C, Conca R, Chung H-R, Hansen G, Potaczek DP, et al. Timing of Blood Sample Processing Affects the Transcriptomic and Epigenomic Profiles in CD4+ T-cells of Atopic Subjects. Cells. 2022; 11(19):2958. https://doi.org/10.3390/cells11192958

Chicago/Turabian Style

Alhamdan, Fahd, Kristina Laubhahn, Christine Happle, Anika Habener, Adan C. Jirmo, Clemens Thölken, Raffaele Conca, Ho-Ryun Chung, Gesine Hansen, Daniel P. Potaczek, and et al. 2022. "Timing of Blood Sample Processing Affects the Transcriptomic and Epigenomic Profiles in CD4+ T-cells of Atopic Subjects" Cells 11, no. 19: 2958. https://doi.org/10.3390/cells11192958

APA Style

Alhamdan, F., Laubhahn, K., Happle, C., Habener, A., Jirmo, A. C., Thölken, C., Conca, R., Chung, H. -R., Hansen, G., Potaczek, D. P., Schaub, B., Grychtol, R., & Garn, H. (2022). Timing of Blood Sample Processing Affects the Transcriptomic and Epigenomic Profiles in CD4+ T-cells of Atopic Subjects. Cells, 11(19), 2958. https://doi.org/10.3390/cells11192958

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop