A Whole-Genome Sequencing Study Implicates GRAMD1B in Multiple Sclerosis Susceptibility.

While the role of common genetic variants in multiple sclerosis (MS) has been elucidated in large genome-wide association studies, the contribution of rare variants to the disease remains unclear. Herein, a whole-genome sequencing study in four affected and four healthy relatives of a consanguineous Italian family identified a novel missense c.1801T > C (p.S601P) variant in the GRAMD1B gene that is shared within MS cases and resides under a linkage peak (LOD: 2.194). Sequencing GRAMD1B in 91 familial MS cases revealed two additional rare missense and two splice-site variants, two of which (rs755488531 and rs769527838) were not found in 1000 Italian healthy controls. Functional studies demonstrated that GRAMD1B, a gene with unknown function in the central nervous system (CNS), is expressed by several cell types, including astrocytes, microglia and neurons as well as by peripheral monocytes and macrophages. Notably, GRAMD1B was downregulated in vessel-associated astrocytes of active MS lesions in autopsied brains and by inflammatory stimuli in peripheral monocytes, suggesting a possible role in the modulation of inflammatory response and disease pathophysiology.

Relationship between clone metrics and clinical outcome in clonal cytopenia.

Clonal cytopenia of undetermined significance (CCUS) is associated with an increased risk of developing a myeloid neoplasm with myelodysplasia (MN). To identify the features of the mutant clone(s) that is associated with clinical phenotype and progression, we studied the following cohorts of individuals: 311 patients with idiopathic cytopenia of undetermined significance (ICUS), 532 community-dwelling individuals without hematologic phenotype (n = 355) or with unexplained anemia (n = 177), and 592 patients with overt MN. Ninety-two of 311 (30%) patients with ICUS carried a somatic genetic lesion that signaled CCUS. Clonal hematopoiesis (CH) was detected in 19.7% and 27.7% of nonanemic and anemic community-dwelling individuals, respectively. Different mutation patterns and variant allele frequencies (VAFs) (clone metrics parameters) were observed in the conditions studied. Recurrent mutation patterns exhibited different VAFs associated with marrow dysplasia (0.17-0.48), indicating variable clinical expressivity of mutant clones. Unsupervised clustering analysis based on mutation profiles identified 2 major clusters, characterized by isolated DNMT3A mutations (CH-like cluster) or combinatorial mutation patterns (MN-like cluster), and showing different overall survival (HR, 1.8). In patients with CCUS, the 2 clusters had different risk of progression to MN (HR, 2.7). Within the MN-like cluster, distinct subsets with different risk of progression to MN were identified based on clone metrics. These findings unveil marked variability in the clinical expressivity of myeloid driver genes and underline the limitations of morphologic dysplasia for clinical staging of mutant hematopoietic clones. Clone metrics appears to be critical for informing clinical decision-making in patients with clonal cytopenia.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color.

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

A bird's-eye view of Italian genomic variation through whole-genome sequencing.

The genomic variation of the Italian peninsula populations is currently under characterised: the only Italian whole-genome reference is represented by the Tuscans from the 1000 Genome Project. To address this issue, we sequenced a total of 947 Italian samples from three different geographical areas. First, we defined a new Italian Genome Reference Panel (IGRP1.0) for imputation, which improved imputation accuracy, especially for rare variants, and we tested it by GWAS analysis on red blood traits. Furthermore, we extended the catalogue of genetic variation investigating the level of population structure, the pattern of natural selection, the distribution of deleterious variants and occurrence of human knockouts (HKOs). Overall the results demonstrate a high level of genomic differentiation between cohorts, different signatures of natural selection and a distinctive distribution of deleterious variants and HKOs, confirming the necessity of distinct genome references for the Italian population.

Big Data in Medicine, the Present and Hopefully the Future.

The emergence of data coming from different venues, as several "omic" approaches, is providing already compelling evidence that the smart use of this information could provide invaluable information to prevent, diagnose and treat human diseases. However, the most daunting challenges remain ahead, as the explosive accumulation of data from additional perspectives, including social graphs, biosensors, and imaging, promise to deliver crucial information that could be exploited for the improvement of the entire human race, both in developed, and developing countries, optimizing health expenses and reaching also the less fortunate sections of the societies. And yet, formidable challenges remain, that pertain for the most part to the collection of the data, their organization, and most relevantly their integration. Here we provide few, pointed examples to the present relevance of these big data approaches in human health as well potential road maps toward the implementation of broader data collections and analyses.

Publisher Correction: Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability.

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

Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability.

