Relevance of the deletion polymorphisms of the glutathione S-transferases GSTT1 and GSTM1 in pharmacology and toxicology

Although cytosolic glutathione S-transferase (GST) enzymes occupy a key position in biological detoxification processes, two of the most relevant human isoenzymes, GSTT1-1 and GSTM1-1, are genetically deleted (non-functional alleles GSTT1*0 and GSTM1*0) in a high percentage of the human population, with major ethnic differences. The structures of the GSTT and GSTM gene areas explain the underlying genetic processes. GSTT1-1 is highly conserved during evolution and plays a major role in phase-II biotransformation of a number of drugs and industrial chemicals, e.g. cytostatic drugs, hydrocarbons and halogenated hydrocarbons. GSTM1-1 is particularly relevant in the deactivation of carcinogenic intermediates of polycyclic aromatic hydrocarbons. Several lines of evidence suggest that hGSTT1-1 and/or hGSTM1-1 play a role in the deactivation of reactive oxygen species that are likely to be involved in cellular processes of inflammation, ageing and degenerative diseases. There is cumulating evidence that combinations of the GSTM1*0 state with other genetic traits affecting the metabolism of carcinogens (CYP1A1, GSTP1) may predispose the aero-digestive tract and lung, especially in smokers, to a higher risk of cancer. The GSTM1*0 status appears also associated with a modest increase in the risk of bladder cancer, consistent with a GSTM1 interaction with carcinogenic tobacco smoke constituents. Both human GST deletions, although largely counterbalanced by overlapping substrate affinities within the GST superfamily, have consequences when the organism comes into contact with distinct man-made chemicals. This appears relevant in industrial toxicology and in drug metabolism.

Cystine-mediated oxidative defence in Lactobacillus reuteri BR11

Lactobacillus reuteri BR11 possesses an abundant cystine uptake (Cyu) ABC-transporter that was previously found to be involved in a novel mechanism of oxidative defence mediated by cystine. The current study aimed to elucidate this mechanism with a focus on the role of the co-transcribed cystathionine ã-lyase (Cgl). Growth studies of wild-type L. reuteri BR11 and mutants inactivated in cgl and the cystine-binding protein encoding gene cyuC showed that in contrast to the Cyu transporter, whose inactivation led to growth arrest in aerated cultures, Cgl is not crucial for oxidative defence. However, the role of Cgl in oxidative defence became apparent in the presence of severe oxidative damage and cysteine deprivation. Cysteine was found to be protective against oxidative stress, and the action of Cgl in both cysteine biosynthesis and degradation poses a seemingly futile pathway that deprives the intracellular cysteine pool. To further characterise the relationship between Cgl activity and cysteine and their roles in oxidative defence, enzymatic assays were performed on purified Cgl, and intracellular concentrations of cysteine, cystathionine and methionine were determined. Cgl was highly active towards cystine and cystathionine and less active towards cysteine in vitro, suggesting the main function of Cgl to be cysteine biosynthesis. Cysteine was found at high concentrations in the cell, but the levels were not significantly affected by inactivation of cgl or growth under aerobic conditions. It was concluded that both anabolic and catabolic activities of Cgl towards cysteine contribute to oxidative defence, the former by maintaining an intracellular reservoir of thiol analogous to glutathione, and the latter by producing H2S which is readily secreted, thus creating a reducing extracellular environment. The significance of the Cyu transporter to the physiology of L. reuteri BR11 prompted a phylogenetic study to determine its presence in bacteria. Orthologs of the Cyu transporter that are closest matches to the Cyu transporter are only limited to several species of Lactobacillus and Leuconostoc. Outside the Lactobacillales order, the closest matching orthologs belong to Proteobacteria, and there are more orthologs in Proteobacteria than non-Lactobacillales Firmicutes, suggesting that the Cyu transporter locus was present in the ancestor of the Proteobacteria and Firmicutes, and over evolutionary time has been lost or diverged in many Firmicutes. The clustering of the Cyu transporter locus with a gene encoding a Cgl family protein is even rarer. It was only found in L. reuteri, Lactobacillus vaginalis, Weissella paramesenteroides, the Lactobacillus casei group, and several Campylobacter sp. An accompanying phylogenetic study of L. reuteri BR11 using multi-locus sequence analysis showed that L. reuteri BR11 had diverged from more than 100 strains of L. reuteri isolated from various hosts and geographical locations. However, comparison with other Lactobacillus species supported the current classification of BR11 as L. reuteri. The most closely related species to L. reuteri is L. vaginalis or Lactobacillus antri, depending on the housekeeping gene used for analysis. The close evolutionary relationship of L. vaginalis to L. reuteri and the high degree of sequence identity between the cgl-cyuABC loci in both species suggest that the Cyu system is highly likely to perform similar functions in L. vaginalis. In search of other genes that function in oxidative defence, a number of mutants which were inactivated in genes that confer increased resistance to oxidative stress in other bacteria were constructed. The genes targeted were ahpC (peroxidase component of the alkyl hydroperoxide reductase system), tpx (thiol peroxidase), osmC (osmotically induced protein C), mntH (Mn2+/Fe2+ transporter), gshA (ã-glutamylcysteine synthetase) and msrA (methionine sulfoxide reductase). The ahpC and mntH mutants had slightly lower minimum inhibitory concentrations of organic peroxides, suggesting these genes might be involved in resistance to organic peroxides in L. reuteri. However, none of the mutants exhibited growth defects in aerated cultures, in stark contrast to the cyuC mutant. This may be due to compensatory functions of other genes, a hypothesis which cannot be tested until a robust protocol for constructing markerless multiple gene deletion mutants in L. reuteri is developed. These results highlight the importance of the Cyu transporter in oxidative defence and provide a foundation for extending the research of this system in other bacteria.

