Circulating immune cells are frequently phenotyped by flow cytometry starting from frozen samples. However, cryopreservation can affect marker expression and cell recovery. To understand which antigens are detectable and reliable after sample cryopreservation, we compared 438 antibodies measured on B, T, and NK-enriched cells and monocytes in frozen lympho-monocytes and blood with the corresponding fresh blood samples. Cryopreservation affected the expression of 283 markers in lympho-monocytes and 262 in blood, modifying them by more than 20% with respect to fresh blood. Thus, it is essential to carefully evaluate antibody performance when working with frozen samples. To maximize the usability of our results, make them publicly accessible and ready to visualize, we created a catalogue of marker expression variability before and after freezing, namely FlowCLOc. This catalogue simplifies flow cytometry panel design, reducing time-consuming preliminary tests to select the most appropriate and specific markers for staining both frozen and fresh samples.

Neopterin is a pro-inflammatory molecule upregulated in several diseases; however, its role in pathophysiology is unclear and its genetic regulation is unexplored. We observed that neopterin levels increase during senescence (P-value = 1.88×10-13, beta = 0.96) and positively correlate with age-related neurodegeneration and inflammation markers. The heritability estimation of neopterin variation was 35%. We then conducted a genome-wide association study on 999 Sardinians, identifying two signals in the GTP cyclohydrolase (GCH1) gene that were suggestively associated with neopterin levels. The first signal, led by rs140884539-C (P-value = 7.05×10-08, beta = 0.59), was in strong linkage disequilibrium with variants associated with predisposition to rheumatoid arthritis, decrease in dopamine, increased levels of GCH1 transcript, dopamine metabolites, and galectin-3. The second signal, represented by rs12323905-T (P-value = 8.17×10-08, beta = 0.30), colocalised with GCH1 splicing and Parkinson’s disease signals. Transcriptome analysis of 605 Sardinians showed that the Parkinson’s disease-predisposing variant was significantly associated with an increase in a shorter and inactive form of GCH1, whose presence is predicted to reduce the GCH1 decamer stability. The GCH1 homo-decamer regulates neopterin and tetrahydrobiopterin production, a cofactor required for the synthesis of dopamine and serotonin. Our data motivate experimental work to test whether modulating GCH1 expression or isoform ratio alters dopaminergic function in Parkinson’s disease models.

BackgroundThe immunoglobulin level following vaccination against SARS-CoV-2 results from a multifaceted immunological process involving cells and molecules that determine its efficacy and protect us against severe infection outcomes. The observed heterogeneity in the immune response to vaccination is partly attributable to host genetic factors; however, this genetic contribution has only been partially explored.MethodsTo elucidate the mechanisms underlying the antibody response elicited by the Pfizer-BioNTech vaccine, we conducted a genome-wide association study on anti-spike immunoglobulin G levels measured in 1,968 Italian individuals who received two doses of the vaccine, selected from a larger cohort of 7,169 volunteers characterized for 8 million genetic variants. Sex, age, body mass index, smoking habit, and time elapsed between vaccine administration and blood draw were accounted as covariates in the linear regression model.ResultsWe identified a novel signal of association on chromosome 3 in the intronic region of the CMTM8 gene and confirmed one previously identified at the HLA locus close to the HLA-B gene. The lead SNP in the CMTM8 gene, rs7643677 (p-value = 2.095×10-8), is associated with anti-S IgG levels and with the expression level of CD66b on granulocytic/polymorphonuclear myeloid-derived suppressor cells.ConclusionsThese findings support a role for CMTM8 in regulating the suppressive activity of specific immune cells, and suggest a potential interplay among genetic, humoral, and cellular mechanisms underlying the immune response to SARS-CoV-2 vaccination.

