Person:

Hsu, Yi-Hsiang

Background: Genome-wide association (GWA) studies that use population-based association approaches may identify spurious associations in the presence of population admixture. In this paper, we propose a novel three-stage approach that is computationally efficient and robust to population admixture and more powerful than the family-based association test (FBAT) for GWA studies with family data. We propose a three-stage approach for GWA studies with family data. The first stage is to perform linear regression ignoring phenotypic correlations among family members. SNPs with a first stage p-value below a liberal cut-off (e.g. 0.1) are then analyzed in the second stage that employs a linear mixed effects (LME) model that accounts for within family correlations. Next, SNPs that reach genome-wide significance (e.g. (10^{-6}) for 34,625 genotyped SNPs in this paper) are analyzed in the third stage using FBAT, with correction of multiple testing only for SNPs that enter the third stage. Simulations are performed to evaluate type I error and power of the proposed method compared to LME adjusting for 10 principal components (PC) of the genotype data. We also apply the three-stage approach to the GWA analyses of uric acid in Framingham Heart Study's SNP Health Association Resource (SHARe) project. Results: Our simulations show that whether or not population admixture is present, the three-stage approach has no inflated type I error. In terms of power, using LME adjusting PC is only slightly more powerful than the three-stage approach. When applied to the GWA analyses of uric acid in the SHARe project of FHS, the three-stage approach successfully identified and confirmed three SNPs previously reported as genome-wide significant signals. Conclusions: For GWA analyses of quantitative traits with family data, our three-stage approach provides another appealing solution to population admixture, in addition to LME adjusting for genetic PC.

Background: Ion channel transient receptor potential membrane melastatin 6 and 7 (TRPM6 and TRPM7) play a central role in magnesium homeostasis, which is critical for maintaining glucose and insulin metabolism. However, it is unclear whether common genetic variation in TRPM6 and TRPM7 contributes to risk of type 2 diabetes. Methods: We conducted a nested case-control study in the Women's Health Study. During a median of 10 years of follow-up, 359 incident diabetes cases were diagnosed and matched by age and ethnicity with 359 controls. We analyzed 20 haplotype-tagging single nucleotide polymorphisms (SNPs) in TRPM6 and 5 common SNPs in TRPM7 for their association with diabetes risk. Results: Overall, there was no robust and significant association between any single SNP and diabetes risk. Neither was there any evidence of association between common TRPM6 and TRPM7 haplotypes and diabetes risk. Our haplotype analyses suggested a significant risk of type 2 diabetes among carriers of both the rare alleles from two non-synomous SNPs in TRPM6 (Val1393Ile in exon 26 [rs3750425] and Lys1584Glu in exon 27 [rs2274924]) when their magnesium intake was lower than 250 mg per day. Compared with non-carriers, women who were carriers of the haplotype 1393Ile-1584Glu had an increased risk of type 2 diabetes (OR, 4.92, 95% CI, 1.05–23.0) only when they had low magnesium intake (less than 250 mg/day). Conclusion: Our results provide suggestive evidence that two common non-synonymous TRPM6 coding region variants, Ile1393Val and Lys1584Glu polymorphisms, might confer susceptibility to type 2 diabetes in women with low magnesium intake. Further replication in large-scale studies is warranted.

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

Osteoporosis is a complex disorder and commonly leads to fractures in elderly persons. Genome-wide association studies (GWAS) have become an unbiased approach to identify variations in the genome that potentially affect health. However, the genetic variants identified so far only explain a small proportion of the heritability for complex traits. Due to the modest genetic effect size and inadequate power, true association signals may not be revealed based on a stringent genome-wide significance threshold. Here, we take advantage of SNP and transcript arrays and integrate GWAS and expression signature profiling relevant to the skeletal system in cellular and animal models to prioritize the discovery of novel candidate genes for osteoporosis-related traits, including bone mineral density (BMD) at the lumbar spine (LS) and femoral neck (FN), as well as geometric indices of the hip (femoral neck-shaft angle, NSA; femoral neck length, NL; and narrow-neck width, NW). A two-stage meta-analysis of GWAS from 7,633 Caucasian women and 3,657 men, revealed three novel loci associated with osteoporosis-related traits, including chromosome 1p13.2 (RAP1A, p = 3.661028), 2q11.2 (TBC1D8), and 18q11.2 (OSBPL1A), and confirmed a previously reported region near TNFRSF11B/OPG gene. We also prioritized 16 suggestive genome-wide significant candidate genes based on their potential involvement in skeletal metabolism. Among them, 3 candidate genes were associated with BMD in women. Notably, 2 out of these 3 genes (GPR177, p = 2.6610213; SOX6, p = 6.4610210) associated with BMD in women have been successfully replicated in a large-scale meta-analysis of BMD, but none of the non-prioritized candidates (associated with BMD) did. Our results support the concept of our prioritization strategy. In the absence of direct biological support for identified genes, we highlighted the efficiency of subsequent functional characterization using publicly available expression profiling relevant to the skeletal system in cellular or whole animal models to prioritize candidate genes for further functional validation.