Evidence for u.v. induction of CDKN2 mutations in melanoma cell lines

The CDKN2 gene, encoding the cyclin dependent kinase inhibitor p16, is a tumour suppressor gene involved in melanoma and maps to chromosome band 9p22. Mutations or interstitial deletions of this gene have been found both in the germline of familial melanoma cases and somatically in melanoma cell lines. Previous mutation analyses of melanoma cell lines have indicated a high frequency of C:G to T:A transitions, with all of these mutations occurring at dipyrimidine sites. Including three melanoma cell lines carrying tandem CC to TT mutations, the spectrum of mutations so far reported indicates a possible role for u.v. radiation in the mutagenesis of this gene in some tumours. To further examine this hypothesis we have characterised mutations of the CDKN2 gene in 30 melanoma cell lines. Nineteen lines carried complete or partial homozygous deletions of the gene. Of the remaining cell lines, eight were shown by direct sequencing of PCR products from exon 1 and exon 2 to carry a total of nine different mutations of CDKN2. Two cell lines carried tandem CC to TT mutations and a high rate of C:G to T:A transitions was observed. This study provides further evidence for the role of u.v. light in the genesis of melanoma, with one target being the CDKN2 tumour suppressor gene.

CDKN2A/p16 is inactivated in most melanoma cell lines

The CDKN2A gene maps to chromosome 9p21-22 and is responsible for melanoma susceptibility in some families. Its product, p16, binds specifically to CDK4 and CDK6 in vitro and in vivo, inhibiting their kinase activity. CDKN2A is homozygously deleted or mutated in a large proportion of tumor cell lines and some primary tumors, including melanomas. The aim of this study was to investigate the involvement of CDKN2A and elucidate the mechanisms of p16 inactivation in a panel of 60 cell lines derived from sporadic melanomas. Twenty-six (43%) of the melanoma lines were homozygously deleted for CDKN2A, and an additional 15 (25%) lines carried missense, nonsense, or frameshift mutations. All but one of the latter group were shown by microsatellite analysis to be hemizygous for the region of 9p surrounding CDKN2A. p16 was detected by Western blotting in only five of the cell lines carrying mutations. Immunoprecipitation of p16 in these lines, followed by Western blotting to detect the coprecipitation of CDK4 and CDK6, revealed that p16 was functionally compromised in all cell lines but the one that carried a heterozygous CDKN2A mutation. In the remaining 19 lines that carried wild-type CDKN2A alleles, Western blot analysis and immunoprecipitation indicated that 11 cell lines expressed a wild-type protein. Northern blotting was performed on the remaining eight cell lines and revealed that one cell line carried an aberrantly sized RNA transcript, and two other cell lines failed to express RNA. The promoter was found to be methylated in five cell lines that expressed CDKN2A transcript but not p16. Presumably, the message seen by Northern blotting in these cell lines is the result of cross-hybridization of the total cDNA probe with the exon 1beta transcript. Microsatellite analysis revealed that the majority of these cell lines were hemi/homozygous for the region surrounding CDKN2A, indicating that the wild-type allele had been lost. In the 11 cell lines that expressed functional p16, microsatellite analysis revealed loss of heterozygosity at the markers immediately surrounding CDKN2A in five cases, and the previously characterized R24C mutation of CDK4 was identified in one of the remaining 6 lines. These data indicate that 55 of 60 (92%) melanoma cell lines demonstrated some aberration of CDKN2A or CDK4, thus suggesting that this pathway is a primary genetic target in melanoma development.

