Gut microbiota accelerates obesity in peri-/post-menopausal women via Bacteroides fragilis and acetic acid.

OBJECTIVE: Many animal experiments and epidemiological studies have shown that the gut microbiota (GM) plays an important role in the development of obesity, but the specific biological mechanism involved in the pathogenesis of disease remain unknown. We aimed to examine the relationships and functional mechanisms of GM on obesity in peri- and post-menopausal women. METHODS: We recruited 499 Chinese peri- and post-menopausal women and performed comprehensive analyses of the gut microbiome, targeted metabolomics for short-chain fatty acids in serum, and host whole-genome sequencing by various association analysis methods. RESULTS: Through constrained linear regression analysis, we found that an elevated abundance of Bacteroides fragilis (B. fragilis) was associated with obesity. We also found that serum levels of acetic acid were negatively associated with obesity, and that B. fragilis was negatively associated with serum acetic acid levels by partial Spearman correlation analysis. Mendelian randomization analysis indicated that B. fragilis increases the risk of obesity and may causally down-regulate acetic acid levels. CONCLUSIONS: We found the gut with B. fragilis may accelerate obesity, in part, by suppressing acetic acid levels. Therefore, B. fragilis and acetic acid may represent important therapeutic targets for obesity intervention in peri- and post-menopausal women.

Identification of IDUA and WNT16 Phosphorylation-Related Non-Synonymous Polymorphisms for Bone Mineral Density in Meta-Analyses of Genome-Wide Association Studies.

Protein phosphorylation regulates a wide variety of cellular processes. Thus, we hypothesize that single-nucleotide polymorphisms (SNPs) that may modulate protein phosphorylation could affect osteoporosis risk. Based on a previous conventional genome-wide association (GWA) study, we conducted a three-stage meta-analysis targeting phosphorylation-related SNPs (phosSNPs) for femoral neck (FN)-bone mineral density (BMD), total hip (HIP)-BMD, and lumbar spine (LS)-BMD phenotypes. In stage 1, 9593 phosSNPs were meta-analyzed in 11,140 individuals of various ancestries. Genome-wide significance (GWS) and suggestive significance were defined by α = 5.21 × 10(-6) (0.05/9593) and 1.00 × 10(-4), respectively. In stage 2, nine stage 1-discovered phosSNPs (based on α = 1.00 × 10(-4)) were in silico meta-analyzed in Dutch, Korean, and Australian cohorts. In stage 3, four phosSNPs that replicated in stage 2 (based on α = 5.56 × 10(-3), 0.05/9) were de novo genotyped in two independent cohorts. IDUA rs3755955 and rs6831280, and WNT16 rs2707466 were associated with BMD phenotypes in each respective stage, and in three stages combined, achieving GWS for both FN-BMD (p = 8.36 × 10(-10), p = 5.26 × 10(-10), and p = 3.01 × 10(-10), respectively) and HIP-BMD (p = 3.26 × 10(-6), p = 1.97 × 10(-6), and p = 1.63 × 10(-12), respectively). Although in vitro studies demonstrated no differences in expressions of wild-type and mutant forms of IDUA and WNT16B proteins, in silico analyses predicts that WNT16 rs2707466 directly abolishes a phosphorylation site, which could cause a deleterious effect on WNT16 protein, and that IDUA phosSNPs rs3755955 and rs6831280 could exert indirect effects on nearby phosphorylation sites. Further studies will be required to determine the detailed and specific molecular effects of these BMD-associated non-synonymous variants.

Identification of a novel FGFRL1 MicroRNA target site polymorphism for bone mineral density in meta-analyses of genome-wide association studies.

