Intestinal Blastocystis is linked to healthier diets and more favorable cardiometabolic outcomes in 56,989 individuals from 32 countries.

Diet impacts human health, influencing body adiposity and the risk of developing cardiometabolic diseases. The gut microbiome is a key player in the diet-health axis, but while its bacterial fraction is widely studied, the role of micro-eukaryotes, including Blastocystis, is underexplored. We performed a global-scale analysis on 56,989 metagenomes and showed that human Blastocystis exhibits distinct prevalence patterns linked to geography, lifestyle, and dietary habits. Blastocystis presence defined a specific bacterial signature and was positively associated with more favorable cardiometabolic profiles and negatively with obesity (p < 1e-16) and disorders linked to altered gut ecology (p < 1e-8). In a diet intervention study involving 1,124 individuals, improvements in dietary quality were linked to weight loss and increases in Blastocystis prevalence (p = 0.003) and abundance (p < 1e-7). Our findings suggest a potentially beneficial role for Blastocystis, which may help explain personalized host responses to diet and downstream disease etiopathogenesis.

Genetic variants regulating immune cell levels in health and disease.

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

Multiple sclerosis disease-modifying therapies and COVID-19 vaccines: a practical review and meta-analysis.

Studies among people with multiple sclerosis (pwMS) receiving disease-modifying therapies (DMTs) have provided adequate evidence for an appraisal of COVID-19 vaccination policies among them. To synthesise the available evidence addressing the effect of MS DMTs on COVID-19 vaccines' immunogenicity and effectiveness, following the Cochrane guidelines, we systematically reviewed all observational studies available in MEDLINE, Scopus, Web of Science, MedRxiv and Google Scholar from January 2021 to January 2022 and extracted their relevant data. Immunogenicity data were then synthesised in a quantitative, and other data in a qualitative manner. Evidence from 28 studies suggests extensively lower B-cell responses in sphingosine-1-phosphate receptor modulator (S1PRM) treated and anti-CD20 (aCD20) treated, and lower T-cell responses in interferon-treated, S1PRM-treated and cladribine-treated pwMS-although most T cell evidence currently comprises of low or very low certainty. With every 10-week increase in aCD20-to-vaccine period, a 1.94-fold (95% CI 1.57 to 2.41, p<0.00001) increase in the odds of seroconversion was observed. Furthermore, the evidence points out that B-cell-depleting therapies may accelerate postvaccination humoral waning, and boosters' immunogenicity is predictable with the same factors affecting the initial vaccination cycle. Four real-world studies further indicate that the comparative incidence/severity of breakthrough COVID-19 has been higher among the pwMS treated with S1PRM and aCD20-unlike the ones treated with other DMTs. S1PRM and aCD20 therapies were the only DMTs reducing the real-world effectiveness of COVID-19 vaccination among pwMS. Hence, it could be concluded that optimisation of humoral immunogenicity and ensuring its durability are the necessities of an effective COVID-19 vaccination policy among pwMS who receive DMTs.

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

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

Overexpression of the Cytokine BAFF and Autoimmunity Risk.

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

Mitogenome Diversity in Sardinians: A Genetic Window onto an Island's Past.

Sardinians are "outliers" in the European genetic landscape and, according to paleogenomic nuclear data, the closest to early European Neolithic farmers. To learn more about their genetic ancestry, we analyzed 3,491 modern and 21 ancient mitogenomes from Sardinia. We observed that 78.4% of modern mitogenomes cluster into 89 haplogroups that most likely arose in situ. For each Sardinian-specific haplogroup (SSH), we also identified the upstream node in the phylogeny, from which non-Sardinian mitogenomes radiate. This provided minimum and maximum time estimates for the presence of each SSH on the island. In agreement with demographic evidence, almost all SSHs coalesce in the post-Nuragic, Nuragic and Neolithic-Copper Age periods. For some rare SSHs, however, we could not dismiss the possibility that they might have been on the island prior to the Neolithic, a scenario that would be in agreement with archeological evidence of a Mesolithic occupation of Sardinia.

Population genomic analysis of ancient and modern genomes yields new insights into the genetic ancestry of the Tyrolean Iceman and the genetic structure of Europe.