CDKN2A is not the principal target of deletions on the short arm of chromosome 9 in neuroendocrine (Merkel cell) carcinoma of the skin

The majority of small-cell lung cancers (SCLCs) express p16 but not pRb. Given our previous study showing loss of pRb in Merkel cell carcinoma (MCC)/neuroendocrine carcinoma of the skin and the clinicopathological similarities between SCLC and MCC, we wished to determine if this was also the case in MCC. Twenty-nine MCC specimens from 23 patients were examined for deletions at 10 loci on 9p and 1 on 9q. No loss of heterozygosity (LOH) was seen in 9 patients including 2 for which tumour and cell line DNAs were examined. Four patients had LOH for all informative loci on 9p. Ten tumours showed more limited regions of loss on 9p, and from these 2 common regions of deletion were determined. Half of all informative cases had LOH at D9S168, the most telomeric marker examined, and 3 specimens showed loss of only D9S168. A second region (IFNA-D9S126) showed LOH in 10 (44%) cases, and case MCC26 showed LOH for only D9S126, implicating genes centromeric of the CDKN2A locus. No mutations in the coding regions of p16 were seen in 7 cell lines tested, and reactivity to anti-p16 antibody was seen in all 11 tumour specimens examined and in 6 of 7 cell lines from 6 patients. Furthermore, all cell lines examined reacted with anti-p14(ARF) antibody. These results suggest that neither transcript of the CDKN2A locus is the target of deletions on 9p in MCC and imply the existence of tumour-suppressor genes mapping both centromeric and telomeric of this locus.

PTEN inactivation is rare in melanoma tumours but occurs frequently in melanoma cell lines

Deletions detected in cytogenetic and loss of heterozygosity (LOH) studies indicate that at least one tumour suppressor gene maps to the long arm of chromosome 10. Previous deletion mapping studies have observed LOH on 10q in about 30% of melanomas analysed. The PTEN gene, mapping to chromosome band 10q23.3, encodes a protein with both lipid and protein phosphatase activity. Somatic mutations and deletions in have been detected in a variety of cell lines and tumours, including melanoma samples. We performed mutation analyses and extensive allelic loss studies to investigate the role this gene plays in melanoma pathogenesis. We found that a total of 34 out of 57 (60%) melanoma cell lines carried hemizygous deletions of chromosome 10q encompassing the PTEN locus. A further three cell lines carried smaller deletions excluding PTEN. Inactivation of both PTEN alleles by exon-specific homozygous deletion or mutation was observed in 13 out of 57 (23%) melanoma cell lines. The mutation spectrum observed does not indicate an important role for ultraviolet radiation in the genesis of these mutations, and evidence from three cell lines supports the acquisition of PTEN aberrations in culture. Ten out of 49 (20%) matched melanoma tumour/normal samples harboured hemizygous deletions of either the whole chromosome or most of the long arm. Mutations within were detected in only one of the 10 tumours demonstrating LOH at 10q23 that were analysed. These results suggest that PTEN inactivation may be important for the propagation of melanoma cells in culture, and that another chromosome 10 tumour suppressor gene may be important for melanoma pathogenesis.

The null allele of GSTM1 does not affect susceptibility to solar keratoses in the Australian white population

