Nicola Pirastu | Università degli Studi di Trieste (original) (raw)

Papers by Nicola Pirastu

Nature, 2015

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darw... more Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one se...

Human Heredity, 2008

length. Haplotypes, block boundaries, and patterns of LD are similar in both populations, suggest... more length. Haplotypes, block boundaries, and patterns of LD are similar in both populations, suggesting a shared distribution of recombination hotspots. Conclusions: About 15% fewer haplotype tagging SNPs need to be genotyped in the isolate, and possibly 70% fewer if selecting SNPs evenly spaced on the metric LD map.

PLoS ONE, 2009

To better design association studies for complex traits in isolated populations it's important to... more To better design association studies for complex traits in isolated populations it's important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not ''a priori'' be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of subpopulations and stratification within genetic isolates.

PLoS ONE, 2014

Taste is one of the main factors determining food choices. Differences in PROP bitter taste perce... more Taste is one of the main factors determining food choices. Differences in PROP bitter taste perception have been implicated in individual differences in food preferences and selection. The present study examined associations between, PROP phenotypes, self-reported food liking and TAS2R38 polymorphisms, the major gene implicated in PROP bitterness, in six different populations of the Caucasus and Central Asia, located along the ancient Silk Road. Differences in the distribution of PROP phenotypes across populations were detected, with a higher frequency of super tasters in Tajikistan (31.3%) and Armenia (39.0%) and a higher frequency of non tasters in Georgia (50.9%). While no relationships were observed between PROP phenotypes and food liking using standard statistical tests, we used an approach based on comparison of distance matrices derived from these data. The first matrix compared the food liking ratings of each population to all others pairwise using the Kruskal-Wallis test (at p,0.00063), and the second one compared the distribution of PROP phenotypes across all populations in a similar manner calculating the chi-square statistic as a distance measure. A strong correlation between the two matrices was found (Mantel test: r = 0.67, p-value = 0.03), suggesting that the pattern of food liking across populations was closely related to the distribution of PROP phenotypes. This same relationship was not observed when TAS2R38 genotypes were substituted for PROP phenotypes in this analysis. Our data suggest that a population-based approach utilizing distance matrices is a useful technique for detecting PROP-related differences in food liking and can be applied to other taste phenotypes. Citation: Robino A, Mezzavilla M, Pirastu N, Dognini M, Tepper BJ, et al. (2014) A Population-Based Approach to Study the Impact of PROP Perception on Food Liking in Populations along the Silk Road. PLoS ONE 9(3): e91716.

PLoS genetics, 2014

Many existing cohorts contain a range of relatedness between genotyped individuals, either by des... more Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

PloS one, 2014

Taste is one of the main factors determining food choices. Differences in PROP bitter taste perce... more Taste is one of the main factors determining food choices. Differences in PROP bitter taste perception have been implicated in individual differences in food preferences and selection. The present study examined associations between, PROP phenotypes, self-reported food liking and TAS2R38 polymorphisms, the major gene implicated in PROP bitterness, in six different populations of the Caucasus and Central Asia, located along the ancient Silk Road. Differences in the distribution of PROP phenotypes across populations were detected, with a higher frequency of super tasters in Tajikistan (31.3%) and Armenia (39.0%) and a higher frequency of non tasters in Georgia (50.9%). While no relationships were observed between PROP phenotypes and food liking using standard statistical tests, we used an approach based on comparison of distance matrices derived from these data. The first matrix compared the food liking ratings of each population to all others pairwise using the Kruskal-Wallis test (at p,0.00063), and the second one compared the distribution of PROP phenotypes across all populations in a similar manner calculating the chi-square statistic as a distance measure. A strong correlation between the two matrices was found (Mantel test: r = 0.67, p-value = 0.03), suggesting that the pattern of food liking across populations was closely related to the distribution of PROP phenotypes. This same relationship was not observed when TAS2R38 genotypes were substituted for PROP phenotypes in this analysis. Our data suggest that a population-based approach utilizing distance matrices is a useful technique for detecting PROP-related differences in food liking and can be applied to other taste phenotypes. Citation: Robino A, Mezzavilla M, Pirastu N, Dognini M, Tepper BJ, et al. (2014) A Population-Based Approach to Study the Impact of PROP Perception on Food Liking in Populations along the Silk Road. PLoS ONE 9(3): e91716.

