From musk to body odor: Decoding olfaction through genetic variation

The olfactory system combines input from multiple receptor types to represent odor information, but there are few explicit examples relating olfactory receptor (OR) activity patterns to odor perception. To uncover these relationships, we performed genome-wide scans on odor-perception phenotypes for ten odors in 1000 Han Chinese and validated results for six of these odors in an ethnically diverse population (n = 364). In both populations, consistent with previous studies, we replicated three previously reported associations (β-ionone/OR5A1, androstenone/OR7D4, cis-3-hexen-1-ol/OR2J3 LD-band), but not for odors containing aldehydes, suggesting that olfactory phenotype/genotype studies are robust across populations. Two novel associations between an OR and odor perception contribute to our understanding of olfactory coding. First, we found a SNP in OR51B2 that associated with trans-3-methyl-2-hexenoic acid, a key component of human underarm odor. Second, we found two linked SNPs associated with the musk Galaxolide in a novel musk receptor, OR4D6, which is also the first human OR shown to drive specific anosmia to a musk compound. We noticed that SNPs detected for odor intensity were enriched with amino acid substitutions, implying functional changes of odor receptors. Furthermore, we also found that the derived alleles of the SNPs tend to be associated with reduced odor intensity, supporting the hypothesis that the primate olfactory gene repertoire has degenerated over time. This study provides information about coding for human body odor, and gives us insight into broader mechanisms of olfactory coding, such as how differential OR activation can converge on a similar percept.     

Mutual cross-adaptation of the volatile steroid androstenone and a non-steroid perceptual analog

Self-and cross-adaptation are believed to result from stimulationof the same olfactory sensory channels. These adaptation phenomenawere studied after exposures to 5-androst-16-en-3-one (androstenone)and a synthetic perceptual analog (DMCMC), viz. a racemic mixtureof the isomers 4(R)-(4',4'-dimethyl-cyclohexyl)-2(R)-methylcyclohexanoneand 4(S)-(4',4'-dimethylcyclohexyl)-2(S)-methylcyclohexanone.In Experiment 1, six subjects very sensitive to androstenonereceived four randomized sequences of six concentrations offour ordants (androstenone, DMCMC, amyl acetate, and Galaxolide*;plus blanks) before and following adaptation to either androstenoneor DMCMC. Exposure to each odorant resulted in self-adaptation.Measures of stimulus intensity and identification thresholdrevealed reciprocal cross-adaptation between androstenone andDMCMC, but no cross-adaptation to amyl acetate or Galaxolide.The degree of cross-adaptation was asymmetric; adaptation toDMCMC resulted in more complete adaptation to androstenone thanvice versa. This asymmetry was apparently due to intensity differences;when stimuli were matched for intensity, the asymmetry disappeared(Experiment 2). These results demonstrate cross-adaptation forqualitatively similar, but not dissimilar, odors and suggestthat androstenone and its perceptual analog DMCMC share thesame sensory channels.     

GWAS of follicular lymphoma reveals allelic heterogeneity at 6p21.32 and suggests shared genetic susceptibility with diffuse large B-cell lymphoma

Non-Hodgkin lymphoma (NHL) represents a diverse group of hematological malignancies, of which follicular lymphoma (FL) is a prevalent subtype. A previous genome-wide association study has established a marker, rs10484561 in the human leukocyte antigen (HLA) class II region on 6p21.32 associated with increased FL risk. Here, in a three-stage genome-wide association study, starting with a genome-wide scan of 379 FL cases and 791 controls followed by validation in 1,049 cases and 5,790 controls, we identified a second independent FL-associated locus on 6p21.32, rs2647012 (OR(combined) = 0.64, P(combined) = 2 × 10(-21)) located 962 bp away from rs10484561 (r(2)<0.1 in controls). After mutual adjustment, the associations at the two SNPs remained genome-wide significant (rs2647012:OR(adjusted) = 0.70, P(adjusted) = 4 × 10(-12); rs10484561:OR(adjusted) = 1.64, P(adjusted) = 5 × 10(-15)). Haplotype and coalescence analyses indicated that rs2647012 arose on an evolutionarily distinct haplotype from that of rs10484561 and tags a novel allele with an opposite (protective) effect on FL risk. Moreover, in a follow-up analysis of the top 6 FL-associated SNPs in 4,449 cases of other NHL subtypes, rs10484561 was associated with risk of diffuse large B-cell lymphoma (OR(combined) = 1.36, P(combined) = 1.4 × 10(-7)). Our results reveal the presence of allelic heterogeneity within the HLA class II region influencing FL susceptibility and indicate a possible shared genetic etiology with diffuse large B-cell lymphoma. These findings suggest that the HLA class II region plays a complex yet important role in NHL.     

