Cilia are important organelles formed by cell membrane protrusions; however, little is known about their regulation by membrane lipids. We characterize a novel activation mechanism for glycogen synthase kinase-3 (GSK3) by the sphingolipids phytoceramide and ceramide that is critical for ciliogenesis in Chlamydomonas and murine ependymal cells, respectively. We show for the first time that Chlamydomonas expresses serine palmitoyl transferase (SPT), the first enzyme in (phyto)ceramide biosynthesis. Inhibition of SPT in Chlamydomonas by myriocin led to loss of flagella and reduced tubulin acetylation, which was prevented by supplementation with the precursor dihydrosphingosine. Immunocytochemistry showed that (phyto)ceramide was colocalized with phospho-Tyr-216-GSK3 (pYGSK3) at the base and tip of Chlamydomonas flagella and motile cilia in ependymal cells. The (phyto)ceramide distribution was consistent with that of a bifunctional ceramide analogue UV cross-linked and visualized by click-chemistry-mediated fluorescent labeling. Ceramide depletion, by myriocin or neutral sphingomyelinase deficiency (fro/fro mouse), led to GSK3 dephosphorylation and defective flagella and cilia. Motile cilia were rescued and pYGSK3 localization restored by incubation of fro/fro ependymal cells with exogenous C24:1 ceramide, which directly bound to pYGSK3. Our findings suggest that (phyto)ceramide-mediated translocation of pYGSK into flagella and cilia is an evolutionarily conserved mechanism fundamental to the regulation of ciliogenesis.

In the International Visible Trait Genetics (VisiGen) Consortium, we investigated the genetics of human skin color by combining a series of genome-wide association studies (GWAS) in a total of 17,262 Europeans with functional follow-up of discovered loci. Our GWAS provide the first genome-wide significant evidence for chromosome 20q11.22 harboring the ASIP gene being explicitly associated with skin color in Europeans. In addition, genomic loci at 5p13.2 (SLC45A2), 6p25.3 (IRF4), 15q13.1 (HERC2/OCA2), and 16q24.3 (MC1R) were confirmed to be involved in skin coloration in Europeans. In follow-up gene expression and regulation studies of 22 genes in 20q11.22, we highlighted two novel genes EIF2S2 and GSS, serving as competing functional candidates in this region and providing future research lines. A genetically inferred skin color score obtained from the 9 top-associated SNPs from 9 genes in 940 worldwide samples (HGDP-CEPH) showed a clear gradual pattern in Western Eurasians similar to the distribution of physical skin color, suggesting the used 9 SNPs as suitable markers for DNA prediction of skin color in Europeans and neighboring populations, relevant in future forensic and anthropological investigations.

Inter-individual variation in facial shape is one of the most noticeable phenotypes in humans, and it is clearly under genetic regulation; however, almost nothing is known about the genetic basis of normal human facial morphology. We therefore conducted a genome-wide association study for facial shape phenotypes in multiple discovery and replication cohorts, considering almost ten thousand individuals of European descent from several countries. Phenotyping of facial shape features was based on landmark data obtained from three-dimensional head magnetic resonance images (MRIs) and two-dimensional portrait images. We identified five independent genetic loci associated with different facial phenotypes, suggesting the involvement of five candidate genes--PRDM16, PAX3, TP63, C5orf50, and COL17A1--in the determination of the human face. Three of them have been implicated previously in vertebrate craniofacial development and disease, and the remaining two genes potentially represent novel players in the molecular networks governing facial development. Our finding at PAX3 influencing the position of the nasion replicates a recent GWAS of facial features. In addition to the reported GWA findings, we established links between common DNA variants previously associated with NSCL/P at 2p21, 8q24, 13q31, and 17q22 and normal facial-shape variations based on a candidate gene approach. Overall our study implies that DNA variants in genes essential for craniofacial development contribute with relatively small effect size to the spectrum of normal variation in human facial morphology. This observation has important consequences for future studies aiming to identify more genes involved in the human facial morphology, as well as for potential applications of DNA prediction of facial shape such as in future forensic applications.

Individuals with exceptional longevity and their offspring have significantly larger high-density lipoprotein concentrations (HDL-C) particle sizes due to the increased homozygosity for the I405V variant in the cholesteryl ester transfer protein (CETP) gene. In this study, we investigate the association of CETP and HDL-C further to identify novel, independent CETP variants associated with HDL-C in humans.

