Genetic variants regulating immune cell levels in health and disease.

The complex network of specialized cells and molecules in the immune system has evolved to defend against pathogens, but inadvertent immune system attacks on "self" result in autoimmune disease. Both genetic regulation of immune cell levels and their relationships with autoimmunity are largely undetermined. Here, we report genetic contributions to quantitative levels of 95 cell types encompassing 272 immune traits, in a cohort of 1,629 individuals from four clustered Sardinian villages. We first estimated trait heritability, showing that it can be substantial, accounting for up to 87% of the variance (mean 41%). Next, by assessing ∼8.2 million variants that we identified and confirmed in an extended set of 2,870 individuals, 23 independent variants at 13 loci associated with at least one trait. Notably, variants at three loci (HLA, IL2RA, and SH2B3/ATXN2) overlap with known autoimmune disease associations. These results connect specific cellular phenotypes to specific genetic variants, helping to explicate their involvement in disease.

The longitudinal dynamics and natural history of clonal haematopoiesis.

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis1-4. The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited5,6. Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2P95H). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection.

GWAS of genetic factors affecting white blood cell morphological parameters in Sardinians uncovers influence of chromosome 11 innate immunity gene cluster on eosinophil morphology.

Few genome-wide association studies (GWAS) analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of 12 morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (P-value = 8.33 x 10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

Autoimmunity Is a Significant Feature of Idiopathic Pulmonary Arterial Hypertension.

Rationale: Autoimmunity is believed to play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear whether this is causative or a bystander of disease and if it carries any prognostic or treatment significance. Objectives: To study autoimmunity in IPAH using a large cross-sectional cohort. Methods: Assessment of the circulating immune cell phenotype was undertaken using flow cytometry, and the profile of serum immunoglobulins was generated using a standardized multiplex array of 19 clinically validated autoantibodies in 473 cases and 946 control subjects. Additional glutathione S-transferase fusion array and ELISA data were used to identify a serum autoantibody to BMPR2 (bone morphogenetic protein receptor type 2). Clustering analyses and clinical correlations were used to determine associations between immunogenicity and clinical outcomes. Measurements and Main Results: Flow cytometric immune profiling demonstrates that IPAH is associated with an altered humoral immune response in addition to raised IgG3. Multiplexed autoantibodies were significantly raised in IPAH, and clustering demonstrated three distinct clusters: "high autoantibody," "low autoantibody," and a small "intermediate" cluster exhibiting high concentrations of ribonucleic protein complex. The high-autoantibody cluster had worse hemodynamics but improved survival. A small subset of patients demonstrated immunoglobulin reactivity to BMPR2. Conclusions: This study establishes aberrant immune regulation and presence of autoantibodies as key features in the profile of a significant proportion of patients with IPAH and is associated with clinical outcomes.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets.

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

Genetic variant within CDK6 regulates immune response to palbociclib treatment.

Everyone carries a set of genetic variants that contribute to regulation of the levels of blood cells, with unknown clinical impact. One of them, rs445 within the cell-cycle checkpoint gene CDK6, reduces the levels of myeloid cell types including granulocytes. We treated CD3+ T cells and whole blood with palbociclib in 41 individuals, who were stratified by genotype for analyses. In T cells we assessed cell cycle and apoptosis, whereas in whole blood, apoptosis in activated (CD11b+), unactivated (CD11b-) granulocytes, cytotoxic (CD8 + CD4-), and helper (CD8-CD4+) T cells. We find that rs445 modulates the immune response of CD8+ T cells. It also increases the level of apoptotic CD11b + activated granulocytes after palbociclib treatment, which, in synergy with neutropenia, may affect drug related adverse events. These results suggest that the effect of palbociclib treatment may depend on underlying genetically encoded individual immune response as well as the direct response to the drug.

Immunologic resilience and COVID-19 survival advantage.

BACKGROUND: The risk of severe coronavirus disease 2019 (COVID-19) varies significantly among persons of similar age and is higher in males. Age-independent, sex-biased differences in susceptibility to severe COVID-19 may be ascribable to deficits in a sexually dimorphic protective attribute that we termed immunologic resilience (IR). OBJECTIVE: We sought to examine whether deficits in IR that antedate or are induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection independently predict COVID-19 mortality. METHODS: IR levels were quantified with 2 novel metrics: immune health grades (IHG-I [best] to IHG-IV) to gauge CD8+ and CD4+ T-cell count equilibrium, and blood gene expression signatures. IR metrics were examined in a prospective COVID-19 cohort (n = 522); primary outcome was 30-day mortality. Associations of IR metrics with outcomes in non-COVID-19 cohorts (n = 13,461) provided the framework for linking pre-COVID-19 IR status to IR during COVID-19, as well as to COVID-19 outcomes. RESULTS: IHG-I, tracking high-grade equilibrium between CD8+ and CD4+ T-cell counts, was the most common grade (73%) among healthy adults, particularly in females. SARS-CoV-2 infection was associated with underrepresentation of IHG-I (21%) versus overrepresentation (77%) of IHG-II or IHG-IV, especially in males versus females (P < .01). Presentation with IHG-I was associated with 88% lower mortality, after controlling for age and sex; reduced risk of hospitalization and respiratory failure; lower plasma IL-6 levels; rapid clearance of nasopharyngeal SARS-CoV-2 burden; and gene expression signatures correlating with survival that signify immunocompetence and controlled inflammation. In non-COVID-19 cohorts, IR-preserving metrics were associated with resistance to progressive influenza or HIV infection, as well as lower 9-year mortality in the Framingham Heart Study, especially in females. CONCLUSIONS: Preservation of immunocompetence with controlled inflammation during antigenic challenges is a hallmark of IR and associates with longevity and AIDS resistance. Independent of age, a male-biased proclivity to degrade IR before and/or during SARS-CoV-2 infection predisposes to severe COVID-19.

Quantifying the Detrimental Effects of Multiple Freeze/Thaw Cycles on Primary Human Lymphocyte Survival and Function.

The use of cryopreserved peripheral blood mononuclear cells is common in biological research. It is widely accepted that primary cells are rendered unusable by several freezing cycles, although this practice might be very helpful when the biological material is valuable and its re-collection is impractical. To determine the extent to which primary cells undergoing repeated freezing cycles are comparable to one another and to fresh samples, we evaluated overall lymphocyte viability, their proliferation and cytokine production capabilities, as well as the levels of 27 cell subtypes in ten human peripheral blood mononuclear cells frozen for five years and repeatedly thawed. As expected, we observed a progressive increase in cell death percentages on three rounds of thawing, but the frequency of the main lymphocyte subsets was stable across the three thawings. Nevertheless, we observed a significant reduction of B cell frequency in frozen samples compared to fresh ones. On repeated thawings and subsequent conventional stimulation, lymphocyte proliferation significantly decreased, and IL-10, IL-6, GM-CSF, IFN-gamma, and IL-8 showed a trend to lower values.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity.

BACKGROUND: Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm. OBJECTIVE: The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D). METHODS: We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians. RESULTS: We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10-5, OR = 0.82. CONCLUSION: Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

