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

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

In the present study, we characterized the distinctive signatures of the gut microbiota (GM) from overweight/obese patients (OB), and normal-weight controls (NW), both of Sardinian origin. Fecal bacterial composition of 46 OB patients (BMI = 36.6 ± 6.0; F/M = 40/6) was analyzed and compared to that of 46 NW subjects (BMI = 21.6 ± 2.1; F/M = 41/5), matched for sex, age and smoking status, by using 16S rRNA gene sequencing on MiSeq Illumina platform. The gut microbial community of OB patients exhibited a significant decrease in the relative abundance of several Bacteroidetes taxa (i.e. Flavobacteriaceae, Porphyromonadaceae, Sphingobacteriaceae, Flavobacterium, Rikenella spp., Pedobacter spp., Parabacteroides spp., Bacteroides spp.) when compared to NW; instead, several Firmicutes taxa were significantly increased in the same subjects (Lachnospiraceae, Gemellaceae, Paenibacillaceae, Streptococcaceae, Thermicanaceae, Gemella, Mitsuokella, Streptococcus, Acidaminococcus spp., Eubacterium spp., Ruminococcus spp., Megamonas spp., Streptococcus, Thermicanus, Megasphaera spp. and Veillonella spp.). Correlation analysis indicated that body fatness and waist circumference negatively correlated with Bacteroidetes taxa, while Firmicutes taxa positively correlated with body fat and negatively with muscle mass and/or physical activity level. Furthermore, the relative abundance of several bacterial taxa belonging to Enterobacteriaceae family, known to exhibit endotoxic activity, was increased in the OB group compared to NW. The results extend our knowledge on the GM profiles in Italian OB, identifying novel taxa linking obesity and intestine.

Genetic variation within the male-specific portion of the Y chromosome (MSY) can clarify the origins of contemporary populations, but previous studies were hampered by partial genetic information. Population sequencing of 1204 Sardinian males identified 11,763 MSY single-nucleotide polymorphisms, 6751 of which have not previously been observed. We constructed a MSY phylogenetic tree containing all main haplogroups found in Europe, along with many Sardinian-specific lineage clusters within each haplogroup. The tree was calibrated with archaeological data from the initial expansion of the Sardinian population ~7700 years ago. The ages of nodes highlight different genetic strata in Sardinia and reveal the presumptive timing of coalescence with other human populations. We calculate a putative age for coalescence of ~180,000 to 200,000 years ago, which is consistent with previous mitochondrial DNA–based estimates.

Parkinson’s disease is a neurodegenerative disorder characterized by the accumulation of intracellular aggregates of misfolded alpha-synuclein along the cerebral axis. Several studies report the association between intestinal dysbiosis and Parkinson’s disease, although a cause-effect relationship remains to be established. Herein, the gut microbiota composition of 64 Italian patients with Parkinson’s disease and 51 controls was determined using a next-generation sequencing approach. A real metagenomics shape based on gas chromatography-mass spectrometry was also investigated. The most significant changes within the Parkinson’s disease group highlighted a reduction in bacterial taxa, which are linked to anti-inflammatory/neuroprotective effects, particularly in the Lachnospiraceae family and key members, such as Butyrivibrio, Pseudobutyrivibrio, Coprococcus, and Blautia. The direct evaluation of fecal metabolites revealed changes in several classes of metabolites. Changes were seen in lipids (linoleic acid, oleic acid, succinic acid, and sebacic acid), vitamins (pantothenic acid and nicotinic acid), amino acids (isoleucine, leucine, phenylalanine, glutamic acid, and pyroglutamic acid) and other organic compounds (cadaverine, ethanolamine, and hydroxy propionic acid). Most modified metabolites strongly correlated with the abundance of members belonging to the Lachnospiraceae family, suggesting that these gut bacteria correlate with altered metabolism rates in Parkinson’s disease. IMPORTANCE To our knowledge, this is one of the few studies thus far that correlates the composition of the gut microbiota with the direct analysis of fecal metabolites in patients with Parkinson’s disease. Overall, our data highlight microbiota modifications correlated with numerous fecal metabolites. This suggests that Parkinson’s disease is associated with gut dysregulation that involves a synergistic relationship between gut microbes and several bacterial metabolites favoring altered homeostasis. Interestingly, a reduction of short-chain fatty acid (SCFA)-producing bacteria influenced the shape of the metabolomics profile, affecting several metabolites with potential protective effects in the Parkinson group. On the other hand, the extensive impact that intestinal dysbiosis has at the level of numerous metabolic pathways could encourage the identification of specific biomarkers for the diagnosis and treatment of Parkinson’s disease, also in light of the effect that specific drugs have on the composition of the intestinal microbiota.