Maristella Steri

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.

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.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.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.Large-scale human genetic and biomarker data are consistent with a causal association between IL6R-related pathways and coronary heart disease.British Heart Foundation; UK Medical Research Council; UK National Institute of Health Research, Cambridge Biomedical Research Centre; BUPA Foundation.

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.

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

We report ∼17.6 million genetic variants from whole-genome sequencing of 2,120 Sardinians; 22% are absent from previous sequencing-based compilations and are enriched for predicted functional consequences. Furthermore, ∼76,000 variants common in our sample (frequency >5%) are rare elsewhere (<0.5% in the 1000 Genomes Project). We assessed the impact of these variants on circulating lipid levels and five inflammatory biomarkers. We observe 14 signals, including 2 major new loci, for lipid levels and 19 signals, including 2 new loci, for inflammatory markers. The new associations would have been missed in analyses based on 1000 Genomes Project data, underlining the advantages of large-scale sequencing in this founder population.

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.

Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10(-8)) loci, some including known iron-related genes (HFE, SLC40A1, TF, TFR2, TFRC, TMPRSS6) and others novel (ABO, ARNTL, FADS2, NAT2, TEX14). SNPs at ARNTL, TF, and TFR2 affect iron markers in HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease.

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

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.

Estimates of the heritability of plasma fibrinogen concentration, an established predictor of cardiovascular disease, range from 34% to 50%. Genetic variants so far identified by genome-wide association studies explain only a small proportion (<2%) of its variation.We conducted a meta-analysis of 28 genome-wide association studies including >90 000 subjects of European ancestry, the first genome-wide association meta-analysis of fibrinogen levels in 7 studies in blacks totaling 8289 samples, and a genome-wide association study in Hispanics totaling 1366 samples. Evaluation for association of single-nucleotide polymorphisms with clinical outcomes included a total of 40 695 cases and 85 582 controls for coronary artery disease, 4752 cases and 24 030 controls for stroke, and 3208 cases and 46 167 controls for venous thromboembolism. Overall, we identified 24 genome-wide significant (P<5×10(-8)) independent signals in 23 loci, including 15 novel associations, together accounting for 3.7% of plasma fibrinogen variation. Gene-set enrichment analysis highlighted key roles in fibrinogen regulation for the 3 structural fibrinogen genes and pathways related to inflammation, adipocytokines, and thyrotrophin-releasing hormone signaling. Whereas lead single-nucleotide polymorphisms in a few loci were significantly associated with coronary artery disease, the combined effect of all 24 fibrinogen-associated lead single-nucleotide polymorphisms was not significant for coronary artery disease, stroke, or venous thromboembolism.We identify 23 robustly associated fibrinogen loci, 15 of which are new. Clinical outcome analysis of these loci does not support a causal relationship between circulating levels of fibrinogen and coronary artery disease, stroke, or venous thromboembolism.

Genome Wide Association Studies (GWAS) have mapped thousands of genetic variants associated with complex disease risk and regulating quantitative traits, thus exploiting an unprecedented high‐resolution genetic characterization of the human genome. A small fraction (3.7%) of the identified associations is located in untranslated regions (UTRs), and the molecular mechanism has been elucidated for few of them. Genetic variations at UTRs may modify regulatory elements affecting the interaction of the UTRs with proteins and microRNAs. The overall functional consequences include modulation of messenger RNA (mRNA) transcription, secondary structure, stability, localization, translation, and access to regulators like microRNAs (miRNAs) and RNA‐binding proteins (RBPs). Alterations of these regulatory mechanisms are known to modify molecular pathways and cellular processes, potentially leading to disease processes. Here, we analyze some examples of genetic risk variants mapping in the UTR regulatory elements. We describe a recently identified genetic variant localized in the 3′UTR of the TNFSF13B gene, associated with autoimmunity risk and responsible of an increased stability and translation of TNFSF13B mRNA. We discuss how the correct use and interpretation of public GWAS repositories could lead to a better understanding of etiopathogenetic mechanisms and the generation of robust biological hypothesis as starting point for further functional studies. This article is categorized under: RNA Structure and Dynamics &gt; RNA Structure, Dynamics and Chemistry RNA Evolution and Genomics &gt; Computational Analyses of RNA RNA in Disease and Development &gt; RNA in Disease