CXCR3 Identifies Human Naive CD8+ T Cells with Enhanced Effector Differentiation Potential.

In mice, the ability of naive T (TN) cells to mount an effector response correlates with TCR sensitivity for self-derived Ags, which can be quantified indirectly by measuring surface expression levels of CD5. Equivalent findings have not been reported previously in humans. We identified two discrete subsets of human CD8+ TN cells, defined by the absence or presence of the chemokine receptor CXCR3. The more abundant CXCR3+ TN cell subset displayed an effector-like transcriptional profile and expressed TCRs with physicochemical characteristics indicative of enhanced interactions with peptide-HLA class I Ags. Moreover, CXCR3+ TN cells frequently produced IL-2 and TNF in response to nonspecific activation directly ex vivo and differentiated readily into Ag-specific effector cells in vitro. Comparative analyses further revealed that human CXCR3+ TN cells were transcriptionally equivalent to murine CXCR3+ TN cells, which expressed high levels of CD5. These findings provide support for the notion that effector differentiation is shaped by heterogeneity in the preimmune repertoire of human CD8+ T cells.

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

BACKGROUND: Genomewide association studies of autoimmune diseases have mapped hundreds of susceptibility regions in the genome. However, only for a few association signals has the causal gene been identified, and for even fewer have the causal variant and underlying mechanism been defined. Coincident associations of DNA variants affecting both the risk of autoimmune disease and quantitative immune variables provide an informative route to explore disease mechanisms and drug-targetable pathways. METHODS: Using case-control samples from Sardinia, Italy, we performed a genomewide association study in multiple sclerosis followed by TNFSF13B locus-specific association testing in systemic lupus erythematosus (SLE). Extensive phenotyping of quantitative immune variables, sequence-based fine mapping, cross-population and cross-phenotype analyses, and gene-expression studies were used to identify the causal variant and elucidate its mechanism of action. Signatures of positive selection were also investigated. RESULTS: A variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), was associated with multiple sclerosis as well as SLE. The disease-risk allele was also associated with up-regulated humoral immunity through increased levels of soluble BAFF, B lymphocytes, and immunoglobulins. The causal variant was identified: an insertion-deletion variant, GCTGT→A (in which A is the risk allele), yielded a shorter transcript that escaped microRNA inhibition and increased production of soluble BAFF, which in turn up-regulated humoral immunity. Population genetic signatures indicated that this autoimmunity variant has been evolutionarily advantageous, most likely by augmenting resistance to malaria. CONCLUSIONS: A TNFSF13B variant was associated with multiple sclerosis and SLE, and its effects were clarified at the population, cellular, and molecular levels. (Funded by the Italian Foundation for Multiple Sclerosis and others.).

