Disruption of P2RY5, an orphan G protein-coupled receptor, underlies autosomal recessive woolly hair

The genetic determinants of hair texture in humans are largely unknown. Several human syndromes exist in which woolly hair comprises a part of the phenotype; however, simple autosomal recessive inheritance of isolated woolly hair has only rarely been reported. To identify a gene involved in controlling hair texture, we performed genetic linkage analysis in six families of Pakistani origin with autosomal recessive woolly hair (ARWH; OMIM 278150). All six families showed linkage to chromosome 13q14.2-14.3 (Z = 17.97). In all cases, we discovered pathogenic mutations in P2RY5, which encodes a G protein-coupled receptor and is a nested gene residing within intron 17 of the retinoblastoma 1 (RB1) gene. P2RY5 is expressed in both Henle's and Huxley's layers of the inner root sheath of the hair follicle. Our findings indicate that disruption of P2RY5 underlies ARWH and, more broadly, uncover a new gene involved in determining hair texture in humans.     

Genome-wide association study and mouse model identify interaction between RET and EDNRB pathways in Hirschsprung disease

Genetic studies of Hirschsprung disease, a common congenital malformation, have identified eight genes with mutations that can be associated with this condition. Mutations at individual loci are, however, neither necessary nor sufficient to cause clinical disease. We conducted a genome-wide association study in 43 Mennonite family trios using 2,083 microsatellites and single-nucleotide polymorphisms and a new multipoint linkage disequilibrium method that searches for association arising from common ancestry. We identified susceptibility loci at 10q11, 13q22 and 16q23; the gene at 13q22 is EDNRB, encoding a G protein-coupled receptor (GPCR) and the gene at 10q11 is RET, encoding a receptor tyrosine kinase (RTK). Statistically significant joint transmission of RET and EDNRB alleles in affected individuals and non-complementation of aganglionosis in mouse intercrosses between Ret null and the Ednrb hypomorphic piebald allele are suggestive of epistasis between EDNRB and RET. Thus, genetic interaction between mutations in RET and EDNRB is an underlying mechanism for this complex disorder.     

Molecular evolution of the human interleukin-8 receptor gene cluster.

Interleukin-8 (IL-8) is the prototype for a family of at least eight neutrophil chemoattractants whose genes map to human chromosome 4q13-q21. Two human IL-8 receptors, IL8RA and IL8RB, are known from cDNA cloning; IL8RA is a promiscuous receptor for at least two other related ligands, GRO alpha and NAP-2. We now report cloning of the genes for IL8RA, IL8RB and a recently inactivated pseudogene of receptor A (IL8RAP). These form a cluster of only three genes in the superfamily of G protein-coupled receptors (GPCRs) and map to 2q34-q35. The coevolutionary diversity displayed by the IL-8 ligand-receptor complex--ligand promiscuity for IL-8, receptor promiscuity for IL8RA, gene duplication for both ligands and receptors and gene extinction in the case of IL8RAP--is unprecedented for the GPCR superfamily.     

Dominant and recessive deafness caused by mutations of a novel gene, TMC1, required for cochlear hair-cell function

Positional cloning of hereditary deafness genes is a direct approach to identify molecules and mechanisms underlying auditory function. Here we report a locus for dominant deafness, DFNA36, which maps to human chromosome 9q13-21 in a region overlapping the DFNB7/B11 locus for recessive deafness. We identified eight mutations in a new gene, transmembrane cochlear-expressed gene 1 (TMC1), in a DFNA36 family and eleven DFNB7/B11 families. We detected a 1.6-kb genomic deletion encompassing exon 14 of Tmc1 in the recessive deafness (dn) mouse mutant, which lacks auditory responses and has hair-cell degeneration. TMC1 and TMC2 on chromosome 20p13 are members of a gene family predicted to encode transmembrane proteins. Tmc1 mRNA is expressed in hair cells of the postnatal mouse cochlea and vestibular end organs and is required for normal function of cochlear hair cells.     