PloS one, 2014

Coffee, one of the most popular beverages in the world, contains many different physiologically a... more Coffee, one of the most popular beverages in the world, contains many different physiologically active compounds with a potential impact on people's health. Despite the recent attention given to the genetic basis of its consumption, very little has been done in understanding genes influencing coffee preference among different individuals. Given its markedly bitter taste, we decided to verify if bitter receptor genes (TAS2Rs) variants affect coffee liking. In this light, 4066 people from different parts of Europe and Central Asia filled in a field questionnaire on coffee liking. They have been consequently recruited and included in the study. Eighty-eight SNPs covering the 25 TAS2R genes were selected from the available imputed ones and used to run association analysis for coffee liking. A significant association was detected with three SNP: one synonymous and two functional variants (W35S and H212R) on the TAS2R43 gene. Both variants have been shown to greatly reduce in vitro protein activity. Surprisingly the wild type allele, which corresponds to the functional form of the protein, is associated to higher liking of coffee. Since the hTAS2R43 receptor is sensible to caffeine, we verified if the detected variants produced differences in caffeine bitter perception on a subsample of people coming from the FVG cohort. We found a significant association between differences in caffeine perception and the H212R variant but not with the W35S, which suggests that the effect of the TAS2R43 gene on coffee liking is mediated by caffeine and in particular by the H212R variant. No other significant association was found with other TAS2R genes. In conclusion, the present study opens new perspectives in the understanding of coffee liking. Further studies are needed to clarify the role of the TAS2R43 gene in coffee hedonics and to identify which other genes and pathways are involved in its genetics.

PLoS genetics, 2012

Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by... more Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (x 2 = 83.89, df = 1; p = 5.2610 220 ). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Human molecular genetics, 2012

The male-to-female sex ratio at birth is constant across world populations with an average of 1.0... more The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency…

PloS one, 2013

Although several studies demonstrated that platelet count is higher in women, decreases with age,... more Although several studies demonstrated that platelet count is higher in women, decreases with age, and is influenced by genetic background, most clinical laboratories still use the reference interval 150-400610 9 platelets/L for all subjects. The present study was to identify age-and sex-specific reference intervals for platelet count.