Solar keratoses (SKs) are induced by exposure to UV radiation and are capable of undergoing transformation to squamous cell carcinoma (SCC).1 The two main factors influencing the occurrence of SK are the sensitivity of the skin to sunlight and the total duration of solar exposure. These factors are responsible for the high incidence of SK in Australia. Although the influence of genetic factors is not defined, there is evidence that the gene encoding the enzyme, glutathione S-transferase, may be implicated in cancer predisposition and therefore SK. Glutathione S-transferase Mu-1 (GSTM1) is an isoenzyme involved in the detoxification of carcinogens. The GSTM1 protein is completely absent in approximately 50% of white persons. This absence is caused by a homozygous gene deletion on chromosome 1p resulting in a null genotype.2 Katoh3 showed that the frequency of the GSTM1 null genotype was significantly higher in 85 patients with urothelial cancer (61.2%; p < 0.05), suggesting that the null genotype may increase cancer susceptibility. This finding was supported by Lafuente et al.4 who found evidence that persons who lack the GSTM1 gene have approximately twice the chance of experiencing malignant melanoma. Further research in the United Kingdom found that patients with two or more skin tumors of different types, basal cell carcinoma (BCC) and SCC, had a significantly higher frequency of GSTM1 null genotypes than controls (71%; p = 0.033). However the GSTM1 genotype in patients with only SCC was not excessive in this population.5 Persons residing in northern Australia have the highest incidence of nonmelanoma skin cancer (SCC and BCC) in the world6 and receive far greater solar exposure than persons residing in the United Kingdom. It is possible that the GSTM1 null genotype may affect susceptibility to SK, which may act as SCC precursors, in Australians exposed to these high levels of solar radiation.

Cyclooxygenase-2-linked attenuation of hypoxia-induced pulmonary hypertension and intravascular thrombosis

Exogenous prostacyclin is effective in reducing pulmonary vascular resistance in some forms of human pulmonary hypertension (PH). To explore whether endogenous prostaglandins played a similar role in pulmonary hypertension, we examined the effect of deleting cyclooxygenase (COX)-gene isoforms in a chronic hypoxia model of PH. Pulmonary hypertension, examined by direct measurement of right ventricular end systolic pressure (RVESP), right ventricular hypertrophy (n = 8), and hematocrit (n = 3), was induced by 3 weeks of hypobarichypoxia in wild-type and COX-knockout (KO) mice. RVESP was increased in wild-type hypoxic mice compared with normoxic controls (24.4 ± 1.4 versus 13.8 ± 1.9 mm Hg; n = 8; p < 0.05). COX-2 KO mice showed a greater increase in RVESP following hypoxia (36.8 ± 2.7 mm Hg; p < 0.05). Urinary thromboxane (TX)B2 excretion increased following hypoxia (44.6 ± 11.1 versus 14.7 ± 1.8 ng/ml; n = 6; p < 0.05), an effect that was exacerbated by COX-2 gene disruption (54.5 ± 10.8 ng/ml; n = 6). In contrast, the increase in 6-keto-prostacyclin1α excretion following hypoxia was reduced by COX-2 gene disruption (29 ± 3 versus 52 ± 4.6 ng/ml; p < 0.01). Tail cut bleed times were lower following hypoxia, and there was evidence of intravascular thrombosis in lung vessels that was exacerbated by disruption of COX-2 and reduced by deletion of COX-1. The TXA2/endoperoxide receptor antagonist ifetroban (50 mg/kg/day) offset the effect of deleting the COX-2 gene, attenuating the hypoxia-induced rise in RVESP and intravascular thrombosis. COX-2 gene deletion exacerbates pulmonary hypertension, enhances sensitivity to TXA2, and induces intravascular thrombosis in response to hypoxia. The data provide evidence that endogenous prostaglandins modulate the pulmonary response to hypoxia. Copyright © 2008 by The American Society for Pharmacology and Experimental Therapeutics.

The role of the molecular footprint of EGFR in tailoring treatment decisions in NSCLC