WNT16 influences bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk

We aimed to identify genetic variants associated with cortical bone thickness (CBT) and bone mineral density (BMD) by performing two separate genome-wide association study (GWAS) meta-analyses for CBT in 3 cohorts comprising 5,878 European subjects and for BMD in 5 cohorts comprising 5,672 individuals. We then assessed selected single-nucleotide polymorphisms (SNPs) for osteoporotic fracture in 2,023 cases and 3,740 controls. Association with CBT and forearm BMD was tested for ～2.5 million SNPs in each cohort separately, and results were meta-analyzed using fixed effect meta-analysis. We identified a missense SNP (Thr>Ile; rs2707466) located in the WNT16 gene (7q31), associated with CBT (effect size of -0.11 standard deviations [SD] per C allele, P = 6.2 × 10(-9)). This SNP, as well as another nonsynonymous SNP rs2908004 (Gly>Arg), also had genome-wide significant association with forearm BMD (-0.14 SD per C allele, P = 2.3 × 10(-12), and -0.16 SD per G allele, P = 1.2 × 10(-15), respectively). Four genome-wide significant SNPs arising from BMD meta-analysis were tested for association with forearm fracture. SNP rs7776725 in FAM3C, a gene adjacent to WNT16, was associated with a genome-wide significant increased risk of forearm fracture (OR = 1.33, P = 7.3 × 10(-9)), with genome-wide suggestive signals from the two missense variants in WNT16 (rs2908004: OR = 1.22, P = 4.9 × 10(-6) and rs2707466: OR = 1.22, P = 7.2 × 10(-6)). We next generated a homozygous mouse with targeted disruption of Wnt16. Female Wnt16(-/-) mice had 27% (P<0.001) thinner cortical bones at the femur midshaft, and bone strength measures were reduced between 43%-61% (6.5 × 10(-13)相似文献   

CDKN2B polymorphism is associated with primary open-angle glaucoma (POAG) in the Afro-Caribbean population of Barbados, West Indies

The purpose of this study was to confirm previously reported associations of common variants in or near CDC7/TGFBR3, ZP4, SRBD1, ELOVL5, CAV1/CAV2, TLR4, CDKN2B, CDKN2B-AS1, ATOH7, PLXDC2, TMTC2, SIX1, and CARD10, with primary open angle glaucoma (POAG) in the Afro-Caribbean population of Barbados, West Indies. A total of 437 unrelated subjects from the Barbados Family Study of Open Angle Glaucoma (BFSG), including 272 with POAG and 165 unaffected individuals were included in this study. Eighteen SNPs were genotyped by using the multiplex SNaPshot method. Allelic, genotypic and model-based (dominant, recessive, and additive) associations of the SNPs with POAG were analyzed using Chi-squared tests and logistic regression. SNP rs1063192 (near CDKN2B) was found to be significantly associated with POAG (allelic P = 0.0008, genotypic P = 0.0029), and the minor allele C of rs1063192 was protective against POAG (OR = 0.39; 95%CI = 0.22-0.69). Suggestive association was also noted for rs7916697 (near ATHO7, allelic P = 0.0096, genotypic P = 0.01) with the minor allele being protective (OR = 0.67; 95% CI = 0.50-0.91), although this finding did not withstand correction for multiple testing. However, a significant interactive effect on POAG risk was identified between rs1063192 and rs7916697 (P-interaction = 2.80 × 10(-5)). Individuals with the rs1063192 protective genotype CC or CT and also rs7916697 genotypes GG or GA show a significantly decreased risk of POAG (OR = 0.17, 95%CI: 0.07-0.41). Our study confirms the significant association between SNP rs1063192 (CDKN2B, previously shown to influence vertical cup-to-disc ratio and POAG at 9p21) and POAG in the Afro-Caribbean population of Barbados. The minor allele of rs1063192 interacts with that of rs7916697 (ATOH7)) to reduce POAG risk. Our results also suggest that rs1063912 is a common protective variant for POAG in populations of African as well as European descent.     