Nature, 2012

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth... more Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P , 10 28 , which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Nature genetics, 2013

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthriti... more Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >40,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (8 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

BMC genetics, 2011

Background: Association studies consist in identifying the genetic variants which are related to ... more Background: Association studies consist in identifying the genetic variants which are related to a specific disease through the use of statistical multiple hypothesis testing or segregation analysis in pedigrees. This type of studies has been very successful in the case of Mendelian monogenic disorders while it has been less successful in identifying genetic variants related to complex diseases where the insurgence depends on the interactions between different genes and the environment. The current technology allows to genotype more than a million of markers and this number has been rapidly increasing in the last years with the imputation based on templates sets and whole genome sequencing. This type of data introduces a great amount of noise in the statistical analysis and usually requires a great number of samples. Current methods seldom take into account gene-gene and gene-environment interactions which are fundamental especially in complex diseases. In this paper we propose to use a non-parametric additive model to detect the genetic variants related to diseases which accounts for interactions of unknown order. Although this is not new to the current literature, we show that in an isolated population, where the most related subjects share also most of their genetic code, the use of additive models may be improved if the available genealogical tree is taken into account. Specifically, we form a sample of cases and controls with the highest inbreeding by means of the Hungarian method, and estimate the set of genes/environmental variables, associated with the disease, by means of Random Forest. Results: We have evidence, from statistical theory, simulations and two applications, that we build a suitable procedure to eliminate stratification between cases and controls and that it also has enough precision in identifying genetic variants responsible for a disease. This procedure has been successfully used for the betathalassemia, which is a well known Mendelian disease, and also to the common asthma where we have identified candidate genes that underlie to the susceptibility of the asthma. Some of such candidate genes have been also found related to common asthma in the current literature. Conclusions: The data analysis approach, based on selecting the most related cases and controls along with the Random Forest model, is a powerful tool for detecting genetic variants associated to a disease in isolated populations. Moreover, this method provides also a prediction model that has accuracy in estimating the unknown disease status and that can be generally used to build kit tests for a wide class of Mendelian diseases.

Journal of food science, 2012

Abstract: Food preferences are the main factor driving food intake and choice. There are good re... more Abstract: Food preferences are the main factor driving food intake and choice. There are good reasons to suspect some genetic influence on food acceptance, not least because genetic factors are implicated in a number of factors that are likely to be related to food choice. In addition, some food dislikes show themselves early in life, before there is any evidence for aversive experiences. Although taste has been widely studied in regards of pure tastes such as bitter or sweet perception, the relationship between taste-related genes and food preferences has seldom been explored. In this work we investigated relationship of 37 taste-related genes with food preferences. The study was carried out during a scientific expedition through Caucasus and Central Asia (Silk Road) analyzing more than 400 samples from 5 different countries. A food preference questionnaire was administered to each participant and a DNA sample was obtained. Other information, such as age, sex, life style and anthropometrical measures, were also collected. We found significant associations with variants of: (1) TAS1R2 [Correction added after initial online publication on 27 Aug 2012. TAS1R3 was changed to TAS1R2.] gene and liking of Vodka (P= 1.6 × 10−3), white wine (P= 4.0 × 10−4) and lamb meat (P= 1.6 × 10−3); (2) PCLB2 gene and preference for Hot Tea (P= 8.0 × 10−4); (3) TPRV1 gene and beet liking (P= 3.8 × 10−5); and (4) ITPR3 gene and liking of both lamb meat (5.8 × 10−4) and sheep cheese (8.9×10−4). These findings give a new insight on a better understanding, of genetic factors influencing food preferences which is critical to the development of effective dietary interventions, especially for people that may be genetically not predisposed for liking specific nutrients.

PLoS genetics, 2013

Intracerebral hemorrhage (ICH) is the stroke subtype with the worst prognosis and has no establis... more Intracerebral hemorrhage (ICH) is the stroke subtype with the worst prognosis and has no established acute treatment. ICH is classified as lobar or nonlobar based on the location of ruptured blood vessels within the brain. These different locations also signal different underlying vascular pathologies. Heritability estimates indicate a substantial genetic contribution to risk of ICH in both locations. We report a genome-wide association study of this condition that meta-analyzed data from six studies that enrolled individuals of European ancestry. Case subjects were ascertained by neurologists blinded to genotype data and classified as lobar or nonlobar based on brain computed tomography. ICH-free control subjects were sampled from ambulatory clinics or random digit dialing. Replication of signals identified in the discovery cohort with p < 1 3 10 À6 was pursued in an independent multiethnic sample utilizing both direct and genome-wide genotyping. The discovery phase included a case cohort of 1,545 individuals (664 lobar and 881 nonlobar cases) and a control cohort of 1,481 individuals and identified two susceptibility loci: for lobar ICH, chromosomal region 12q21.1 (rs11179580, odds ratio [OR] ¼ 1.56, p ¼ 7.0 3 10 À8 ); and for nonlobar ICH, chromosomal region 1q22 (rs2984613, OR ¼ 1.44, p ¼ 1.6 3 10 À8 ).

PloS one, 2013

The red blood cell related traits are highly heritable but their genetics are poorly defined. Onl... more The red blood cell related traits are highly heritable but their genetics are poorly defined. Only 5-10% of the total observed variance is explained by the genetic loci found to date, suggesting that additional loci should be searched using approaches alternative to large meta analysis. GWAS (Genome Wide Association Study) for red blood cell traits in a founder population cohort from Northern Italy identified a new locus for mean corpuscular hemoglobin concentration (MCHC) in the TAF3 gene. The association was replicated in two cohorts (rs1887582, P = 4.25E-09). TAF3 encodes a transcription cofactor that participates in core promoter recognition complex, and is involved in zebrafish and mouse erythropoiesis. We show here that TAF3 is required for transcription of the SPTA1 gene, encoding alpha spectrin, one of the proteins that link the plasma membrane to the actin cytoskeleton. Mutations in SPTA1 are responsible for hereditary spherocytosis, a monogenic disorder of MCHC, as well as for the normal MCHC level. Based on our results, we propose that TAF3 is required for normal erythropoiesis in human and that it might have a role in controlling the ratio between hemoglobin (Hb) and cell volume and in the dynamics of RBC maturation in healthy individuals. Finally, TAF3 represents a potential candidate or a modifier gene for disorders of red cell membrane.

Journal of the International AIDS Society, 2014

Introduction: With the aim of searching genetic factors associated with the response to an immune... more Introduction: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment.

PLOS Genetics, 2011

White blood cell (WBC) count is a common clinical measure from complete blood count assays, and i... more White blood cell (WBC) count is a common clinical measure from complete blood count assays, and it varies widely among healthy individuals. Total WBC count and its constituent subtypes have been shown to be moderately heritable, with the heritability estimates varying across cell types. We studied 19,509 subjects from seven cohorts in a discovery analysis, and 11,823 subjects from ten cohorts for replication analyses, to determine genetic factors influencing variability within the normal hematological range for total WBC count and five WBC subtype measures. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We identified and replicated ten associations with total WBC count and five WBC subtypes at seven different genomic loci (total WBC count-6p21 in the HLA region, 17q21 near ORMDL3, and CSF3; neutrophil count-17q21; basophil count-3p21 near RPN1 and C3orf27; lymphocyte count-6p21, 19p13 at EPS15L1; monocyte count-2q31 at ITGA4, 3q21, 8q24 an intergenic region, 9q31 near EDG2), including three previously reported associations and seven novel associations. To investigate functional relationships among variants contributing to variability in the six WBC traits, we utilized gene expression-and pathways-based analyses. We implemented gene-clustering algorithms to evaluate functional connectivity among implicated loci and showed functional relationships across cell types. Gene expression data from whole blood was utilized to show that significant biological consequences can be extracted from our genome-wide analyses, with effect estimates for significant loci from the meta-analyses being highly corellated with the proximal gene expression. In addition, collaborative efforts between the groups contributing to this study and related studies conducted by the COGENT and RIKEN groups allowed for the examination of effect homogeneity for genome-wide significant associations across populations of diverse ancestral backgrounds. (SKG) . These authors contributed equally to this work.

Nature Genetics, 2011

Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood... more Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 × 10(-8) to P = 2.3 × 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.

Nature, 2015

Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darw... more Homozygosity has long been associated with rare, often devastating, Mendelian disorders, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one se...

Human Heredity, 2008

length. Haplotypes, block boundaries, and patterns of LD are similar in both populations, suggest... more length. Haplotypes, block boundaries, and patterns of LD are similar in both populations, suggesting a shared distribution of recombination hotspots. Conclusions: About 15% fewer haplotype tagging SNPs need to be genotyped in the isolate, and possibly 70% fewer if selecting SNPs evenly spaced on the metric LD map.

PLoS ONE, 2009

To better design association studies for complex traits in isolated populations it's important to... more To better design association studies for complex traits in isolated populations it's important to understand how history and isolation moulded the genetic features of different communities. Population isolates should not ''a priori'' be considered homogeneous, even if the communities are not distant and part of a small region. We studied a particular area of Sardinia called Ogliastra, characterized by the presence of several distinct villages that display different history, immigration events and population size. Cultural and geographic isolation characterized the history of these communities. We determined LD parameters in 8 villages and defined population structure through high density SNPs (about 360 K) on 360 unrelated people (45 selected samples from each village). These isolates showed differences in LD values and LD map length. Five of these villages show high LD values probably due to their reduced population size and extreme isolation. High genetic differentiation among villages was detected. Moreover population structure analysis revealed a high correlation between genetic and geographic distances. Our study indicates that history, geography and biodemography have influenced the genetic features of Ogliastra communities producing differences in LD and population structure. All these data demonstrate that we can consider each village an isolate with specific characteristics. We suggest that, in order to optimize the study design of complex traits, a thorough characterization of genetic features is useful to identify the presence of subpopulations and stratification within genetic isolates.

PLoS ONE, 2014

Taste is one of the main factors determining food choices. Differences in PROP bitter taste perce... more Taste is one of the main factors determining food choices. Differences in PROP bitter taste perception have been implicated in individual differences in food preferences and selection. The present study examined associations between, PROP phenotypes, self-reported food liking and TAS2R38 polymorphisms, the major gene implicated in PROP bitterness, in six different populations of the Caucasus and Central Asia, located along the ancient Silk Road. Differences in the distribution of PROP phenotypes across populations were detected, with a higher frequency of super tasters in Tajikistan (31.3%) and Armenia (39.0%) and a higher frequency of non tasters in Georgia (50.9%). While no relationships were observed between PROP phenotypes and food liking using standard statistical tests, we used an approach based on comparison of distance matrices derived from these data. The first matrix compared the food liking ratings of each population to all others pairwise using the Kruskal-Wallis test (at p,0.00063), and the second one compared the distribution of PROP phenotypes across all populations in a similar manner calculating the chi-square statistic as a distance measure. A strong correlation between the two matrices was found (Mantel test: r = 0.67, p-value = 0.03), suggesting that the pattern of food liking across populations was closely related to the distribution of PROP phenotypes. This same relationship was not observed when TAS2R38 genotypes were substituted for PROP phenotypes in this analysis. Our data suggest that a population-based approach utilizing distance matrices is a useful technique for detecting PROP-related differences in food liking and can be applied to other taste phenotypes. Citation: Robino A, Mezzavilla M, Pirastu N, Dognini M, Tepper BJ, et al. (2014) A Population-Based Approach to Study the Impact of PROP Perception on Food Liking in Populations along the Silk Road. PLoS ONE 9(3): e91716.

PLoS genetics, 2014

Many existing cohorts contain a range of relatedness between genotyped individuals, either by des... more Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

PloS one, 2014

Taste is one of the main factors determining food choices. Differences in PROP bitter taste perce... more Taste is one of the main factors determining food choices. Differences in PROP bitter taste perception have been implicated in individual differences in food preferences and selection. The present study examined associations between, PROP phenotypes, self-reported food liking and TAS2R38 polymorphisms, the major gene implicated in PROP bitterness, in six different populations of the Caucasus and Central Asia, located along the ancient Silk Road. Differences in the distribution of PROP phenotypes across populations were detected, with a higher frequency of super tasters in Tajikistan (31.3%) and Armenia (39.0%) and a higher frequency of non tasters in Georgia (50.9%). While no relationships were observed between PROP phenotypes and food liking using standard statistical tests, we used an approach based on comparison of distance matrices derived from these data. The first matrix compared the food liking ratings of each population to all others pairwise using the Kruskal-Wallis test (at p,0.00063), and the second one compared the distribution of PROP phenotypes across all populations in a similar manner calculating the chi-square statistic as a distance measure. A strong correlation between the two matrices was found (Mantel test: r = 0.67, p-value = 0.03), suggesting that the pattern of food liking across populations was closely related to the distribution of PROP phenotypes. This same relationship was not observed when TAS2R38 genotypes were substituted for PROP phenotypes in this analysis. Our data suggest that a population-based approach utilizing distance matrices is a useful technique for detecting PROP-related differences in food liking and can be applied to other taste phenotypes. Citation: Robino A, Mezzavilla M, Pirastu N, Dognini M, Tepper BJ, et al. (2014) A Population-Based Approach to Study the Impact of PROP Perception on Food Liking in Populations along the Silk Road. PLoS ONE 9(3): e91716.

PloS one, 2014

Coffee, one of the most popular beverages in the world, contains many different physiologically a... more Coffee, one of the most popular beverages in the world, contains many different physiologically active compounds with a potential impact on people's health. Despite the recent attention given to the genetic basis of its consumption, very little has been done in understanding genes influencing coffee preference among different individuals. Given its markedly bitter taste, we decided to verify if bitter receptor genes (TAS2Rs) variants affect coffee liking. In this light, 4066 people from different parts of Europe and Central Asia filled in a field questionnaire on coffee liking. They have been consequently recruited and included in the study. Eighty-eight SNPs covering the 25 TAS2R genes were selected from the available imputed ones and used to run association analysis for coffee liking. A significant association was detected with three SNP: one synonymous and two functional variants (W35S and H212R) on the TAS2R43 gene. Both variants have been shown to greatly reduce in vitro protein activity. Surprisingly the wild type allele, which corresponds to the functional form of the protein, is associated to higher liking of coffee. Since the hTAS2R43 receptor is sensible to caffeine, we verified if the detected variants produced differences in caffeine bitter perception on a subsample of people coming from the FVG cohort. We found a significant association between differences in caffeine perception and the H212R variant but not with the W35S, which suggests that the effect of the TAS2R43 gene on coffee liking is mediated by caffeine and in particular by the H212R variant. No other significant association was found with other TAS2R genes. In conclusion, the present study opens new perspectives in the understanding of coffee liking. Further studies are needed to clarify the role of the TAS2R43 gene in coffee hedonics and to identify which other genes and pathways are involved in its genetics.

PLoS genetics, 2012

Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by... more Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (x 2 = 83.89, df = 1; p = 5.2610 220 ). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

Human molecular genetics, 2012

The male-to-female sex ratio at birth is constant across world populations with an average of 1.0... more The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency…

PloS one, 2013

Although several studies demonstrated that platelet count is higher in women, decreases with age,... more Although several studies demonstrated that platelet count is higher in women, decreases with age, and is influenced by genetic background, most clinical laboratories still use the reference interval 150-400610 9 platelets/L for all subjects. The present study was to identify age-and sex-specific reference intervals for platelet count.

Nature, 2012

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth... more Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P , 10 28 , which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Nature genetics, 2013

Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthriti... more Elevated serum urate concentrations can cause gout, a prevalent and painful inflammatory arthritis. By combining data from >40,000 individuals of European ancestry within the Global Urate Genetics Consortium (GUGC), we identified and replicated 28 genome-wide significant loci in association with serum urate concentrations (8 new regions in or near TRIM46, INHBB, SFMBT1, TMEM171, VEGFA, BAZ1B, PRKAG2, STC1, HNF4G, A1CF, ATXN2, UBE2Q2, IGF1R, NFAT5, MAF, HLF, ACVR1B-ACVRL1 and B3GNT4). Associations for many of the loci were of similar magnitude in individuals of non-European ancestry. We further characterized these loci for associations with gout, transcript expression and the fractional excretion of urate. Network analyses implicate the inhibins-activins signaling pathways and glucose metabolism in systemic urate control. New candidate genes for serum urate concentration highlight the importance of metabolic control of urate production and excretion, which may have implications for the treatment and prevention of gout.

BMC genetics, 2011

Background: Association studies consist in identifying the genetic variants which are related to ... more Background: Association studies consist in identifying the genetic variants which are related to a specific disease through the use of statistical multiple hypothesis testing or segregation analysis in pedigrees. This type of studies has been very successful in the case of Mendelian monogenic disorders while it has been less successful in identifying genetic variants related to complex diseases where the insurgence depends on the interactions between different genes and the environment. The current technology allows to genotype more than a million of markers and this number has been rapidly increasing in the last years with the imputation based on templates sets and whole genome sequencing. This type of data introduces a great amount of noise in the statistical analysis and usually requires a great number of samples. Current methods seldom take into account gene-gene and gene-environment interactions which are fundamental especially in complex diseases. In this paper we propose to use a non-parametric additive model to detect the genetic variants related to diseases which accounts for interactions of unknown order. Although this is not new to the current literature, we show that in an isolated population, where the most related subjects share also most of their genetic code, the use of additive models may be improved if the available genealogical tree is taken into account. Specifically, we form a sample of cases and controls with the highest inbreeding by means of the Hungarian method, and estimate the set of genes/environmental variables, associated with the disease, by means of Random Forest. Results: We have evidence, from statistical theory, simulations and two applications, that we build a suitable procedure to eliminate stratification between cases and controls and that it also has enough precision in identifying genetic variants responsible for a disease. This procedure has been successfully used for the betathalassemia, which is a well known Mendelian disease, and also to the common asthma where we have identified candidate genes that underlie to the susceptibility of the asthma. Some of such candidate genes have been also found related to common asthma in the current literature. Conclusions: The data analysis approach, based on selecting the most related cases and controls along with the Random Forest model, is a powerful tool for detecting genetic variants associated to a disease in isolated populations. Moreover, this method provides also a prediction model that has accuracy in estimating the unknown disease status and that can be generally used to build kit tests for a wide class of Mendelian diseases.

Journal of food science, 2012

Abstract: Food preferences are the main factor driving food intake and choice. There are good re... more Abstract: Food preferences are the main factor driving food intake and choice. There are good reasons to suspect some genetic influence on food acceptance, not least because genetic factors are implicated in a number of factors that are likely to be related to food choice. In addition, some food dislikes show themselves early in life, before there is any evidence for aversive experiences. Although taste has been widely studied in regards of pure tastes such as bitter or sweet perception, the relationship between taste-related genes and food preferences has seldom been explored. In this work we investigated relationship of 37 taste-related genes with food preferences. The study was carried out during a scientific expedition through Caucasus and Central Asia (Silk Road) analyzing more than 400 samples from 5 different countries. A food preference questionnaire was administered to each participant and a DNA sample was obtained. Other information, such as age, sex, life style and anthropometrical measures, were also collected. We found significant associations with variants of: (1) TAS1R2 [Correction added after initial online publication on 27 Aug 2012. TAS1R3 was changed to TAS1R2.] gene and liking of Vodka (P= 1.6 × 10−3), white wine (P= 4.0 × 10−4) and lamb meat (P= 1.6 × 10−3); (2) PCLB2 gene and preference for Hot Tea (P= 8.0 × 10−4); (3) TPRV1 gene and beet liking (P= 3.8 × 10−5); and (4) ITPR3 gene and liking of both lamb meat (5.8 × 10−4) and sheep cheese (8.9×10−4). These findings give a new insight on a better understanding, of genetic factors influencing food preferences which is critical to the development of effective dietary interventions, especially for people that may be genetically not predisposed for liking specific nutrients.

PLoS genetics, 2013

Intracerebral hemorrhage (ICH) is the stroke subtype with the worst prognosis and has no establis... more Intracerebral hemorrhage (ICH) is the stroke subtype with the worst prognosis and has no established acute treatment. ICH is classified as lobar or nonlobar based on the location of ruptured blood vessels within the brain. These different locations also signal different underlying vascular pathologies. Heritability estimates indicate a substantial genetic contribution to risk of ICH in both locations. We report a genome-wide association study of this condition that meta-analyzed data from six studies that enrolled individuals of European ancestry. Case subjects were ascertained by neurologists blinded to genotype data and classified as lobar or nonlobar based on brain computed tomography. ICH-free control subjects were sampled from ambulatory clinics or random digit dialing. Replication of signals identified in the discovery cohort with p < 1 3 10 À6 was pursued in an independent multiethnic sample utilizing both direct and genome-wide genotyping. The discovery phase included a case cohort of 1,545 individuals (664 lobar and 881 nonlobar cases) and a control cohort of 1,481 individuals and identified two susceptibility loci: for lobar ICH, chromosomal region 12q21.1 (rs11179580, odds ratio [OR] ¼ 1.56, p ¼ 7.0 3 10 À8 ); and for nonlobar ICH, chromosomal region 1q22 (rs2984613, OR ¼ 1.44, p ¼ 1.6 3 10 À8 ).

PloS one, 2013

The red blood cell related traits are highly heritable but their genetics are poorly defined. Onl... more The red blood cell related traits are highly heritable but their genetics are poorly defined. Only 5-10% of the total observed variance is explained by the genetic loci found to date, suggesting that additional loci should be searched using approaches alternative to large meta analysis. GWAS (Genome Wide Association Study) for red blood cell traits in a founder population cohort from Northern Italy identified a new locus for mean corpuscular hemoglobin concentration (MCHC) in the TAF3 gene. The association was replicated in two cohorts (rs1887582, P = 4.25E-09). TAF3 encodes a transcription cofactor that participates in core promoter recognition complex, and is involved in zebrafish and mouse erythropoiesis. We show here that TAF3 is required for transcription of the SPTA1 gene, encoding alpha spectrin, one of the proteins that link the plasma membrane to the actin cytoskeleton. Mutations in SPTA1 are responsible for hereditary spherocytosis, a monogenic disorder of MCHC, as well as for the normal MCHC level. Based on our results, we propose that TAF3 is required for normal erythropoiesis in human and that it might have a role in controlling the ratio between hemoglobin (Hb) and cell volume and in the dynamics of RBC maturation in healthy individuals. Finally, TAF3 represents a potential candidate or a modifier gene for disorders of red cell membrane.

Journal of the International AIDS Society, 2014

Introduction: With the aim of searching genetic factors associated with the response to an immune... more Introduction: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment.

PLOS Genetics, 2011

White blood cell (WBC) count is a common clinical measure from complete blood count assays, and i... more White blood cell (WBC) count is a common clinical measure from complete blood count assays, and it varies widely among healthy individuals. Total WBC count and its constituent subtypes have been shown to be moderately heritable, with the heritability estimates varying across cell types. We studied 19,509 subjects from seven cohorts in a discovery analysis, and 11,823 subjects from ten cohorts for replication analyses, to determine genetic factors influencing variability within the normal hematological range for total WBC count and five WBC subtype measures. Cohort specific data was supplied by the CHARGE, HeamGen, and INGI consortia, as well as independent collaborative studies. We identified and replicated ten associations with total WBC count and five WBC subtypes at seven different genomic loci (total WBC count-6p21 in the HLA region, 17q21 near ORMDL3, and CSF3; neutrophil count-17q21; basophil count-3p21 near RPN1 and C3orf27; lymphocyte count-6p21, 19p13 at EPS15L1; monocyte count-2q31 at ITGA4, 3q21, 8q24 an intergenic region, 9q31 near EDG2), including three previously reported associations and seven novel associations. To investigate functional relationships among variants contributing to variability in the six WBC traits, we utilized gene expression-and pathways-based analyses. We implemented gene-clustering algorithms to evaluate functional connectivity among implicated loci and showed functional relationships across cell types. Gene expression data from whole blood was utilized to show that significant biological consequences can be extracted from our genome-wide analyses, with effect estimates for significant loci from the meta-analyses being highly corellated with the proximal gene expression. In addition, collaborative efforts between the groups contributing to this study and related studies conducted by the COGENT and RIKEN groups allowed for the examination of effect homogeneity for genome-wide significant associations across populations of diverse ancestral backgrounds. (SKG) . These authors contributed equally to this work.

Nature Genetics, 2011

Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood... more Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 × 10(-8) to P = 2.3 × 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.