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.

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 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor two conditionally distinct, non-coding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR=0.180; MAFEUR=0.008). 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.

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.

A large-scale screening identified in USH2A gene the P3272L founder pathogenic variant explaining familial Usher syndrome in Sardinia, Italy.

BACKGROUND: Usher syndrome (USH) encompasses a group of disorders characterized by congenital sensorineural hearing loss (SNHL) and retinitis pigmentosa (RP). We described the clinical findings, natural history, and molecular analyses of USH patients identified during a large-scale screening to identify quantitative traits related to ocular disorders in the SardiNIA project cohort. METHODS: We identified 3 USH-affected families out of a cohort of 6,148 healthy subjects. 9 subjects presented a pathological phenotype, with SNHL and RP. All patients and their family members underwent a complete ophthalmic examination including best-corrected visual acuity, slit-lamp biomicroscopy, fundoscopy, fundus autofluorescence, spectral-domain optical coherence tomography, and electrophysiological testing. Audiological evaluation was performed with a clinical audiometer. Genotyping was performed using several arrays integrated with whole genome sequence data providing approximately 22 million markers equally distributed for each subject analyzed. Molecular diagnostics focused on analysis of the following candidate genes: MYO7A, USH1C, CDH23, PCDH15, USH1G, CIB2, USH2A, GPR98, DFNB31, CLRN1, and PDZD7. RESULTS: A single missense causal variant in USH2A gene was identified in homozygous status in all patients and in heterozygous status in unaffected parents. The presence of multiple homozygous patients with the same phenotypic severity of the syndromic form suggests that the Sardinian USH phenotype is the result of a founder effect on a specific pathogenic variant related haplotype. The frequency of heterozygotes in general Sardinian population is 1.89. Additionally, to provide new insights into the structure of usherin and the pathological mechanisms caused by small pathogenic in-frame variants, like p.Pro3272Leu, molecular dynamics simulations of native and mutant protein-protein and protein-ligand complexes were performed that predicted a destabilization of the protein with a decrease in the free energy change. CONCLUSIONS: Our results suggest that our approach is effective for the genetic diagnosis of USH. Based on the heterozygous frequency, targeted screening of this variant in the general population and in families at risk or with familial USH can be suggested. This can lead to more accurate molecular diagnosis, better genetic counseling, and improved molecular epidemiology data that are critical for future intervention plans. TRIAL REGISTRATION: We did not perform any health-related interventions for the participants.

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.

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.

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.

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.).

MTHFR, XRCC1 and OGG1 genetic polymorphisms in breast cancer: a case-control study in a population from North Sardinia.

BACKGROUND: Despite conflicting results, considerable evidence suggests the association between single nucleotide polymorphisms in MTHFR, XRCC1 and OGG1 genes and, risk of developing breast cancer. Here a case-control study is reported, including 135 breat cancer patients and 112 healthy women, all representative of Northern Sardinian population. METHODS: Polymerase chain reaction/restriction fragment length polymorphism method was used to determine the genotypes of five polymorphisms: MTHFR C677T (rs1801133) and A1298C (rs1801131), XRCC1 Arg194Trp (rs1799782) and Arg399Gln (rs25487) and OGG1 Ser326Cys (rs1052133). Allelic, genotypic and haplotype association analyses with disease risk and clinicopathological parameters were performed. RESULTS: A nominally significant association with breast cancer risk was observed for MTHFR C677T polymorphism heterozygous genotype in the codominant model (OR: 0.57, 95% CI: 0.32-1.00, p = 0.049) and for Cys/Cys genotype of the OGG1 Ser326Cys polymorphism in the recessive model (OR: 0.23, 95% CI: 0.05-1.11, p = 0.0465). No significant differences were found at genotype-level for A1298C polymorphism of the MTHFR gene and Arg194Trp and Arg399Gln of the XRCC1 gene. Furthermore, the OGG1 and XRCC1 rs25487 polymorphisms were nominally associated with PgR, Her2 status and with sporadic breast cancer, respectively. CONCLUSIONS: Based on genetic characteristics of individuals included in this study, results suggest that MTHFR CT and OGG1 Cys/Cys genotypes have a protective effect that may have an influence on breast cancer risk in a representative Northern Sardinian population.