Hair color is one of the most recognizable visual traits in European populations and is under strong genetic control. Here we report the results of a genome-wide association study meta-analysis of almost 300,000 participants of European descent. We identified 123 autosomal and one X-chromosome loci significantly associated with hair color; all but 13 are novel. Collectively, single-nucleotide polymorphisms associated with hair color within these loci explain 34.6% of red hair, 24.8% of blond hair, and 26.1% of black hair heritability in the study populations. These results confirm the polygenic nature of complex phenotypes and improve our understanding of melanin pigment metabolism in humans.

Reconstructing the genetic history of Italians: new insights from a male (Y-chromosome) perspective.

BACKGROUND: Due to its central and strategic position in Europe and in the Mediterranean Basin, the Italian Peninsula played a pivotal role in the first peopling of the European continent and has been a crossroad of peoples and cultures since then. AIM: This study aims to gain more information on the genetic structure of modern Italian populations and to shed light on the migration/expansion events that led to their formation. SUBJECTS AND METHODS: High resolution Y-chromosome variation analysis in 817 unrelated males from 10 informative areas of Italy was performed. Haplogroup frequencies and microsatellite haplotypes were used, together with available data from the literature, to evaluate Mediterranean and European inputs and date their arrivals. RESULTS: Fifty-three distinct Y-chromosome lineages were identified. Their distribution is in general agreement with geography, southern populations being more differentiated than northern ones. CONCLUSIONS: A complex genetic structure reflecting the multifaceted peopling pattern of the Peninsula emerged: southern populations show high similarity with those from the Middle East and Southern Balkans, while those from Northern Italy are close to populations of North-Western Europe and the Northern Balkans. Interestingly, the population of Volterra, an ancient town of Etruscan origin in Tuscany, displays a unique Y-chromosomal genetic structure.

Genome-Wide Meta-Analysis Unravels Interactions between Magnesium Homeostasis and Metabolic Phenotypes.

Magnesium (Mg2+) homeostasis is critical for metabolism. However, the genetic determinants of the renal handling of Mg2+, which is crucial for Mg2+ homeostasis, and the potential influence on metabolic traits in the general population are unknown. We obtained plasma and urine parameters from 9099 individuals from seven cohorts, and conducted a genome-wide meta-analysis of Mg2+ homeostasis. We identified two loci associated with urinary magnesium (uMg), rs3824347 (P=4.4×10-13) near TRPM6, which encodes an epithelial Mg2+ channel, and rs35929 (P=2.1×10-11), a variant of ARL15, which encodes a GTP-binding protein. Together, these loci account for 2.3% of the variation in 24-hour uMg excretion. In human kidney cells, ARL15 regulated TRPM6-mediated currents. In zebrafish, dietary Mg2+ regulated the expression of the highly conserved ARL15 ortholog arl15b, and arl15b knockdown resulted in renal Mg2+ wasting and metabolic disturbances. Finally, ARL15 rs35929 modified the association of uMg with fasting insulin and fat mass in a general population. In conclusion, this combined observational and experimental approach uncovered a gene-environment interaction linking Mg2+ deficiency to insulin resistance and obesity.

Fertility Preservation in Endometriosis Patients: Anti-Müllerian Hormone Is a Reliable Marker of the Ovarian Follicle Density.

OBJECTIVE: To analyze the ovarian reserve via measurement of follicular density and anti-Müllerian hormone (AMH) in endometriosis patients participating to a clinical program of cortical ovarian cryopreservation. DESIGN: Retrospective analysis of serum AMH levels and prospective investigation of ovarian follicle number. SETTING: University Hospital. PATIENTS: Two hundred and two women with endometriosis and 400 controls. INTERVENTIONS: Blood samples and ovarian biopsies. MAIN OUTCOME MEASURES: Correlation of serum AMH levels and the number of non-growing follicles in the biopsied cortical tissues in endometriosis and control subjects, including age, type of AMH kit, and the laboratory performing the analysis as covariates. RESULTS: AMH levels were shown to decrease with age in untreated endometriosis patients (P < 1.0 × 10-5) but they were significantly lower in endometriosis compared to controls only in patients over 36 years old (P = 2.7 × 10-4). The AMH decrease was faster in endometriosis compared to controls (beta = 0.27, P = 4.0 × 10-4). Primordial follicle number decreased with the reduction of AMH levels in both cases and controls (beta = 0.3; P = 0.04). CONCLUSION: AMH is a reliable marker of ovarian reserve in endometriosis patients, and it can predict follicular density in women undergoing ovarian tissue cryopreservation.

Enrichment of low-frequency functional variants revealed by whole-genome sequencing of multiple isolated European populations.

