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.

Whole genome analysis of plasma fibrinogen reveals population-differentiated genetic regulators with putative liver roles.

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation. Key Points: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.

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.

Polygenic risk score and biochemical/environmental variables predict a low-risk profile of age-related macular degeneration in Sardinia.

PURPOSE: To ascertain the prevalence and clinical and genetic features of age-related macular degeneration (AMD) in subjects living in the Lanusei valley, Central Sardinia, Italy, involved in a study on ageing (SardiNIA project). METHODS: A total of 814 volunteers aged ≥ 50 years, randomly selected from the SardiNIA project dataset, were included. A color fundus (CF) photograph of the 30° central retina of each eye was obtained and graded according to the Age-Related Eye Disease Study system. Life-style choices were investigated using standardized questionnaires. The concentrations of several inflammatory biomarkers (i.e., complement component, fibrinogen, and C-reactive protein) were measured. Polygenic risk score (PRS) was calculated and compared with results obtained from a European cohort. RESULTS: A total of 756 subjects had gradable CF photographs for AMD detection. In 91.3%, no signs of AMD were observed. The prevalence rates of early and late AMDs were 6.9% and 0.6%, respectively. A total of 85% of subjects were physically active; only 13.5% were current smokers. Low concentrations of complement component, fibrinogen, and C-reactive protein were found. We calculated the polygenic risk scores (PRS) using 40 AMD markers distributed on several candidate genes in Europeans and Sardinians. The mean PRS value was significantly lower in Sardinians than in the Europeans (0.21 vs. 0.248, respectively, p = 1.18 × 10-77). CONCLUSIONS: In our cohort, most subjects showed no sign of any AMD type and late AMD was a condition rarely observed. Results of genetic, biochemical, and life-style investigation support the hypothesis that Sardinia population may present of a peculiar background with a protective effect against AMD development.

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.

Cross-sectional analysis of the humoral response after SARS-CoV-2 vaccination in Sardinian multiple sclerosis patients, a follow-up study.

Monitoring immune responses to SARS-CoV-2 vaccination and its clinical efficacy over time in Multiple Sclerosis (MS) patients treated with disease-modifying therapies (DMTs) help to establish the optimal strategies to ensure adequate COVID-19 protection without compromising disease control offered by DMTs. Following our previous observations on the humoral response one month after two doses of BNT162b2 vaccine (T1) in MS patients differently treated, here we present a cross-sectional and longitudinal follow-up analysis six months following vaccination (T2, n=662) and one month following the first booster (T3, n=185). Consistent with results at T1, humoral responses were decreased in MS patients treated with fingolimod and anti-CD20 therapies compared with untreated patients also at the time points considered here (T2 and T3). Interestingly, a strong upregulation one month after the booster was observed in patients under every DMTs analyzed, including those treated with fingolimod and anti-CD20 therapies. Although patients taking these latter therapies had a higher rate of COVID-19 infection five months after the first booster, only mild symptoms that did not require hospitalization were reported for all the DMTs analyzed here. Based on these findings we anticipate that additional vaccine booster shots will likely further improve immune responses and COVID-19 protection in MS patients treated with any DMT.

Impact on breast cancer susceptibility and clinicopathological traits of common genetic polymorphisms in TP53, MDM2 and ATM genes in Sardinian women.

Common variants of genes involved in DNA damage correction [tumor protein p53 (TP53), murine double 2 homolog oncoprotein (MDM2) and ataxia-telengiectasia mutated (ATM)] may serve a role in cancer predisposition. The purpose of the present study was to investigate the association of five variants in these genes with breast cancer risk and clinicopathological traits in a cohort of 261 women from northern Sardinia. Polymorphic variants in TP53 (rs17878362, rs1042522 and rs1625895), MDM2 (rs2279744) and ATM (rs1799757) were determined by PCR and TaqMan single nucleotide polymorphism assay in patients with breast cancer (n=136) and healthy controls (n=125). Association with clinicopathological (e.g., age at diagnosis, lymph node involvement, clinical stage) and lifestyle factors (e.g., smoking status, alcohol intake, contraceptive use) was also evaluated. TP53 rs17878362 and rs1625895 polymorphisms were associated with decreased risk of BC diagnosis in patients older than 50 years (codominant and recessive models) and post-menopause (recessive model). Furthermore, there was a significant association between lymph node status (positive vs. negative) and ATM rs1799757-delT in dominant and additive models and between MDM2 rs2279744-allele and use of oral contraceptives. This analysis suggested that TP53 rs17878362 and rs1625895 may affect age of onset of breast cancer and ATM rs1799757 and MDM2 rs2279744 may be associated with lymph node status and prolonged use of oral contraceptives, respectively.

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.

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.