Background: The association between ischemic stroke and 2 single nucleotide polymorphisms (SNPs) on chromosome 12p13, rs12425791 and rs11833579 appears inconsistent across different samples. These SNPs are close to the ninjurin2 gene which may alter the risk of stroke by affecting brain response to ischemic injury. The purpose of this study was to investigate the association between these two SNPs and ischemic stroke risk, as well as prognostic outcomes in a Taiwanese sample. Methods We examined the relations of these two SNPs to the odds of new-onset ischemic stroke, ischemic stroke subtypes, and to the one year risk of stroke-related death or recurrent stroke following initial stroke in a case-control study. A total of 765 consecutive patients who had first-ever ischemic stroke were compared to 977 stroke-free, age-matched controls. SNPs were genotyped by Taqman fluorescent allelic discrimination assay. The association between ischemic stroke and SNPs were analyzed by multivariate logistic regression. Cox proportional hazard model was used to assess the effect of individual SNPs on stroke-related mortality or recurrent stroke. Results: There was no significant association between SNP rs12425791 and rs11833579 and ischemic stroke after multiple testing corrections. However, the marginal significant association was observed between SNP rs12425791 and large artery atherosclerosis under recessive model (OR, 2.30; 95%CI, 1.22-4.34; q-value = 0.062). Among the 765 ischemic stroke patients, 59 died or developed a recurrent stroke. After adjustment for age, sex, vascular risk factors and baseline stroke severity, Cox proportional hazard analysis indicated that the hazard ratios were 2.76 (95%CI, 1.34-5.68; q-value, 0.02) and 2.15 (95%CI, 1.15-4.02; q-value, 0.03) for individuals with homozygous variant allele of rs12425791 and rs11833579, respectively. Conclusions: This is a precedent study that found genetic variants of rs12425791 and rs11833579 on chromosome 12p13 are independent predictors of stroke-related mortality or stroke recurrence in patients with incident ischemic stroke in Taiwan. Further study is needed to explore the details of the physiological function and the molecular mechanisms underlying the association of this genetic locus with ischemic stroke.

The length of female reproductive lifespan is associated with multiple adverse outcomes, including breast cancer, cardiovascular disease and infertility. The biological processes that govern the timing of the beginning and end of reproductive life are not well understood. Genetic variants are known to contribute to ∼50% of the variation in both age at menarche and menopause, but to date the known genes explain <15% of the genetic component. We have used genome-wide association in a bivariate meta-analysis of both traits to identify genes involved in determining reproductive lifespan. We observed significant genetic correlation between the two traits using genome-wide complex trait analysis. However, we found no robust statistical evidence for individual variants with an effect on both traits. A novel association with age at menopause was detected for a variant rs1800932 in the mismatch repair gene MSH6 (P = 1.9 × 10−9), which was also associated with altered expression levels of MSH6 mRNA in multiple tissues. This study contributes to the growing evidence that DNA repair processes play a key role in ovarian ageing and could be an important therapeutic target for infertility.

Osteoporosis is a systemic skeletal disease characterized by a decreased bone mineral density that results in an increased risk of fragility fractures. Previous studies indicated that genetic factors are involved in the pathogenesis of osteoporosis. Polymorphisms of the FONG (FTCDNL1) gene (rs7605378) were reported to be associated with the risk of osteoporosis in a Japanese population. To assess whether polymorphisms of the FTCDNL1 gene contribute to the susceptibility and severity of osteoporosis in a Taiwanese population, 326 osteoporosis patients and 595 controls of a Taiwanese population were included in this study. Our results indicated that rs10203122 was significantly associated with osteoporosis susceptibility among female. Our findings provide evidence that rs10203122 in FTCDNL1 is associated with a susceptibility to osteoporosis.

