GWAS of human bitter taste perception identifies new loci and reveals additional complexity of bitter taste genetics.

Human perception of bitterness displays pronounced interindividual variation. This phenotypic variation is mirrored by equally pronounced genetic variation in the family of bitter taste receptor genes. To better understand the effects of common genetic variations on human bitter taste perception, we conducted a genome-wide association study on a discovery panel of 504 subjects and a validation panel of 104 subjects from the general population of São Paulo in Brazil. Correction for general taste-sensitivity allowed us to identify a SNP in the cluster of bitter taste receptors on chr12 (10.88- 11.24 Mb, build 36.1) significantly associated (best SNP: rs2708377, P = 5.31 × 10(-13), r(2) = 8.9%, β = -0.12, s.e. = 0.016) with the perceived bitterness of caffeine. This association overlaps with-but is statistically distinct from-the previously identified SNP rs10772420 influencing the perception of quinine bitterness that falls in the same bitter taste cluster. We replicated this association to quinine perception (P = 4.97 × 10(-37), r(2) = 23.2%, β = 0.25, s.e. = 0.020) and additionally found the effect of this genetic locus to be concentration specific with a strong impact on the perception of low, but no impact on the perception of high concentrations of quinine. Our study, thus, furthers our understanding of the complex genetic architecture of bitter taste perception.

Association of eGFR-Related Loci Identified by GWAS with Incident CKD and ESRD.

Family studies suggest a genetic component to the etiology of chronic kidney disease (CKD) and end stage renal disease (ESRD). Previously, we identified 16 loci for eGFR in genome-wide association studies, but the associations of these single nucleotide polymorphisms (SNPs) for incident CKD or ESRD are unknown. We thus investigated the association of these loci with incident CKD in 26,308 individuals of European ancestry free of CKD at baseline drawn from eight population-based cohorts followed for a median of 7.2 years (including 2,122 incident CKD cases defined as eGFR <60ml/min/1.73m(2) at follow-up) and with ESRD in four case-control studies in subjects of European ancestry (3,775 cases, 4,577 controls). SNPs at 11 of the 16 loci (UMOD, PRKAG2, ANXA9, DAB2, SHROOM3, DACH1, STC1, SLC34A1, ALMS1/NAT8, UBE2Q2, and GCKR) were associated with incident CKD; p-values ranged from p = 4.1e-9 in UMOD to p = 0.03 in GCKR. After adjusting for baseline eGFR, six of these loci remained significantly associated with incident CKD (UMOD, PRKAG2, ANXA9, DAB2, DACH1, and STC1). SNPs in UMOD (OR = 0.92, p = 0.04) and GCKR (OR = 0.93, p = 0.03) were nominally associated with ESRD. In summary, the majority of eGFR-related loci are either associated or show a strong trend towards association with incident CKD, but have modest associations with ESRD in individuals of European descent. Additional work is required to characterize the association of genetic determinants of CKD and ESRD at different stages of disease progression.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment.

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Polygenic determinants of white matter volume derived from GWAS lack reproducibility in a replicate sample

A recent publication reported an exciting polygenic effect of schizophrenia (SCZ) risk variants, identified by a large genome-wide association study (GWAS), on total brain and white matter volumes in schizophrenic patients and, even more prominently, in healthy subjects. The aim of the present work was to replicate and then potentially extend these findings. According to the original publication, polygenic risk scores using single nucleotide polymorphism (SNP) information of SCZ GWAS (polygenic SCZ risk scores; PSS) were calculated in 122 healthy subjects, enrolled in a structural magnetic resonance imaging (MRI) study. These scores were computed based on P-values and odds ratios available through the Psychiatric GWAS Consortium. In addition, polygenic white matter scores (PWM) were calculated, using the respective SNP subset in the original publication. None of the polygenic scores, either PSS or PWM, were found to be associated with total brain, white matter or gray matter volume in our replicate sample. Minor differences between the original and the present study that might have contributed to lack of reproducibility (but unlikely explain it fully), are number of subjects, ethnicity, age distribution, array technology, SNP imputation quality and MRI scanner type. In contrast to the original publication, our results do not reveal the slightest signal of association of the described sets of GWAS-identified SCZ risk variants with brain volumes in adults. Caution is indicated in interpreting studies building on polygenic risk scores without replication sample.

Using Prior Information from the Medical Literature in GWAS of Oral Cancer Identifies Novel Susceptibility Variant on Chromosome 4-the AdAPT Method