Kidney size in relation to ageing, gender, renal function, birthweight and chronic kidney disease risk factors in a general population.

BACKGROUND: The relationship of kidney size to ageing, kidney function and kidney disease risk factors is not fully understood. METHODS: Ultrasound length and parenchymal kidney volume were determined from a population-based sample of 3972 Sardinians (age range 18-100 years). We then identified the subset of 2256 'healthy' subjects to define age- and sex-specific reference ranges (2.5-97.5 percentile) of kidney volume. Logistic regression (accounting for family clustering) was used to identify the clinical characteristics associated with abnormally large kidneys or abnormally small kidneys. RESULTS: In the healthy subset, kidney volume and length increased up to the fourth to fifth decade of life followed by a progressive decrease in men, whereas there was a gradual kidney volume decrease throughout the lifespan of women. In the whole sample, independent predictors of lower kidney volume (<2.5 percentile for age and sex) were male sex, low body mass index, short height, low waist:hip ratio and high serum creatinine (SCr); the independent predictors of larger kidney volume (>97.5 percentile for age and sex) were younger age, female sex, diabetes, obesity, high height, high waist:hip ratio and lower SCr. Estimated heritability for kidney volume was 15%, and for length 27%; kidney volume correlated strongly with birthweight. CONCLUSIONS: Overall, in a general healthy population, kidney measures declined with age differently in men and women. The determinants of kidney parenchymal volume include genetic factors and modifiable clinical factors.

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.

A meta-analysis of 120 246 individuals identifies 18 new loci for fibrinogen concentration.

Genome-wide association studies have previously identified 23 genetic loci associated with circulating fibrinogen concentration. These studies used HapMap imputation and did not examine the X-chromosome. 1000 Genomes imputation provides better coverage of uncommon variants, and includes indels. We conducted a genome-wide association analysis of 34 studies imputed to the 1000 Genomes Project reference panel and including ∼120 000 participants of European ancestry (95 806 participants with data on the X-chromosome). Approximately 10.7 million single-nucleotide polymorphisms and 1.2 million indels were examined. We identified 41 genome-wide significant fibrinogen loci; of which, 18 were newly identified. There were no genome-wide significant signals on the X-chromosome. The lead variants of five significant loci were indels. We further identified six additional independent signals, including three rare variants, at two previously characterized loci: FGB and IRF1. Together the 41 loci explain 3% of the variance in plasma fibrinogen concentration.

Adipocytokine correlations with thyroid function and autoimmunity in euthyroid sardinians.

OBJECTIVE: Cytokines release by adipocytes could interact with TSH secretion. We evaluated the relationship between adipocytokines and TSH. We further tested for association of cytokines and thyroid autoimmunity. METHODS: We conducted a cross-sectional study in a community-based sample including 5385 individuals (2964 female) with TSH within the reference range. Subjects who reported taking thyroid medications or drugs that alter thyroid function were excluded. TSH, FT4, adiponectin, leptin, antibody against thyroperoxidase and against thyroglobulin were measured. Linear and logistic regression models were used to test for association. RESULTS: Females had higher adiponectin and leptin level and increased frequency of thyroid antibodies. In multiple regression analysis TSH was directly associated with leptin (ß = 0.003, p = 0.001) and the presence of circulating antibody against thyroperoxidase (ß = 0.315, p < 0.001), but negatively associated with age (ß = -0.012, p < 0.001) and FT4 (ß = -0.359, p < 0.001). Adiponectin, the presence of antibody against thyroglobulin and smoking habit were not associated with TSH levels (p = 0.223, p = 0.174 and p = 0.788, respectively). Logistic regression analysis revealed that higher adiponectin levels were associated with thyroid autoimmunity. CONCLUSIONS: Leptin is positively associated with TSH levels in euthyroid individuals, suggesting an effect of the adipokine on TSH secretion. Our results support the hypothesis that the leptin and pituitary-thyroid axis might interact in the context of energy homeostasis. The effect of adiponectin on thyroid autoimmunity will require more studies.