BACKGROUND:Growth differentiation factor 15 (GDF-15) levels are emerging as a candidate biomarker of aging. The present study aimed to: (1) characterize the association of GDF-15 with the continuum of arterial stiffening, assessed as carotid-femoral pulse wave velocity, as age increases; (2) determine the predictive role of serum GDF-15 levels on mortality; and (3) identify genetic determinants of serum GDF-15 levels. METHODS AND RESULTS:Serum levels of GDF-15 and established cardiovascular risk factors, including pulse wave velocity, were assessed in a large (4736 individual) Sardinian population. Serum levels of GDF-15, which can be reliably measured repeatedly over time, increase with age; are associated with a stiffer aorta; "mediate" a large proportion of the age-associated increase in arterial stiffness; pose risks because of their association with greater mortality; and are significantly associated with the variant rs11549407, which causes thalassemia major in homozygosity. CONCLUSIONS:Because of its consistent ability to predict functional and clinical outcomes, including all-cause mortality, we conclude that GDF-15 serum levels serve as a robust biomarker for the continuum from health to the emergence of clinical disease during aging and, subsequently, to the likelihood of mortality.

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

IntroductionDisease-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.MethodsHere, 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 DiscussionMS 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.

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.

During vertebrate somitogenesis, various transcriptional factors function coordinately to determine the position of the somite boundary. Previously, we reported on the signaling crosstalk that occurs between two major transcription factors involved in somitogenesis, Tbx6 and mespb/mesp2. These factors synergistically activated the expression of a downstream gene, bowline/Ripply2, which is essential for precise formation of the somite boundary. However, the molecular mechanism underlying this synergistic effect remains unclear. In this report, we found that the Tbx6 and mespb proteins interacted physically with each other. Pulldown assays with various deletion mutants of these proteins identified the essential domains for this physical interaction. Finally, we found that interference with the physical interaction by a dominant-negative form of mespb, mespbΔDBD, abrogated the expression of the bowline gene during Xenopus somitogenesis. These results indicate that the appropriate expression of bowline/Ripply2 is regulated by a direct interaction between the Tbx6 and mespb proteins during Xenopus somitogenesis.

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

Valeria Orrù§, Maristella Steri§, Carlo Sidore§, Michele Marongiu, Valentina Serra, Stefania Olla, Gabriella Sole, Sandra Lai, Mariano Dei, Antonella Mulas, Francesca Virdis, Maria Grazia Piras, Monia Lobina, Mara Marongiu, Maristella Pitzalis, Francesca Deidda, Annalisa Loizedda, Stefano Onano, Magdalena Zoledziewska, Marcella Devoto, Matteo Floris, Mauro Pala, David Schlessinger, Edoardo Fiorillo, Francesco Cucca.

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic 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.)

Identifying functional non-coding variants can enhance genome interpretation and inform novel genetic risk factors. We used whole genomes and peripheral white blood cell transcriptomes from 624 Sardinian individuals to identify non-coding variants that contribute to population, family, and individual differences in transcript abundance. We identified 21,183 independent expression quantitative trait loci (eQTLs) and 6,768 independent splicing quantitative trait loci (sQTLs) influencing 73 and 41% of all tested genes. When we compared Sardinian eQTLs to those previously identified in Europe, we identified differentiated eQTLs at genes involved in malarial resistance and multiple sclerosis, reflecting the long-term epidemiological history of the island’s population. Taking advantage of pedigree data for the population sample, we identify segregating patterns of outlier gene expression and allelic imbalance in 61 Sardinian trios. We identified 809 expression outliers (median z-score of 2.97) averaging 13.3 genes with outlier expression per individual. We then connected these outlier expression events to rare non-coding variants. Our results provide new insight into the effects of non-coding variants and their relationship to population history, traits and individual genetic risk.