Disordered gambling is a moderately heritable trait, but the underlying genetic basis is largely unknown. We performed a genome-wide association study (GWAS) for disordered gambling using a quantitative factor score in 1312 twins from 894 Australian families. Association was conducted for 2 381 914 single-nucleotide polymorphisms (SNPs) using the family-based association test in Merlin followed by gene and pathway enrichment analyses. Although no SNP reached genome-wide significance, six achieved P-values < 1 × 10(-5) with variants in three genes (MT1X, ATXN1 and VLDLR) implicated in disordered gambling. Secondary case-control analyses found two SNPs on chromosome 9 (rs1106076 and rs12305135 near VLDLR) and rs10812227 near FZD10 on chromosome 12 to be significantly associated with lifetime Diagnostic and Statistical Manual of Mental Disorders, fourth edition pathological gambling and South Oaks Gambling Screen classified probable pathological gambling status. Furthermore, several addiction-related pathways were enriched for SNPs associated with disordered gambling. Finally, gene-based analysis of 24 candidate genes for dopamine agonist-induced gambling in individuals with Parkinson's disease suggested an enrichment of SNPs associated with disordered gambling. We report the first GWAS of disordered gambling. While further replication is required, the identification of susceptibility loci and biological pathways will be important in characterizing the biological mechanisms that underpin disordered gambling.

The perceived taste intensities of quinine HCl, caffeine, sucrose octaacetate (SOA) and propylthiouracil (PROP) solutions were examined in 1457 twins and their siblings. Previous heritability modeling of these bitter stimuli indicated a common genetic factor for quinine, caffeine and SOA (22-28%), as well as separate specific genetic factors for PROP (72%) and quinine (15%). To identify the genes involved, we performed a genome-wide association study with the same sample as the modeling analysis, genotyped for approximately 610,000 single-nucleotide polymorphisms (SNPs). For caffeine and SOA, no SNP association reached a genome-wide statistical criterion. For PROP, the peak association was within TAS2R38 (rs713598, A49P, P = 1.6 × 10(-104)), which accounted for 45.9% of the trait variance. For quinine, the peak association was centered in a region that contains bitter receptor as well as salivary protein genes and explained 5.8% of the trait variance (TAS2R19, rs10772420, R299C, P = 1.8 × 10(-15)). We confirmed this association in a replication sample of twins of similar ancestry (P = 0.00001). The specific genetic factor for the perceived intensity of PROP was identified as the gene previously implicated in this trait (TAS2R38). For quinine, one or more bitter receptor or salivary proline-rich protein genes on chromosome 12 have alleles which affect its perception but tight linkage among very similar genes precludes the identification of a single causal genetic variant.

Chronic kidney disease (CKD) affects ~10% of the global population, with considerable ethnic differences in prevalence and aetiology. We assemble genome-wide association studies of estimated glomerular filtration rate (eGFR), a measure of kidney function that defines CKD, in 312,468 individuals of diverse ancestry. We identify 127 distinct association signals with homogeneous effects on eGFR across ancestries and enrichment in genomic annotations including kidney-specific histone modifications. Fine-mapping reveals 40 high-confidence variants driving eGFR associations and highlights putative causal genes with cell-type specific expression in glomerulus, and in proximal and distal nephron. Mendelian randomisation supports causal effects of eGFR on overall and cause-specific CKD, kidney stone formation, diastolic blood pressure and hypertension. These results define novel molecular mechanisms and putative causal genes for eGFR, offering insight into clinical outcomes and routes to CKD treatment development.

C-reactive protein (CRP) is a sensitive biomarker of chronic low-grade inflammation and is associated with multiple complex diseases. The genetic determinants of chronic inflammation remain largely unknown, and the causal role of CRP in several clinical outcomes is debated. We performed two genome-wide association studies (GWASs), on HapMap and 1000 Genomes imputed data, of circulating amounts of CRP by using data from 88 studies comprising 204,402 European individuals. Additionally, we performed in silico functional analyses and Mendelian randomization analyses with several clinical outcomes. The GWAS meta-analyses of CRP revealed 58 distinct genetic loci (p < 5 × 10-8). After adjustment for body mass index in the regression analysis, the associations at all except three loci remained. The lead variants at the distinct loci explained up to 7.0% of the variance in circulating amounts of CRP. We identified 66 gene sets that were organized in two substantially correlated clusters, one mainly composed of immune pathways and the other characterized by metabolic pathways in the liver. Mendelian randomization analyses revealed a causal protective effect of CRP on schizophrenia and a risk-increasing effect on bipolar disorder. Our findings provide further insights into the biology of inflammation and could lead to interventions for treating inflammation and its clinical consequences.

The skin's tendency to sunburn rather than tan is a major risk factor for skin cancer. Here we report a large genome-wide association study of ease of skin tanning in 176,678 subjects of European ancestry. We identify significant association with tanning ability at 20 loci. We confirm previously identified associations at six of these loci, and report 14 novel loci, of which ten have never been associated with pigmentation-related phenotypes. Our results also suggest that variants at the AHR/AGR3 locus, previously associated with cutaneous malignant melanoma the underlying mechanism of which is poorly understood, might act on disease risk through modulation of tanning ability.