Genome sequencing of the 5,300-year-old mummy of the Tyrolean Iceman, found in 1991 on a glacier near the border of Italy and Austria, has yielded new insights into his origin and relationship to modern European populations. A key finding of that study was an apparent recent common ancestry with individuals from Sardinia, based largely on the Y chromosome haplogroup and common autosomal SNP variation. Here, we compiled and analyzed genomic datasets from both modern and ancient Europeans, including genome sequence data from over 400 Sardinians and two ancient Thracians from Bulgaria, to investigate this result in greater detail and determine its implications for the genetic structure of Neolithic Europe. Using whole-genome sequencing data, we confirm that the Iceman is, indeed, most closely related to Sardinians. Furthermore, we show that this relationship extends to other individuals from cultural contexts associated with the spread of agriculture during the Neolithic transition, in contrast to individuals from a hunter-gatherer context. We hypothesize that this genetic affinity of ancient samples from different parts of Europe with Sardinians represents a common genetic component that was geographically widespread across Europe during the Neolithic, likely related to migrations and population expansions associated with the spread of agriculture.

Genome-wide association study of susceptibility loci for breast cancer in Sardinian population.

BACKGROUND: Despite progress in identifying genes associated with breast cancer, many more risk loci exist. Genome-wide association analyses in genetically-homogeneous populations, such as that of Sardinia (Italy), could represent an additional approach to detect low penetrance alleles. METHODS: We performed a genome-wide association study comparing 1431 Sardinian patients with non-familial, BRCA1/2-mutation-negative breast cancer to 2171 healthy Sardinian blood donors. DNA was genotyped using GeneChip Human Mapping 500 K Arrays or Genome-Wide Human SNP Arrays 6.0. To increase genomic coverage, genotypes of additional SNPs were imputed using data from HapMap Phase II. After quality control filtering of genotype data, 1367 cases (9 men) and 1658 controls (1156 men) were analyzed on a total of 2,067,645 SNPs. RESULTS: Overall, 33 genomic regions (67 candidate SNPs) were associated with breast cancer risk at the p < 0(-6) level. Twenty of these regions contained defined genes, including one already associated with breast cancer risk: TOX3. With a lower threshold for preliminary significance to p < 10(-5), we identified 11 additional SNPs in FGFR2, a well-established breast cancer-associated gene. Ten candidate SNPs were selected, excluding those already associated with breast cancer, for technical validation as well as replication in 1668 samples from the same population. Only SNP rs345299, located in intron 1 of VAV3, remained suggestively associated (p-value, 1.16 x 10(-5)), but it did not associate with breast cancer risk in pooled data from two large, mixed-population cohorts. CONCLUSIONS: This study indicated the role of TOX3 and FGFR2 as breast cancer susceptibility genes in BRCA1/2-wild-type breast cancer patients from Sardinian population.

The burden of multiple sclerosis variants in continental Italians and Sardinians.

BACKGROUND: Recent studies identified > 100 non-HLA (human leukocyte antigen) multiple sclerosis (MS) susceptibility variants in Northern European populations, but their role in Southern Europeans is largely unexplored. OBJECTIVE: We aimed to investigate the cumulative impact of those variants in two Mediterranean populations: Continental Italians and Sardinians. METHODS: We calculated four weighted Genetic Risk Scores (wGRS), using up to 102 non-HLA MS risk variants and 5 HLA MS susceptibility markers in 1691 patients and 2194 controls from continental Italy; and 2861 patients and 3034 controls from Sardinia. We then assessed the differences between populations using Nagelkerke's R(2) and the area under the Receiver Operating Characteristic (ROC) curves. RESULTS: As expected, the genetic burden (mean wGRS value) was significantly higher in MS patients than in controls, in both populations. Of note, the burden was significantly higher in Sardinians. Conversely, the proportion of variability explained and the predictive power were significantly higher in continental Italians. Notably, within the Sardinian patients, we also observed a significantly higher burden of non-HLA variants in individuals who do not carry HLA risk alleles. CONCLUSIONS: The observed differences in MS genetic burden between the two Mediterranean populations highlight the need for more genetic studies in South Europeans, to further expand the knowledge of MS genetics.

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

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

Intermittent and periodic fasting - Evidence and perspectives in multiple sclerosis.