Background: Genome-wide association studies (GWAS) have identified multiple genetic loci for C-reactive protein (CRP) and lipids, of which some overlap. We aimed to identify genetic pleiotropy among CRP and lipids in order to better understand the shared biology of chronic inflammation and lipid metabolism. Results: In a bivariate GWAS, we combined summary statistics of published GWAS on CRP (n = 66,185) and lipids, including LDL-cholesterol, HDL-cholesterol, triglycerides, and total cholesterol (n = 100,184), using an empirical weighted linear-combined test statistic. We sought replication for novel CRP associations in an independent sample of 17,743 genotyped individuals, and performed in silico replication of novel lipid variants in 93,982 individuals. Fifty potentially pleiotropic SNPs were identified among CRP and lipids: 21 for LDL-cholesterol and CRP, 20 for HDL-cholesterol and CRP, 21 for triglycerides, and CRP and 20 for total cholesterol and CRP. We identified and significantly replicated three novel SNPs for CRP in or near CTSB/FDFT1 (rs10435719, Preplication: 2.6 × 10−5), STAG1/PCCB (rs7621025, Preplication: 1.4 × 10−3) and FTO (rs1558902, Preplication: 2.7 × 10−5). Seven pleiotropic lipid loci were replicated in the independent set of MetaboChip samples of the Global Lipids Genetics Consortium. Annotating the effect of replicated CRP SNPs to the expression of nearby genes, we observed an effect of rs10435719 on gene expression of FDFT1, and an effect of rs7621025 on PCCB. Conclusions: Our large scale combined GWAS analysis identified numerous pleiotropic loci for CRP and lipids providing further insight in the genetic interrelation between lipids and inflammation. In addition, we provide evidence for FDFT1, PCCB and FTO to be associated with CRP levels. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2712-4) contains supplementary material, which is available to authorized users.

Background: Molds can cause respiratory symptoms and asthma. We sought to use isolated peripheral blood mononuclear cells (PBMCs) to understand changes in cytokine and chemokine levels in response to mold and mycotoxin exposures and to link these levels with respiratory symptoms in humans. We did this by utilizing an ex vivo assay approach to differentiate mold-exposed patients and unexposed controls. While circulating plasma chemokine and cytokine levels from these two groups might be similar, we hypothesized that by challenging their isolated white blood cells with mold or mold extracts, we would see a differential chemokine and cytokine release. Methods and Findings: Peripheral blood mononuclear cells (PBMCs) were isolated from blood from 33 patients with a history of mold exposures and from 17 controls. Cultured PBMCs were incubated with the most prominent Stachybotrys chartarum mycotoxin, satratoxin G, or with aqueous mold extract, ionomycin, or media, each with or without PMA. Additional PBMCs were exposed to spores of Aspergillus niger, Cladosporium herbarum and Penicillium chrysogenum. After 18 hours, cytokines and chemokines released into the culture medium were measured by multiplex assay. Clinical histories, physical examinations and pulmonary function tests were also conducted. After ex vivo PBMC exposures to molds or mycotoxins, the chemokine and cytokine profiles from patients with a history of mold exposure were significantly different from those of unexposed controls. In contrast, biomarker profiles from cells exposed to media alone showed no difference between the patients and controls. Conclusions: These findings demonstrate that chronic mold exposures induced changes in inflammatory and immune system responses to specific mold and mycotoxin challenges. These responses can differentiate mold-exposed patients from unexposed controls. This strategy may be a powerful approach to document immune system responsiveness to molds and other inflammation-inducing environmental agents.

Imputation using the 1000 Genomes haplotype reference panel has been widely adapted to estimate genotypes in genome wide association studies. To evaluate imputation quality with a relatively larger reference panel and a reference panel composed of different ethnic populations, we conducted imputations in the Framingham Heart Study and the North Chinese Study using a combined reference panel from the 1000 Genomes (N = 1,092) and UK10K (N = 3,781) projects. For rare variants with 0.01% < MAF ≤ 0.5%, imputation in the Framingham Heart Study with the combined reference panel increased well-imputed genotypes (with imputation quality score ≥0.4) from 62.9% to 76.1% when compared to imputation with the 1000 Genomes. For the North Chinese samples, imputation of rare variants with 0.01% < MAF ≤ 0.5% with the combined reference panel increased well-imputed genotypes by from 49.8% to 61.8%. The predominant European ancestry of the UK10K and the combined reference panels may explain why there was less of an increase in imputation success in the North Chinese samples. Our results underscore the importance and potential of larger reference panels to impute rare variants, while recognizing that increasing ethnic specific variants in reference panels may result in better imputation for genotypes in some ethnic groups.