The genetic features of isolated populations can boost power in complex-trait association studies, and an in-depth understanding of how their genetic variation has been shaped by their demographic history can help leverage these advantageous characteristics. Here, we perform a comprehensive investigation using 3,059 newly generated low-depth whole-genome sequences from eight European isolates and two matched general populations, together with published data from the 1000 Genomes Project and UK10K. Sequencing data give deeper and richer insights into population demography and genetic characteristics than genotype-chip data, distinguishing related populations more effectively and allowing their functional variants to be studied more fully. We demonstrate relaxation of purifying selection in the isolates, leading to enrichment of rare and low-frequency functional variants, using novel statistics, DVxy and SVxy. We also develop an isolation-index (Isx) that predicts the overall level of such key genetic characteristics and can thus help guide population choice in future complex-trait association studies.

Clinical significance of somatic mutation in unexplained blood cytopenia.

Unexplained blood cytopenias, in particular anemia, are often found in older persons. The relationship between these cytopenias and myeloid neoplasms like myelodysplastic syndromes is currently poorly defined. We studied a prospective cohort of patients with unexplained cytopenia with the aim to estimate the predictive value of somatic mutations for identifying subjects with, or at risk of, developing a myeloid neoplasm. The study included a learning cohort of 683 consecutive patients investigated for unexplained cytopenia, and a validation cohort of 190 patients referred for suspected myeloid neoplasm. Using granulocyte DNA, we looked for somatic mutations in 40 genes that are recurrently mutated in myeloid malignancies. Overall, 435/683 patients carried a somatic mutation in at least 1 of these genes. Carrying a somatic mutation with a variant allele frequency ≥0.10, or carrying 2 or more mutations, had a positive predictive value for diagnosis of myeloid neoplasm equal to 0.86 and 0.88, respectively. Spliceosome gene mutations and comutation patterns involving TET2, DNMT3A, or ASXL1 had positive predictive values for myeloid neoplasm ranging from 0.86 to 1.0. Within subjects with inconclusive diagnostic findings, carrying 1 or more somatic mutations was associated with a high probability of developing a myeloid neoplasm during follow-up (hazard ratio = 13.9, P < .001). The predictive values of mutation analysis were confirmed in the independent validation cohort. The findings of this study indicate that mutation analysis on peripheral blood granulocytes may significantly improve the current diagnostic approach to unexplained cytopenia and more generally the diagnostic accuracy of myeloid neoplasms.

COL6A5 variants in familial neuropathic chronic itch.

Itch is thought to represent the peculiar response to stimuli conveyed by somatosensory pathways shared with pain through the activation of specific neurons and receptors. It can occur in association with dermatological, systemic and neurological diseases, or be the side effect of certain drugs. However, some patients suffer from chronic idiopathic itch that is frequently ascribed to psychological distress and for which no biomarker is available to date. We investigated three multigenerational families, one of which diagnosed with joint hypermobility syndrome/Ehlers-Danlos syndrome hypermobility type (JHS/EDS-HT), characterized by idiopathic chronic itch with predominantly proximal distribution. Skin biopsy was performed in all eight affected members and revealed in six of them reduced intraepidermal nerve fibre density consistent with small fibre neuropathy. Whole exome sequencing identified two COL6A5 rare variants co-segregating with chronic itch in eight affected members and absent in non-affected members, and in one unrelated sporadic patient with type 1 painless diabetic neuropathy and chronic itch. Two families and the diabetic patient carried the nonsense c.6814G>T (p.Glu2272*) variant and another family carried the missense c.6486G>C (p.Arg2162Ser) variant. Both variants were predicted as likely pathogenic by in silico analyses. The two variants were rare (minor allele frequency < 0.1%) in 6271 healthy controls and absent in 77 small fibre neuropathy and 167 JHS/EDS-HT patients without itch. Null-allele test on cDNA from patients' fibroblasts of both families carrying the nonsense variant demonstrated functional haploinsufficiency due to activation of nonsense mediated RNA decay. Immunofluorescence microscopy and western blotting revealed marked disorganization and reduced COL6A5 synthesis, respectively. Indirect immunofluorescence showed reduced COL6A5 expression in the skin of patients carrying the nonsense variant. Treatment with gabapentinoids provided satisfactory itch relief in the patients carrying the mutations. Our findings first revealed an association between COL6A5 gene and familiar chronic itch, suggesting a new contributor to the pathogenesis of neuropathic itch and identifying a new candidate therapeutic target.

MCM8 and MCM9 Nucleotide Variants in Women With Primary Ovarian Insufficiency.