MicroRNAs (miRNAs) are critical post-transcriptional regulators. Based on a previous genome-wide association (GWA) scan, we conducted a polymorphism in microRNA target sites (poly-miRTS)-centric multistage meta-analysis for lumbar spine (LS)-, total hip (HIP)- and femoral neck (FN)-bone mineral density (BMD). In stage I, 41 102 poly-miRTSs were meta-analyzed in seven cohorts with a genome-wide significance (GWS) α = 0.05/41 102 = 1.22 × 10(-6). By applying α = 5 × 10(-5) (suggestive significance), 11 poly-miRTSs were selected, with FGFRL1 rs4647940 and PRR5 rs3213550 as top signals for FN-BMD (P = 7.67 × 10(-6) and 1.58 × 10(-5)) in gender-combined sample. In stage II in silico replication (two cohorts), FGFRL1 rs4647940 was the only signal marginally replicated for FN-BMD (P = 5.08 × 10(-3)) at α = 0.10/11 = 9.09 × 10(-3). PRR5 rs3213550 was also selected based on biological significance. In stage III de novo genotyping replication (two cohorts), FGFRL1 rs4647940 was the only signal significantly replicated for FN-BMD (P = 7.55 × 10(-6)) at α = 0.05/2 = 0.025 in gender-combined sample. Aggregating three stages, FGFRL1 rs4647940 was the single stage I-discovered and stages II- and III-replicated signal attaining GWS for FN-BMD (P = 8.87 × 10(-12)). Dual-luciferase reporter assays demonstrated that FGFRL1 3' untranslated region harboring rs4647940 appears to be hsa-miR-140-5p's target site. In a zebrafish microinjection experiment, dre-miR-140-5p is shown to exert a dramatic impact on craniofacial skeleton formation. Taken together, we provided functional evidence for a novel FGFRL1 poly-miRTS rs4647940 in a previously known 4p16.3 locus, and experimental and clinical genetics studies have shown both FGFRL1 and hsa-miR-140-5p are important for bone formation.

Multistage genome-wide association meta-analyses identified two new loci for bone mineral density.

Aiming to identify novel genetic variants and to confirm previously identified genetic variants associated with bone mineral density (BMD), we conducted a three-stage genome-wide association (GWA) meta-analysis in 27 061 study subjects. Stage 1 meta-analyzed seven GWA samples and 11 140 subjects for BMDs at the lumbar spine, hip and femoral neck, followed by a Stage 2 in silico replication of 33 SNPs in 9258 subjects, and by a Stage 3 de novo validation of three SNPs in 6663 subjects. Combining evidence from all the stages, we have identified two novel loci that have not been reported previously at the genome-wide significance (GWS; 5.0 × 10(-8)) level: 14q24.2 (rs227425, P-value 3.98 × 10(-13), SMOC1) in the combined sample of males and females and 21q22.13 (rs170183, P-value 4.15 × 10(-9), CLDN14) in the female-specific sample. The two newly identified SNPs were also significant in the GEnetic Factors for OSteoporosis consortium (GEFOS, n = 32 960) summary results. We have also independently confirmed 13 previously reported loci at the GWS level: 1p36.12 (ZBTB40), 1p31.3 (GPR177), 4p16.3 (FGFRL1), 4q22.1 (MEPE), 5q14.3 (MEF2C), 6q25.1 (C6orf97, ESR1), 7q21.3 (FLJ42280, SHFM1), 7q31.31 (FAM3C, WNT16), 8q24.12 (TNFRSF11B), 11p15.3 (SOX6), 11q13.4 (LRP5), 13q14.11 (AKAP11) and 16q24 (FOXL1). Gene expression analysis in osteogenic cells implied potential functional association of the two candidate genes (SMOC1 and CLDN14) in bone metabolism. Our findings independently confirm previously identified biological pathways underlying bone metabolism and contribute to the discovery of novel pathways, thus providing valuable insights into the intervention and treatment of osteoporosis.

Meta-analysis of genome-wide association data identifies novel susceptibility loci for obesity.

Obesity is a major public health problem with strong genetic determination. Multiple genetic variants have been implicated for obesity by conducting genome-wide association (GWA) studies, primarily focused on body mass index (BMI). Fat body mass (FBM) is phenotypically more homogeneous than BMI and is more appropriate for obesity research; however, relatively few studies have been conducted on FBM. Aiming to identify variants associated with obesity, we carried out meta-analyses of seven GWA studies for BMI-related traits including FBM, and followed these analyses by de novo replication. The discovery cohorts consisted of 21 969 individuals from diverse ethnic populations and a total of over 4 million genotyped or imputed SNPs. The de novo replication cohorts consisted of 6663 subjects from two independent samples. To complement individual SNP-based association analyses, we also carried out gene-based GWA analyses in which all variations within a gene were considered jointly. Individual SNP-based association analyses identified a novel locus 1q21 [rs2230061, CTSS (Cathepsin S)] that was associated with FBM after the adjustment of lean body mass (LBM) (P = 3.57 × 10(-8)) at the genome-wide significance level. Gene-based association analyses identified a novel gene NLK (nemo-like kinase) in 17q11 that was significantly associated with FBM adjusted by LBM. In addition, we confirmed three previously reported obesity susceptibility loci: 16q12 [rs62033400, P = 1.97 × 10(-14), FTO (fat mass and obesity associated)], 18q22 [rs6567160, P = 8.09 × 10(-19), MC4R (melanocortin 4 receptor)] and 2p25 [rs939583, P = 1.07 × 10(-7), TMEM18 (transmembrane protein 18)]. We also found that rs6567160 may exert pleiotropic effects to both FBM and LBM. Our results provide additional insights into the molecular genetic basis of obesity and may provide future targets for effective prevention and therapeutic intervention.

