The plasma level of the tumor necrosis factor-alpha receptor 2 (TNFR2) is associated with obesity phenotypes. However, the genetic polymorphisms for such an association have rarely been explored and are generally unknown. In this study, by employing a large sample of 1,873 subjects from 405 Caucasian nuclear families, we explored the association of 12 SNPs of the TNFR2 gene and obesity-related phenotypes, including body mass index (BMI), fat mass, and percentage fat mass (PFM). The within-family quantitative transmission disequilibrium test, which is robust to sample stratification, was implemented to evaluate the association of TNFR2 gene with obesity phenotypes. Evidence of association was obtained at SNP9 (rs5746059) with fat mass (P = 0.0002), BMI (P = 0.002), and PFM (P = 0.0006). The contribution of this polymorphism to the variation of fat mass and PFM was 6.24 and 7.82%, respectively. Individuals carrying allele A at the SNP9 site had a 4.6% higher fat mass and a 2.5% increased PFM compared to noncarriers. The results remained significant even after correction for multiple testing. Evidence of association between the TNFR2 gene and obesity phenotypes are also found in 700 independent Chinese Han and 1,000 random Caucasians samples. The results suggest that the TNFR2 gene polymorphisms contribute to the variation of obesity phenotypes.
Alcohol dependence (AD) is a complex disorder characterized by impaired control over drinking. It is determined by both genetic and environmental factors. The recent approach of genome-wide association study (GWAS) is a powerful tool for identifying complex disease-associated susceptibility alleles, however, a few GWASs have been conducted for AD, and their results are largely inconsistent. The present study aimed to screen the loci associated with alcohol-related phenotypes using GWAS technology.
A genome-wide association study with the behavior of regular alcohol drinking and alcohol consumption was performed to identify susceptibility genes associated with AD, using the Affymetrix 500K SNP array in an initial sample consisting of 904 unrelated Caucasian subjects. Then, the initial results in GWAS were replicated in three independent samples: 1972 Caucasians in 593 nuclear families, 761 unrelated Caucasian subjects, and 2955 unrelated Chinese Hans.
Several genes were associated with the alcohol-related phenotypes at the genome-wide significance level, with the ankyrin repeat domain 7 gene (ANKRD7) showing the strongest statistical evidence for regular alcohol drinking and suggestive statistical evidence for alcohol consumption. In addition, certain haplotypes within the ANKRD7 and cytokine-like1 (CYTL1) genes were significantly associated with regular drinking behavior, such as one ANKRD7 block composed of the SNPs rs6466686-rs4295599-rs12531086 (P = 6.51×10–8). The association of alcohol consumption was successfully replicated with rs4295599 in ANKRD7 gene in independent Caucasian nuclear families and independent unrelated Chinese Hans, and with rs16836497 in CYTL1 gene in independent unrelated Caucasians. Meta-analyses based on both the GWAS and replication samples further supported the observed significant associations between the ANKRD7 or CYTL1 gene and alcohol consumption.
The evidence suggests that ANKRD7 and CYTL1 genes may play an important role in the variance in AD risk.
alcohol dependence; ANKRD7; CYTL1; genome-wide association study
Lean body mass (LBM) is a heritable trait predicting a series of health problems, such as osteoporotic fracture and sarcopenia. We aim to identify sequence variants associated with LBM by a genome-wide association study (GWAS) of copy number variants (CNVs). We genotyped genome-wide CNVs of 1627 individuals of the Chinese population with Affymetrix SNP6.0 genotyping platform, which comprised of 9 40 000 copy number probes. We then performed a GWAS of CNVs with lean mass at seven sites: left and right arms, left and right legs, total of limb, trunk and whole body. We identified a CNV that is associated with LBM variation at the genome-wide significance level (CNV2073, Bonferroni corrected P-value 0.002 at right arm). CNV2073 locates at chromosome 15q13.3, which has been implicated as a candidate region for LBM by our previous linkage studies. The nearest gene, gremlin1, has a key role in the regulation of skeletal muscle formation and repair. Our results suggest that the gremlin1 gene is a potentially important gene for LBM variation. Our findings also show the utility and efficacy of CNV as genetic markers in association studies.
association; copy number variation; gremlin1 gene; lean body mass; 15q13.3
Prostate cancer, one of the most lethal forms of urinary system cancer, remains resistant to currently available treatments. Therefore, novel mechanism and target-based approaches are needed for the management of this neoplasm. PI3K/AKT signaling pathway activation correlates with human prostate cancer progression and metastasis. However, the role of mTOR in prostate cancer is not well-established. Here, we demonstrate that mTOR is over-expressed in both clinical tissue specimens and cultured human prostate cancer cells when compared to normal prostate tissues, respectively. Further, mTOR gene knockdown via lentivirus mediated mTOR specific shRNA resulted in a significant decrease in the viability and growth of prostate cancer cells without affecting normal human prostate cells. In addition, mTOR inhibition resulted in a significant i) decrease in 4EBP1, S6K, PI3K and AKT protein, ii) increase in PARP protein of prostate cancer cells. Most importantly, mTOR inhibition triggered apoptosis and suppressed pancreatic carcinoma growth in vivo in a mouse xenograft model. We suggest that targeting of mTOR may be a viable approach for the treatment of prostate cancer.
mTOR; prostatic carcinoma; apoptosis
In the present study, serum samples from 650 goats were collected from five counties between May and June 2012 and antibodies to Toxoplasma gondii were detected by indirect haemagglutination assay; 58 (9%) had antibodies to T. gondii with antibody titres of 1:64 to 1:1024. Seropositive samples were distributed in all five counties: seroprevalences in Kuandian county (15%, 21/139, 95% confidence interval [CI] 9–21%) were statistically different from the four other counties (Gaizhou, Huanren, Xiuyan and Liaoyang), and the seroprevalence difference between Xiuyan county (12%, 15/127, 95% CI 6–17%) and two other counties (Huanren, Liaoyang) was significantly different (P < 0.05). To our knowledge, this is the first report of the seroprevalence of T. gondii exposure in Liaoning cashmere goat in China. Our results indicated that Liaoning cashmere goat could be a potential reservoir for the transmission of T. gondii in Liaoning Province.
Toxoplasma gondii; seroprevalence; Liaoning cashmere goat; indirect haemagglutination assay
Objective: To construct a recombinant lentivirus vector driven by the PSMA promoter carrying mTOR-shRNA, and to obtain the effect on the mTOR gene silencing in human prostate cancer xenografts. Methods: The complimentary oligos of small interference RNA (siRNA) with hairpin structures targeting the mTOR gene and a negative control were synthesized, then ligated with pLV-PSMA-promoter vector and sequenced. The recombinant vectors were then transfected with viral packaging mix into 293T cells, viral supernatant was harvested to determine the titer. Prostate cancer cells infected by virus were harvested and the expression of mTOR (LV-PSMA-shmTOR), target proteins and cell growth were detected by reverse transcription-PCR (RT-PCR), Western blot and MTT separately. In established tumors derived from human prostate cancer cells, concentrated LV-PSMA-shmTOR lentivirus was injected intravenously in the tail vein of C4-2b tumor bearing female severe combined immunodeficient (SCID) mice. Tumor volume and immunohistochemistry was assessed. Results: Sequencing data showed that the constructed plasmids contained the correct sequences of mTOR siRNA transcript templates. A vector producing cell line 293T was established, and the titer for transfection was obtained. RT-PCR, Western blot and MTT analyses demonstrated that mTOR shRNA expression construct could suppress the expression of mTOR and inhibit the prostate cancer cell growth, specially. The tumor growth was suppressed in nude mouse. Conclusion: A PSMA driven lentivirus mediated siRNA targeting mTOR gene was successfully constructed, which decreased the expression of mTOR and induced the prostate cancer cell growth in vitro and in vivo. It has set up a research platform for the gene therapy of tumors which take mTOR as the target in the prostate cancer field.
mTOR; PSMA; prostatic carcinoma; apoptosis
Prostate cancer represents the leading cause of male death across the world. A recent genome-wide association study (GWAS) identified five novel susceptibility loci for prostate cancer in the Japanese population. This study is to replicate and fine map the potential association of these five loci with prostate cancer in the Chinese Han population.
In Phase I of the study, we tested the five single nucleotide polymorphisms (SNPs) which showed the strongest association evidence in the original GWAS in Japanese. The study sample consists of 1,169 Chinese Hans, comprising 483 patients and 686 healthy controls. Then in phase II, flanking SNPs of the successfully replicated SNPs in Phase I were genotyped and tested for association with prostate cancer to fine map those significant association signals.
We successfully replicated the association of rs13385191 (located in the C2orf43 gene, P = 8.60×10−5), rs12653946 (P = 1.33×10−6), rs1983891 (FOXP4, P = 6.22×10−5), and rs339331 (GPRC6A/RFX6, P = 1.42×10−5) with prostate cancer. The most significant odds ratio (OR) was recorded as 1.41 (95% confidence interval 1.18–1.68) for rs12653946. Rs9600079 did not show significant association (P = 8.07×10−2) with prostate cancer in this study. The Phase II study refined these association signals, and identified several SNPs showing more significant association with prostate cancer than the very SNPs tested in Phase I.
Our results provide further support for association of the C2orf43, FOXP4, GPRC6A and RFX6 genes with prostate cancer in Eastern Asian populations. This study also characterized the novel loci reported in the original GWAS with more details. Further work is still required to determine the functional variations and finally clarify the underlying biological mechanisms.
Osteoporotic hip fracture (HF) is a serious global public health problem associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of key measurable risk factors for HF, independent of bone mineral density (BMD). Hip BS is highly genetically determined, but genetic factors underlying BS variation are still poorly defined. Here, we performed an initial genome-wide copy number variation (CNV) association analysis for hip BS in 1,627 Chinese Han subjects using Affymetrix GeneChip Human Mapping SNP 6.0 Array and a follow-up replicate study in 2,286 unrelated US Caucasians sample. We found that a copy number polymorphism (CNP267) located at chromosome 2q12.2 was significantly associated with hip BS in both initial Chinese and replicate Caucasian samples with p values of 4.73E-03 and 5.66E-03, respectively. An important candidate gene, four and a half LIM domains 2 (FHL2), was detected at the downstream of CNP267, which plays important roles in bone metabolism by binding to several bone formation regulator, such as insulin-like growth factor-binding protein 5 (IGFBP-5) and androgen receptor (AR). Our findings suggest that CNP267 region may be associated with hip BS which might influence the FHL2 gene downstream.
Hip OF carries the highest morbidity and mortality. Previous studies revealed that individual genes/loci in the Tumor Necrosis Factor (TNF) -Related Apoptosis-Inducing Ligand (TRAIL) pathway were associated with bone metabolism. This study aims to verify the potential association between hip OF and TRAIL pathway.
Using genome-wide genotype data from Affymetrix 500 K SNP arrays, we performed novel pathway-based association analyses for hip OF in 700 elderly Chinese Han subjects (350 with hip OF and 350 healthy matched controls).
The TRAIL pathway achieved a significant p value (p = 0.01) for association with hip OF. Among the 38 genes in the TRAIL pathway, seven genes achieved nominally significant association with hip OF (p<0.05); the TNFSF10 (TRAIL) gene obtained the most significant p value (p = 1.70×10−4). SNPs (rs719126, rs6533015, rs9594738, rs1805034, rs11160706) from five genes (CFLAR, NFKB1, TNFSF11, TNFRSF11A, TRAF3) of the pathway had minor alleles that appear to be protective to hip OF. SNPs (rs6445063 and rs4259415) from two genes (TNFSF10 and TNFRSF10B) of the pathway had minor alleles (A) that are associated with an increased risk of hip OF, with the ORs (odds ratios) of 16.51 (95%CI:3.83–71.24) and 1.37 (95%CI:1.08–1.74), respectively.
Our study supports the potential role of the TRAIL pathway in the pathogenesis of hip OF in Chinese Han population. Further functional study of this pathway will be pursued to determine the mechanism by which it confers risk to hip OF.
Wrist fracture is not only one of the most common osteoporotic fractures but also a predictor of future fractures at other sites. Wrist bone mineral density (BMD) is an important determinant of wrist fracture risk, with high heritability. Specific genes underlying wrist BMD variation are largely unknown. Most published genome-wide association studies (GWASs) have focused only on a few top-ranking single-nucleotide polymorphisms (SNPs)/genes and considered each of the identified SNPs/genes independently. To identify biologic pathways important to wrist BMD variation, we used a novel pathway-based analysis approach in our GWAS of wrist ultradistal radius (UD) BMD, examining approximately 500,000 SNPs genome-wide from 984 unrelated whites. A total of 963 biologic pathways/gene sets were analyzed. We identified the regulation-of-autophagy (ROA) pathway that achieved the most significant result (p = .005, qfdr = 0.043, pfwer = 0.016) for association with UD BMD. The ROA pathway also showed significant association with arm BMD in the Framingham Heart Study sample containing 2187 subjects, which further confirmed our findings in the discovery cohort. Earlier studies indicated that during endochondral ossification, autophagy occurs prior to apoptosis of hypertrophic chondrocytes, and it also has been shown that some genes in the ROA pathway (e.g., INFG) may play important roles in osteoblastogenesis or osteoclastogenesis. Our study supports the potential role of the ROA pathway in human wrist BMD variation and osteoporosis. Further functional evaluation of this pathway to determine the mechanism by which it regulates wrist BMD should be pursued to provide new insights into the pathogenesis of wrist osteoporosis. © 2010 American Society for Bone and Mineral Research.
osteoporosis; bone mineral density; genome-wide association; regulation of autophagy; whites
The objective of this study was to identify genetic variants that are associated with adult leisure-time exercise behavior using genome-wide association in two independent samples.
Exercise behavior was measured in 1,772 unrelated Dutch and 978 unrelated American adults with detailed questions about type, frequency and duration of exercise. Individuals were classified into regular exercisers or non-exercisers using a threshold of 4 METhours (metabolic equivalents*hours per week). Regular exercisers were further divided into 5 categories of METhours, ranging from moderate (>=4 METhours) to highly vigorous (>=40 METhours) exercisers. Genome-wide association analyses with a total of 470,719 SNPs were conducted in both samples independently using regression-based techniques in SNPtest, including sex, age and BMI as covariates.
SNPs located in SGIP1, DNASE2B, PRSS16, ERCC2 and PAPSS2 were associated with exercise participation (combined p-value between 0.0004 and 4.5*10-6 with the same direction of allelic effects in both samples). Associations of candidate genes based on existing literature were replicated for the LEPR gene in the American sample (rs12405556, p=0.0005) and for the CYP19A1 gene in the Dutch sample (rs2470158, 0.0098).
Two genes (SGIP1 and LEPR) are expressed in the hypothalamus and involved in the regulation of energy homeostasis. Their effects were independent of BMI, suggesting a direct role of hypothalamic factors in the drive to exercise.
Physical activity; sports participation; genetics; genotype imputation; energy homeostasis
Poor femoral neck bone geometry at the femur is an important risk factor for hip fracture. We conducted a genome-wide association study (GWAS) of femoral neck bone geometry, examining approximately 379,000 eligible single-nucleotide polymorphisms (SNPs) in 1000 Caucasians. A common genetic variant, rs7430431 in the receptor transporting protein 3 (RTP3) gene, was identified in strong association with the buckling ratio (BR, P = 1.6 × 10−7), an index of bone structural instability, and with femoral cortical thickness (CT, P = 1.9 × 10−6). The RTP3 gene is located in 3p21.31, a region that we found to be linked with CT (LOD = 2.19, P = 6.0 × 10−4) in 3998 individuals from 434 pedigrees. The replication analyses in 1488 independent Caucasians and 2118 Chinese confirmed the association of rs7430431 to BR and CT (combined P = 7.0 × 10−3 for BR and P = 1.4 × 10−2 for CT). In addition, 350 hip fracture patients and 350 healthy control individuals were genotyped to assess the association of the RTP3 gene with the risk of hip fracture. Significant association between a nearby common SNP, rs10514713 of the RTP3 gene, and hip fracture (P = 1.0 × 10−3) was found. Our observations suggest that RTP3 may be a novel candidate gene for femoral neck bone geometry. © 2010 American Society for Bone and Mineral Research
genome-wide association; femoral neck bone geometry; bone fracture; RTP3
Recent success of genome-wide association studies (GWASs) on human height variation emphasized the effects of individual loci or genes. In this study, we used a developed pathway-based approach to further test biological pathways for potential association with stature, by examining ∼370 000 single-nucleotide polymorphisms (SNPs) across the human genome in 618 unrelated elder Han Chinese. A total of 626 biological pathways annotated by any of the three major public pathway databases (KEGG, BioCarta and Ambion GeneAssist Pathway Atlas) were tested. The regulation-of-autophagy (ROA) (nominal P=0.012) pathway was marginally significantly associated with human stature after our family wise error rate multiple-testing correction. We also used 1000 random recruited US whites for further replication. Interestingly, the ROA pathway presented the strongest signals in whites for height variation (nominal P=0.002). The results correspond to biological roles of the ROA pathway in human long bone development and growth. Our findings also implied that multiple-genetic factors may work jointly as a functional unit (pathway), and the traditional GWASs could have missed important genetic information imbedded in those less significant markers.
autophagy; GWAS; height; pathway; stature
Bone mineral density (BMD) and femoral neck cross-sectional geometric parameters (FNCSGPs) are under strong genetic control. DNA copy number variation (CNV) is an important source of genetic diversity for human diseases. This study aims to identify CNVs associated with BMD and FNCSGPs.
Genome-wide CNV association analyses were conducted in 1,000 unrelated Caucasian subjects for BMD at the spine, hip, femoral neck, and for three FNCSGPs - cortical thickness (CT), cross section area (CSA), and buckling ratio (BR). BMD was measured by dual energy X-ray absorptiometry (DEXA). CT, CSA, and BR were estimated using DEXA measurements. Affymetrix 500K arrays and copy number analysis tool was used to identify CNVs.
A CNV in VPS13B gene was significantly associated with spine, hip and FN BMDs, and CT, CSA, and BR (p < 0.05). Compared to subjects with 2 copies of the CNV, carriers of one copy had an average of 14.6%, 12.4%, and 13.6% higher spine, hip, and FN BMD, 20.0% thicker CT, 10.6% larger CSA, and 12.4% lower BR. Thus a decrease of the CNV consistently produced stronger bone, thereby reducing osteoporotic fracture risk.
VPS13B gene, via affecting BMD and FNCSGPs, is a novel osteoporosis risk gene
copy number variation; bone mineral density; bone geometry; osteoporosis
Femoral neck compression strength index (fCSI), a novel phenotypic parameter that integrates bone density, bone size, and body size, has significant potential to improve hip fracture risk assessment. The genetic factors underlying variations in fCSI, however, remain largely unknown. Given the important roles of the receptor activator of the nuclear factor-κB ligand/receptor activator of the nuclear factor-κB/osteoprotegerin (RANKL/RANK/OPG) pathway in the regulation of bone remodeling, we tested the associations between RANKL/RANK/OPG polymorphisms and variations in fCSI as well as its components (femoral neck bone mineral density [fBMD], femoral neck width [FNW], and weight). This was accomplished with a sample comprising 1873 subjects from 405 Caucasian nuclear families. Of the 37 total SNPs studied in these three genes, 3 SNPs, namely, rs12585014, rs7988338, and rs2148073, of RANKL were significantly associated with fCSI (P = 0.0007, 0.0007, and 0.0005, respectively) after conservative Bonferroni correction. Moreover, the three SNPs were approximately in complete linkage disequilibrium. Haplotype-based association tests corroborated the single-SNP results since haplotype 1 of block 1 of the RANKL gene achieved an even more significant association with fCSI (P = 0.0003) than any of the individual SNPs. However, we did not detect any significant associations of these genes with fBMD, FNW, or weight. In summary, our findings suggest that the RANKL gene may play an important role in variation in fCSI, independent of fBMD and non-fBMD components.
Femoral neck compression strength index; Femoral neck bone mineral density; Femoral neck width; RANKL/RANK/OPG gene; Quantitative transmission disequilibrium test
Human stature, as an important physical index in clinical practice and a usual covariate in gene mapping of complex disorders, is a highly heritable complex trait. To identify specific genes underlying stature, a genome-wide association study was performed in 1000 unrelated homogeneous Caucasian subjects using Affymetrix 500K arrays. A group of seven contiguous markers in the region of SBF2 gene (Set-binding factor 2) are associated with stature, significantly so at the genome-wide level after false discovery rate (FDR) correction (FDR q = 0.034–0.042). Three SNPs in another SNP group in the Filamin B (FLNB) gene were also associated with stature, significantly so with FDR q = 0.042–0.048. In follow-up independent replication studies, rs10734652 in the SBF2 gene was significantly (P = 0.036) and suggestively (P = 0.07) associated with stature in Caucasian families and 1306 unrelated Caucasian subjects, respectively, and rs9834312 in the FLNB gene was also associated with stature in such two independent Caucasian populations (P = 0.008 in unrelated sample and P = 0.049 in family sample). Particularly, additional significant replication association signals were detected in Chinese, an ethnic population different from Caucasian, between rs9834312 and stature in 619 unrelated northern Chinese subjects (P = 0.017), as well as between rs10734652 and stature in 2953 unrelated southern Chinese subjects (P = 0.048). This study also provides additional replication evidence for some of the already published stature loci. These results, together with the known functional relevance of the SBF2 and FLNB genes to skeletal linear growth and bone formation, support that two regions containing FLNB and SBF2 genes are two novel loci underlying stature variation.
Cigarette smoking is the leading preventable cause of death in the US. Although smoking behavior has a significant genetic determination, the specific genes and associated mechanisms underlying smoking behavior are largely unknown. Here, we performed a genome-wide association study on smoking behavior in 840 Caucasians, including 417 males and 423 females, in which we examined ∼380,000 SNPs. We found that a cluster of nine SNPs upstream from the IL15 gene were associated with smoking status in males, with the most significant SNP, rs4956302, achieving a p value (8.80×10−8) of genome-wide significance. Another SNP, rs17354547, that is highly conserved across multiple species, achieved a p value of 5.65×10−5. These two SNPs, together with two additional SNPs (rs1402812 and rs4956396) were selected from the above nine SNPs for replication in an African-American sample containing 1,251 subjects, including 412 males and 839 females. The SNP rs17354547 was successfully replicated in the male subgroup of the replication sample; it was associated with smoking quantity (SQ), the Heaviness of Smoking Index (HSI) and the Fagerstrom Test for Nicotine Dependence (FTND), with p values of 0.031, 0.0046 and 0.019, respectively. In addition, a haplotype formed by rs17354547, rs1402812 and rs4956396 was also associated with SQ, HSI and FTND, achieving p values of 0.039, 0.0093 and 0.0093, respectively. To further confirm our findings, we performed an in silico replication study of the nine SNPs in a Framingham Heart Study sample containing 7,623 Caucasians from 1,731 families, among which, 3,491 subjects are males and 4,132 are females. Again, male-specific association with smoking status was observed, for which seven of the nine SNPs achieved significant p values (p<0.05) and two achieved marginally significant p values (p<0.10) in males. Several of the nine SNPs, including the highly conserved one across species, rs17354547, are located at potential transcription factor binding sites, suggesting transcription regulation as a possible function for these SNPs. Through this function, the SNPs may modulate gene expression of IL15, a key cytokine regulating immune function. As the immune system has long been recognized to influence drug addiction behavior, our association findings suggest a novel mechanism for smoking addiction involving immune modulation via the IL15 pathway.
smoking; nicotine addiction; IL15; genomewide association; genetics
Current genome-wide association studies (GWAS) are normally implemented in a univariate framework and analyze different phenotypes in isolation. This univariate approach ignores the potential genetic correlation between important disease traits. Hence this approach is difficult to detect pleiotropic genes, which may exist for obesity and osteoporosis, two common diseases of major public health importance that are closely correlated genetically.
To identify such pleiotropic genes and the key mechanistic links between the two diseases, we here performed the first bivariate GWAS of obesity and osteoporosis. We searched for genes underlying co-variation of the obesity phenotype, body mass index (BMI), with the osteoporosis risk phenotype, hip bone mineral density (BMD), scanning ∼380,000 SNPs in 1,000 unrelated homogeneous Caucasians, including 499 males and 501 females. We identified in the male subjects two SNPs in intron 1 of the SOX6 (SRY-box 6) gene, rs297325 and rs4756846, which were bivariately associated with both BMI and hip BMD, achieving p values of 6.82×10−7 and 1.47×10−6, respectively. The two SNPs ranked at the top in significance for bivariate association with BMI and hip BMD in the male subjects among all the ∼380,000 SNPs examined genome-wide. The two SNPs were replicated in a Framingham Heart Study (FHS) cohort containing 3,355 Caucasians (1,370 males and 1,985 females) from 975 families. In the FHS male subjects, the two SNPs achieved p values of 0.03 and 0.02, respectively, for bivariate association with BMI and femoral neck BMD. Interestingly, SOX6 was previously found to be essential to both cartilage formation/chondrogenesis and obesity-related insulin resistance, suggesting the gene's dual role in both bone and fat.
Our findings, together with the prior biological evidence, suggest the SOX6 gene's importance in co-regulation of obesity and osteoporosis.
In Caucasian, several studies have identified some common variants associated with human stature variation. However, no such study was performed in Chinese, which is the largest population in the world and evidently differs from Caucasian in genetic background. To identify common or ethnic specific genes for stature in Chinese, an initial GWAS and follow-up replication study were performed. Our initial GWAS study found that a group of 13 contiguous SNPs, which span a region of ∼150 kb containing two neighboring genes, zinc finger protein (ZNP) 510 and ZNP782, achieved strong signals for association with stature, with P values ranging from 9.71 × 10−5 to 3.11 × 10−6. After false discovery rate correction for multiple testing, 9 of the 13 SNPs remain significant (FDR q = 0.036–0.046). The follow-up replication study in an independent 2,953 unrelated southern Chinese confirmed the association of rs10816533 with stature (P = 0.029). All the13 SNPs were in consistently strong linkage disequilibrium (D′ > 0.99) and formed a single perfect haplotype block. The minor allele frequencies for the 13 contiguous SNPs have evidently ethnic difference, which range from 0.21 to 0.33 in Chinese but have as low as ∼0.017 reported in dbSNP database in Caucasian. The present results suggest that the genomic region containing the ZNP510 and ZNP782 genes is an ethnic specific locus associated with stature variation in Chinese.
Obesity is a major public health problem with strong genetic determination; however, the genetic factors underlying obesity are largely unknown. In this study, we performed a genome-wide association scan for obesity by examining approximately 500 000 single-nucleotide polymorphisms (SNPs) in a sample of 1000 unrelated US Caucasians. We identified a novel gene, CTNNBL1, which has multiple SNPs associated with body mass index (BMI) and fat mass. The most significant SNP, rs6013029, achieved experiment-wise P-values of 2.69 × 10−7 for BMI and of 4.99 × 10−8 for fat mass, respectively. The SNP rs6013029 minor allele T confers an average increase in BMI and fat mass of 2.67 kg/m2 and 5.96 kg, respectively, compared with the alternative allele G. We further genotyped the five most significant CTNNBL1 SNPs in a French case–control sample comprising 896 class III obese adults (BMI ≥ 40 kg/m2) and 2916 lean adults (BMI < 25 kg/m2). All five SNPs showed consistent associations with obesity (8.83 × 10−3 < P < 6.96 × 10−4). Those subjects who were homozygous for the rs6013029 T allele had 1.42-fold increased odds of obesity compared with those without the T allele. The protein structure of CTNNBL1 is homologous to β-catenin, a family of proteins containing armadillo repeats, suggesting similar biological functions. β-Catenin is involved in the Wnt/β-catenin-signaling pathway which appears to contribute to maintaining the undifferentiated state of pre-adipocytes by inhibiting adipogenic gene expression. Our study hence suggests a novel mechanism for the development of obesity, where CTNNBL1 may play an important role. Our study also provided supportive evidence for previously identified associations between obesity and INSIG2 and PFKP, but not FTO.
For females, menarche is a most significant physiological event. Age at menarche (AAM) is a trait with high genetic determination and is associated with major complex diseases in women. However, specific genes for AAM variation are largely unknown. To identify genetic factors underlying AAM variation, a genome-wide association study (GWAS) examining about 380,000 SNPs was conducted in 477 Caucasian women. A follow-up replication study was performed to validate our major GWAS findings using two independent Caucasian cohorts with 854 siblings and 762 unrelated subjects, respectively, and one Chinese cohort of 1,387 unrelated subjects—all females. Our GWAS identified a novel gene, SPOCK (Sparc/Osteonectin, CWCV, and Kazal-like domains proteoglycan), which had seven SNPs associated with AAM with genome-wide false discovery rate (FDR) q<0.05. Six most significant SNPs of the gene were selected for validation in three independent replication cohorts. All of the six SNPs were replicated in at least one cohort. In particular, SNPs rs13357391 and rs1859345 were replicated both within and across different ethnic groups in all three cohorts, with p values of 5.09×10−3 and 4.37×10−3, respectively, in the Chinese cohort and combined p values (obtained by Fisher's method) of 5.19×10−5 and 1.02×10−4, respectively, in all three replication cohorts. Interestingly, SPOCK can inhibit activation of MMP-2 (matrix metalloproteinase-2), a key factor promoting endometrial menstrual breakdown and onset of menstrual bleeding. Our findings, together with the functional relevance, strongly supported that the SPOCK gene underlies variation of AAM.
Menarche is a physical milestone in a woman's life. Age at menarche (AAM) is related to many common female health problems. AAM is mainly determined by genetic factors. However, the specific genes and the associated mechanisms underlying AAM are largely unknown. Here, taking advantage of the most recent technological advances in the field of human genetics, we identified multiple genetic variants in a gene, SPOCK, which are associated with AAM variation in a group of Caucasian women. This association was subsequently confirmed not only in two independent groups of Caucasian women but also across ethnic boundaries in one group of Chinese women. In addition, SPOCK has a function in regulating a key factor involved in menstrual cycles, MMP-2, which provides further support to our findings. Our study provides a solid basis for further investigation of the gene, which may help to reveal the underlying mechanisms for the timing of menarche and for AAM's relationship with women's health in general.
Osteoporosis, the most prevalent metabolic bone disease among older people, increases risk for low trauma hip fractures (HF) that are associated with high morbidity and mortality. Hip bone size (BS) has been identified as one of the key measurable risk factors for HF. Although hip BS is highly genetically determined, genetic factors underlying the trait are still poorly defined. Here, we performed the first genome-wide association study (GWAS) of hip BS interrogating ∼380,000 SNPs on the Affymetrix platform in 1,000 homogeneous unrelated Caucasian subjects, including 501 females and 499 males. We identified a gene, PLCL1 (phospholipase c-like 1), that had four SNPs associated with hip BS at, or approaching, a genome-wide significance level in our female subjects; the most significant SNP, rs7595412, achieved a p value of 3.72×10−7. The gene's importance to hip BS was replicated using the Illumina genotyping platform in an independent UK cohort containing 1,216 Caucasian females. Two SNPs of the PLCL1 gene, rs892515 and rs9789480, surrounded by the four SNPs identified in our GWAS, achieved p values of 8.62×10−3 and 2.44×10−3, respectively, for association with hip BS. Imputation analyses on our GWAS and the UK samples further confirmed the replication signals; eight SNPs of the gene achieved combined imputed p values<10−5 in the two samples. The PLCL1 gene's relevance to HF was also observed in a Chinese sample containing 403 females, including 266 with HF and 177 control subjects. A SNP of the PLCL1 gene, rs3771362 that is only ∼0.6 kb apart from the most significant SNP detected in our GWAS (rs7595412), achieved a p value of 7.66×10−3 (odds ratio = 0.26) for association with HF. Additional biological support for the role of PLCL1 in BS comes from previous demonstrations that the PLCL1 protein inhibits IP3 (inositol 1,4,5-trisphosphate)-mediated calcium signaling, an important pathway regulating mechanical sensing of bone cells. Our findings suggest that PLCL1 is a novel gene associated with variation in hip BS, and provide new insights into the pathogenesis of HF.