Genome-wide association meta-analysis identifies 48 risk variants and highlights the role of the stria vascularis in hearing loss.

Hearing loss is one of the top contributors to years lived with disability and is a risk factor for dementia. Molecular evidence on the cellular origins of hearing loss in humans is growing. Here, we performed a genome-wide association meta-analysis of clinically diagnosed and self-reported hearing impairment on 723,266 individuals and identified 48 significant loci, 10 of which are novel. A large proportion of associations comprised missense variants, half of which lie within known familial hearing loss loci. We used single-cell RNA-sequencing data from mouse cochlea and brain and mapped common-variant genomic results to spindle, root, and basal cells from the stria vascularis, a structure in the cochlea necessary for normal hearing. Our findings indicate the importance of the stria vascularis in the mechanism of hearing impairment, providing future paths for developing targets for therapeutic intervention in hearing loss.

Using genetic variation to disentangle the complex relationship between food intake and health outcomes.

Diet is considered as one of the most important modifiable factors influencing human health, but efforts to identify foods or dietary patterns associated with health outcomes often suffer from biases, confounding, and reverse causation. Applying Mendelian randomization in this context may provide evidence to strengthen causality in nutrition research. To this end, we first identified 283 genetic markers associated with dietary intake in 445,779 UK Biobank participants. We then converted these associations into direct genetic effects on food exposures by adjusting them for effects mediated via other traits. The SNPs which did not show evidence of mediation were then used for MR, assessing the association between genetically predicted food choices and other risk factors, health outcomes. We show that using all associated SNPs without omitting those which show evidence of mediation, leads to biases in downstream analyses (genetic correlations, causal inference), similar to those present in observational studies. However, MR analyses using SNPs which have only a direct effect on the exposure on food exposures provided unequivocal evidence of causal associations between specific eating patterns and obesity, blood lipid status, and several other risk factors and health outcomes.

TAS2R38 Genotype Does Not Affect SARS-CoV-2 Infection in Primary Ciliary Dyskinesia.

Several chronic respiratory diseases could be risk factors for acquiring SARS-CoV-2 infection: among them, Primary Ciliary Dyskinesia (PCD) is a rare (about 1:10.000) inherited ciliopathy (MIM 242650) characterized by recurrent upper and lower respiratory tract infections due to a dysfunction of the respiratory cilia. In this study, we aimed to investigate whether PCD subjects are more susceptible to infection by SARS-CoV-2 and whether some polymorphisms of the TAS2R38 bitter taste receptor correlate with an increased prevalence of SARS-CoV-2 infection and severity of symptoms. Patients answered several questions about possible SARS-CoV-2 infection, experienced symptoms, and vaccinations; in the case of infection, they also filled out a SNOT-22 questionnaire and ARTIQ. Forty PCD adult patients (mean age, 36.6 ± 16.7 years; 23 females, 17 males) participated in this study, out of which 30% had tested positive for COVID-19 during the last four years; most of them reported a mildly symptomatic disease. We found no differences in age or sex, but a statistically significant difference (p = 0.03) was observed in body mass index (BMI), which was higher in the COVID-acquired group (23.2 ± 3.3 vs. 20.1 ± 4.1 kg/m2). Genotyping for TAS2R38 polymorphisms showed a prevalence of 28.6% PAV/PAV, 48.6% PAV/AVI, and 22.8% AVI/AVI individuals in our cohort. In contrast to our hypothesis, we did not observe a protective role of the PAV allele towards SARS-CoV-2 infection. Conclusions: Our findings suggest that subjects with PCD may not be at increased risk of severe outcomes from COVID-19 and the TAS2R38 bitter taste receptor genotype does not affect SARS-CoV-2 infection.

Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin.

BACKGROUND: Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown. METHODS: We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing. RESULTS: Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion. CONCLUSIONS: Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