A genome-wide association scan on the levels of markers of inflammation in Sardinians reveals associations that underpin its complex regulation.

Identifying the genes that influence levels of pro-inflammatory molecules can help to elucidate the mechanisms underlying this process. We first conducted a two-stage genome-wide association scan (GWAS) for the key inflammatory biomarkers Interleukin-6 (IL-6), the general measure of inflammation erythrocyte sedimentation rate (ESR), monocyte chemotactic protein-1 (MCP-1), and high-sensitivity C-reactive protein (hsCRP) in a large cohort of individuals from the founder population of Sardinia. By analysing 731,213 autosomal or X chromosome SNPs and an additional ∼1.9 million imputed variants in 4,694 individuals, we identified several SNPs associated with the selected quantitative trait loci (QTLs) and replicated all the top signals in an independent sample of 1,392 individuals from the same population. Next, to increase power to detect and resolve associations, we further genotyped the whole cohort (6,145 individuals) for 293,875 variants included on the ImmunoChip and MetaboChip custom arrays. Overall, our combined approach led to the identification of 9 genome-wide significant novel independent signals-5 of which were identified only with the custom arrays-and provided confirmatory evidence for an additional 7. Novel signals include: for IL-6, in the ABO gene (rs657152, p = 2.13×10(-29)); for ESR, at the HBB (rs4910472, p = 2.31×10(-11)) and UCN119B/SPPL3 (rs11829037, p = 8.91×10(-10)) loci; for MCP-1, near its receptor CCR2 (rs17141006, p = 7.53×10(-13)) and in CADM3 (rs3026968, p = 7.63×10(-13)); for hsCRP, within the CRP gene (rs3093077, p = 5.73×10(-21)), near DARC (rs3845624, p = 1.43×10(-10)), UNC119B/SPPL3 (rs11829037, p = 1.50×10(-14)), and ICOSLG/AIRE (rs113459440, p = 1.54×10(-08)) loci. Confirmatory evidence was found for IL-6 in the IL-6R gene (rs4129267); for ESR at CR1 (rs12567990) and TMEM57 (rs10903129); for MCP-1 at DARC (rs12075); and for hsCRP at CRP (rs1205), HNF1A (rs225918), and APOC-I (rs4420638). Our results improve the current knowledge of genetic variants underlying inflammation and provide novel clues for the understanding of the molecular mechanisms regulating this complex process.

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 role of genetically predicted serum iron levels on neurodegenerative and cardiovascular traits.

Iron is an essential mineral that supports numerous biological functions. Studies have reported associations between iron dysregulation and certain cardiovascular and neurodegenerative diseases, but the direction of influence is not clear. Our goal was to use computational approaches to better understand the role of genetically predicted iron levels on disease risk. We meta-analyzed genome-wide association study summary statistics for serum iron levels from two cohorts and two previous meta-analyses. We then obtained summary statistics from 11 neurodegenerative, cerebrovascular, cardiovascular or lipid traits to assess global and regional genetic correlation between iron levels and these traits. We used two-sample Mendelian randomization (MR) to estimate causal effects. Sex-stratified analyses were also carried out to identify effects potentially differing by sex. Overall, we identified three significant global correlations between iron levels and (i) coronary heart disease, (ii) triglycerides, and (iii) high-density lipoprotein (HDL) cholesterol levels. A total of 194 genomic regions had significant (after correction for multiple testing) local correlations between iron levels and the 11 tested traits. MR analysis revealed two potential causal relationships, between genetically predicted iron levels and (i) total cholesterol or (ii) non-HDL cholesterol. Sex-stratified analyses suggested a potential protective effect of iron levels on Parkinson's disease risk in females, but not in males. Our results will contribute to a better understanding of the genetic basis underlying iron in cardiovascular and neurological health in aging, and to the eventual identification of new preventive interventions or therapeutic avenues for diseases which affect women and men worldwide.