Effect of Different Disease-Modifying Therapies on Humoral Response to BNT162b2 Vaccine in Sardinian Multiple Sclerosis Patients.

Objectives: Vaccination against COVID-19 is highly recommended to patients affected by multiple sclerosis (MS); however, the impact of MS disease-modifying therapies (DMTs) on the immune response following vaccination has been only partially investigated. Here, we aimed to elucidate the effect of DMTs on the humoral immune response to mRNA-based anti-SARS-CoV-2 vaccines in MS patients. Methods: We obtained sera from 912 Sardinian MS patients and 63 healthy controls 30 days after the second dose of BNT162b2 vaccine and tested them for SARS-CoV-2 response using anti-Spike (S) protein-based serology. Previous SARS-CoV-2 infection was assessed by anti-Nucleocapsid (N) serology. Patients were either untreated or undergoing treatment with a total of 13 different DMTs. Differences between treatment groups comprised of at least 10 patients were assessed by generalized linear mixed-effects model. Demographic and clinical data and smoking status were analyzed as additional factors potentially influencing humoral immunity from COVID-19 vaccine. Results: MS patients treated with natalizumab, teriflunomide, azathioprine, fingolimod, ocrelizumab, and rituximab showed significantly lower humoral responses compared to untreated patients. We did not observe a statistically significant difference in response between patients treated with the other drugs (dimethyl fumarate, interferon, alemtuzumab and glatiramer acetate) and untreated patients. In addition, older age, male sex and active smoking were significantly associated with lower antibody titers against SARS-CoV-2. MS patients previously infected with SARS-CoV-2 had significantly higher humoral responses to vaccine than uninfected patients. Conclusion: Humoral response to BNT162b2 is significantly influenced by the specific DMTs followed by patients, as well as by other factors such as previous SARS-CoV-2 infection, age, sex, and smoking status. These results are important to inform targeted strategies to prevent clinically relevant COVID-19 in MS patients.

Combining Human Genetics of Multiple Sclerosis with Oxidative Stress Phenotype for Drug Repositioning.

In multiple sclerosis (MS), oxidative stress (OS) is implicated in the neurodegenerative processes that occur from the beginning of the disease. Unchecked OS initiates a vicious circle caused by its crosstalk with inflammation, leading to demyelination, axonal damage and neuronal loss. The failure of MS antioxidant therapies relying on the use of endogenous and natural compounds drives the application of novel approaches to assess target relevance to the disease prior to preclinical testing of new drug candidates. To identify drugs that can act as regulators of intracellular oxidative homeostasis, we applied an in silico approach that links genome-wide MS associations and molecular quantitative trait loci (QTLs) to proteins of the OS pathway. We found 10 drugs with both central nervous system and oral bioavailability, targeting five out of the 21 top-scoring hits, including arginine methyltransferase (CARM1), which was first linked to MS. In particular, the direction of brain expression QTLs for CARM1 and protein kinase MAPK1 enabled us to select BIIB021 and PEITC drugs with the required target modulation. Our study highlights OS-related molecules regulated by functional MS variants that could be targeted by existing drugs as a supplement to the approved disease-modifying treatments.

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.

Application of Genetic Studies to Flow Cytometry Data and Its Impact on Therapeutic Intervention for Autoimmune Disease.