SNP rs6265 regulates protein phosphorylation and osteoblast differentiation and influences BMD in humans.

Bone mineral density (BMD) is a major index for diagnosing osteoporosis. PhosSNPs are nonsynonymous SNPs that affect protein phosphorylation. The relevance and significance of phosSNPs to BMD and osteoporosis is unknown. This study aimed to identify and characterize phosSNPs significant for BMD in humans. We conducted a pilot genomewide phosSNP association study for BMD in three independent population samples, involving ∼5000 unrelated individuals. We identified and replicated three phosSNPs associated with both spine BMD and hip BMD in Caucasians. Association with hip BMD for one of these phosSNPs, ie, rs6265 (major/minor allele: G/A) in BDNF gene, was also suggested in Chinese. Consistently in both ethnicities, individuals carrying the AA genotype have significantly lower hip BMD than carriers of the GA and GG genotypes. Through in vitro molecular and cellular studies, we found that compared to osteoblastic cells transfected with wild-type BDNF-Val66 (encoded with allele G at rs6265), transfection of variant BDNF-Met66 (encoded with allele A at rs6265) significantly decreased BDNF protein phosphorylation (at amino acid residue T62), expression of osteoblastic genes (OPN, BMP2, and ALP), and osteoblastic activity. The findings are consistent with and explain our prior observations in general human populations. We conclude that phosSNP rs6265, by regulating BDNF protein phosphorylation and osteoblast differentiation, influences hip BMD in humans. This study represents our first endeavor to dissect the functions of phosSNPs in bone, which might stimulate extended large-scale studies on bone or similar studies on other human complex traits and diseases.

Bivariate genome-wide association analyses identified genes with pleiotropic effects for femoral neck bone geometry and age at menarche.

Femoral neck geometric parameters (FNGPs), which include cortical thickness (CT), periosteal diameter (W), buckling ratio (BR), cross-sectional area (CSA), and section modulus (Z), contribute to bone strength and may predict hip fracture risk. Age at menarche (AAM) is an important risk factor for osteoporosis and bone fractures in women. Some FNGPs are genetically correlated with AAM. In this study, we performed a bivariate genome-wide association study (GWAS) to identify new candidate genes responsible for both FNGPs and AAM. In the discovery stage, we tested 760,794 SNPs in 1,728 unrelated Caucasian subject, followed by replication analyses in independent samples of US Caucasians (with 501 subjects) and Chinese (with 826 subjects). We found six SNPs that were associated with FNGPs and AAM. These SNPs are located in three genes (i.e. NRCAM, IDS and LOC148145), suggesting these three genes may co-regulate FNGPs and AAM. Our findings may help improve the understanding of genetic architecture and pathophysiological mechanisms underlying both osteoporosis and AAM.

ANKRD7 and CYTL1 are novel risk genes for alcohol drinking behavior.

BACKGROUND: 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. METHODS: 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. RESULTS: 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. CONCLUSION: The evidence suggests that ANKRD7 and CYTL1 genes may play an important role in the variance in AD risk.

Genome-wide association study of copy number variants suggests LTBP1 and FGD4 are important for alcohol drinking.