Cooperative translational control of polymorphic BAFF by NF90 and miR-15a.

Polymorphisms in untranslated regions (UTRs) of disease-associated mRNAs can alter protein production. We recently identified a genetic variant in the 3'UTR of the TNFSF13B gene, encoding the cytokine BAFF (B-cell-activating factor), that generates an alternative polyadenylation site yielding a shorter, more actively translated variant, BAFF-var mRNA. Accordingly, individuals bearing the TNFSF13B variant had higher circulating BAFF and elevated risk of developing autoimmune diseases. Here, we investigated the molecular mechanisms controlling the enhanced translation of BAFF-var mRNA. We identified nuclear factor 90 (NF90, also known as ILF3) as an RNA-binding protein that bound preferentially the wild-type (BAFF-WT mRNA) but not BAFF-var mRNA in human monocytic leukemia THP-1 cells. NF90 selectively suppressed BAFF translation by recruiting miR-15a to the 3'UTR of BAFF-WT mRNA. Our results uncover a paradigm whereby an autoimmunity-causing BAFF polymorphism prevents NF90-mediated recruitment of microRNAs to suppress BAFF translation, raising the levels of disease-associated BAFF.

CD8+CD28-CD127loCD39+ regulatory T-cell expansion: A new possible pathogenic mechanism for HIV infection?