In recent years, systematic genome-wide association studies of quantitative immune cell traits, represented by circulating levels of cell subtypes established by flow cytometry, have revealed numerous association signals, a large fraction of which overlap perfectly with genetic signals associated with autoimmune diseases. By identifying further overlaps with association signals influencing gene expression and cell surface protein levels, it has also been possible, in several cases, to identify causal genes and infer candidate proteins affecting immune cell traits linked to autoimmune disease risk. Overall, these results provide a more detailed picture of how genetic variation affects the human immune system and autoimmune disease risk. They also highlight druggable proteins in the pathogenesis of autoimmune diseases; predict the efficacy and side effects of existing therapies; provide new indications for use for some of them; and optimize the research and development of new, more effective and safer treatments for autoimmune diseases. Here we review the genetic-driven approach that couples systematic multi-parametric flow cytometry with high-resolution genetics and transcriptomics to identify endophenotypes of autoimmune diseases for the development of new therapies.

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.

Evolutionarily Selected Overexpression of the Cytokine BAFF Enhances Mucosal Immune Response Against P. falciparum.

We have previously shown that a variant of the TNFSF13B gene that we called BAFF-var increases the production of the cytokine BAFF, upregulating humoral immunity and increasing the risk for certain autoimmune diseases. In addition, genetic population signatures revealed that BAFF-var was evolutionarily advantageous, most likely by increasing resistance to malaria infection, which is a prime candidate for selective pressure. To evaluate whether the increased soluble BAFF (sBAFF) production confers protection, we experimentally assessed the role of BAFF-var in response to malaria antigens. Lysates of erythrocytes infected with Plasmodium falciparum (iRBCs) or left uninfected (uRBCs, control) were used to treat peripheral blood mononuclear cells (PBMCs) with distinct BAFF genotypes. The PBMCs purified from BAFF-var donors and treated with iRBCs showed different levels of specific cells, immunoglobulins, and cytokines as compared with BAFF-WT. In particular, a relevant differential effect on mucosal immunity B subpopulations have been observed. These findings point to specific immune cells and molecules through which the evolutionary selected BAFF-var may have improved fitness during P. falciparum infection.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy.

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.

Long Noncoding RNAs and Circular RNAs in Autoimmune Diseases.

Immune responses are essential for the clearance of pathogens and the repair of injured tissues; however, if these responses are not properly controlled, autoimmune diseases can occur. Autoimmune diseases (ADs) are a family of disorders characterized by the body's immune response being directed against its own tissues, with consequent chronic inflammation and tissue damage. Despite enormous efforts to identify new drug targets and develop new therapies to prevent and ameliorate AD symptoms, no definitive solutions are available today. Additionally, while substantial progress has been made in drug development for some ADs, most treatments only ameliorate symptoms and, in general, ADs are still incurable. Hundreds of genetic loci have been identified and associated with ADs by genome-wide association studies. However, the whole list of molecular factors that contribute to AD pathogenesis is still unknown. Noncoding (nc)RNAs, such as microRNAs, circular (circ)RNAs, and long noncoding (lnc)RNAs, regulate gene expression at different levels in various diseases, including ADs, and serve as potential drug targets as well as biomarkers for disease progression and response to therapy. In this review, we will focus on the potential roles and genetic regulation of ncRNA in four autoimmune diseases-systemic lupus erythematosus, rheumatoid arthritis, multiple sclerosis, and type 1 diabetes mellitus.

Delta-Globin Gene Expression Is Enhanced in vivo by Interferon Type I.

Beta hemoglobinopathies are widely spread monogenic lethal diseases. Delta-globin gene activation has been proposed as a possible approach for curing these pathologies. The therapeutic potential of delta-globin, the non-alpha component of Hemoglobin A2 (α2δ2; HbA2), has been demonstrated in a mouse model of beta thalassemia, while its anti-sickling effect, comparable to that of gamma globin, was established some time ago. Here we show that the delta-globin mRNA level is considerably increased in a Deoxyribonuclease II-alpha knockout mouse model in which type 1 interferon (interferon beta, IFNb) is activated. IFNb activation in the fetal liver improves the delta-globin mRNA level, while the beta-globin mRNA level is significantly reduced. In addition, we show that HbA2 is significantly increased in patients with multiple sclerosis under type 1 interferon treatment. Our results represent a proof of principle that delta-globin expression can be enhanced through the use of molecules. This observation is potentially interesting in view of a pharmacological approach able to increase the HbA2 level.