Mutations in the V2 vasopressin receptor gene are associated with X-linked nephrogenic diabetes insipidus.

X-linked nephrogenic diabetes insipidus (NDI) is a rare disorder in which the kidney is insensitive to the antidiuretic hormone, vasopressin. It has been proposed that the kidney-specific V2 vasopressin receptor, a G protein-coupled receptor, is defective in this disorder as both the disease and the receptor map to Xq28. We report six unique mutations in the V2 receptor gene of five unrelated NDI patients, with one patient having two mutations. The most severely affected patient has a nonsense mutation which would terminate the protein in transmembrane domain III. Other mutations include three missense mutations, a frameshift and one small in-frame deletion. These results represent one of the first examples of recessive mutations affecting a G protein-coupled receptor.     

Genome-wide association study identifies five loci associated with susceptibility to pancreatic cancer in Chinese populations

Pancreatic cancer has the lowest survival rate among human cancers, and there are no effective markers for its screening and early diagnosis. To identify genetic susceptibility markers for this cancer, we carried out a genome-wide association study on 981 individuals with pancreatic cancer (cases) and 1,991 cancer-free controls of Chinese descent using 666,141 autosomal SNPs. Promising associations were replicated in an additional 2,603 pancreatic cancer cases and 2,877 controls recruited from 25 hospitals in 16 provinces or cities in China. We identified five new susceptibility loci at chromosomes 21q21.3, 5p13.1, 21q22.3, 22q13.32 and 10q26.11 (P = 2.24 × 10(-13) to P = 4.18 × 10(-10)) in addition to 13q22.1 previously reported in populations of European ancestry. These results advance our understanding of the development of pancreatic cancer and highlight potential targets for the prevention or treatment of this cancer.     

A genome-wide association study identifies two new lung cancer susceptibility loci at 13q12.12 and 22q12.2 in Han Chinese

Lung cancer is the leading cause of cancer-related deaths worldwide. To identify genetic factors that modify the risk of lung cancer in individuals of Chinese ancestry, we performed a genome-wide association scan in 5,408 subjects (2,331 individuals with lung cancer (cases) and 3,077 controls) followed by a two-stage validation among 12,722 subjects (6,313 cases and 6,409 controls). The combined analyses identified six well-replicated SNPs with independent effects and significant lung cancer associations (P < 5.0 × 10(-8)) located in TP63 (rs4488809 at 3q28, P = 7.2 × 10(-26)), TERT-CLPTM1L (rs465498 and rs2736100 at 5p15.33, P = 1.2 × 10(-20) and P = 1.0 × 10(-27), respectively), MIPEP-TNFRSF19 (rs753955 at 13q12.12, P = 1.5 × 10(-12)) and MTMR3-HORMAD2-LIF (rs17728461 and rs36600 at 22q12.2, P = 1.1 × 10(-11) and P = 6.2 × 10(-13), respectively). Two of these loci (13q12.12 and 22q12.2) were newly identified in the Chinese population. These results suggest that genetic variants in 3q28, 5p15.33, 13q12.12 and 22q12.2 may contribute to the susceptibility of lung cancer in Han Chinese.     

Seven prostate cancer susceptibility loci identified by a multi-stage genome-wide association study

Prostate cancer (PrCa) is the most frequently diagnosed male cancer in developed countries. We conducted a multi-stage genome-wide association study for PrCa and previously reported the results of the first two stages, which identified 16 PrCa susceptibility loci. We report here the results of stage 3, in which we evaluated 1,536 SNPs in 4,574 individuals with prostate cancer (cases) and 4,164 controls. We followed up ten new association signals through genotyping in 51,311 samples in 30 studies from the Prostate Cancer Association Group to Investigate Cancer Associated Alterations in the Genome (PRACTICAL) consortium. In addition to replicating previously reported loci, we identified seven new prostate cancer susceptibility loci on chromosomes 2p11, 3q23, 3q26, 5p12, 6p21, 12q13 and Xq12 (P = 4.0 × 10(-8) to P = 2.7 × 10(-24)). We also identified a SNP in TERT more strongly associated with PrCa than that previously reported. More than 40 PrCa susceptibility loci, explaining ～25% of the familial risk in this disease, have now been identified.     