Objective: To assess the frequency of variants, including biallelic pathogenic variants, in minichromosome maintenance 8 (MCM8) and minichromosome maintenance 9 (MCM9), other genes related to MCM8-MCM9, and DNA damage repair (DDR) pathway in participants with primary ovarian insufficiency (POI). Design: MCM8, MCM9, and genes encoding DDR proteins that have been implicated in reproductive aging were sequenced among POI participants. Setting: Academic research institution. Participants: All were diagnosed with POI prior to age 40 years and presented with elevated follicle-stimulating hormone levels. Interventions: None. Main Outcome Measures: We identified nucleotide variants in MCM8, MCM9, and genes thought to be involved in the DNA damage response pathway and/or implicated in reproductive aging. Results: MCM8 was sequenced in 155 POI participants, whereas MCM9 was sequenced in 151 participants. Three of 155 (2%) participants carried possibly damaging heterozygous variants in MCM8, whereas 7 of 151 (5%) individuals carried possibly damaging heterozygous variants in MCM9. One participant carried a novel homozygous variant, c.1651C>T, p.Gln551*, in MCM9, which is predicted to introduce a premature stop codon in exon 9. Biallelic damaging heterozygous variants in both MCM8 and MCM9 were identified in 1 participant. Of a total of 10 participants carrying damaging heterozygous variants in either MCM8 or MCM9, 2 individuals carried heterozygous damaging variants in genes associated with either MCM8 or MCM9 or the DDR pathway. Conclusions: We identified a significant number of potentially damaging and novel variants in MCM8 and MCM9 among participants with POI and examined multiallelic association with variants in DDR and MCM8-MCM9 interactome genes.

Common variants in CLDN14 are associated with differential excretion of magnesium over calcium in urine.

The nature and importance of genetic factors regulating the differential handling of Ca2+ and Mg2+ by the renal tubule in the general population are poorly defined. We conducted a genome-wide meta-analysis of urinary magnesium-to-calcium ratio to identify associated common genetic variants. We included 9320 adults of European descent from four genetic isolates and three urban cohorts. Urinary magnesium and calcium concentrations were measured centrally in spot urine, and each study conducted linear regression analysis of urinary magnesium-to-calcium ratio on ~2.5 million single-nucleotide polymorphisms (SNPs) using an additive model. We investigated, in mouse, the renal expression profile of the top candidate gene and its variation upon changes in dietary magnesium. The genome-wide analysis evidenced a top locus (rs172639, p = 1.7 × 10-12), encompassing CLDN14, the gene coding for claudin-14, that was genome-wide significant when using urinary magnesium-to-calcium ratio, but not either one taken separately. In mouse, claudin-14 is expressed in the distal nephron segments specifically handling magnesium, and its expression is regulated by chronic changes in dietary magnesium content. A genome-wide approach identified common variants in the CLDN14 gene exerting a robust influence on the differential excretion of Mg2+ over Ca2+ in urine. These data highlight the power of urinary electrolyte ratios to unravel genetic determinants of renal tubular function. Coupled with mouse experiments, these results support a major role for claudin-14, a gene associated with kidney stones, in the differential paracellular handling of divalent cations by the renal tubule.

Meta-GWAS and Meta-Analysis of Exome Array Studies Do Not Reveal Genetic Determinants of Serum Hepcidin.

Serum hepcidin concentration is regulated by iron status, inflammation, erythropoiesis and numerous other factors, but underlying processes are incompletely understood. We studied the association of common and rare single nucleotide variants (SNVs) with serum hepcidin in one Italian study and two large Dutch population-based studies. We genotyped common SNVs with genome-wide association study (GWAS) arrays and subsequently performed imputation using the 1000 Genomes reference panel. Cohort-specific GWAS were performed for log-transformed serum hepcidin, adjusted for age and gender, and results were combined in a fixed-effects meta-analysis (total N 6,096). Six top SNVs (p<5x10-6) were genotyped in 3,821 additional samples, but associations were not replicated. Furthermore, we meta-analyzed cohort-specific exome array association results of rare SNVs with serum hepcidin that were available for two of the three cohorts (total N 3,226), but no exome-wide significant signal (p<1.4x10-6) was identified. Gene-based meta-analyses revealed 19 genes that showed significant association with hepcidin. Our results suggest the absence of common SNVs and rare exonic SNVs explaining a large proportion of phenotypic variation in serum hepcidin. We recommend extension of our study once additional substantial cohorts with hepcidin measurements, GWAS and/or exome array data become available in order to increase power to identify variants that explain a smaller proportion of hepcidin variation. In addition, we encourage follow-up of the potentially interesting genes that resulted from the gene-based analysis of low-frequency and rare variants.

