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.

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

Phosphodiesterase 8B gene variants are associated with serum TSH levels and thyroid function.

Thyroid-stimulating hormone (TSH) controls thyroid growth and hormone secretion through binding to its G protein-coupled receptor (TSHR) and production of cyclic AMP (cAMP). Serum TSH is a sensitive indicator of thyroid function, and overt abnormalities in thyroid function lead to common endocrine disorders affecting approximately 10% of individuals over a life span. By genotyping 362,129 SNPs in 4,300 Sardinians, we identified a strong association (p = 1.3 x 10(-11)) between alleles of rs4704397 and circulating TSH levels; each additional copy of the minor A allele was associated with an increase of 0.13 muIU/ml in TSH. The single-nucleotide polymorphism (SNP) is located in intron 1 of PDE8B, encoding a high-affinity cAMP-specific phosphodiesterase. The association was replicated in 4,158 individuals, including additional Sardinians and two genetically distant cohorts from Tuscany and the Old Order Amish (overall p value = 1.9 x 10(-20)). In addition to association of TSH levels with SNPs in PDE8B, our genome scan provided evidence for association with PDE10A and several biologically interesting candidates in a focused analysis of 24 genes. In particular, we found evidence for association of TSH levels with SNPs in the THRB (rs1505287, p = 7.3 x 10(-5)), GNAQ (rs10512065, p = 2.0 x 10(-4)), TG (rs2252696, p = 2.2 x 10(-3)), POU1F1 (rs1976324, p = 3.9 x 10(-3)), PDE4D (rs27178, p = 8.3 x 10(-3)), and TSHR (rs4903957, p = 8.6 x 10(-3)) loci. Overall, the results suggest a primary effect of PDE8B variants on cAMP levels in the thyroid. This would affect production of T4 and T3 and feedback to alter TSH release by the pituitary. PDE8B may thus provide a candidate target for the treatment of thyroid dysfunction.

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.

The value of some Corsican sub-populations for genetic association studies.

BACKGROUND: Genetic isolates with a history of a small founder population, long-lasting isolation and population bottlenecks represent exceptional resources in the identification of disease genes. In these populations the disease allele reveals Linkage Disequilibrium (LD) with markers over significant genetic intervals, therefore facilitating disease locus identification. In a previous study we examined the LD extension on the Xq13 region in three Corsican sub-populations from the inner mountainous region of the island. On the basis of those previous results we have proposed a multistep procedure to carry out studies aimed at the identification of genes involved in complex diseases in Corsica. A prerequisite to carry out the proposed multi-step procedure was the presence of different degrees of LD on the island and a common genetic derivation of the different Corsican sub-populations. In order to evaluate the existence of these conditions in the present paper we extended the analysis to the Corsican coastal populations. METHODS: Samples were analyzed using seven dinucleotide microsatellite markers on chromosome Xq13-21: DXS983, DXS986, DXS8092, DXS8082, DXS1225, DXS8037 and DXS995 spanning approximately 4.0 cM (13.3 Mb). We have also investigated the distribution of the DXS1225-DXS8082 haplotype which has been recently proposed as a good marker of population genetic history due to its low recombination rate. RESULTS: the results obtained indicate a decrease of LD on the island from the central mountainous toward the coastal sub-populations. In addition the analysis of the DXS1225-DXS8082 haplotype revealed: 1) the presence of a particular haplotype with high frequency; 2) the derivation from a common genetic pool of the sub-populations examined in the present study. CONCLUSION: These results indicate the Corsican sub-populations useful for the fine mapping of genes contributing to complex diseases.

A mutation in the TMPRSS6 gene, encoding a transmembrane serine protease that suppresses hepcidin production, in familial iron deficiency anemia refractory to oral iron.