Implications of disease-modifying therapies for multiple sclerosis on immune cells and response to COVID-19 vaccination.

Introduction: Disease-modifying therapies (DMTs) have been shown to improve disease outcomes in multiple sclerosis (MS) patients. They may also impair the immune response to vaccines, including the SARS-CoV-2 vaccine. However, available data on both the intrinsic immune effects of DMTs and their influence on cellular response to the SARS-CoV-2 vaccine are still incomplete. Methods: Here, we evaluated the immune cell effects of 3 DMTs on the response to mRNA SARS-CoV-2 vaccination by comparing MS patients treated with one specific therapy (fingolimod, dimethyl fumarate, or natalizumab) with both healthy controls and untreated patients. We profiled 23 B-cell traits, 57 T-cell traits, and 10 cytokines, both at basal level and after stimulation with a pool of SARS-CoV-2 spike peptides, in 79 MS patients, treated with DMTs or untreated, and 32 healthy controls. Measurements were made before vaccination and at three time points after immunization. Results and Discussion: MS patients treated with fingolimod showed the strongest immune cell dysregulation characterized by a reduction in all measured lymphocyte cell classes; the patients also had increased immune cell activation at baseline, accompanied by reduced specific immune cell response to the SARS-CoV-2 vaccine. Also, anti-spike specific B cells progressively increased over the three time points after vaccination, even when antibodies measured from the same samples instead showed a decline. Our findings demonstrate that repeated booster vaccinations in MS patients are crucial to overcoming the immune cell impairment caused by DMTs and achieving an immune response to the SARS-CoV-2 vaccine comparable to that of healthy controls.

Interleukin-6 receptor pathways in coronary heart disease: a collaborative meta-analysis of 82 studies.

BACKGROUND: Persistent inflammation has been proposed to contribute to various stages in the pathogenesis of cardiovascular disease. Interleukin-6 receptor (IL6R) signalling propagates downstream inflammation cascades. To assess whether this pathway is causally relevant to coronary heart disease, we studied a functional genetic variant known to affect IL6R signalling. METHODS: In a collaborative meta-analysis, we studied Asp358Ala (rs2228145) in IL6R in relation to a panel of conventional risk factors and inflammation biomarkers in 125,222 participants. We also compared the frequency of Asp358Ala in 51,441 patients with coronary heart disease and in 136,226 controls. To gain insight into possible mechanisms, we assessed Asp358Ala in relation to localised gene expression and to postlipopolysaccharide stimulation of interleukin 6. FINDINGS: The minor allele frequency of Asp358Ala was 39%. Asp358Ala was not associated with lipid concentrations, blood pressure, adiposity, dysglycaemia, or smoking (p value for association per minor allele ≥0·04 for each). By contrast, for every copy of 358Ala inherited, mean concentration of IL6R increased by 34·3% (95% CI 30·4-38·2) and of interleukin 6 by 14·6% (10·7-18·4), and mean concentration of C-reactive protein was reduced by 7·5% (5·9-9·1) and of fibrinogen by 1·0% (0·7-1·3). For every copy of 358Ala inherited, risk of coronary heart disease was reduced by 3·4% (1·8-5·0). Asp358Ala was not related to IL6R mRNA levels or interleukin-6 production in monocytes. INTERPRETATION: Large-scale human genetic and biomarker data are consistent with a causal association between IL6R-related pathways and coronary heart disease. FUNDING: British Heart Foundation; UK Medical Research Council; UK National Institute of Health Research, Cambridge Biomedical Research Centre; BUPA Foundation.