Genome-wide association study of treatment refractory schizophrenia in Han Chinese

We report the first genome-wide association study of a joint analysis using 795 Han Chinese individuals with treatment-refractory schizophrenia (TRS) and 806 controls. Three loci showed suggestive significant association with TRS were identified. These loci include: rs10218843 (P = 3.04 × 10(-7)) and rs11265461 (P = 1.94 × 10(-7)) are adjacent to signaling lymphocytic activation molecule family member 1 (SLAMF1); rs4699030 (P = 1.94 × 10(-6)) and rs230529 (P = 1.74 × 10(-7)) are located in the gene nuclear factor of kappa light polypeptide gene enhancer in B-cells 1 (NFKB1); and rs13049286 (P = 3.05 × 10(-5)) and rs3827219 (P = 1.66 × 10(-5)) fall in receptor-interacting serine/threonine-protein kinase 4 (RIPK4). One isolated single nucleotide polymorphism (SNP), rs739617 (P = 3.87 × 10(-5)) was also identified to be associated with TRS. The -94delATTG allele (rs28362691) located in the promoter region of NFKB1 was identified by resequencing and was found to associate with TRS (P = 4.85 × 10(-6)). The promoter assay demonstrated that the -94delATTG allele had a significant lower promoter activity than the -94insATTG allele in the SH-SY5Y cells. This study suggests that rs28362691 in NFKB1 might be involved in the development of TRS.     

Genetic variation of an odorant receptor OR7D4 and sensory perception of cooked meat containing androstenone

Although odour perception impacts food preferences, the effect of genotypic variation of odorant receptors (ORs) on the sensory perception of food is unclear. Human OR7D4 responds to androstenone, and genotypic variation in OR7D4 predicts variation in the perception of androstenone. Since androstenone is naturally present in meat derived from male pigs, we asked whether OR7D4 genotype correlates with either the ability to detect androstenone or the evaluation of cooked pork tainted with varying levels of androstenone within the naturally-occurring range. Consistent with previous findings, subjects with two copies of the functional OR7D4 RT variant were more sensitive to androstenone than subjects carrying a non-functional OR7D4 WM variant. When pork containing varying levels of androstenone was cooked and tested by sniffing and tasting, subjects with two copies of the RT variant tended to rate the androstenone-containing meat as less favourable than subjects carrying the WM variant. Our data is consistent with the idea that OR7D4 genotype predicts the sensory perception of meat containing androstenone and that genetic variation in an odorant receptor can alter food preferences.     

Polymorphisms in the NPY2R gene show significant associations with BMI that are additive to FTO, MC4R, and NPFFR2 gene effects