Mutations in the gene encoding pejvakin, a newly identified protein of the afferent auditory pathway, cause DFNB59 auditory neuropathy

Auditory neuropathy is a particular type of hearing impairment in which neural transmission of the auditory signal is impaired, while cochlear outer hair cells remain functional. Here we report on DFNB59, a newly identified gene on chromosome 2q31.1-q31.3 mutated in four families segregating autosomal recessive auditory neuropathy. DFNB59 encodes pejvakin, a 352-residue protein. Pejvakin is a paralog of DFNA5, a protein of unknown function also involved in deafness. By immunohistofluorescence, pejvakin is detected in the cell bodies of neurons of the afferent auditory pathway. Furthermore, Dfnb59 knock-in mice, homozygous for the R183W variant identified in one DFNB59 family, show abnormal auditory brainstem responses indicative of neuronal dysfunction along the auditory pathway. Unlike previously described sensorineural deafness genes, all of which underlie cochlear cell pathologies, DFNB59 is the first human gene implicated in nonsyndromic deafness due to a neuronal defect.     

Functional screening of an asthma QTL in YAC transgenic mice.

Many quantitative trait loci (QTLs) contributing to genetically complex conditions have been discovered, but few causative genes have been identified. This is mainly due to the large size of QTLs and the subtle connection between genotype and quantitative phenotype associated with these conditions. Transgenic mice have been successfully used to analyse well-characterized genes suspected of contributing to quantitative traits. Although this approach is powerful for examining one gene at a time, it can be impractical for surveying the large genomic intervals containing many genes that are typically associated with QTLs. To screen for genes contributing to an asthma QTL mapped to human chromosome 5q3 (refs 6,7), we characterized a panel of large-insert 5q31 transgenics based on studies demonstrating that altering gene dosage frequently affects quantitative phenotypes normally influenced by that gene. This panel of human YAC transgenics, propagating a 1-Mb interval of chromosome 5q31 containing 6 cytokine genes and 17 partially characterized genes, was screened for quantitative changes in several asthma-associated phenotypes. Multiple independent transgenic lines with altered IgE response to antigen treatment shared a 180-kb region containing 5 genes, including those encoding human interleukin 4 (IL4) and interleukin 13 (IL13 ), which induce IgE class switching in B cells. Further analysis of these mice and mice transgenic for mouse Il4 and Il13 demonstrated that moderate changes in Il4 and Il13 expression affect asthma-associated phenotypes in vivo. This functional screen of large-insert transgenics enabled us to identify genes that influence the QTL phenotype in vivo.     

Tas1r3, encoding a new candidate taste receptor, is allelic to the sweet responsiveness locus Sac

The ability to taste the sweetness of carbohydrate-rich foodstuffs has a critical role in the nutritional status of humans. Although several components of bitter transduction pathways have been identified, the receptors and other sweet transduction elements remain unknown. The Sac locus in mouse, mapped to the distal end of chromosome 4 (refs. 7-9), is the major determinant of differences between sweet-sensitive and -insensitive strains of mice in their responsiveness to saccharin, sucrose and other sweeteners. To identify the human Sac locus, we searched for candidate genes within a region of approximately one million base pairs of the sequenced human genome syntenous to the region of Sac in mouse. From this search, we identified a likely candidate: T1R3, a previously unknown G protein-coupled receptor (GPCR) and the only GPCR in this region. Mouse Tas1r3 (encoding T1r3) maps to within 20,000 bp of the marker closest to Sac (ref. 9) and, like human TAS1R3, is expressed selectively in taste receptor cells. By comparing the sequence of Tas1r3 from several independently derived strains of mice, we identified a specific polymorphism that assorts between taster and non-taster strains. According to models of its structure, T1r3 from non-tasters is predicted to have an extra amino-terminal glycosylation site that, if used, would interfere with dimerization.     