Alcohol dependence (AD) is a complex disorder characterized by psychiatric and physiological dependence on alcohol. AD is reflected by regular alcohol drinking, which is highly inheritable. In this study, to identify susceptibility genes associated with alcohol drinking, we performed a genome-wide association study of copy number variants (CNVs) in 2,286 Caucasian subjects with Affymetrix SNP6.0 genotyping array. We replicated our findings in 1,627 Chinese subjects with the same genotyping array. We identified two CNVs, CNV207 (combined p-value 1.91E-03) and CNV1836 (combined p-value 3.05E-03) that were associated with alcohol drinking. CNV207 and CNV1836 are located at the downstream of genes LTBP1 (870 kb) and FGD4 (400 kb), respectively. LTBP1, by interacting TGFB1, may down-regulate enzymes directly participating in alcohol metabolism. FGD4 plays a role in clustering and trafficking GABA(A) receptor and subsequently influence alcohol drinking through activating CDC42. Our results provide suggestive evidence that the newly identified CNV regions and relevant genes may contribute to the genetic mechanism of alcohol dependence.

Genome-wide association study of copy number variation identified gremlin1 as a candidate gene for lean body mass.

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.

An integrative study ascertained SOD2 as a susceptibility gene for osteoporosis in Chinese.

Osteoporosis is characterized by low BMD and has strong genetic determination. However, specific genetic variants influencing BMD and contributing to the pathogenesis of osteoporosis are largely uncharacterized. Current genetic studies in bone, which are aimed at identification of osteoporosis risk genes, are focused mostly on DNA, RNA, or the protein level individually, lacking integrative evidence from the three levels of genetic information flow to confidently ascertain the significance of genes for osteoporosis. Our previous proteomics study discovered that superoxide dismutase 2 (SOD2) in circulating monocytes (CMCs, ie, potential osteoclast precursors) was significantly upregulated at protein level in vivo in Chinese with low versus high hip BMD. Herein, at mRNA level, we found that SOD2 gene expression also was upregulated in CMCs (p < 0.05) in Chinese with low versus high hip BMD. At the DNA level, in 1627 unrelated Chinese subjects, we identified eight single-nucleotide polymorphisms (SNPs) at the SOD2 gene locus that were suggestively associated with hip BMD (peak signal at rs11968525, p = 0.048). Among the eight SNPs, three SNPs (rs7754103, rs7754295, and rs2053949) were associated with the SOD2 mRNA expression level (p < 0.05), suggesting that they are expression quantitative trait loci (eQTLs) regulating SOD2 gene expression. In conclusion, this integrative evidence from DNA, RNA, and protein levels support SOD2 as a susceptibility gene for osteoporosis.

Pathway-based association analyses identified TRAIL pathway for osteoporotic fractures.

INTRODUCTION: 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. METHODS: 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). RESULTS: 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. CONCLUSIONS: 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.

The regulation-of-autophagy pathway may influence Chinese stature variation: evidence from elder adults.

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 approximately 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.

Genome-wide association study identifies ALDH7A1 as a novel susceptibility gene for osteoporosis.

Osteoporosis is a major public health problem. It is mainly characterized by low bone mineral density (BMD) and/or low-trauma osteoporotic fractures (OF), both of which have strong genetic determination. The specific genes influencing these phenotypic traits, however, are largely unknown. Using the Affymetrix 500K array set, we performed a case-control genome-wide association study (GWAS) in 700 elderly Chinese Han subjects (350 with hip OF and 350 healthy matched controls). A follow-up replication study was conducted to validate our major GWAS findings in an independent Chinese sample containing 390 cases with hip OF and 516 controls. We found that a SNP, rs13182402 within the ALDH7A1 gene on chromosome 5q31, was strongly associated with OF with evidence combined GWAS and replication studies (P = 2.08x10(-9), odds ratio = 2.25). In order to explore the target risk factors and potential mechanism underlying hip OF risk, we further examined this candidate SNP's relevance to hip BMD both in Chinese and Caucasian populations involving 9,962 additional subjects. This SNP was confirmed as consistently associated with hip BMD even across ethnic boundaries, in both Chinese and Caucasians (combined P = 6.39x10(-6)), further attesting to its potential effect on osteoporosis. ALDH7A1 degrades and detoxifies acetaldehyde, which inhibits osteoblast proliferation and results in decreased bone formation. Our findings may provide new insights into the pathogenesis of osteoporosis.

Gene expression profiling in monocytes and SNP association suggest the importance of the STAT1 gene for osteoporosis in both Chinese and Caucasians.