BACKGROUND: Hepcidin plays a key role in body iron metabolism by preventing the release of iron from macrophages and intestinal cells. Defective hepcidin synthesis causes iron loading, while overproduction results in defective reticuloendothelial iron release and iron absorption. DESIGN AND METHODS: We studied a Sardinian family in which microcytic anemia due to defective iron absorption and utilization is inherited as a recessive character. Five members showed iron deficiency anemia that was not responsive to oral iron and only partially responsive to parenteral iron administration. To investigate the involvement of known genes implicated in iron metabolism we carried out linkage analysis with microsatellite markers mapping close to these genes. Afterwards, a genome-wide search was performed. RESULTS: No linkage was found between the phenotype of the patients and several known human genes involved in iron metabolism (DMT1, TF, TFRC, ZIRTL, HAMP, HJV). Genome-wide scanning by microsatellites and single nucleotide polymorphisms showed a multipoint LOD score of 5.6 on chromosome 22q12.3-13.1, where the matriptase-2 (also known as transmembrane protease, serine 6 or TMPRSS6) gene is located. Its murine counterpart (Tmprss6) has recently been found to be an essential component of a pathway that detects iron deficiency and suppresses hepcidin production. Sequencing analysis of TMPRSS6 revealed a homozygous causal mutation, predicting a splicing error and a truncated TMPRSS6 protein in affected members. Homozygous subjects had inappropriately elevated levels of serum and urinary hepcidin. CONCLUSIONS: The findings of this study suggest that the observed TMPRSS6 mutation leads to overproduction of hepcidin and, in turn, to defective iron absorption and utilization. More generally, they confirm in humans the inhibitory effect of matriptase-2 on hepcidin synthesis already demonstrated in mice.

Genetic isolates in Corsica (France): linkage disequilibrium extension analysis on the Xq13 region.

Genetic isolates with a history of a small founder population, long-lasting isolation and population bottlenecks represent exceptional resources in the identification of genes involved in the pathogenesis of multifactorial diseases. In these populations, the disease allele reveals linkage disequilibrium (LD) with markers over significant genetic intervals, therefore facilitating disease locus identification. This study has been designed to examine the background LD extension in some subpopulations of Corsica. Our interest in the island of Corsica is due to its geographical and genetic proximity to the other Mediterranean island of Sardinia. Sardinian isolates in which the extension of the background LD is particularly high have been recently identified and are now the object of studies aimed at the mapping of genes involved in complex diseases. Recent evidence has highlighted that the genetic proximity between the populations of Corsica and Sardinia is particularly true for the internal conservative populations. Given these considerations, Sardinia and Corsica may represent a unique system to carry out parallel association studies whose results could be validated by comparison. In the present study, we have analyzed the LD extension on the Xq13 genomic region in three subpopulations of Corsica: Corte, Niolo and Bozio, all located in the mountainous north-center of the island. Our results show a strong degree of LD over long distance for the population of Bozio and to a less extent for the population of Niolo. Their LD extent is comparable to or higher than that reported for other isolates.

Serum IgE reactivity profiling in an asthma affected cohort.

BACKGROUND: Epidemiological evidence indicates that atopic asthma correlates with high serum IgE levels though the contribution of allergen specific IgE to the pathogenesis and the severity of the disease is still unclear. METHODS: We developed a microarray immunoassay containing 103 allergens to study the IgE reactivity profiles of 485 asthmatic and 342 non-asthmatic individuals belonging to families whose members have a documented history of asthma and atopy. We employed k-means clustering, to investigate whether a particular IgE reactivity profile correlated with asthma and other atopic conditions such as rhinitis, conjunctivitis and eczema. RESULTS: Both case-control and parent-to-siblings analyses demonstrated that while the presence of specific IgE against individual allergens correlated poorly with pathological conditions, particular reactivity profiles were significantly associated with asthma (p<10E-09). An artificial neural network (ANN)-based algorithm, calibrated with the profile reactivity data, correctly classified as asthmatic or non-asthmatic 78% of the individual examined. Multivariate statistical analysis demonstrated that the familiar relationships of the study population did not affect the observed correlations. CONCLUSIONS: These findings indicate that asthma is a higher-order phenomenon related to patterns of IgE reactivity rather than to single antibody reactions. This notion sheds new light on the pathogenesis of the disease and can be readily employed to distinguish asthmatic and non-asthmatic individuals on the basis of their serum reactivity profile.

[Heterozonal electrocardiographic changes in ST-elevation myocardial infarction]. / Coinvolgimento del miocardio remoto nel corso di infarto miocardico acuto con sopralivellamento del tratto ST: le alterazioni ischemiche eterozonali all'elettrocardiogramma.