A genome-wide association study identifies two new risk loci for Graves' disease

Graves' disease is a common autoimmune disorder characterized by thyroid stimulating hormone receptor autoantibodies (TRAb) and hyperthyroidism. To investigate the genetic architecture of Graves' disease, we conducted a genome-wide association study in 1,536 individuals with Graves' disease (cases) and 1,516 controls. We further evaluated a group of associated SNPs in a second set of 3,994 cases and 3,510 controls. We confirmed four previously reported loci (in the major histocompatibility complex, TSHR, CTLA4 and FCRL3) and identified two new susceptibility loci (the RNASET2-FGFR1OP-CCR6 region at 6q27 (P(combined) = 6.85 × 10(-10) for rs9355610) and an intergenic region at 4p14 (P(combined) = 1.08 × 10(-13) for rs6832151)). These newly associated SNPs were correlated with the expression levels of RNASET2 at 6q27, of CHRNA9 and of a previously uncharacterized gene at 4p14, respectively. Moreover, we identified strong associations of TSHR and major histocompatibility complex class II variants with persistently TRAb-positive Graves' disease.     

Rheumatoid arthritis susceptibility loci at chromosomes 10p15, 12q13 and 22q13

The WTCCC study identified 49 SNPs putatively associated with rheumatoid arthritis at P = 1 x 10(-4) - 1 x 10(-5) (tier 3). Here we show that three of these SNPs, mapping to chromosome 10p15 (rs4750316), 12q13 (rs1678542) and 22q13 (rs3218253), are also associated (trend P = 4 x 10(-5), P = 4 x 10(-4) and P = 4 x 10(-4), respectively) in a validation study of 4,106 individuals with rheumatoid arthritis and an expanded reference group of 11,238 subjects, confirming them as true susceptibility loci in individuals of European ancestry.     

Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture

Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10(-8)). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10(-4), Bonferroni corrected), of which six reached P < 5 × 10(-8), including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.     

The transmembrane protein meckelin (MKS3) is mutated in Meckel-Gruber syndrome and the wpk rat

Meckel-Gruber syndrome is a severe autosomal, recessively inherited disorder characterized by bilateral renal cystic dysplasia, developmental defects of the central nervous system (most commonly occipital encephalocele), hepatic ductal dysplasia and cysts and polydactyly. MKS is genetically heterogeneous, with three loci mapped: MKS1, 17q21-24 (ref. 4); MKS2, 11q13 (ref. 5) and MKS3 (ref. 6). We have refined MKS3 mapping to a 12.67-Mb interval (8q21.13-q22.1) that is syntenic to the Wpk locus in rat, which is a model with polycystic kidney disease, agenesis of the corpus callosum and hydrocephalus. Positional cloning of the Wpk gene suggested a MKS3 candidate gene, TMEM67, for which we identified pathogenic mutations for five MKS3-linked consanguineous families. MKS3 is a previously uncharacterized, evolutionarily conserved gene that is expressed at moderate levels in fetal brain, liver and kidney but has widespread, low levels of expression. It encodes a 995-amino acid seven-transmembrane receptor protein of unknown function that we have called meckelin.     

A genome-wide association study identifies colorectal cancer susceptibility loci on chromosomes 10p14 and 8q23.3

To identify colorectal cancer (CRC) susceptibility alleles, we conducted a genome-wide association study. In phase 1, we genotyped 550,163 tagSNPs in 940 familial colorectal tumor cases (627 CRC, 313 high-risk adenoma) and 965 controls. In phase 2, we genotyped 42,708 selected SNPs in 2,873 CRC cases and 2,871 controls. In phase 3, we evaluated 11 SNPs showing association at P < 10(-4) in a joint analysis of phases 1 and 2 in 4,287 CRC cases and 3,743 controls. Two SNPs were taken forward to phase 4 genotyping (10,731 CRC cases and 10,961 controls from eight centers). In addition to the previously reported 8q24, 15q13 and 18q21 CRC risk loci, we identified two previously unreported associations: rs10795668, located at 10p14 (P = 2.5 x 10(-13) overall; P = 6.9 x 10(-12) replication), and rs16892766, at 8q23.3 (P = 3.3 x 10(-18) overall; P = 9.6 x 10(-17) replication), which tags a plausible causative gene, EIF3H. These data provide further evidence for the 'common-disease common-variant' model of CRC predisposition.     