The majority of patients with non-small-cell lung cancer (NSCLC) present with advanced disease, with targeted therapies providing some improvement in clinical outcomes. The epidermal growth factor receptor (EGFR) tyrosine kinase (TK) plays an important role in the pathogenesis of NSCLC. Tyrosine kinase inhibitors (TKIs), which target the EGFR TK domain, have proven to be an effective treatment strategy; however, patient responses to treatment vary considerably. Therefore, the identification of patients most likely to respond to treatment is essential to optimise the benefit of TKIs. Tumour-associated activating mutations in EGFR can identify patients with NSCLC who are likely to have a good response to TKIs. Nonetheless, the majority of patients relapse within a year of starting treatment. Studies of tumours at relapse have demonstrated expression of a T790M mutation in exon 20 of the EGFR TK domain in approximately 50% of cases. Although conferring resistance to reversible TKIs, these patients may remain sensitive to new-generation irreversible/panerb inhibitors. A number of techniques have been employed for genotypic assessment of tumourassociated DNA to identify EGFR mutations, each of which has advantages and disadvantages. This review presents an overview of the current methodologies used to identify such molecular markers. Recent developments in technology may make the monitoring of changes in patients' tumour genotypes easier in clinical practice, which may enable patients' treatment regimens to be tailored during the course of their disease, potentially leading to improved patient outcomes.

Ran is a potential therapeutic target for cancer cells with molecular changes associated with activation of the PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways

Purpose Cancer cells have been shown to be more susceptible to Ran knockdown than normal cells. We now investigate whether Ran is a potential therapeutic target of cancers with frequently found mutations that lead to higher Ras/MEK/ERK [mitogen-activated protein/extracellular signal-regulated kinase (ERK; MEK)] and phosphoinositide 3-kinase (PI3K)/Akt/mTORC1 activities. Experimental Design Apoptosis was measured by flow cytometry [propidium iodide (PI) and Annexin V staining] and MTT assay in cancer cells grown under different conditions after knockdown of Ran. The correlations between Ran expression and patient survival were examined in breast and lung cancers. Results Cancer cells with their PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways inhibited are less susceptible to Ran silencing-induced apoptosis. K-Ras-mutated, c-Met-amplified, and Pten-deleted cancer cells are also more susceptible to Ran silencing-induced apoptosis than their wild-type counterparts and this effect is reduced by inhibitors of the PI3K/Akt/mTORC1 and MEK/ERK pathways. Overexpression of Ran in clinical specimens is significantly associated with poor patient outcome in both breast and lung cancers. This association is dramatically enhanced in cancers with increased c-Met or osteopontin expression, or with oncogenic mutations of K-Ras or PIK3CA, all of which are mutations that potentially correlate with activation of the PI3K/Akt/mTORC1 and/or Ras/MEK/ERK pathways. Silencing Ran also results in dysregulation of nucleocytoplasmic transport of transcription factors and downregulation of Mcl-1 expression, at the transcriptional level, which are reversed by inhibitors of the PI3K/Akt/mTORC1 and MEK/ERK pathways. Conclusion Ran is a potential therapeutic target for treatment of cancers with mutations/changes of expression in protooncogenes that lead to activation of the PI3K/Akt/mTORC1 and Ras/MEK/ERK pathways. ©2011 AACR.

De Novo Transcriptome Sequence Assembly and Analysis of RNA Silencing Genes of Nicotiana benthamiana

Background: Nicotiana benthamiana has been widely used for transient gene expression assays and as a model plant in the study of plant-microbe interactions, lipid engineering and RNA silencing pathways. Assembling the sequence of its transcriptome provides information that, in conjunction with the genome sequence, will facilitate gaining insight into the plant's capacity for high-level transient transgene expression, generation of mobile gene silencing signals, and hyper-susceptibility to viral infection. Methodology/Results: RNA-seq libraries from 9 different tissues were deep sequenced and assembled, de novo, into a representation of the transcriptome. The assembly, of16GB of sequence, yielded 237,340 contigs, clustering into 119,014 transcripts (unigenes). Between 80 and 85% of reads from all tissues could be mapped back to the full transcriptome. Approximately 63% of the unigenes exhibited a match to the Solgenomics tomato predicted proteins database. Approximately 94% of the Solgenomics N. benthamiana unigene set (16,024 sequences) matched our unigene set (119,014 sequences). Using homology searches we identified 31 homologues that are involved in RNAi-associated pathways in Arabidopsis thaliana, and show that they possess the domains characteristic of these proteins. Of these genes, the RNA dependent RNA polymerase gene, Rdr1, is transcribed but has a 72 nt insertion in exon1 that would cause premature termination of translation. Dicer-like 3 (DCL3) appears to lack both the DEAD helicase motif and second dsRNA binding motif, and DCL2 and AGO4b have unexpectedly high levels of transcription. Conclusions: The assembled and annotated representation of the transcriptome and list of RNAi-associated sequences are accessible at www.benthgenome.com alongside a draft genome assembly. These genomic resources will be very useful for further study of the developmental, metabolic and defense pathways of N. benthamiana and in understanding the mechanisms behind the features which have made it such a well-used model plant. © 2013 Nakasugi et al.