A Genome-Wide Association Study in isolated populations reveals new genes associated to common food likings.

Food preferences are the first factor driving food choice and thus nutrition. They involve numerous different senses such as taste and olfaction as well as various other factors such as personal experiences and hedonistic aspects. Although it is clear that several of these have a genetic basis, up to now studies have focused mostly on the effects of polymorphisms of taste receptor genes. Therefore, we have carried out one of the first large scale (4611 individuals) GWAS on food likings assessed for 20 specific food likings belonging to 4 different categories (vegetables, fatty, dairy and bitter). A two-step meta-analysis using three different isolated populations from Italy for the discovery step and two populations from The Netherlands and Central Asia for replication, revealed 15 independent genome-wide significant loci (p < 5 × 10(-8)) for 12 different foods. None of the identified genes coded for either taste or olfactory receptors suggesting that genetics impacts in determining food likings in a much broader way than simple differences in taste perception. These results represent a further step in uncovering the genes that underlie liking of common foods that in the end will greatly help understanding the genetics of human nutrition in general.

Six Novel Loci Associated with Circulating VEGF Levels Identified by a Meta-analysis of Genome-Wide Association Studies.

Vascular endothelial growth factor (VEGF) is an angiogenic and neurotrophic factor, secreted by endothelial cells, known to impact various physiological and disease processes from cancer to cardiovascular disease and to be pharmacologically modifiable. We sought to identify novel loci associated with circulating VEGF levels through a genome-wide association meta-analysis combining data from European-ancestry individuals and using a dense variant map from 1000 genomes imputation panel. Six discovery cohorts including 13,312 samples were analyzed, followed by in-silico and de-novo replication studies including an additional 2,800 individuals. A total of 10 genome-wide significant variants were identified at 7 loci. Four were novel loci (5q14.3, 10q21.3, 16q24.2 and 18q22.3) and the leading variants at these loci were rs114694170 (MEF2C, P = 6.79 x 10(-13)), rs74506613 (JMJD1C, P = 1.17 x 10(-19)), rs4782371 (ZFPM1, P = 1.59 x 10(-9)) and rs2639990 (ZADH2, P = 1.72 x 10(-8)), respectively. We also identified two new independent variants (rs34528081, VEGFA, P = 1.52 x 10(-18); rs7043199, VLDLR-AS1, P = 5.12 x 10(-14)) at the 3 previously identified loci and strengthened the evidence for the four previously identified SNPs (rs6921438, LOC100132354, P = 7.39 x 10(-1467); rs1740073, C6orf223, P = 2.34 x 10(-17); rs6993770, ZFPM2, P = 2.44 x 10(-60); rs2375981, KCNV2, P = 1.48 x 10(-100)). These variants collectively explained up to 52% of the VEGF phenotypic variance. We explored biological links between genes in the associated loci using Ingenuity Pathway Analysis that emphasized their roles in embryonic development and function. Gene set enrichment analysis identified the ERK5 pathway as enriched in genes containing VEGF associated variants. eQTL analysis showed, in three of the identified regions, variants acting as both cis and trans eQTLs for multiple genes. Most of these genes, as well as some of those in the associated loci, were involved in platelet biogenesis and functionality, suggesting the importance of this process in regulation of VEGF levels. This work also provided new insights into the involvement of genes implicated in various angiogenesis related pathologies in determining circulating VEGF levels. The understanding of the molecular mechanisms by which the identified genes affect circulating VEGF levels could be important in the development of novel VEGF-related therapies for such diseases.

Genome sequencing elucidates Sardinian genetic architecture and augments association analyses for lipid and blood inflammatory markers.

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.

Multicohort analysis of the maternal age effect on recombination.

Several studies have reported that the number of crossovers increases with maternal age in humans, but others have found the opposite. Resolving the true effect has implications for understanding the maternal age effect on aneuploidies. Here, we revisit this question in the largest sample to date using single nucleotide polymorphism (SNP)-chip data, comprising over 6,000 meioses from nine cohorts. We develop and fit a hierarchical model to allow for differences between cohorts and between mothers. We estimate that over 10 years, the expected number of maternal crossovers increases by 2.1% (95% credible interval (0.98%, 3.3%)). Our results are not consistent with the larger positive and negative effects previously reported in smaller cohorts. We see heterogeneity between cohorts that is likely due to chance effects in smaller samples, or possibly to confounders, emphasizing that care should be taken when interpreting results from any specific cohort about the effect of maternal age on recombination.