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.

Gorlin syndrome (GS) is inherited in an autosomal dominant pattern with high-penetrance and is characterized by a range of developmental anomalies and increased risk of developing basal cell carcinoma and medulloblastoma. Between 50% and 85% of patients with GS harbor germ line mutations in the only susceptibility gene identified to date, PTCH1, a key component in the Sonic Hedgehog signaling pathway. Another component in this pathway, SUFU, is known to be involved in susceptibility to medulloblastoma but has never been reported in GS patients to date. We have identified the known c.1022 + 1G>A SUFU germ line splicing mutation in a family that was PTCH1-negative and who had signs and symptoms of GS, including medulloblastoma. This is the first report of a germ line SUFU mutation associated with GS.

Rett syndrome is a neurological disorder affecting predominantly females with regression loss of speech and purposeful hand use, after a few months of almost normal development. Postnatal microcephaly, hand dispraxia, stereotypic 'hand-washing' activities, ataxia, and abnormal breathing are among its most characteristic features. Another aspect of this disorder is growth failure. The preserved speech variant (PSV) shares with Rett syndrome the same course and the stereotypic hand-washing activities but it differs in that patients typically recover some degree of speech and hand use and usually do not show growth failure. Progressive scoliosis, epilepsy and other minor handicaps, usually present in Rett syndrome, are rare in the preserved speech variant. Here we explore the spectrum of mutations affecting the MECP2 gene in a group of 25 classic Rett syndrome girls and in three patients with the preserved speech variant. Among the Rett syndrome group, two novel mutational hot spots (R270X and R294X), four novel mutations, two novel small deletions, as well as the previously reported 806delG, R168X and R255X mutations, were identified in 20/25 patients. Of note, among the preserved speech variants, two patients carry deletions of 41 bp and 44 bp each, which are strikingly similar to those observed in classic Rett syndrome. Our results confirm the presence of mutational hot spots in MECP2, broaden the spectrum of mutations, pinpoint additional mutational hot spots and establish that the preserved speech variant is indeed allelic of the classic form. Phenotype variability is only partially dependent on the kind of MECP2 mutation and other mechanisms such as skewed X-inactivation, and/or modifier gene effects should be investigated to explain the variable recovery in speech and hand use.

Gnathodiaphyseal dysplasia (GDD) is an autosomal dominant syndrome characterized by frequent bone fractures at a young age, bowing of tubular bones and cemento-osseus lesions of the jawbones. Anoctamin 5 (ANO5) belongs to the anoctamin protein family that includes calcium-activated chloride channels. However, recent data together with our own experiments reported here add weight to the hypothesis that ANO5 may not function as calcium-activated chloride channel. By sequencing the entire ANO5 gene coding region and untranslated regions in a large Italian GDD family, we found a novel missense mutation causing the p.Thr513Ile substitution. The mutation segregates with the disease in the family and has never been described in any database as a polymorphism. To date, only two mutations on the same cysteine residue at position 356 of ANO5 amino-acid sequence have been described in GDD families. As ANO5 has also been found to be mutated in two different forms of muscular dystrophy, the finding of this third mutation in GDD adds clues to the role of ANO5 in these disorders.