Background: Genome-wide association studies (GWAS) require large sample sizes to obtain adequate statistical power, but it may be possible to increase the power by incorporating complementary data. In this study we investigated the feasibility of automatically retrieving information from the medical literature and leveraging this information in GWAS. Methods: We developed a method that searches through PubMed abstracts for pre-assigned keywords and key concepts, and uses this information to assign prior probabilities of association for each single nucleotide polymorphism (SNP) with the phenotype of interest - the Adjusting Association Priors with Text (AdAPT) method. Association results from a GWAS can subsequently be ranked in the context of these priors using the Bayes False Discovery Probability (BFDP) framework. We initially tested AdAPT by comparing rankings of known susceptibility alleles in a previous lung cancer GWAS, and subsequently applied it in a two-phase GWAS of oral cancer. Results: Known lung cancer susceptibility SNPs were consistently ranked higher by AdAPT BFDPs than by p-values. In the oral cancer GWAS, we sought to replicate the top five SNPs as ranked by AdAPT BFDPs, of which rs991316, located in the ADH gene region of 4q23, displayed a statistically significant association with oral cancer risk in the replication phase (per-rare-allele log additive p-value [p(trend)] = 2.5 x 10(-3)). The combined OR for having one additional rare allele was 0.83 (95% CI: 0.76-0.90), and this association was independent of previously identified susceptibility SNPs that are associated with overall UADT cancer in this gene region. We also investigated if rs991316 was associated with other cancers of the upper aerodigestive tract (UADT), but no additional association signal was found. Conclusion: This study highlights the potential utility of systematically incorporating prior knowledge from the medical literature in genome-wide analyses using the AdAPT methodology. AdAPT is available online (url: http://services.gate.ac.uk/lld/gwas/service/config).

Association of GWAS loci with PD in China

Genome-wide association studies (GWAS) have identified numerous single-nucleotide polymorphisms (SNPs) at four loci (SNCA, PARK16, LRRK2, BST1) that can modulate the risk of Parkinson's disease (PD). The strength of these associations has yet to be clarified in Mainland China. Ethnic specific effect is an important consideration in GWAS analysis. Using a case-control methodology, we genotyped multiple SNPs at these four loci to investigate their association with risk of PD in Mainland China. A total of 1,146 study subjects comprising 636 patients with PD and 510 unrelated healthy controls were recruited. The minor alleles at SNPs rs894278, rs1994090, rs2046932, rs4698412, and rs7304279 were found to be significantly higher in cases than in controls, while the minor alleles were found to significantly reduce the risk of developing PD at SNPs rs823128, rs823156, rs6532194, rs1191532, and rs16856139. These associations remained after taking into considerations the effects of age and gender. We showed that multiple SNPs at LRRK2 and SNCA increase risk of PD, while PARK16 SNPs are associated with a lower risk of PD in China. Our study findings will contribute to further research using GWAS-linked data and research on ethnic specific effect of common variants.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

THE EXPRESSION PROFILE OF A GWAS DETECTED GENE (NINJ2) IN THE BRAIN OF A MOUSE MODEL OF ISCHEMIC STROKE

Stroke is the third leading cause of death and a major debilitating disease in the United States. Multiple factors, including genetic factors, contribute to the development of the disease. Genome-wide association studies (GWAS) have contributed to the identification of genetic loci influencing risk for complex diseases, such as stroke. In 2010, a GWAS of incident stroke was performed in four large prospective cohorts from the USA and Europe and identified an association of two Single Nucleotide Polymorphisms (SNPs) on chromosome 12p13 with a greater risk of ischemic stroke in individuals of European and African-American ancestry. These SNPs are located 11 Kb upstream of the nerve injury-induced gene 2, Ninjurin2 (NINJ2), suggesting that this gene may be involved in stroke pathogenesis. NINJ2 is a cell adhesion molecule induced in the distal glial cells from a sciatic-nerve injury at 7-days after injury. In an effort to ascribe a possible role of NINJ2 in stroke, we have assessed changes in the level of gene and protein expression of NINJ2 following a time-course from a transiently induced middle cerebral artery ischemic stroke in mice brains. We report an increase in the gene expression of NINJ2 in the ischemic and peri-infarct (ipsilateral) cortical tissues at 7 and 14-days after stroke. We also report an increase in the protein expression of NINJ2 in the cortex of both the ipsilateral and contralateral cortical tissues at the same time-points. We conclude that the expression of NINJ2 is regulated by an ischemic stroke in the cortex and is consistent with NINJ2 being involved in the recovery time-points of stroke.

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

Identificação de regiões genômicas por GWAS em trigo associadas à tolerância à seca - primeiros resultados.

2016

GWAS of a maize Phosphorus-Starvation Tolerance 1 homolog expression in a tropical diverse panel.

2016

Gwas analysis for the identification of molecular basis responsible for yield related traits and disease resistance in durum wheat

The PhD thesis was developed in the framework of Innovar H2020 project. This project aimed at using genomics, transcriptomics and phenotyping techniques to update varietal registration procedure used in Europe for Value of Cultivation and Use (VCU) and Distinctiness Uniformity and Stability (DUS) protocols. The phenotypic and genotypic diversity of a durum wheat panel were assessed for different agronomic traits, connected with wheat development, disease resistance and spike fertility. A panel of 253 durum wheat varieties was characterized for VCU and DUS traits and genotyped with Illumina 90K SNP Chip array (Wang et al., 2014). GWAS analysis was performed, detecting strong QTLs confirmed also by literature review. Candidate genes were identified for each trait and molecular markers will be developed to be used for marker assisted selection in breeding programs. As for disease resistance, the panel was evaluated for resistance to Soil-Borne-Cereal-Mosaic-Virus (SBCMV). A major QTL, sbm2, was detected on chromosome 2B responsible for durum wheat resistance (Maccaferri et al., 2011). The sbm2 interval was explored by fine mapping on segregant population using KASP markers and by RNASeq analysis, detecting candidate genes involved in plant-pathogen reaction. As regards yield related traits, detailed analysis was performed on the GNI-2A QTL (Milner et al., 2016), responsible for increased number spike fertility. Fine mapping analysis was performed on durum panel identifying hox2 a strong candidate gene, codifying for transcription factor protein. The gene is paralogue of GNI-1 (Sakuma et al., 2019), and it has a 4 kbp deletion responsible for increased number of florets per spikelet. To conclude, the herein reported thesis shows a complete characterization of agronomic and disease resistance traits in modern durum wheat varieties. The results obtained will augment available information for each variety, identifying informative molecular markers for breeding purposes and QTLs/candidate genes responsible for different agronomic traits.

An empirical evaluation of imputation accuracy for association statistics reveals increased type-I error rates in genome-wide associations

Background: Genome wide association studies (GWAS) are becoming the approach of choice to identify genetic determinants of complex phenotypes and common diseases. The astonishing amount of generated data and the use of distinct genotyping platforms with variable genomic coverage are still analytical challenges. Imputation algorithms combine directly genotyped markers information with haplotypic structure for the population of interest for the inference of a badly genotyped or missing marker and are considered a near zero cost approach to allow the comparison and combination of data generated in different studies. Several reports stated that imputed markers have an overall acceptable accuracy but no published report has performed a pair wise comparison of imputed and empiric association statistics of a complete set of GWAS markers. Results: In this report we identified a total of 73 imputed markers that yielded a nominally statistically significant association at P < 10(-5) for type 2 Diabetes Mellitus and compared them with results obtained based on empirical allelic frequencies. Interestingly, despite their overall high correlation, association statistics based on imputed frequencies were discordant in 35 of the 73 (47%) associated markers, considerably inflating the type I error rate of imputed markers. We comprehensively tested several quality thresholds, the haplotypic structure underlying imputed markers and the use of flanking markers as predictors of inaccurate association statistics derived from imputed markers. Conclusions: Our results suggest that association statistics from imputed markers showing specific MAF (Minor Allele Frequencies) range, located in weak linkage disequilibrium blocks or strongly deviating from local patterns of association are prone to have inflated false positive association signals. The present study highlights the potential of imputation procedures and proposes simple procedures for selecting the best imputed markers for follow-up genotyping studies.

Genetic Background of Patients from a University Medical Center in Manhattan: Implications for Personalized Medicine

Background: The rapid progress currently being made in genomic science has created interest in potential clinical applications; however, formal translational research has been limited thus far. Studies of population genetics have demonstrated substantial variation in allele frequencies and haplotype structure at loci of medical relevance and the genetic background of patient cohorts may often be complex. Methods and Findings: To describe the heterogeneity in an unselected clinical sample we used the Affymetrix 6.0 gene array chip to genotype self-identified European Americans (N = 326), African Americans (N = 324) and Hispanics (N = 327) from the medical practice of Mount Sinai Medical Center in Manhattan, NY. Additional data from US minority groups and Brazil were used for external comparison. Substantial variation in ancestral origin was observed for both African Americans and Hispanics; data from the latter group overlapped with both Mexican Americans and Brazilians in the external data sets. A pooled analysis of the African Americans and Hispanics from NY demonstrated a broad continuum of ancestral origin making classification by race/ethnicity uninformative. Selected loci harboring variants associated with medical traits and drug response confirmed substantial within-and between-group heterogeneity. Conclusion: As a consequence of these complementary levels of heterogeneity group labels offered no guidance at the individual level. These findings demonstrate the complexity involved in clinical translation of the results from genome-wide association studies and suggest that in the genomic era conventional racial/ethnic labels are of little value.

A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p <= 5 x 10(-7)). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1 x 10(-8)) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2 x 10(-8)) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5 x 10(-8); rs1229984-ADH1B, p = 7 x 10(-9); and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.