Osteoporosis is characterized mainly by low bone mineral density (BMD). Many cytokines and chemokines have been related with bone metabolism. Monocytes in the immune system are important sources of cytokines and chemokines for bone metabolism. However, no study has investigated in vivo expression of a large number of various factors simultaneously in human monocytes underlying osteoporosis. This study explored the in vivo expression pattern of general cytokines, chemokines, and their receptor genes in human monocytes and validated the significant genes by qRT-PCR and genetic association analyses. Expression profilings were performed in monocyte samples from 26 Chinese and 20 Caucasian premenopausal women with discordant BMD. Genome-wide association analysis with BMD variation was conducted in 1000 unrelated Caucasians. We selected 168 cytokines, chemokines, osteoclast-related factors, and their receptor genes for analyses. Significantly, the signal transducer and activator of transcription 1 (STAT1) gene was upregulated in the low versus the high BMD groups in both Chinese and Caucasians. We also revealed a significant association of the STAT1 gene with BMD variation in the 1000 Caucasians. Thus we conclude that the STAT1 gene is important in human circulating monocytes in the etiology of osteoporosis.

Genome-wide association and replication studies identified TRHR as an important gene for lean body mass.

Low lean body mass (LBM) is related to a series of health problems, such as osteoporotic fracture and sarcopenia. Here we report a genome-wide association (GWA) study on LBM variation, by using Affymetrix 500K single-nucleotide polymorphism (SNP) arrays. In the GWA scan, we tested 379,319 eligible SNPs in 1,000 unrelated US whites and found that two SNPs, rs16892496 (p = 7.55 x 10(-8)) and rs7832552 (p = 7.58 x 10(-8)), within the thyrotropin-releasing hormone receptor (TRHR) gene were significantly associated with LBM. Subjects carrying unfavorable genotypes at rs16892496 and rs7832552 had, on average, 2.70 and 2.55 kg lower LBM, respectively, compared to those with alternative genotypes. We replicated the significant associations in three independent samples: (1) 1488 unrelated US whites, (2) 2955 Chinese unrelated subjects, and (3) 593 nuclear families comprising 1972 US whites. Meta-analyses of the GWA scan and the replication studies yielded p values of 5.53 x 10(-9) for rs16892496 and 3.88 x 10(-10) for rs7832552. In addition, we found significant interactions between rs16892496 and polymorphisms of several other genes involved in the hypothalamic-pituitary-thyroid and the growth hormone-insulin-like growth factor-I axes. Results of this study, together with the functional relevance of TRHR in muscle metabolism, support the TRHR gene as an important gene for LBM variation.

Genome-wide association analyses identify SPOCK as a key novel gene underlying age at menarche.

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 x 10(-3) and 4.37 x 10(-3), respectively, in the Chinese cohort and combined p values (obtained by Fisher's method) of 5.19 x 10(-5) and 1.02 x 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.

An in vivo genome wide gene expression study of circulating monocytes suggested GBP1, STAT1 and CXCL10 as novel risk genes for the differentiation of peak bone mass.

Peak bone mass (PBM) is an important determinant of osteoporosis. Circulating monocytes serve as early progenitors of osteoclasts and produce important molecules for bone metabolism. To search for genes functionally important for PBM variation, we performed a whole genome gene differential expression study of circulating monocytes in human premenopausal subjects with extremely low (N=12) vs. high (N=14) PBM. We used Affymetrix HG-U133 plus2.0 GeneChip arrays. We identified 70 differential expression probe sets (p<0.01) corresponding to 49 unique genes. After false discovery rate adjustment, three genes [STAT1, signal transducer and activator of transcription 1; GBP1, guanylate binding protein 1; CXCL10, Chemokine (C-X-C motif) ligand 10] expressed significantly differentially (p<0.05). The RT-PCR results independently confirmed the significantly differential expression of GBP1 gene, and the differential expression trend of STAT1. Functional analyses suggested that the three genes are associated with the osteoclastogenic processes of proliferation, migration, differentiation, migration, chemotaxis, adhesion. Therefore, we may tentatively hypothesize that the three genes may potentially contribute to differential osteoclastogenesis, which may in the end lead to differential PBM. Our results indicate that the GBP1, STAT1 and CXCL10 may be novel risk genes for the differentiation of PBM at the monocyte stage.

Genome-wide association scan for stature in Chinese: evidence for ethnic specific loci.

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 approximately 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 x 10(-5) to 3.11 x 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 the 13 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 approximately 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.

Genome-wide association study identifies two novel loci containing FLNB and SBF2 genes underlying stature variation.

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.