Mutant frizzled-4 disrupts retinal angiogenesis in familial exudative vitreoretinopathy

Familial exudative vitreoretinopathy (FEVR) is a hereditary ocular disorder characterized by a failure of peripheral retinal vascularization. Loci associated with FEVR map to 11q13-q23 (EVR1; OMIM 133780, ref. 1), Xp11.4 (EVR2; OMIM 305390, ref. 2) and 11p13-12 (EVR3; OMIM 605750, ref. 3). Here we have confirmed linkage to the 11q13-23 locus for autosomal dominant FEVR in one large multigenerational family and refined the disease locus to a genomic region spanning 1.55 Mb. Mutations in FZD4, encoding the putative Wnt receptor frizzled-4, segregated completely with affected individuals in the family and were detected in affected individuals from an additional unrelated family, but not in normal controls. FZD genes encode Wnt receptors, which are implicated in development and carcinogenesis. Injection of wildtype and mutated FZD4 into Xenopus laevis embryos revealed that wildtype, but not mutant, frizzled-4 activated calcium/calmodulin-dependent protein kinase II (CAMKII) and protein kinase C (PKC), components of the Wnt/Ca(2+) signaling pathway. In one of the mutants, altered subcellular trafficking led to defective signaling. These findings support a function for frizzled-4 in retinal angiogenesis and establish the first association between a Wnt receptor and human disease.     

Variants of ENPP1 are associated with childhood and adult obesity and increase the risk of glucose intolerance and type 2 diabetes

We identified a locus on chromosome 6q16.3-q24.2 (ref. 1) associated with childhood obesity that includes 2.4 Mb common to eight genome scans for type 2 diabetes (T2D) or obesity. Analysis of the gene ENPP1 (also called PC-1), a candidate for insulin resistance, in 6,147 subjects showed association between a three-allele risk haplotype (K121Q, IVS20delT-11 and A-->G+1044TGA; QdelTG) and childhood obesity (odds ratio (OR) = 1.69, P = 0.0006), morbid or moderate obesity in adults (OR = 1.50, P = 0.006 or OR = 1.37, P = 0.02, respectively) and T2D (OR = 1.56, P = 0.00002). The Genotype IBD Sharing Test suggested that this obesity-associated ENPP1 risk haplotype contributes to the observed chromosome 6q linkage with childhood obesity. The haplotype confers a higher risk of glucose intolerance and T2D to obese children and their parents and associates with increased serum levels of soluble ENPP1 protein in children. Expression of a long ENPP1 mRNA isoform, which includes the obesity-associated A-->G+1044TGA SNP, was specific for pancreatic islet beta cells, adipocytes and liver. These findings suggest that several variants of ENPP1 have a primary role in mediating insulin resistance and in the development of both obesity and T2D, suggesting that an underlying molecular mechanism is common to both conditions.     

A genome-wide association study of nonsynonymous SNPs identifies a type 1 diabetes locus in the interferon-induced helicase (IFIH1) region

In this study we report convincing statistical support for a sixth type 1 diabetes (T1D) locus in the innate immunity viral RNA receptor gene region IFIH1 (also known as mda-5 or Helicard) on chromosome 2q24.3. We found the association in an interim analysis of a genome-wide nonsynonymous SNP (nsSNP) scan, and we validated it in a case-control collection and replicated it in an independent family collection. In 4,253 cases, 5,842 controls and 2,134 parent-child trio genotypes, the risk ratio for the minor allele of the nsSNP rs1990760 A --> G (A946T) was 0.86 (95% confidence interval = 0.82-0.90) at P = 1.42 x 10(-10).