BACKGROUND: It has been observed that in patients with ST-elevation myocardial infarction (STEMI) the presence of ST-segment depression (ST) in heterozonal electrocardiographic leads (remote STI) worsens the patient's prognosis. The aim of this study was to observe in an unselected population with a first STEMI the incidence of remote STI and the risk factors related to this condition. METHODS: We evaluated retrospectively 350 patients with a first STEMI; we excluded from our analysis 139 patients because no data about their coronary anatomy was available. ST-segment depression in the heterozonal myocardium was considered significant if > 0.1 mV at 60 ms from the J point, in at least two electrocardiographic leads. RESULTS: Patients were classified according to the presence (group 1, 117 patients) or absence (group II, 94 patients) of remote ST. The two groups did not differ for age, sex and coronary anatomy. In group I we found more heterozonal wall motion abnormalities than group II (32 vs. 18%, p = 0.018). In this group there was a higher incidence of smokers (56 vs. 33%, p = 0.025) and less patients were treated with statins when the STEMI occurred (6 vs. 14%, p = 0.047). Patients with remote ST had higher total cholesterol (214.6 +/- 48.9 vs. 192.3 +/- 29.8 mg/dl, p < 0.001) and low-density lipoprotein cholesterol (138.7 +/- 40.7 vs. 123.2 +/- 22.9 mg/dl, p < 0.0001) levels. Conclusions. In patients with STEMI the presence of remote ST is rather frequent, and seems to indicate a real heterozonal ischemia, independently of an epicardial coronary stenosis of the non-infarct-related artery. Remote ST. is associated with a higher incidence of risk factors related to microcirculatory dysfunction, such as cigarette smoking, a worse lipid profile and less protective factors, such as the use of statins prior to acute myocardial infarction.

IRAK-M is involved in the pathogenesis of early-onset persistent asthma.

Asthma is a multifactorial disease influenced by genetic and environmental factors. In the past decade, several loci and >100 genes have been found to be associated with the disease in at least one population. Among these loci, region 12q13-24 has been implicated in asthma etiology in multiple populations, suggesting that it harbors one or more asthma susceptibility genes. We performed linkage and association analyses by transmission/disequilibrium test and case-control analysis in the candidate region 12q13-24, using the Sardinian founder population, in which limited heterogeneity of pathogenetic alleles for monogenic and complex disorders as well as of environmental conditions should facilitate the study of multifactorial traits. We analyzed our cohort, using a cutoff age of 13 years at asthma onset, and detected significant linkage to a portion of 12q13-24. We identified IRAK-M as the gene contributing to the linkage and showed that it is associated with early-onset persistent asthma. We defined protective and predisposing SNP haplotypes and replicated associations in an outbred Italian population. Sequence analysis in patients found mutations, including inactivating lesions, in the IRAK-M coding region. Immunohistochemistry of lung biopsies showed that IRAK-M is highly expressed in epithelial cells. We report that IRAK-M is involved in the pathogenesis of early-onset persistent asthma. IRAK-M, a negative regulator of the Toll-like receptor/IL-1R pathways, is a master regulator of NF- kappa B and inflammation. Our data suggest a mechanistic link between hyperactivation of the innate immune system and chronic airway inflammation and indicate IRAK-M as a potential target for therapeutic intervention against asthma.

A locus for familial skewed X chromosome inactivation maps to chromosome Xq25 in a family with a female manifesting Lowe syndrome.

In mammals, X-linked gene products can be dosage compensated between males and females by inactivation of one of the two X chromosomes in the developing female embryos. X inactivation choice is usually random in embryo mammals, but several mechanisms can influence the choice determining skewed X inactivation. As a consequence, females heterozygous for X-linked recessive disease can manifest the full phenotype. Herein, we report a family with extremely skewed X inactivation that produced the full phenotype of Lowe syndrome, a recessive X-linked disease, in a female. The X chromosome inactivation studies detected an extremely skewed inactivation pattern with a ratio of 100:0 in the propositus as well as in five out of seven unaffected female relatives in four generations. The OCRL1 "de novo" mutation resides in the active paternally inherited X chromosome. X chromosome haplotype analysis suggests the presence of a locus for the familial skewed X inactivation in chromosome Xq25 most likely controlling X chromosome choice in X inactivation or cell proliferation. The description of this case adds Lowe syndrome to the list of X-linked disorders which may manifest the full phenotype in females because of the skewed X inactivation.