Neuropeptide Y (NPY) is an appetite hormone that acts centrally to control feeding behavior. The 5' and exon 2 regions of NPY2R, one of five NPY receptor genes, have been weakly and inconsistently implicated with obesity. With the ATG start site of the gene at the beginning of exon 2, single-nucleotide polymorphisms (SNPs) across intron 1 may show stronger associations with obesity than expected. Two 5' SNPs, three intron 1 SNPs, and one synonymous exon 2 SNP were genotyped on 2,985 white Utah subjects. Previously associated FTO, NPY, NPY1R, MC4R, PPARGC1A, OR7D4, and four NPFFR2 SNPs were also genotyped and related to BMI. One NPY2R 5' SNP (rs12649641, P = 0.008), an exon 2 SNP (rs2880415, P = 0.009), and an intron 1 SNP (rs17376826, P = 7 × 10(-6)) were each significantly associated with BMI. All three SNPs, plus FTO (rs9939609, P = 1.5 × 10(-6)) and two NPFFR2 SNPs (rs4129733, P = 3.7 × 10(-13) and rs11940196, 4.2 × 10(-10)) remained significant in a multiple regression additive model. Diplotypes using the estimated haplotypes of NPY2R, NPFFR2, and MC4R were significantly associated with BMI (P = 1.0 × 10(-10), 3.2 × 10(-8), and 1.1 × 10(-4), respectively). Haplotypes of NPY2R, NPFFR2, and MC4R, plus the FTO SNP, explained 9.6% of the BMI variance. SNP effect sizes per allele for the four genes ranged from 0.8 to 3.5 kg/m(2). We conclude that haplotypes containing the rs17376826 SNP in intron 1 of NPY2R have strong associations with BMI, some NPFFR2 haplotypes are strongly protective against or increase risk of obesity, and both NPY2R and NPFFR2 play important roles in obesity predisposition independent of FTO and MC4R.     

Mitochondrial C150T polymorphism increases the risk of cervical cancer and HPV infection

During a survey of control region (D-loop) sequence variances in 142 cervical cancer (CC) patients and 136 controls, all Chinese women, including both HPV-positive (human papillomavirus) and HPV-negative subjects, we determined that the C150T polymorphism increased the CC risk in a case-control study (OR=3.0, 95% CI=1.8-5.0, P<0.05). HPV-positive individuals were more likely to carry the C150T polymorphism than HPV-negative controls (OR=5.8, 95% CI=2.6-13.2, P=2.3×10(-5)). HPV-positive CC patients were more likely to carry the C150T polymorphism than HPV-negative controls (OR=4.9, 95% CI=2.6-9.3, P=9.9×10(-7)). In all subjects, an increased risk of HPV infection was also associated with the C150T polymorphism (OR=4.5, 95% CI=2.5-8.1, P=6.6×10(-7)). However, no significant difference in the frequency of other alleles was found at the variable sites in D146, D152, D310 and D514. These results indicated that the C150T polymorphism increased the risk of HPV infection and CC progression. Additionally, we assessed the association of mtDNA copy number with CC risk or the C150T polymorphism in 45 CC patients and 43 controls. There was no significant association of mtDNA copy number with CC risk or the C150T polymorphism. To the best of our knowledge, this is the first report to suggest that mtDNA C150T polymorphism was positively associated with HPV infection and subsequent CC risk among Chinese women.     

Association of common variants in TNFRSF13B, TNFSF13, and ANXA3 with serum levels of non-albumin protein and immunoglobulin isotypes in Japanese

We performed a genome-wide association study (GWAS) on levels of serum total protein (TP), albumin (ALB), and non-albumin protein (NAP). We analyzed SNPs on autosomal chromosomes using data from 9,103 Japanese individuals, followed by a replication study of 1,600 additional individuals. We confirmed the previously- reported association of GCKR on chromosome 2p23.3 with serum ALB (rs1260326, P(meta) = 3.1 × 10(-9)), and additionally identified the significant genome-wide association of rs4985726 in TNFRSF13B on 17p11.2 with both TP and NAP (P(meta) = 1.2 × 10(-14) and 7.1 × 10(-24), respectively). For NAP, rs3803800 and rs11552708 in TNFSF13 on 17p13.1 (P(meta) = 7.2 × 10(-15) and 7.5 × 10(-10), respectively) as well as rs10007186 on 4q21.2 near ANXA3 (P(meta) = 1.3 × 10(-9)) also indicated significant associations. Interestingly, TNFRSF13B and TNFSF13 encode a tumor necrosis factor (TNF) receptor and its ligand, which together constitute an important receptor-ligand axis for B-cell homeostasis and immunoglobulin production. Furthermore, three SNPs, rs4985726, rs3803800, and rs11552708 in TNFRSF13B and TNFSF13, were indicated to be associated with serum levels of IgG (P<2.3 × 10(-3)) and IgM (P<0.018), while rs3803800 and rs11552708 were associated with IgA (P<0.013). Rs10007186 in 4q21.2 was associated with serum levels of IgA (P = 0.036), IgM (P = 0.019), and IgE (P = 4.9 × 10(-4)). Our results should add interesting knowledge about the regulation of major serum components.     