Directional dominance on stature and cognition in diverse human populations.

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 second, general cognitive ability and educational attainment (P < 1 × 10(-300), 2.1 × 10(-6), 2.5 × 10(-10) and 1.8 × 10(-10), respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months' less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Estimation of metabolic syndrome heritability in three large populations including full pedigree and genomic information.

Metabolic syndrome is a complex human disorder characterized by a cluster of conditions (increased blood pressure, hyperglycemia, excessive body fat around the waist, and abnormal cholesterol or triglyceride levels). Any of these conditions increases the risk of serious disorders such as diabetes or cardiovascular disease. Currently, the degree of genetic regulation of this syndrome is under debate and partially unknown. The principal aim of this study was to estimate the genetic component and the common environmental effects in different populations using full pedigree and genomic information. We used three large populations (Gubbio, ARIC, and Ogliastra cohorts) to estimate the heritability of metabolic syndrome. Due to both pedigree and genotyped data, different approaches were applied to summarize relatedness conditions. Linear mixed models (LLM) using average information restricted maximum likelihood (AIREML) algorithm were applied to partition the variances and estimate heritability (h2) and common sib-household effect (c2). Globally, results obtained from pedigree information showed a significant heritability (h2: 0.286 and 0.271 in Gubbio and Ogliastra, respectively), whereas a lower, but still significant heritability was found using SNPs data ([Formula: see text]: 0.167 and 0.254 in ARIC and Ogliastra). The remaining heritability between h2 and [Formula: see text] ranged between 0.031 and 0.237. Finally, the common environmental c2 in Gubbio and Ogliastra were also significant accounting for about 11% of the phenotypic variance. Availability of different kinds of populations and data helped us to better understand what happened when heritability of metabolic syndrome is estimated and account for different possible confounding. Furthermore, the opportunity of comparing different results provided more precise and less biased estimation of heritability.

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.

Investigation of the link between PROP taste perception and vegetables consumption using FAOSTAT data.

In this work we investigated, in populations located in Central Asia, the relationship between PROP taste perception and vegetables liking and consumption using FAOSTAT dataset. Collected data were analysed using distance matrices, Mantel test and Pearson correlation. Populations showing similar ability in tasting PROP bitterness are more similar as respect to vegetable consumption (r = 0.63, p-value = .05). Moreover, a significant negative correlation was found between the percentage of Non Taster (NT) in different countries and the percentage of vegetable consumption (r = -0.87, p-value = .02), while a significant positive correlation emerged between the percentage of Super Taster (ST) and the percentage of vegetable liking (r = 0.87, p-value = .02). In our work we showed that differences in bitter perception among populations contributes to differences in vegetable liking and vegetable consumption. More in detail, populations with higher percentage of ST consume more vegetables than population where the majority of individuals are NT.

A genome-wide association study of corneal astigmatism: The CREAM Consortium.

Purpose: To identify genes and genetic markers associated with corneal astigmatism. Methods: A meta-analysis of genome-wide association studies (GWASs) of corneal astigmatism undertaken for 14 European ancestry (n=22,250) and 8 Asian ancestry (n=9,120) cohorts was performed by the Consortium for Refractive Error and Myopia. Cases were defined as having >0.75 diopters of corneal astigmatism. Subsequent gene-based and gene-set analyses of the meta-analyzed results of European ancestry cohorts were performed using VEGAS2 and MAGMA software. Additionally, estimates of single nucleotide polymorphism (SNP)-based heritability for corneal and refractive astigmatism and the spherical equivalent were calculated for Europeans using LD score regression. Results: The meta-analysis of all cohorts identified a genome-wide significant locus near the platelet-derived growth factor receptor alpha (PDGFRA) gene: top SNP: rs7673984, odds ratio=1.12 (95% CI:1.08-1.16), p=5.55×10-9. No other genome-wide significant loci were identified in the combined analysis or European/Asian ancestry-specific analyses. Gene-based analysis identified three novel candidate genes for corneal astigmatism in Europeans-claudin-7 (CLDN7), acid phosphatase 2, lysosomal (ACP2), and TNF alpha-induced protein 8 like 3 (TNFAIP8L3). Conclusions: In addition to replicating a previously identified genome-wide significant locus for corneal astigmatism near the PDGFRA gene, gene-based analysis identified three novel candidate genes, CLDN7, ACP2, and TNFAIP8L3, that warrant further investigation to understand their role in the pathogenesis of corneal astigmatism. The much lower number of genetic variants and genes demonstrating an association with corneal astigmatism compared to published spherical equivalent GWAS analyses suggest a greater influence of rare genetic variants, non-additive genetic effects, or environmental factors in the development of astigmatism.

NFAT5 and SLC4A10 Loci Associate with Plasma Osmolality.

Disorders of water balance, an excess or deficit of total body water relative to body electrolyte content, are common and ascertained by plasma hypo- or hypernatremia, respectively. We performed a two-stage genome-wide association study meta-analysis on plasma sodium concentration in 45,889 individuals of European descent (stage 1 discovery) and 17,637 additional individuals of European descent (stage 2 replication), and a transethnic meta-analysis of replicated single-nucleotide polymorphisms in 79,506 individuals (63,526 individuals of European descent, 8765 individuals of Asian Indian descent, and 7215 individuals of African descent). In stage 1, we identified eight loci associated with plasma sodium concentration at P<5.0 × 10-6 Of these, rs9980 at NFAT5 replicated in stage 2 meta-analysis (P=3.1 × 10-5), with combined stages 1 and 2 genome-wide significance of P=5.6 × 10-10 Transethnic meta-analysis further supported the association at rs9980 (P=5.9 × 10-12). Additionally, rs16846053 at SLC4A10 showed nominally, but not genome-wide, significant association in combined stages 1 and 2 meta-analysis (P=6.7 × 10-8). NFAT5 encodes a ubiquitously expressed transcription factor that coordinates the intracellular response to hypertonic stress but was not previously implicated in the regulation of systemic water balance. SLC4A10 encodes a sodium bicarbonate transporter with a brain-restricted expression pattern, and variant rs16846053 affects a putative intronic NFAT5 DNA binding motif. The lead variants for NFAT5 and SLC4A10 are cis expression quantitative trait loci in tissues of the central nervous system and relevant to transcriptional regulation. Thus, genetic variation in NFAT5 and SLC4A10 expression and function in the central nervous system may affect the regulation of systemic water balance.

Salt-inducible kinase 3, SIK3, is a new gene associated with hearing.

Hearing function is known to be heritable, but few significant and reproducible associations of genetic variants have been identified to date in the adult population. In this study, genome-wide association results of hearing function from the G-EAR consortium and TwinsUK were used for meta-analysis. Hearing ability in eight population samples of Northern and Southern European ancestry (n = 4591) and the Silk Road (n = 348) was measured using pure-tone audiometry and summarized using principal component (PC) analysis. Genome-wide association analyses for PC1-3 were conducted separately in each sample assuming an additive model adjusted for age, sex and relatedness of subjects. Meta-analysis was performed using 2.3 million single-nucleotide polymorphisms (SNPs) tested against each of the three PCs of hearing ability in 4939 individuals. A single SNP lying in intron 6 of the salt-inducible kinase 3 (SIK3) gene was found to be associated with hearing PC2 (P = 3.7×10(-8)) and further supported by whole-genome sequence in a subset. To determine the relevance of this gene in the ear, expression of the Sik3 protein was studied in mouse cochlea of different ages. Sik3 was expressed in murine hair cells during early development and in cells of the spiral ganglion during early development and adulthood. Our results suggest a developmental role of Sik3 in hearing and may be required for the maintenance of adult auditory function.

Application of a new method for GWAS in a related case/control sample with known pedigree structure: identification of new loci for nephrolithiasis.

In contrast to large GWA studies based on thousands of individuals and large meta-analyses combining GWAS results, we analyzed a small case/control sample for uric acid nephrolithiasis. Our cohort of closely related individuals is derived from a small, genetically isolated village in Sardinia, with well-characterized genealogical data linking the extant population up to the 16(th) century. It is expected that the number of risk alleles involved in complex disorders is smaller in isolated founder populations than in more diverse populations, and the power to detect association with complex traits may be increased when related, homogeneous affected individuals are selected, as they are more likely to be enriched with and share specific risk variants than are unrelated, affected individuals from the general population. When related individuals are included in an association study, correlations among relatives must be accurately taken into account to ensure validity of the results. A recently proposed association method uses an empirical genotypic covariance matrix estimated from genome-screen data to allow for additional population structure and cryptic relatedness that may not be captured by the genealogical data. We apply the method to our data, and we also investigate the properties of the method, as well as other association methods, in our highly inbred population, as previous applications were to outbred samples. The more promising regions identified in our initial study in the genetic isolate were then further investigated in an independent sample collected from the Italian population. Among the loci that showed association in this study, we observed evidence of a possible involvement of the region encompassing the gene LRRC16A, already associated to serum uric acid levels in a large meta-analysis of 14 GWAS, suggesting that this locus might lead a pathway for uric acid metabolism that may be involved in gout as well as in nephrolithiasis.

Clenching the Strings of Bruxism Etiopathogenesis: Association Analyses on Genetics and Environmental Risk Factors in a Deeply Characterized Italian Cohort.

Bruxism is a worldwide oral health problem. Although there is a consensus about its multifactorial nature, its precise etiopathogenetic mechanisms are unclear. This study, taking advantage of a deeply characterized cohort of 769 individuals (aged 6-89 years) coming from Northern Italy's genetically isolated populations, aims to epidemiologically describe environmental risk factors for bruxism development and identify genes potentially involved through a Genome-Wide Association Study (GWAS) approach. Logistic mixed models adjusted for age and sex were performed to evaluate associations between bruxism and possible risk factors, e.g., anxiety, smoking, and alcohol and caffeine intake. A case-control GWAS (135 cases, 523 controls), adjusted for age, sex, and anxiety, was conducted to identify new candidate genes. The GTEx data analysis was performed to evaluate the identified gene expression in human body tissues. Statistical analyses determined anxiety as a bruxism risk factor (OR = 2.54; 95% CI: 1.20-5.38; p-value = 0.015), and GWAS highlighted three novel genes potentially associated with bruxism: NLGN1 (topSNP = rs2046718; p-value = 2.63 × 10-7), RIMBP2 (topSNP = rs571497947; p-value = 4.68 × 10-7), and LHFP (topSNP = rs2324342; p-value = 7.47 × 10-6). The GTEx data analysis showed their expression in brain tissues. Overall, this work provided a deeper understanding of bruxism etiopathogenesis with the long-term perspective of developing personalized therapeutic approaches for improving affected individuals' quality of life.

The Bittersweet Symphony of COVID-19: Associations between TAS1Rs and TAS2R38 Genetic Variations and COVID-19 Symptoms.

The innate immune system is crucial in fighting SARS-CoV-2 infection, which is responsible for coronavirus disease 2019 (COVID-19). Therefore, deepening our understanding of the underlying immune response mechanisms is fundamental for the development of novel therapeutic strategies. The role of extra-oral bitter (TAS2Rs) and sweet (TAS1Rs) taste receptors in immune response regulation has yet to be fully understood. However, a few studies have investigated the association between taste receptor genes and COVID-19 symptom severity, with controversial results. Therefore, this study aims to deepen the relationship between COVID-19 symptom presence/severity and TAS1R and TAS2R38 (TAS2Rs member) genetic variations in a cohort of 196 COVID-19 patients. Statistical analyses detected significant associations between rs307355 of the TAS1R3 gene and the following COVID-19-related symptoms: chest pain and shortness of breath. Specifically, homozygous C/C patients are exposed to an increased risk of manifesting severe forms of chest pain (OR 8.11, 95% CI 2.26-51.99) and shortness of breath (OR 4.83, 95% CI 1.71-17.32) in comparison with T/C carriers. Finally, no significant associations between the TAS2R38 haplotype and the presence/severity of COVID-19 symptoms were detected. This study, taking advantage of a clinically and genetically characterised cohort of COVID-19 patients, revealed TAS1R3 gene involvement in determining COVID-19 symptom severity independently of TAS2R38 activity, thus providing novel insights into the role of TAS1Rs in regulating the immune response to viral infections.

Genome-Wide Interaction Analyses of Serum Calcium on Ventricular Repolarization Time in 125 393 Participants.

BACKGROUND: Ventricular repolarization time (ECG QT and JT intervals) is associated with malignant arrhythmia. Genome-wide association studies have identified 230 independent loci for QT and JT; however, 50% of their heritability remains unexplained. Previous work supports a causal effect of lower serum calcium concentrations on longer ventricular repolarization time. We hypothesized calcium interactions with QT and JT variant associations could explain a proportion of the missing heritability. METHODS AND RESULTS: We performed genome-wide calcium interaction analyses for QT and JT intervals. Participants were stratified by their calcium level relative to the study distribution (top or bottom 20%). We performed a 2-stage analysis (genome-wide discovery [N=62 532] and replication [N=59 861] of lead variants) and a single-stage genome-wide meta-analysis (N=122 393, [European ancestry N=117 581, African ancestry N=4812]). We also calculated 2-degrees of freedom joint main and interaction and 1-degree of freedom interaction P values. In 2-stage and single-stage analyses, 50 and 98 independent loci, respectively, were associated with either QT or JT intervals (2-degrees of freedom joint main and interaction P value <5×10-8). No lead variant had a significant interaction result after correcting for multiple testing and sensitivity analyses provided similar findings. Two loci in the single-stage meta-analysis were not reported previously (SPPL2B and RFX6). CONCLUSIONS: We have found limited support for an interaction effect of serum calcium on QT and JT variant associations despite sample sizes with suitable power to detect relevant effects. Therefore, such effects are unlikely to explain a meaningful proportion of the heritability of QT and JT, and factors including rare variation and other environmental interactions need to be considered.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum.

Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Regulator of G-Protein Signalling 9: A New Candidate Gene for Sweet Food Liking?

Genetics plays an important role in individual differences in food liking, which influences food choices and health. Sweet food liking is a complex trait and has been associated with increased body mass index (BMI) and related comorbidities. This genome-wide association study (GWAS) aimed to investigate the genetics of sweet food liking using two adult discovery cohorts (n = 1109, n = 373) and an independent replication cohort (n = 1073). In addition, we tested the association of our strongest result on parameters related to behaviour (food adventurousness (FA) and reward dependence (RD) and health status (BMI and blood glucose). The results demonstrate a novel strong association between the Regulator of G-Protein Signalling 9 (RGS9I) gene, strongest single nucleotide polymorphism (SNP) rs58931966 (p-value 7.05 × 10-9 in the combined sample of discovery and replication), and sweet food liking, with the minor allele (A) being associated with a decreased sweet food liking. We also found that the A allele of the rs58931966 SNP was associated with decreased FA and RD, and increased BMI and blood glucose (p-values < 0.05). Differences were highlighted in sex-specific analysis on BMI and glucose. Our results highlight a novel genetic association with food liking and are indicative of genetic variation influencing the psychological-biological drivers of food preference. If confirmed in other studies, such genetic associations could allow a greater understanding of chronic disease management from both a habitual dietary intake and reward-related perspective.

Genome-wide association studies on Northern Italy isolated populations provide further support concerning genetic susceptibility for major depressive disorder.

OBJECTIVES: Major depressive disorder (MDD) is a psychiatric disorder with pathogenesis influenced by both genetic and environmental factors. To date, the molecular-level understanding of its aetiology remains unclear. Thus, we aimed to identify genetic variants and susceptibility genes for MDD with a genome-wide association study (GWAS) approach. METHODS: We performed a meta-analysis of GWASs and a gene-based analysis on two Northern Italy isolated populations (cases/controls n = 166/472 and 33/320), followed by replication and polygenic risk score (PRS) analyses in Italian independent samples (cases n = 464, controls n = 339). RESULTS: We identified two novel MDD-associated genes, KCNQ5 (lead SNP rs867262, p = 3.82 × 10-9) and CTNNA2 (rs6729523, p = 1.25 × 10-8). The gene-based analysis revealed another six genes (p < 2.703 × 10-6): GRM7, CTNT4, SNRK, SRGAP3, TRAPPC9, and FHIT. No replication of the genome-wide significant SNPs was found in the independent cohort, even if 14 SNPs around CTNNA2 showed association with MDD and related phenotypes at the nominal level of p (<0.05). Furthermore, the PRS model developed in the discovery cohort discriminated cases and controls in the replication cohort. CONCLUSIONS: Our work suggests new possible genes associated with MDD, and the PRS analysis confirms the polygenic nature of this disorder. Future studies are required to better understand the role of these findings in MDD.

Identifying missing pieces in color vision defects: a genome-wide association study in Silk Road populations.

Introduction: Color vision defects (CVDs) are conditions characterized by the alteration of normal trichromatic vision. CVDs can arise as the result of alterations in three genes (OPN1LW, OPN1MW, OPN1SW) or as a combination of genetic predisposition and environmental factors. To date, apart from Mendelian CVDs forms, nothing is known about multifactorial CVDs forms. Materials and Methods: Five hundred and twenty individuals from Silk Road isolated communities were genotyped and phenotypically characterized for CVDs using the Farnsworth D-15 color test. The CVDs traits Deutan-Protan (DP) and Tritan (TR) were analysed. Genome Wide Association Study for both traits was performed, and results were corrected with a False Discovery Rate linkage-based approach (FDR-p). Gene expression of final candidates was investigated using a published human eye dataset, and pathway analysis was performed. Results: Concerning DP, three genes: PIWIL4 (FDR-p: 9.01*10-9), MBD2 (FDR-p: 4.97*10-8) and NTN1 (FDR-p: 4.98*10-8), stood out as promising candidates. PIWIL4 is involved in the preservation of Retinal Pigmented Epithelium (RPE) homeostasis while MBD2 and NTN1 are both involved in visual signal transmission. With regards to TR, four genes: VPS54 (FDR-p: 4.09*10-9), IQGAP (FDR-p: 6,52*10-10), NMB (FDR-p: 8.34*10-11), and MC5R (FDR-p: 2.10*10-8), were considered promising candidates. VPS54 is reported to be associated with Retinitis pigmentosa; IQGAP1 is reported to regulate choroidal vascularization in Age-Related Macular Degeneration; NMB is involved in RPE homeostasis regulation; MC5R is reported to regulate lacrimal gland function. Discussion: Overall, these results provide novel insights regarding a complex phenotype (i.e., CVDs) in an underrepresented population such as Silk Road isolated communities.