Background and aims: Increased carotid artery intima-media thickness (IMT) and the presence of plaques have been shown to be predictors of cardiovascular disease. The cardiovascular risk in patients with overt thyroid diseases is related to increased risk of atherosclerosis, but there has been no clear evidence about subclinical disorders. We have assessed whether subclinical thyroid dysfunction is associated with arterial thickening and plaque.Methods and results: The SardiNIA study is a population-based survey on the Italian island of Sardinia. We reviewed data from 5815 subjects (aged 14-102 years), none of whom had overt hyperthyroidism or hypothyroidism or was taking thyroid medication. Serum thyrotropin (TSH), free thyroxine, together with carotid ultrasound IMT and the presence of common carotid plaques were analysed in all subjects. Possible association of IMT and carotid plaques with thyroid parameters was evaluated by univariate and multivariate analyses. IMT was significantly associated with age, sex, smoking, low density lipoprotein cholesterol (LDL), high density lipoprotein cholesterol, pulse pressure (PP), history of arterial hypertension, diabetes, and previous cardiovascular events (p = 0.001 or lower, R-2 = 0.47). Carotid plaques were predicted by age, sex, LDL, PP, history of diabetes, previous cardiovascular events, and the use of statins (p = 0.029 or lower). Thyroid hormone was not predictive of carotid atherosclerosis when adjusted for confounders.Conclusion: Thyroid hormone is not associated with increased IMT or with the presence of carotid artery plaque. Our data do not support the idea that treating subclinical disorders might help to prevent arterial remodelling or carotid atherosclerosis. (C) 2015 The Italian Society of Diabetology, the Italian Society for the Study of Atherosclerosis, the Italian Society of Human Nutrition, and the Department of Clinical Medicine and Surgery, Federico II University. Published by Elsevier B.V. All rights reserved.

Objectives. Lower levels of serum brain derived neurotrophic factor (BDNF) is one of the best known biomarkers of depression. To identify genetic variants associated with serum BDNF, we tested the Val66Met (rs6265) functional variant and conducted a genome-wide association scan (GWAS). Methods. In a community-based sample (N = 2054; aged 19-101, M = 51, SD = 15) from Sardinia, Italy, we measured serum BDNF concentration and conducted a GWAS. Results. We estimated the heritability of serum BDNF to be 0.48 from sib-pairs. There was no association between serum BDNF and Val66Met in the SardiNIA sample and in a meta-analysis of published studies (k = 13 studies, total n = 4727, P = 0.92). Although no genome-wide significant associations were identified, some evidence of association was found in the BDNF gene (rs11030102, P = 0.001) and at two loci (rs7170215, P = 4.8 x 10 (5) and rs11073742 P = 1.2 x 10 (5)) near and within NTRK3 gene, a neurotrophic tyrosine kinase receptor. Conclusions. Our study and meta-analysis of the literature indicate that the BDNF Val66Met variant is not associated with serum BDNF, but other variants in the BDNF and NTRK3 genes might regulate the level of serum BDNF.

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.

Thyroid hormone is essential for normal metabolism and development, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over their life span. In addition, even mild alterations in thyroid function are associated with weight changes, atrial fibrillation, osteoporosis, and psychiatric disorders. To identify novel variants underlying thyroid function, we performed a large meta-analysis of genome-wide association studies for serum levels of the highly heritable thyroid function markers TSH and FT4, in up to 26,420 and 17,520 euthyroid subjects, respectively. Here we report 26 independent associations, including several novel loci for TSH (PDE10A, VEGFA, IGFBP5, NFIA, SOX9, PRDM11, FGF7, INSR, ABO, MIR1179, NRG1, MBIP, ITPK1, SASH1, GLIS3) and FT4 (LHX3, FOXE1, AADAT, NETO1/FBXO15, LPCAT2/CAPNS2). Notably, only limited overlap was detected between TSH and FT4 associated signals, in spite of the feedback regulation of their circulating levels by the hypothalamic-pituitary-thyroid axis. Five of the reported loci (PDE8B, PDE10A, MAF/LOC440389, NETO1/FBXO15, and LPCAT2/CAPNS2) show strong gender-specific differences, which offer clues for the known sexual dimorphism in thyroid function and related pathologies. Importantly, the TSH-associated loci contribute not only to variation within the normal range, but also to TSH values outside the reference range, suggesting that they may be involved in thyroid dysfunction. Overall, our findings explain, respectively, 5.64% and 2.30% of total TSH and FT4 trait variance, and they improve the current knowledge of the regulation of hypothalamic-pituitary-thyroid axis function and the consequences of genetic variation for hypo- or hyperthyroidism.

Anna Köttgen and colleagues report genome-wide association studies for serum urate in over 140,000 individuals from the Global Urate Genetics Consortium (GUGC). They identify 18 loci newly associated with serum urate concentrations and confirm 10 known loci, characterize their associations with gout and include a network analysis suggesting a role for inhibins-activins pathways in regulating urate homeostasis. Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >140,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (18 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.