We show that in Madin-Darby canine kidney (MDCK) cells, an apical ceramide-enriched compartment (ACEC) at the base of primary cilia is colocalized with Rab11a. Ceramide and Rab11a vesicles isolated by magnetic sorting contain a highly similar profile of proteins (atypical protein kinase C [aPKC], Cdc42, Sec8, Rab11a, and Rab8) and ceramide species, suggesting the presence of a ciliogenic protein complex associated with ceramide at the ACEC. It is intriguing that C16 and C18 ceramide, although less abundant ceramide species in MDCK cells, are highly enriched in ceramide and Rab11a vesicles. Expression of a ceramide-binding but dominant-negative mutant of aPKC suppresses ciliogenesis, indicating that the association of ceramide with aPKC is critical for the formation of this complex. Our results indicate that ciliogenic ceramide is derived from apical sphingomyelin (SM) that is endocytosed and then converted to the ACEC. Consistently, inhibition of acid sphingomyelinase with imipramine disrupts ACEC formation, association of ciliogenic proteins with Rab11a vesicles, and cilium formation. Ciliogenesis is rescued by the histone deacetylase (HDAC) inhibitor trichostatin A, indicating that ceramide promotes tubulin acetylation in cilia. Taken together, our results suggest that the ACEC is a novel compartment in which SM-derived ceramide induces formation of a ciliogenic lipid-protein complex that sustains primary cilia by preventing deacetylation of microtubules.

The use of spoken and written language is a fundamental human capacity. Individual differences in reading- and language-related skills are influenced by genetic variation, with twin-based heritability estimates of 30 to 80% depending on the trait. The genetic architecture is complex, heterogeneous, and multifactorial, but investigations of contributions of single-nucleotide polymorphisms (SNPs) were thus far underpowered. We present a multicohort genome-wide association study (GWAS) of five traits assessed individually using psychometric measures (word reading, nonword reading, spelling, phoneme awareness, and nonword repetition) in samples of 13,633 to 33,959 participants aged 5 to 26 y. We identified genome-wide significant association with word reading (rs11208009, P = 1.098 × 10-8) at a locus that has not been associated with intelligence or educational attainment. All five reading-/language-related traits showed robust SNP heritability, accounting for 13 to 26% of trait variability. Genomic structural equation modeling revealed a shared genetic factor explaining most of the variation in word/nonword reading, spelling, and phoneme awareness, which only partially overlapped with genetic variation contributing to nonword repetition, intelligence, and educational attainment. A multivariate GWAS of word/nonword reading, spelling, and phoneme awareness maximized power for follow-up investigation. Genetic correlation analysis with neuroimaging traits identified an association with the surface area of the banks of the left superior temporal sulcus, a brain region linked to the processing of spoken and written language. Heritability was enriched for genomic elements regulating gene expression in the fetal brain and in chromosomal regions that are depleted of Neanderthal variants. Together, these results provide avenues for deciphering the biological underpinnings of uniquely human traits.

Vitamin D deficiency is a candidate risk factor for a range of adverse health outcomes. In a genome-wide association study of 25 hydroxyvitamin D (25OHD) concentration in 417,580 Europeans we identify 143 independent loci in 112 1-Mb regions, providing insights into the physiology of vitamin D and implicating genes involved in lipid and lipoprotein metabolism, dermal tissue properties, and the sulphonation and glucuronidation of 25OHD. Mendelian randomization models find no robust evidence that 25OHD concentration has causal effects on candidate phenotypes (e.g. BMI, psychiatric disorders), but many phenotypes have (direct or indirect) causal effects on 25OHD concentration, clarifying the epidemiological relationship between 25OHD status and the health outcomes examined in this study.

The human face represents a combined set of highly heritable phenotypes, but knowledge on its genetic architecture remains limited, despite the relevance for various fields. A series of genome-wide association studies on 78 facial shape phenotypes quantified from 3-dimensional facial images of 10,115 Europeans identified 24 genetic loci reaching study-wide suggestive association (p < 5 × 10-8), among which 17 were previously unreported. A follow-up multi-ethnic study in additional 7917 individuals confirmed 10 loci including six unreported ones (padjusted < 2.1 × 10-3). A global map of derived polygenic face scores assembled facial features in major continental groups consistent with anthropological knowledge. Analyses of epigenomic datasets from cranial neural crest cells revealed abundant cis-regulatory activities at the face-associated genetic loci. Luciferase reporter assays in neural crest progenitor cells highlighted enhancer activities of several face-associated DNA variants. These results substantially advance our understanding of the genetic basis underlying human facial variation and provide candidates for future in-vivo functional studies.

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.