Conjugation of carcinogens by 6 class glutathione S-transferases : mechanisms and relevance to variations in human risk

Conjugation of chemicals with glutathione (GSH) can lead to decreased or increased toxicity. A genetic deficiency in the GSH S-transferase μ class gene M1 has been hypothesized to lead to greater risk of lung cancer in smokers. Recently a gene deletion polymorphism involving the human θ enzyme T1 has been described; the enzyme is present in erythrocytes and can be readily assayed. A rat θ class enzyme, 5-5, has structural and catalytic similarity and the protein was expressed in the Salmonella typhimurium tester strain TA1535. Expression of the cDNA vector increased the mutagenicity of ethylene dibromide and several methylene dihalides. Mutations resulting from the known GSH S-transferase substrate 1,2-epoxy-3-(4′nitrophenoxy)propane were decreased in the presence of the transferase. Expression of transferase 5-5 increased mutations when 1,2,3,4-diepoxybutane (butadiene diepoxide), 4-bromo-1,2-epoxybutane, or 1,3-dichloracetone were added. The latter compound is a model for the putative 1,2-dibromo-3-chloropropane oxidation product 1-bromo-3-chloroacetone. These genotoxicity and genotyping assays may be of use in further studies of the roles of GSH S-transferase θ enzymes in bioactivation and detoxication and any changes in risk due to polymorphism.

Species differences in the glutathione transferase GSTT1-1 activity towards the model substrates methyl chloride and dichloromethane in liver and kidney

Glutathione transferase (GST) GSTT1-1 is involved in the biotransformation of several chemicals widely used in industry, such as butadiene and dichloro methane DCM. The polymorphic hGSTT1-1 may well play a role in the development of kidney tumours after high and long-term occupational exposure against trichloroethylene. Although several studies have investigated the association of this polymorphism with malignant diseases little is known about its enzyme activity in potential extrahepatic target tissues. The known theta-specific substrates methyl chloride (MC) dichloromethane and 1,2-epoxy-3-(p-nitrophenoxy)propane (EPNP) were used to assay GSTT1-1 activity in liver and kidney of rats, mice, hamsters and humans differentiating the three phenotypes (non-conjugators, low conjugators, high conjugators) seen in humans. In addition GSTT1-1 activity towards MC and DCM was determined in human erythrocytes. No GSTT1-1 activity was found in any tissue of non-conjugators (NC). In all organs high conjugators (HC) showed twofold higher activity towards MC and DCM than low conjugators (LC). The activity in human samples towards EPNP was too close to the detection limit to differentiate between the three conjugator phenotypes. GSTT1-1 activity towards MC was two to seven-times higher in liver cytosol than in kidney cytosol. The relation for MC between species was identical in both organs: mouse > HC > rat > LC > hamster > NC. In rats, mice and hamsters GSTT1-1 activity in liver cytosol towards DCM was also two to seven-times higher than in the kidney cytosol. In humans this activity was twice as high in kidney cytosol than in liver cytosol. The relation between species was mouse > rat > HC > LC > hamster > NC for liver, but mouse > HC > LC/rat > hamster/NC for kidney cytosol. The importance to heed the specific environment at potential target sites in risk assessment is emphasized by these results.

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<P<5.9×10-4) at both femur and tibia, compared with their wild-type littermates. Natural variation in humans and targeted disruption in mice demonstrate that WNT16 is an important determinant of CBT, BMD, bone strength, and risk of fracture. © 2012 Zheng et al.