We report sequencing-based whole-genome association analyses to evaluate the impact of rare and founder variants on stature in 6,307 individuals on the island of Sardinia. We identify two variants with large effects. One variant, which introduces a stop codon in the GHR gene, is relatively frequent in Sardinia (0.87% versus <0.01% elsewhere) and in the homozygous state causes Laron syndrome involving short stature. We find that this variant reduces height in heterozygotes by an average of 4.2 cm (-0.64 s.d.). The other variant, in the imprinted KCNQ1 gene (minor allele frequency (MAF) = 7.7% in Sardinia versus <1% elsewhere) reduces height by an average of 1.83 cm (-0.31 s.d.) when maternally inherited. Additionally, polygenic scores indicate that known height-decreasing alleles are at systematically higher frequencies in Sardinians than would be expected by genetic drift. The findings are consistent with selection for shorter stature in Sardinia and a suggestive human example of the proposed 'island effect' reducing the size of large mammals.

The utility of genotype imputation in genome-wide association studies is increasing as progressively larger reference panels are improved and expanded through whole-genome sequencing. Developing general guidelines for optimally cost-effective imputation, however, requires evaluation of performance issues that include the relative utility of study-specific compared with general/multipopulation reference panels; genotyping with various array scaffolds; effects of different ethnic backgrounds; and assessment of ranges of allele frequencies. Here we compared the effectiveness of study-specific reference panels to the commonly used 1000 Genomes Project (1000G) reference panels in the isolated Sardinian population and in cohorts of European ancestry including samples from Minnesota (USA). We also examined different combinations of genome-wide and custom arrays for baseline genotypes. In Sardinians, the study-specific reference panel provided better coverage and genotype imputation accuracy than the 1000G panels and other large European panels. In fact, even gene-centered custom arrays (interrogating ~200 000 variants) provided highly informative content across the entire genome. Gain in accuracy was also observed for Minnesotans using the study-specific reference panel, although the increase was smaller than in Sardinians, especially for rare variants. Notably, a combined panel including both study-specific and 1000G reference panels improved imputation accuracy only in the Minnesota sample, and only at rare sites. Finally, we found that when imputation is performed with a study-specific reference panel, cutoffs different from the standard thresholds of MACH-Rsq and IMPUTE-INFO metrics should be used to efficiently filter badly imputed rare variants. This study thus provides general guidelines for researchers planning large-scale genetic studies.

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.

Electrocardiographic PR interval measures atrio-ventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. Our genome-wide association study of over 92,000 European-descent individuals identifies 44 PR interval loci (34 novel). Examination of these loci reveals known and previously not-yet-reported biological processes involved in cardiac atrial electrical activity. Genes in these loci are over-represented in cardiac disease processes including heart block and atrial fibrillation. Variants in over half of the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with missense variants. Six additional loci were identified either by meta-analysis of ~105,000 African and European-descent individuals and/or by pleiotropic analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation. These findings implicate developmental pathways, and identify transcription factors, ion-channel genes, and cell-junction/cell-signaling proteins in atrio-ventricular conduction, identifying potential targets for drug development.

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

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.

The population of the Mediterranean island of Sardinia has made important contributions to genome-wide association studies of complex disease traits and, based on ancient DNA studies of mainland Europe, Sardinia is hypothesized to be a unique refuge for early Neolithic ancestry. To provide new insights on the genetic history of this flagship population, we analyzed 3,514 whole-genome sequenced individuals from Sardinia. Sardinian samples show elevated levels of shared ancestry with Basque individuals, especially samples from the more historically isolated regions of Sardinia. Our analysis also uniquely illuminates how levels of genetic similarity with mainland ancient DNA samples varies subtly across the island. Together, our results indicate that within-island substructure and sex-biased processes have substantially impacted the genetic history of Sardinia. These results give new insight into the demography of ancestral Sardinians and help further the understanding of sharing of disease risk alleles between Sardinia and mainland populations.