Human perception of bitter substances is partially genetically determined. Previously we discovered a single nucleotide polymorphism (SNP) within the cluster of bitter taste receptor genes on chromosome 12 that accounts for 5.8% of the variance in the perceived intensity rating of quinine, and we strengthened the classic association between TAS2R38 genotype and the bitterness of propylthiouracil (PROP). Here we performed a genome-wide association study (GWAS) using a 40% larger sample (n = 1999) together with a bivariate approach to detect previously unidentified common variants with small effects on bitter perception.

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within-family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

Scores on an optimistic-pessimistic personality scale have been associated with mortality, but optimism and pessimism scores are separable traits and it is unclear which has effects on health or longevity. The Life Orientation Test (LOT), containing items for optimism and pessimism, was included in a twin study on health of Australians aged over 50 in 1993-1995. After a mean of 20 years, participants were matched against death information from the Australian National Death Index. 1,068 out of 2,978 participants with useable LOT scores had died. Survival analysis tested for associations between separate optimism and pessimism scores and mortality from any cause, and from cancers, cardiovascular diseases or other known causes. Age-adjusted scores on the pessimism scale were associated with all-cause and cardiovascular mortality (Hazard Ratios per 1 standard deviation unit, 95% confidence intervals and p-values 1.134, 1.065-1.207, 8.85 × 10-5 and 1.196, 1.045-1.368, 0.0093, respectively) but not with cancer deaths. Optimism scores, which were only weakly correlated with pessimism scores (age-adjusted rank correlation = - 0.176), did not show significant associations with overall or cause-specific mortality. Reverse causation (disease causing pessimism) is unlikely because in that case both cardiovascular diseases and cancers would be expected to lead to pessimism.

We analyzed genome-wide association studies (GWASs), including data from 71,638 individuals from four ancestries, for estimated glomerular filtration rate (eGFR), a measure of kidney function used to define chronic kidney disease (CKD). We identified 20 loci attaining genome-wide-significant evidence of association (p < 5 × 10(-8)) with kidney function and highlighted that allelic effects on eGFR at lead SNPs are homogeneous across ancestries. We leveraged differences in the pattern of linkage disequilibrium between diverse populations to fine-map the 20 loci through construction of "credible sets" of variants driving eGFR association signals. Credible variants at the 20 eGFR loci were enriched for DNase I hypersensitivity sites (DHSs) in human kidney cells. DHS credible variants were expression quantitative trait loci for NFATC1 and RGS14 (at the SLC34A1 locus) in multiple tissues. Loss-of-function mutations in ancestral orthologs of both genes in Drosophila melanogaster were associated with altered sensitivity to salt stress. Renal mRNA expression of Nfatc1 and Rgs14 in a salt-sensitive mouse model was also reduced after exposure to a high-salt diet or induced CKD. Our study (1) demonstrates the utility of trans-ethnic fine mapping through integration of GWASs involving diverse populations with genomic annotation from relevant tissues to define molecular mechanisms by which association signals exert their effect and (2) suggests that salt sensitivity might be an important marker for biological processes that affect kidney function and CKD in humans.

The sphingosine-1-phosphate (S1P) transporter Spns2 regulates myocardial precursor migration in zebrafish and lymphocyte trafficking in mice. However, its function in cancer has not been investigated. We show here that ectopic Spns2 expression induced apoptosis and its knockdown enhanced cell migration in non-small cell lung cancer (NSCLC) cells. Metabolically, Spns2 expression increased the extracellular S1P level while its knockdown the intracellular. Pharmacological inhibition of S1P synthesis abolished the augmented cell migration mediated by Spns2 knockdown, indicating that intracellular S1P plays a key role in this process. Cell signaling studies indicated that Spns2 expression impaired GSK-3β and Stat3 mediated pro-survival pathways. Conversely, these pathways were activated by Spns2 knockdown, which explains the increased cell migration since they are also crucial for migration. Alterations of Spns2 were found to affect several enzymes involved in S1P metabolism, including sphingosine kinases, S1P phosphatases, and S1P lyase 1. Genetically, Spns2 mRNA level was found to be reduced in advanced lung cancer (LC) patients as quantified by using a small scale qPCR array. These data show for the first time that Spns2 plays key roles in regulating the cellular functions in NSCLC cells, and that its down-regulation is a potential risk factor for LC.

Lead is an environmental pollutant that causes acute and chronic toxicity. Surveys have related mean blood lead concentrations to exogenous sources, including industrial activity, use of lead-based paints, or traffic density. However, there has been little investigation of individual differences in lead absorption, distribution, or toxicity, or of genetic causes of such variation.