BACKGROUND: HIV-associated immunodeficiency is related to loss of CD4+ T cells. This mechanism does not explain certain manifestations of HIV disease, such as immunodeficiency events in patients with greater than 500 CD4+ T cells/µL. CD8+CD28-CD127loCD39+ T cells are regulatory T (Treg) lymphocytes that are highly concentrated within the tumor microenvironment and never analyzed in the circulation of HIV-infected patients. OBJECTIVES: We sought to analyze the frequency of CD8+CD28-CD127loCD39+ Treg cells in the circulation of HIV-infected patients. METHODS: The frequency of circulating CD8+CD28-CD127loCD39+ Treg cells was analyzed and correlated with viral load and CD4+ T-cell counts/percentages in 93 HIV-1-infected patients subdivided as follows: naive (n = 63), elite controllers (n = 19), long-term nonprogressors (n = 7), and HIV-infected patients affected by tumor (n = 4). The same analyses were performed in HIV-negative patients with cancer (n = 53), hepatitis C virus-infected patients (n = 17), and healthy donors (n = 173). RESULTS: HIV-infected patients had increased circulating levels of functional CD8+CD28-CD127loCD39+ Treg cells. These cells showed antigen specificity against HIV proteins. Their frequency after antiretroviral therapy (ART) correlated with HIV viremia, CD4+ T-cell counts, and immune activation markers, suggesting their pathogenic involvement in AIDS- or non-AIDS-related complications. Their increase after initiation of ART heralded a lack of virologic or clinical response, and hence their monitoring is clinically relevant. CONCLUSION: HIV infection induces remarkable expansion of CD8+CD28-CD127loCD39+ Treg cells, the frequency of which correlates with both clinical disease and signs of chronic immune cell activation. Monitoring their frequency in the circulation is a new marker of response to ART when effects on viremia and clinical response are not met.

Autoimmune-associated lymphoid tyrosine phosphatase is a gain-of-function variant.

A SNP in the gene PTPN22 is associated with type 1 diabetes, rheumatoid arthritis, lupus, Graves thyroiditis, Addison disease and other autoimmune disorders. T cells from carriers of the predisposing allele produce less interleukin-2 upon TCR stimulation, and the encoded phosphatase has higher catalytic activity and is a more potent negative regulator of T lymphocyte activation. We conclude that the autoimmune-predisposing allele is a gain-of-function mutant.

Immune and glial cell alterations in the rat brain after repeated exposure to the synthetic cannabinoid JWH-018.

The use of synthetic cannabinoid receptor agonists (SCRAs) poses major psychiatric risks. We previously showed that repeated exposure to the prototypical SCRA JWH-018 induces alterations in dopamine (DA) transmission, abnormalities in the emotional state, and glial cell activation in the mesocorticolimbic DA circuits of rats. Despite growing evidence suggesting the relationship between substance use disorders (SUD) and neuroinflammation, little is known about the impact of SCRAs on the neuroimmune system. Here, we investigated whether repeated JWH-018 exposure altered neuroimmune signaling, which could be linked with previously reported central effects. Adult male Sprague-Dawley (SD) rats were exposed to JWH-018 (0.25 mg/kg, i.p.) for fourteen consecutive days, and the expression of cytokines, chemokines, and growth factors was measured seven days after treatment discontinuation in the striatum, cortex, and hippocampus. Moreover, microglial (ionized calcium-binding adaptor molecule 1, IBA-1) and astrocyte (glial fibrillary acidic protein, GFAP) activation markers were evaluated in the caudate-putamen (CPu). Repeated JWH-018 exposure induces a perturbation of neuroimmune signaling specifically in the striatum, as shown by increased levels of cytokines [interleukins (IL) -2, -4, -12p70, -13, interferon (IFN) γ], chemokines [macrophage inflammatory protein (MIP) -1α, -3α], and growth factors [macrophage colony-stimulating factor (M-CSF), vascular endothelial growth factor (VEGF)], together with increased IBA-1 and GFAP expression in the CPu. JWH-018 exposure induces persistant brain region-specific immune alterations up to seven days after drug discontinuation, which may contribute to the behavioral and neurochemical dysregulations in striatal areas that play a role in the reward-related processes that are frequently impaired in SUD.

Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset.

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

The PPARα agonist fenofibrate reduces the cytokine imbalance in a maternal immune activation model of schizophrenia.