Multiple sclerosis (MS) is a complex neurological disease characterized by great variability in clinical presentation, including the radiological features, and degree of disability. Both genetics and environment contribute to disease etiopathogenesis. Because MS is more common in Western countries, and diet has been proposed among the etiologic factors. However, based on the several studies published thus far, principally involving small cohorts, there is no described diet-protocol to be applied in clinical practice as a supplement to the standard immunomodulatory treatment of MS. Diet is an easily changeable factor thus the research on the diet importance in MS has been exploded in last years. Starting from the notions that diet can change lifespan and quality of life in general, and its improvement could be one of many contributing factors with effects on disease evolution, this review examines the evidence of the effects of intermittent fasting in a mouse model of MS; the evidence derived from clinical trials; and future perspectives.

Omics approaches to understanding the efficacy and safety of disease-modifying treatments in multiple sclerosis.

From the perspective of precision medicine, the challenge for the future is to improve the accuracy of diagnosis, prognosis, and prediction of therapeutic responses through the identification of biomarkers. In this framework, the omics sciences (genomics, transcriptomics, proteomics, and metabolomics) and their combined use represent innovative approaches for the exploration of the complexity and heterogeneity of multiple sclerosis (MS). This review examines the evidence currently available on the application of omics sciences to MS, analyses the methods, their limitations, the samples used, and their characteristics, with a particular focus on biomarkers associated with the disease state, exposure to disease-modifying treatments (DMTs), and drug efficacies and safety profiles.

TCRBV20S1 polymorphism does not influence the susceptibility to type 1 diabetes and multiple sclerosis in the Sardinian population.

Among the different T-cell receptor (TCR) BV20S1 polymorphisms, nucleotide substitution at position 524 results in the introduction of a stop codon, whose potential functional relevance is still unknown. We have recently showed in Sardinian subjects the most elevated allele frequency ever reported worldwide for this "null allele" (0.44). As this variant generates a gap in the TCR repertoire, this preliminary finding prompted us to further analyze the role of this polymorphism in the susceptibility to type 1 diabetes (T1D) and multiple sclerosis (MS), which are extremely common in this population. With this aim, we evaluated the influence of the TCRBV20S1 polymorphism by assessing it with the transmission disequilibirum test (TDT) in 652 T1D and 616 MS families, without detecting any significant difference. We conclude that the high frequency of this null allele in Sardinia is not directly related to the high incidence of these autoimmune diseases observed in this founder population.

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

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

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

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

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

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

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

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

Genetic loci linked to type 1 diabetes and multiple sclerosis families in Sardinia.

BACKGROUND: The Mediterranean island of Sardinia has a strikingly high incidence of the autoimmune disorders Type 1 Diabetes (T1D) and Multiple Sclerosis (MS). Furthermore, the two diseases tend to be co-inherited in the same individuals and in the same families. These observations suggest that some unknown autoimmunity variant with relevant effect size could be fairly common in this founder population and could be detected using linkage analysis. METHODS: To search for T1D and MS loci as well as any that predispose to both diseases, we performed a whole genome linkage scan, sequentially genotyping 593 microsatellite marker loci in 954 individuals distributed in 175 Sardinian families. In total, 413 patients were studied; 285 with T1D, 116 with MS and 12 with both disorders. Model-free linkage analysis was performed on the genotyped samples using the Kong and Cox logarithm of odds (LOD) score statistic. RESULTS: In T1D, aside from the HLA locus, we found four regions showing a lod-score > or =1; 1p31.1, 6q26, 10q21.2 and 22q11.22. In MS we found three regions showing a lod-score > or =1; 1q42.2, 18p11.21 and 20p12.3. In the combined T1D-MS scan for shared autoimmunity loci, four regions showed a LOD >1, including 6q26, 10q21.2, 20p12.3 and 22q11.22. When we typed more markers in these intervals we obtained suggestive evidence of linkage in the T1D scan at 10q21.2 (LOD = 2.1), in the MS scan at 1q42.2 (LOD = 2.5) and at 18p11.22 (LOD = 2.6). When all T1D and MS families were analysed jointly we obtained suggestive evidence in two regions: at 10q21.1 (LOD score = 2.3) and at 20p12.3 (LOD score = 2.5). CONCLUSION: This suggestive evidence of linkage with T1D, MS and both diseases indicates critical chromosome intervals to be followed up in downstream association studies.

Genomic history of the Sardinian population.

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.