Identification of Odorant-Receptor Interactions by Global Mapping of the Human Odorome

The human olfactory system recognizes a broad spectrum of odorants using approximately 400 different olfactory receptors (hORs). Although significant improvements of heterologous expression systems used to study interactions between ORs and odorant molecules have been made, screening the olfactory repertoire of hORs remains a tremendous challenge. We therefore developed a chemical systems level approach based on protein-protein association network to investigate novel hOR-odorant relationships. Using this new approach, we proposed and validated new bioactivities for odorant molecules and OR2W1, OR51E1 and OR5P3. As it remains largely unknown how human perception of odorants influence or prevent diseases, we also developed an odorant-protein matrix to explore global relationships between chemicals, biological targets and disease susceptibilities. We successfully experimentally demonstrated interactions between odorants and the cannabinoid receptor 1 (CB1) and the peroxisome proliferator-activated receptor gamma (PPARγ). Overall, these results illustrate the potential of integrative systems chemical biology to explore the impact of odorant molecules on human health, i.e. human odorome.     

Identification of a systemic lupus erythematosus susceptibility locus at 11p13 between PDHX and CD44 in a multiethnic study

Systemic lupus erythematosus (SLE) is considered to be the prototypic autoimmune disease, with a complex genetic architecture influenced by environmental factors. We sought to replicate a putative association at 11p13 not yet exceeding genome-wide significance (p < 5 × 10(-8)) identified in a genome-wide association study (GWAS). Our GWA scan identified two intergenic SNPs located between PDHX and CD44 showing suggestive evidence of association with SLE in cases of European descent (rs2732552, p = 0.004, odds ratio [OR] = 0.78; rs387619, p = 0.003, OR = 0.78). The replication cohort consisted of >15,000 subjects, including 3562 SLE cases and 3491 controls of European ancestry, 1527 cases and 1811 controls of African American (AA) descent, and 1265 cases and 1260 controls of Asian origin. We observed robust association at both rs2732552 (p = 9.03 × 10(-8), OR = 0.83) and rs387619 (p = 7.7 × 10(-7), OR = 0.83) in the European samples with p(meta) = 1.82 × 10(-9) for rs2732552. The AA and Asian SLE cases also demonstrated association at rs2732552 (p = 5 × 10(-3), OR = 0.81 and p = 4.3 × 10(-4), OR = 0.80, respectively). A meta-analysis of rs2732552 for all racial and ethnic groups studied produced p(meta) = 2.36 × 10(-13). This locus contains multiple regulatory sites that could potentially affect expression and functions of CD44, a cell-surface glycoprotein influencing immunologic, inflammatory, and oncologic phenotypes, or PDHX, a subunit of the pyruvate dehydrogenase complex.     

An x chromosome association scan of the Norfolk Island genetic isolate provides evidence for a novel migraine susceptibility locus at Xq12