Maternal infections during pregnancy may increase the risk of psychiatric disorders in offspring. We recently demonstrated that activation of peroxisome proliferator-activate receptor-α (PPARα), with the clinically available agonist fenofibrate (FEN), attenuates the neurodevelopmental disturbances induced by maternal immune activation (MIA) in rat offspring. We hypothesized that fenofibrate might reduce MIA-induced cytokine imbalance using a MIA model based on the viral mimetic polyriboinosinic-polyribocytidilic acid [poly (I:C)]. By using the Bio-Plex Multiplex-Immunoassay-System, we measured cytokine/chemokine/growth factor levels in maternal serum and in the fetal brain of rats treated with fenofibrate, at 6 and 24 h after poly (I:C). We found that MIA induced time-dependent changes in the levels of several cytokines/chemokines/colony-stimulating factors (CSFs). Specifically, the maternal serum of the poly (I:C)/control (CTRL) group showed increased levels of (i) proinflammatory chemokine macrophage inflammatory protein 1-alpha (MIP-1α), (ii) tumor necrosis factor-alpha (TNF-α), the monocyte chemoattractant protein-1 (MCP-1), the macrophage (M-CSF) and granulocyte-macrophage colony-stimulating factor (GM-CSF). Conversely, in the fetal brain of the poly (I:C)/CTRL group, interleukin 12p70 and MIP-1α levels were lower than in vehicle (veh)/CTRL group. Notably, MIP-1α, TNF-α, keratinocyte derived chemokine (GRO/KC), GM-CSF, and M-CSF levels were lower in the poly (I:C)/FEN than in poly (I:C)/CTRL rats, suggesting the protective role of the PPARα agonist. PPARα might represent a therapeutic target to attenuate MIA-induced inflammation.

Immune resilience despite inflammatory stress promotes longevity and favorable health outcomes including resistance to infection.

Some people remain healthier throughout life than others but the underlying reasons are poorly understood. Here we hypothesize this advantage is attributable in part to optimal immune resilience (IR), defined as the capacity to preserve and/or rapidly restore immune functions that promote disease resistance (immunocompetence) and control inflammation in infectious diseases as well as other causes of inflammatory stress. We gauge IR levels with two distinct peripheral blood metrics that quantify the balance between (i) CD8+ and CD4+ T-cell levels and (ii) gene expression signatures tracking longevity-associated immunocompetence and mortality-associated inflammation. Profiles of IR metrics in ~48,500 individuals collectively indicate that some persons resist degradation of IR both during aging and when challenged with varied inflammatory stressors. With this resistance, preservation of optimal IR tracked (i) a lower risk of HIV acquisition, AIDS development, symptomatic influenza infection, and recurrent skin cancer; (ii) survival during COVID-19 and sepsis; and (iii) longevity. IR degradation is potentially reversible by decreasing inflammatory stress. Overall, we show that optimal IR is a trait observed across the age spectrum, more common in females, and aligned with a specific immunocompetence-inflammation balance linked to favorable immunity-dependent health outcomes. IR metrics and mechanisms have utility both as biomarkers for measuring immune health and for improving health outcomes.

Comparison of Whole Blood Cryopreservation Methods for Extensive Flow Cytometry Immunophenotyping.

Fresh blood immunophenotyping by flow cytometry, based on the reliable simultaneous detection of several markers in a cell, is the method of choice to study the circulating human immune system. Especially in large and multicenter studies, high sample quality is difficult to achieve, and adequate collection and storage of samples with fine-tuned whole blood cryopreservation is mandatory. Here, we compared the quality of immunophenotypic data obtained from fresh blood with those obtained after five cryopreservation methods by quantifying the levels of 41 immune cell populations. They comprised B and T lymphocyte subsets and their maturation stages, as well as monocytes and granulocytes. Three methods used fixative solutions and two other methods used dimethyl sulfoxide solutions to preserve cell viability. The fixative methods prevented detection of markers critical for identification of B and T cell subsets, including CD27, CXCR3, and CCR6. The other two methods permitted reliable discrimination of most immune-cell populations in thawed samples, though some cell frequencies varied compared to the corresponding fresh sample. Of those two methods, the one preserving blood in media containing dimethyl sulfoxide produced results that were most similar to those with fresh samples.