We report genome-wide association study results for the levels of A1, A2 and fetal hemoglobins, analyzed for the first time concurrently. Integrating high-density array genotyping and whole-genome sequencing in a large general population cohort from Sardinia, we detected 23 associations at 10 loci. Five signals are due to variants at previously undetected loci: MPHOSPH9, PLTP-PCIF1, ZFPM1 (FOG1), NFIX and CCND3. Among the signals at known loci, ten are new lead variants and four are new independent signals. Half of all variants also showed pleiotropic associations with different hemoglobins, which further corroborated some of the detected associations and identified features of coordinated hemoglobin species production.

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

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

Francesco Cucca, Stephen Montgomery and colleagues identify regulatory variants that influence gene expression and splicing using whole-genome and transcriptome sequence data from 624 Sardinians. They find a high frequency of splicing and expression quantitative trait loci near genes involved in malarial resistance and multiple sclerosis. 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.

Abstract 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 (&amp;lt;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 (&amp;gt;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.

In systemic lupus erythematosus (SLE), perturbed immunoregulation underpins a pathogenic imbalance between regulatory and effector CD4+ T-cell activity. However, to date, the characterisation of the CD4+ regulatory T cell (Treg) compartment in SLE has yielded conflicting results. Here we show that patients have an increased frequency of CD4+FOXP3+ cells in circulation owing to a specific expansion of thymically-derived FOXP3+HELIOS+ Tregs with a demethylated FOXP3 Treg-specific demethylated region. We found that the Treg expansion was strongly associated with markers of recent immune activation, including PD-1, plasma concentrations of IL-2 and the type I interferon biomarker soluble SIGLEC-1. Since the expression of the negative T-cell signalling molecule PTPN22 is increased and a marker of poor prognosis in SLE, we tested the influence of its missense risk allele Trp620 (rs2476601C>T) on Treg frequency. Trp620 was reproducibly associated with increased frequencies of thymically-derived Tregs in blood, and increased PD-1 expression on both Tregs and effector T cells (Teffs). Our results support the hypothesis that FOXP3+ Tregs are increased in SLE patients as a consequence of a compensatory mechanism in an attempt to regulate pathogenic autoreactive Teff activity. We suggest that restoration of IL-2-mediated homeostatic regulation of FOXP3+ Tregs by IL-2 administration could prevent disease flares rather than treating at the height of a disease flare. Moreover, stimulation of PD-1 with specific agonists, perhaps in combination with low-dose IL-2, could be an effective therapeutic strategy in autoimmune disease and in other immune disorders.

At present, specialized companies offering genetic testing services without the involvement of clinicians are growing; this development is a direct consequence of the significant decrease in genotyping and sequencing costs. Online companies offer predictions about the risk of developing complex diseases during one’s life course, and they offer suggestions for personal lifestyle. Several companies have been created that provide nutrigenetics services; these companies suggest dietary indications—a central issue in the prevention and etiopathogenesis of specific diseases—based on one’s personal genetic background. Dietary patterns are defined on the basis of a limited set of genetic markers. In this article, we analyze the online nutrigenetics services offered by 45 companies worldwide, to obtain an overall picture of the costs, the types of nutritional traits considered and the level of scientific precision of the services proposed. Our analysis clearly highlights the need for specific guidelines, in order to ensure a set of minimum quality standards for the nutrigenetics services offered to the customer.

An increasing number of genome-wide association (GWA) studies are now using the higher resolution 1000 Genomes Project reference panel (1000G) for imputation, with the expectation that 1000G imputation will lead to the discovery of additional associated loci when compared to HapMap imputation. In order to assess the improvement of 1000G over HapMap imputation in identifying associated loci, we compared the results of GWA studies of circulating fibrinogen based on the two reference panels. Using both HapMap and 1000G imputation we performed a meta-analysis of 22 studies comprising the same 91,953 individuals. We identified six additional signals using 1000G imputation, while 29 loci were associated using both HapMap and 1000G imputation. One locus identified using HapMap imputation was not significant using 1000G imputation. The genome-wide significance threshold of 5×10-8 is based on the number of independent statistical tests using HapMap imputation, and 1000G imputation may lead to further independent tests that should be corrected for. When using a stricter Bonferroni correction for the 1000G GWA study (P-value < 2.5×10-8), the number of loci significant only using HapMap imputation increased to 4 while the number of loci significant only using 1000G decreased to 5. In conclusion, 1000G imputation enabled the identification of 20% more loci than HapMap imputation, although the advantage of 1000G imputation became less clear when a stricter Bonferroni correction was used. More generally, our results provide insights that are applicable to the implementation of other dense reference panels that are under development.

Objective Thyroid hormone influences lipoprotein metabolism. The role of menopausal status in this association has not been extensively studied. The aim of the present study is to evaluate the association between lipid parameters and mild elevations of thyrotropin (TSH), and whether menopause influences this relationship. Study design A cross-sectional study was conducted with a sample of 2,914 women (aged 14–102 years) from the SardiNIA study. Main outcome measures The association of TSH with blood lipid levels was examined using regression analyses, according to menopausal status. Results Postmenopausal women had lower serum TSH concentrations and higher levels of total cholesterol, low-density lipoprotein cholesterol (LDLc), high-density lipoprotein cholesterol (HDLc), and triglycerides than did premenopausal women (p = 0.001 or less for all). In premenopausal women, after adjusting for the confounders age, BMI, smoking, insulin and glycaemia, TSH showed a direct relation to the levels of total cholesterol (β = 0.046, p = 0.010), LDLc (β = 0.044, p = 0.016) and triglycerides (β = 0.085, p < 0.001), but no association with HDLc level. In the postmenopausal group, TSH was directly associated only with triglyceride levels (β = 0.103, p = 0.014). Conclusions The association between mild elevation of TSH and lipid levels is influenced by menopausal status. Further research is needed to clarify this finding.

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

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.

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.

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.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).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.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.Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

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

Plasma fibrinogen is an acute phase protein playing an important role in the blood coagulation cascade having strong associations with smoking, alcohol consumption and body mass index (BMI). Genome-wide association studies (GWAS) have identified a variety of gene regions associated with elevated plasma fibrinogen concentrations. However, little is yet known about how associations between environmental factors and fibrinogen might be modified by genetic variation. Therefore, we conducted large-scale meta-analyses of genome-wide interaction studies to identify possible interactions of genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentration. The present study included 80,607 subjects of European ancestry from 22 studies. Genome-wide interaction analyses were performed separately in each study for about 2.6 million single nucleotide polymorphisms (SNPs) across the 22 autosomal chromosomes. For each SNP and risk factor, we performed a linear regression under an additive genetic model including an interaction term between SNP and risk factor. Interaction estimates were meta-analysed using a fixed-effects model. No genome-wide significant interaction with smoking status, alcohol consumption or BMI was observed in the meta-analyses. The most suggestive interaction was found for smoking and rs10519203, located in the LOC123688 region on chromosome 15, with a p value of 6.2 × 10(-8). This large genome-wide interaction study including 80,607 participants found no strong evidence of interaction between genetic variants and smoking status, alcohol consumption or BMI on fibrinogen concentrations. Further studies are needed to yield deeper insight in the interplay between environmental factors and gene variants on the regulation of fibrinogen concentrations.

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

Abstract The population of the Mediterranean island of Sardinia has made important contributions to genome-wide association studies of traits and diseases. The history of the Sardinian population has also been the focus of much research, and in recent ancient DNA (aDNA) studies, Sardinia has provided unique insight into the peopling of Europe and the spread of agriculture. In this study, we analyze whole-genome sequences of 3,514 Sardinians to address hypotheses regarding the founding of Sardinia and its relation to the peopling of Europe, including examining fine-scale substructure, population size history, and signals of admixture. We find the population of the mountainous Gennargentu region shows elevated genetic isolation with higher levels of ancestry associated with mainland Neolithic farmers and depleted ancestry associated with more recent Bronze Age Steppe migrations on the mainland. Notably, the Gennargentu region also has elevated levels of pre-Neolithic hunter-gatherer ancestry and increased affinity to Basque populations. Further, allele sharing with pre-Neolithic and Neolithic mainland populations is larger on the X chromosome compared to the autosome, providing evidence for a sex-biased demographic history in Sardinia. These results give new insight to the demography of ancestral Sardinians and help further the understanding of sharing of disease risk alleles between Sardinia and mainland populations.

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.

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

Nature Communications 5: Artice number: 4926 (2014); Published 29 October 2014; Updated 30 March 2015. The original version of this Article contained an error in the spelling of the author Jonas Hälldin, which was incorrectly given as Jonas Häldin. This has now been corrected in both the PDF and HTML versions of the Article.

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.

We have previously shown that a variant of the TNFSF13B gene that we called BAFF-var, increases the production of the cytokine BAFF, up-regulating 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 assessed experimentally 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 PBMCs and B cells 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 particularly 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.

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.

Inherited platelets disorders can be severe, especially after trauma or surgical procedures in some monogenic disorders, as in Bernard–Soulier syndrome [BSS; Mendelian Inheritance in Man (MIM) #231200]. BSS is a rare autosomal recessive macrothrombocytopenia (incidence of about one per million). Its hallmark is a defective adhesion of platelets to the sub-endothelium, resulting from quantitative or qualitative defects in the glycoprotein Ib (GPIb)-IX-V complex, a platelet receptor for von Willebrand Factor (VWF), which is composed of four subunits: GPIbα, GPIbβ, GPIX and GPV.1 Laboratory diagnosis is based on prolonged bleeding time, moderate-to-severe thrombocytopenia (platelet count typically ranges from 20 to 100 × 109/l), giant platelets and deficient ristocetin-dependent platelet agglutination.2 Very little is known about the biochemical and clinical features of heterozygous carriers of the mutations causing BSS, and about the impact in general population individuals of variation in genes encoding the GPIb-IX-V complex when present in heterozygosity. In fact, family members with only one mutated allele are generally asymptomatic, with sub-normal platelet count, slightly enlarged platelets and marginally reduced levels of glycoproteins expression. Here, to dissect the impact of genetic variability on platelet count, a sequencing-based whole-genome association study was performed in 6528 volunteers included in the SardiNIA general population cohort.3 Six signals were identified (Table SI, Data S1 for description), including a novel non-synonymous variant [22:19711445:C/T; minor allele frequency (MAF) = 0·0045; P = 1·172 × 10−16], mapping in the second exon (c.C79T, p.P27S) of the GP1BB gene (Figure S1). Completely independent of previously reported associations in the same genomic region (Data S1), p.P27S is Sardinian-specific, being completely missing in large sequencing datasets such as 1000 Genomes Project,4 Genome of the Netherlands (GoNL) project,5 Genome Aggregation Database (gnomAD),6 the Exome Sequencing Project in the National Heart, Lung, and Blood Institute’s (NHLBI’s) Trans-Omics for Precision Medicine (TOPMed) programme.7 No homozygous and 57 carriers for the rare 22:19711445-T allele were found. The platelet count in wild-type homozygous was 242·87 ± 117·05 × 109/l (mean ± 1·96 × SD), whereas in p.P27S carriers was 174·17 ± 91·51 × 109/l, corresponding to a reduction of 70·13 × 109/l for each copy of the minor allele (Fig 1A). With this large effect, the novel founder mutation explains ~1·05% of phenotypic variance for platelet count, representing the largest phenotypic effect among all the independent variants reported to date in the genome-wide association studies (GWAS) Catalogue (Tables SII and SIII).8 Moreover, in a subset of 2000 individuals, whose mean platelet volume was measured, this variant was associated with notably larger platelets (P = 2·13 × 10−10), consistently with evidence of morphologically enlarged platelets in patients with BSS (Fig 1B). To assess platelet functionality, a seven-colour flow cytometry panel (Table SIV) was set up in 24 of 57 p.P27S carriers (42·1%) and in an equal number of matched unrelated controls. Monoclonal antibodies directed against the GPIIb-IIIa complex (CD41a and CD61), and the VWF receptor complex (CD42a and CD42b) were used to investigate the basal receptor expression in resting platelets. The p.P27S carriers showed increased levels of GPIIb (CD41a, +22·36%; P = 1·61 × 10−4; N = 48) and GPIIIa (CD61, +16·20%; P = 6·61 × 104; N = 48) a typical finding in the presence of enlarged platelets (Figure S2). The expression of GPIX and GPIbα glycoproteins and their correct assembly into the GPIb-IX-V complex are known to be impaired by a defective GPIbβ peptide9. Indeed, despite carrying only one mutated allele, the p.P27S heterozygous showed appreciably lower basal expression levels of both GPIX (−24·69%, P = 2·66 × 10−6, N = 46; Fig 1C) and GPIbα (−26·51%, P = 3·66 × 10−8, N = 48; Fig 1D), and consequently less of the entire complex, compared to controls. This is far more than the normal expression levels of GPIX and GPIbα in carriers of other missense mutations in GP1BB, as recently reported.10 Pre-activation and reactivity changes in p.P27S platelets were investigated after exposure to the agonist adenosine diphosphate (ADP). Indeed, activated αIIbβ3 was prominently induced in the p.P27S carriers, as shown by the extent of procaspase-activating compound 1 (PAC-1) binding to resting and activated platelets (+41·94%, P = 4·84 × 10−3, N = 48, Fig 1E). Notably, no variation in the response of platelets after ADP stimulation was recently reported in patients with BSS and carriers.9 Remarkably, platelet reactivity turned out to be differentially regulated: no changes were observed in surface exposure of neo P-selectin (CD62P, +35·22%, P = 0·138, N = 48; Fig 1F) and neo granulophysin (CD63, +1·86%, P = 0·658, N = 46; Fig 1G), markers of granule content release. The unique functional effects of the p.P27S led us to examine its possible consequences on the molecular structure and conformational changes of GPIbβ by molecular modelling analysis based on the X-ray crystal structure.11 Proline–Serine substitution falls in the leucine-rich repeat N-terminal (LRRNT) domain of the 206 amino acid long protein encoded by GP1BB (Fig 2A,B). Proline residues are expected to be disruptive of structure; and indeed, in that highly conserved region and close to cysteine residues involved in the Cys26-Cys32 disulphide bridge, p.P27S could thus modify the stability and consequently the conformation of GPIbβ. To test this hypothesis, we first performed in silico molecular dynamic simulations, observing an increased conformational mobility of the amino acid backbone close to p.P27S (Figure S3), suggesting the instability of the GPIbβ glycoprotein in accordance with the observed reduction in the expression of GPIX and GPIbα. Strikingly, a greater fluctuation of the amino acids in loop 2 of the p.P27S protein was also recorded, as indicated by root-mean-square-fluctuation (Fig 2C). In summary, all typical findings of macrothrombocytopenias (i.e. BSS) were observed in the p.P27S obligate carriers characterised in the present study: low levels of large platelets and low expression of GPIX and GPIbα glycoproteins, as shown by flow cytometry. As one might anticipate, the most severe cases are caused by deletions and nonsense mutations, but some missense mutations are disabling enough to be clinically significant. In one of the reported cases,12 a charge difference is introduced (p.Asn89Asp); in the other,13 as in this case, the Proline residue is replaced (p.Pro27Leu), which is expected to disrupt secondary structure in the protein. That p.P27S influences conformational changes and stability of GPIbβ, in turn affecting GPIb-IX-V complex function, is further clearly supported by the in silico molecular dynamic analyses. Noteworthy, a critical interaction of GPIbβ with GPIX involves N-terminal residues 15 through 32 of GPIbβ, precisely including Proline 27.14 According to Hardy–Weinberg expectation, at least four p.P27S homozygous individuals, most likely with BSS, are expected in Sardinia, but none have been reported to date: this may suggest that BSS is likely underdiagnosed in Sardinia, consistent with other reports.15 Thus, clinicians should be aware of the novel p.P27S mutation in the molecular characterisation of Sardinian-origin patients with a clinical picture of platelet macrocytosis and platelet count of <100 × 109/l. We thank all the volunteers who generously participated in this study; we are grateful to Mr Mario Lovicu and Mr Nazario Olla for the logistic support provided and helpful suggestions. Supported by contracts N01-AG-1-2109 and HHSN271201100005C from the Intramural Research Program of the National Institute on Aging, National Institutes of Health (NIH). The authors declare no competing interests. Fabio Busonero, Maristella Pitzalis and Francesca Deidda collected samples and extracted genomic DNA from blood; Fabio Busonero, Andrea Maschio, Antonella Mulas and Magdalena Zoledziewska performed genotyping; Fabio Busonero and Andrea Maschio, performed sequencing; Valeria Orrù and Edoardo Fiorillo designed flow cytometric panels; Valeria Orrù, Edoardo Fiorillo and Sandra Lai performed cytometric analysis; Maristella Steri, Gabriella Sole, Carlo Sidore, Matteo Floris, Mauro Pala, Paola Forabosco, Marco Masala and Serena Sanna performed statistical analyses; Maristella Steri, Michele Marongiu and Serena Sanna performed bioinformatic analyses; Fabio Busonero and Serena Sanna performed region-specific analysis and selected candidate genes; Stefania Olla performed in silico analyses; Isadora Asunis and Cristian Antonio Caria performed functional evaluation; Francesco Cucca, Goncalo R. Abecasis and David Schlessinger provided funds and supervised the work; Susanna Barella provided clinical support; Fabio Busonero, Maristella Steri and Francesco Cucca wrote the paper; Valeria Orrù, Gabriella Sole, Stefania Olla, Andrea Maschio, Carlo Sidore, Matteo Floris, David Schlessinger and Serena Sanna revised the paper. All authors read the paper and contributed to its final form. The URLs for data presented herein are as follows: Online Mendelian Inheritance in Man, https://www.omim.org/; SardiNIA Project, https://sardinia.irp.nia.nih.gov/; 1000 Genomes Project data repository, ftp://ftp.1000genomes.ebi.ac.uk/; GoNL, Genome of the Netherlands, http://www.nlgenome.nl/; GnomAD, http://gnomad.broadinstitute.org/; Exome Sequencing Project, https://esp.gs.washington.edu/drupal/; NHLBI TOPMed Programme, https://www.nhlbiwgs.org/; GWAS catalogue, https://www.ebi.ac.uk/gwas/. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

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.

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.

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

ABSTRACT Electrocardiographic PR interval measures atrial and atrioventricular depolarization and conduction, and abnormal PR interval is a risk factor for atrial fibrillation and heart block. We performed a genome-wide association study in over 92,000 individuals of European descent and identified 44 loci associated with PR interval (34 novel). Examination of the 44 loci revealed known and novel biological processes involved in cardiac atrial electrical activity, and genes in these loci were highly over-represented in several cardiac disease processes. Nearly half of the 61 independent index variants in the 44 loci were associated with atrial or blood transcript expression levels, or were in high linkage disequilibrium with one or more missense variants. Cardiac regulatory regions of the genome as measured by cardiac DNA hypersensitivity sites were enriched for variants associated with PR interval, compared to non-cardiac regulatory regions. Joint analyses combining PR interval with heart rate, QRS interval, and atrial fibrillation identified additional new pleiotropic loci. The majority of associations discovered in European-descent populations were also present in African-American populations. Meta-analysis examining over 105,000 individuals of African and European descent identified additional novel PR loci. These additional analyses identified another 13 novel loci. Together, these findings underscore the power of GWAS to extend knowledge of the molecular underpinnings of clinical processes.

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.

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.

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

This article is free to read on the publisher's website Variation in body iron is associated with or causes diseases, including anaemia and iron overload. Here, we analyse genetic association data on biochemical markers of iron status from 11 European-population studies, with replication in eight additional cohorts (total up to 48,972 subjects). We find 11 genome-wide-significant (P<5 × 10−8) loci, some including known iron-related genes (​HFE, ​SLC40A1, ​TF, ​TFR2, ​TFRC, ​TMPRSS6) and others novel (​ABO, ​ARNTL, ​FADS2, ​NAT2, ​TEX14). SNPs at ​ARNTL, ​TF, and ​TFR2 affect iron markers in ​HFE C282Y homozygotes at risk for hemochromatosis. There is substantial overlap between our iron loci and loci affecting erythrocyte and lipid phenotypes. These results will facilitate investigation of the roles of iron in disease.

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

Abstract 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 of 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 a 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. And 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.

<h3>Background</h3> Recently, a variant in TNFSF13B, encoding the cytokine and drug target B-cell activating factor (BAFF), has been associated with Systemic Lupus Erythematosus (SLE).¹ The aim of this study was to explore the BAFF-var effect on serologic and clinical features in a cohort of patients affected with SLE. <h3>Methods</h3> Overall, 190 Sardinian patients affected with SLE according to the modified 1997 ACR classification criteria and 256 Sardinian healthy controls were enrolled in this study and genotyped for the BAFF-var. In each patient demographic, serologic and clinical characteristics retrospectively collected at the time of SLE diagnosis and pre-therapy were recorded. Sera from 76 SLE patients, collected before starting therapy and stored at −80°, and 79 controls were used to measure soluble BAFF cytokine (ELISA). <h3>Results</h3> BAFF-var allelic frequency was higher in SLE patients (0.368) than in healthy controls (0.259) and associated with a higher risk of developing SLE (OR: 1.6; 95% CI: 1.2 to 2.2; p=0.0005). Serum BAFF concentration was significantly increased (p=1.61×10–9) in SLE cases (mean 1530 pg/ml; range 328–9327 pg/ml) versus healthy controls (mean 829 pg/ml; range 527–1410 pg/ml). Notably, when we stratified the data according to BAFF-var, the levels of serum BAFF increased in a genotype dependent way (p=0.001). No association with gender or age at SLE onset and BAFF-var was identified. Stratifying SLE manifestations according to ACR classification criteria, no significant correlation with any of the tested manifestations and the BAFF-var genotype was discovered. However, the quantitative levels of anti-dsDNA autoantibodies increased in a BAFF-var genotype dependent way (p=0.004), being higher in patients with BAFF-var homozygosis (88.5 UI/dl, IQR 4.1–491) than in those with wild-BAFF/BAFF-var heterozygosis (48.5 UI/dl, IQR 9.7–197) and wild BAFF homozygosis (29.0 UI/dl, IQR 3.5–116). <h3>Conclusion</h3> BAFF-var is associated with higher risk of SLE in general population and it is associated with increased serum BAFF and anti-dsDNA levels suggesting that it could also impact on SLE phenotype and outcomes. <h3>Reference</h3> . Steri M, et al. Overexpression of the cytokine BAFF and autoimmunity risk. N Engl J Med2017;376(17):1615–1626.

Abstract To investigate the genetic regulation of platelet (PLT) levels we carried out a whole-genome association analysis in 6,528 Sardinians from the general population of the Lanusei valley. We found 6 variants significantly influencing PLT levels, including a novel rare missense mutation (p.Pro27Ser) in the GP1BB protein that is associated with PLT reduction (P=1.17×10 −16 ). This mutation is rare in the SardiNIA population cohort (frequency of 0.45%), even rarer in the rest of the Sardinian island (frequency of 0.16%), and not reported elsewhere. Notably, GP1BB is involved in Bernard-Soulier syndrome (BSS), a rare autosomal recessive bleeding disorder caused by a defect in the platelet GPIb-IX-V protein complex. Consistently, the 57 identified individuals heterozygous for the p.P27S mutation showed mild thrombocytopenia, morphologically enlarged platelets (P=2.13×10 −10 ), and reduced expression of two GPIb-IX-V-complex components: GPIbα (−26.51%, P=3.66×10 −8 ) and GPIX (−24.69%, P=2.66×10 −6 ). Molecular modeling infers a corresponding reduction in the stability of GP1BB. These observations predict that in homozygosity as well as in individuals carrying specific compound heterozygous configurations, this variant likely causes BSS.

Additional file 1: Table S5 A. Association of MTHFR, XRCC1, OGG1 polymorphisms and ER, PgR, Her-2, Ki67 and Lymph Node status in BC patients. B. Association of MTHFR, XRCC1, OGG1 polymorphisms and Age at diagnosis, BMI, Menopause, BC family history status in BC patients.