Migraine is a common and debilitating neurovascular disorder with a complex envirogenomic aetiology. Numerous studies have demonstrated a preponderance of women affected with migraine and previous pedigree linkage studies in our laboratory have identified susceptibility loci on chromosome Xq24-Xq28. In this study we have used the genetic isolate of Norfolk Island to further analyse the X chromosome for migraine susceptibility loci.An association approach was employed to analyse 14,124 SNPs spanning the entire X chromosome. Genotype data from 288 individuals comprising a large core-pedigree, of which 76 were affected with migraine, were analysed. Although no SNP reached chromosome-wide significance (empirical α = 1 × 10(-5)) ranking by P-value revealed two primary clusters of SNPs in the top 25. A 10 SNP cluster represents a novel migraine susceptibility locus at Xq12 whilst a 11 SNP cluster represents a previously identified migraine susceptibility locus at Xq27. The strongest association at Xq12 was seen for rs599958 (OR = 1.75, P = 8.92 × 10(-4)), whilst at Xq27 the strongest association was for rs6525667 (OR = 1.53, P = 1.65 × 10(-4)). Further analysis of SNPs at these loci was performed in 5,122 migraineurs from the Women's Genome Health Study and provided additional evidence for association at the novel Xq12 locus (P<0.05).Overall, this study provides evidence for a novel migraine susceptibility locus on Xq12. The strongest effect SNP (rs102834, joint P = 1.63 × 10(-5)) is located within the 5'UTR of the HEPH gene, which is involved in iron homeostasis in the brain and may represent a novel pathway for involvement in migraine pathogenesis.     

Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.     

Genome-wide association study on bipolar disorder in the Bulgarian population

Bipolar disorder is a severe psychiatric disorder influenced by environmental and genetic factors. Genetic studies have implicated many variants in the disease's etiology but only few have been successfully replicated. We conducted a genome-wide association study (GWAS) on bipolar disorder in the Bulgarian population followed by a replication study of the top 100 single nucleotide polymorphisms (SNPs) showing the smallest P values. The GWAS was performed on 188 bipolar disorder patients and 376 control subjects genotyped on the Illumina 550 platform. The replication study was conducted on 122 patients and 328 controls. Although our study did not show any association P value that achieved genome-wide significance, and none of the top 100 SNPs reached the Bonferroni-corrected P value in the replication study, the plausible involvement of some variants cannot be entirely discarded. Three polymorphisms, rs8099939 [P = 2.12 × 10(-6), odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.43-2.67] in GRIK5, rs6122972 (P = 3.11 × 10(-6), OR = 2.02, 95% CI = 1.46-2.80) in PARD6B and rs2289700 (P = 9.14 × 10(-6), OR = 2.13, 95% CI = 1.53-2.95) in CTSH remained associated at a similar level after Mantel-Haenszel test for combining the results from the genome-wide and replication studies. A modest association was also detected for SNP rs1012053 (GWAS P = 4.50 × 10(-2)) in DGKH, which has already been reported as the most significant variant in a previous genome-wide scan on bipolar disorder. However, further studies using larger datasets are needed to identify variants with smaller effects that contribute to the risk of bipolar disorder.     

Genetic elucidation of human hyperosmia to isovaleric acid

The genetic basis of odorant-specific variations in human olfactory thresholds, and in particular of enhanced odorant sensitivity (hyperosmia), remains largely unknown. Olfactory receptor (OR) segregating pseudogenes, displaying both functional and nonfunctional alleles in humans, are excellent candidates to underlie these differences in olfactory sensitivity. To explore this hypothesis, we examined the association between olfactory detection threshold phenotypes of four odorants and segregating pseudogene genotypes of 43 ORs genome-wide. A strong association signal was observed between the single nucleotide polymorphism variants in OR11H7P and sensitivity to the odorant isovaleric acid. This association was largely due to the low frequency of homozygous pseudogenized genotype in individuals with specific hyperosmia to this odorant, implying a possible functional role of OR11H7P in isovaleric acid detection. This predicted receptor–ligand functional relationship was further verified using the Xenopus oocyte expression system, whereby the intact allele of OR11H7P exhibited a response to isovaleric acid. Notably, we also uncovered another mechanism affecting general olfactory acuity that manifested as a significant inter-odorant threshold concordance, resulting in an overrepresentation of individuals who were hyperosmic to several odorants. An involvement of polymorphisms in other downstream transduction genes is one possible explanation for this observation. Thus, human hyperosmia to isovaleric acid is a complex trait, contributed to by both receptor and other mechanisms in the olfactory signaling pathway.