MGAS: a powerful tool for multivariate gene-based genome-wide association analysis.

MOTIVATION: Standard genome-wide association studies, testing the association between one phenotype and a large number of single nucleotide polymorphisms (SNPs), are limited in two ways: (i) traits are often multivariate, and analysis of composite scores entails loss in statistical power and (ii) gene-based analyses may be preferred, e.g. to decrease the multiple testing problem. RESULTS: Here we present a new method, multivariate gene-based association test by extended Simes procedure (MGAS), that allows gene-based testing of multivariate phenotypes in unrelated individuals. Through extensive simulation, we show that under most trait-generating genotype-phenotype models MGAS has superior statistical power to detect associated genes compared with gene-based analyses of univariate phenotypic composite scores (i.e. GATES, multiple regression), and multivariate analysis of variance (MANOVA). Re-analysis of metabolic data revealed 32 False Discovery Rate controlled genome-wide significant genes, and 12 regions harboring multiple genes; of these 44 regions, 30 were not reported in the original analysis. CONCLUSION: MGAS allows researchers to conduct their multivariate gene-based analyses efficiently, and without the loss of power that is often associated with an incorrectly specified genotype-phenotype models. AVAILABILITY AND IMPLEMENTATION: MGAS is freely available in KGG v3.0 (http://statgenpro.psychiatry.hku.hk/limx/kgg/download.php). Access to the metabolic dataset can be requested at dbGaP (https://dbgap.ncbi.nlm.nih.gov/). The R-simulation code is available from http://ctglab.nl/people/sophie_van_der_sluis. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multivariate Gene-Based Association Test on Family Data in MGAS.

In analyses of unrelated individuals, the program multivariate gene-based association test by extended Simes (MGAS), which facilitates multivariate gene-based association testing, was shown to have correct Type I error rate and superior statistical power compared to other multivariate gene-based approaches. Here we show, through simulation, that MGAS can also be applied to data including genetically related subjects (e.g., family data), by using p value information obtained in Plink or in generalized estimating equations (with the 'exchangeable' working correlation matrix), both of which account for the family structure on a univariate single nucleotide polymorphism-based level by applying a sandwich correction of standard errors. We show that when applied to family-data, MGAS has correct Type I error rate, and given the details of the simulation setup, adequate power. Application of MGAS to seven eye measurement phenotypes showed statistically significant association with two genes that were not discovered in previous univariate analyses of a composite score. We conclude that MGAS is a useful and convenient tool for multivariate gene-based genome-wide association analysis in both unrelated and related individuals.

Predicting mendelian disease-causing non-synonymous single nucleotide variants in exome sequencing studies.

Exome sequencing is becoming a standard tool for mapping Mendelian disease-causing (or pathogenic) non-synonymous single nucleotide variants (nsSNVs). Minor allele frequency (MAF) filtering approach and functional prediction methods are commonly used to identify candidate pathogenic mutations in these studies. Combining multiple functional prediction methods may increase accuracy in prediction. Here, we propose to use a logit model to combine multiple prediction methods and compute an unbiased probability of a rare variant being pathogenic. Also, for the first time we assess the predictive power of seven prediction methods (including SIFT, PolyPhen2, CONDEL, and logit) in predicting pathogenic nsSNVs from other rare variants, which reflects the situation after MAF filtering is done in exome-sequencing studies. We found that a logit model combining all or some original prediction methods outperforms other methods examined, but is unable to discriminate between autosomal dominant and autosomal recessive disease mutations. Finally, based on the predictions of the logit model, we estimate that an individual has around 5% of rare nsSNVs that are pathogenic and carries ~22 pathogenic derived alleles at least, which if made homozygous by consanguineous marriages may lead to recessive diseases.

HYST: a hybrid set-based test for genome-wide association studies, with application to protein-protein interaction-based association analysis.

The extended Simes' test (known as GATES) and scaled chi-square test were proposed to combine a set of dependent genome-wide association signals at multiple single-nucleotide polymorphisms (SNPs) for assessing the overall significance of association at the gene or pathway levels. The two tests use different strategies to combine association p values and can outperform each other when the number of and linkage disequilibrium between SNPs vary. In this paper, we introduce a hybrid set-based test (HYST) combining the two tests for genome-wide association studies (GWASs). We describe how HYST can be used to evaluate statistical significance for association at the protein-protein interaction (PPI) level in order to increase power for detecting disease-susceptibility genes of moderate effect size. Computer simulations demonstrated that HYST had a reasonable type 1 error rate and was generally more powerful than its parents and other alternative tests to detect a PPI pair where both genes are associated with the disease of interest. We applied the method to three complex disease GWAS data sets in the public domain; the method detected a number of highly connected significant PPI pairs involving multiple confirmed disease-susceptibility genes not found in the SNP- and gene-based association analyses. These results indicate that HYST can be effectively used to examine a collection of predefined SNP sets based on prior biological knowledge for revealing additional disease-predisposing genes of modest effects in GWASs.

GATES: a rapid and powerful gene-based association test using extended Simes procedure.

The gene has been proposed as an attractive unit of analysis for association studies, but a simple yet valid, powerful, and sufficiently fast method of evaluating the statistical significance of all genes in large, genome-wide datasets has been lacking. Here we propose the use of an extended Simes test that integrates functional information and association evidence to combine the p values of the single nucleotide polymorphisms within a gene to obtain an overall p value for the association of the entire gene. Our computer simulations demonstrate that this test is more powerful than the SNP-based test, offers effective control of the type 1 error rate regardless of gene size and linkage-disequilibrium pattern among markers, and does not need permutation or simulation to evaluate empirical significance. Its statistical power in simulated data is at least comparable, and often superior, to that of several alternative gene-based tests. When applied to real genome-wide association study (GWAS) datasets on Crohn disease, the test detected more significant genes than SNP-based tests and alternative gene-based tests. The proposed test, implemented in an open-source package, has the potential to identify additional novel disease-susceptibility genes for complex diseases from large GWAS datasets.

A comprehensive framework for prioritizing variants in exome sequencing studies of Mendelian diseases.

Exome sequencing strategy is promising for finding novel mutations of human monogenic disorders. However, pinpointing the casual mutation in a small number of samples is still a big challenge. Here, we propose a three-level filtration and prioritization framework to identify the casual mutation(s) in exome sequencing studies. This efficient and comprehensive framework successfully narrowed down whole exome variants to very small numbers of candidate variants in the proof-of-concept examples. The proposed framework, implemented in a user-friendly software package, named KGGSeq (http://statgenpro.psychiatry.hku.hk/kggseq), will play a very useful role in exome sequencing-based discovery of human Mendelian disease genes.

Evaluating the effective numbers of independent tests and significant p-value thresholds in commercial genotyping arrays and public imputation reference datasets.

Current genome-wide association studies (GWAS) use commercial genotyping microarrays that can assay over a million single nucleotide polymorphisms (SNPs). The number of SNPs is further boosted by advanced statistical genotype-imputation algorithms and large SNP databases for reference human populations. The testing of a huge number of SNPs needs to be taken into account in the interpretation of statistical significance in such genome-wide studies, but this is complicated by the non-independence of SNPs because of linkage disequilibrium (LD). Several previous groups have proposed the use of the effective number of independent markers (M(e)) for the adjustment of multiple testing, but current methods of calculation for M(e) are limited in accuracy or computational speed. Here, we report a more robust and fast method to calculate M(e). Applying this efficient method [implemented in a free software tool named Genetic type 1 error calculator (GEC)], we systematically examined the M(e), and the corresponding p-value thresholds required to control the genome-wide type 1 error rate at 0.05, for 13 Illumina or Affymetrix genotyping arrays, as well as for HapMap Project and 1000 Genomes Project datasets which are widely used in genotype imputation as reference panels. Our results suggested the use of a p-value threshold of ~10(-7) as the criterion for genome-wide significance for early commercial genotyping arrays, but slightly more stringent p-value thresholds ~5 × 10(-8) for current or merged commercial genotyping arrays, ~10(-8) for all common SNPs in the 1000 Genomes Project dataset and ~5 × 10(-8) for the common SNPs only within genes.

Complete genome sequence of the bacterium Methylovorus sp. strain MP688, a high-level producer of pyrroloquinolone quinone.

Methylotrophic bacteria are widespread microbes which can use one carbon compound as their only carbon and energy sources. Here we report the finished, annotated genome sequence of the methylotrophic bacterium Methylovorus sp. strain MP688, which was isolated from soil for high-level production of pyrroloquinolone quinone (PQQ) in our lab.

An NT-3-releasing bioscaffold supports the formation of TrkC-modified neural stem cell-derived neural network tissue with efficacy in repairing spinal cord injury.

The mechanism underlying neurogenesis during embryonic spinal cord development involves a specific ligand/receptor interaction, which may be help guide neuroengineering to boost stem cell-based neural regeneration for the structural and functional repair of spinal cord injury. Herein, we hypothesized that supplying spinal cord defects with an exogenous neural network in the NT-3/fibroin-coated gelatin sponge (NF-GS) scaffold might improve tissue repair efficacy. To test this, we engineered tropomyosin receptor kinase C (TrkC)-modified neural stem cell (NSC)-derived neural network tissue with robust viability within an NF-GS scaffold. When NSCs were genetically modified to overexpress TrkC, the NT-3 receptor, a functional neuronal population dominated the neural network tissue. The pro-regenerative niche allowed the long-term survival and phenotypic maintenance of the donor neural network tissue for up to 8 weeks in the injured spinal cord. Additionally, host nerve fibers regenerated into the graft, making synaptic connections with the donor neurons. Accordingly, motor function recovery was significantly improved in rats with spinal cord injury (SCI) that received TrkC-modified NSC-derived neural network tissue transplantation. Together, the results suggested that transplantation of the neural network tissue formed in the 3D bioactive scaffold may represent a valuable approach to study and develop therapies for SCI.

IGG3: a tool to rapidly integrate large genotype datasets for whole-genome imputation and individual-level meta-analysis.

SUMMARY: There is an urgent and increasing demand for integrating large genotype datasets across genome-wide association studies and HapMap project for whole-genome imputation and individual-level meta-analysis. A new algorithm was developed to efficiently merge raw genotypes across large datasets and implemented in the latest version of IGG, IGG3. In addition, IGG3 can integrate the latest phased and unphased HapMap genotypes and can flexibly generate complete sets of input files for six popular genotype imputation tools. We demonstrated the efficiency of IGG3 by simulation tests, which could rapidly merge genotypes in tens of thousands of large genotype chips (e.g. Affymetrix Genome-Wide Human SNP Array 6.0 and Illumina Human1m-duo) and in HapMap III project on an ordinary desktop computer. AVAILABILITY: (http://bioinfo.hku.hk/iggweb) (version 3.0).

Power of transmission/disequilibrium tests in admixed populations.

The power of transmission/disequilibrium tests (TDTs) for detecting disease susceptibility loci is expected to be influenced by population admixture through its impact on the degree of linkage disequilibrium (LD) between the genetic marker and the DSL. However, few studies have been done to systematically examine this behavior of the TDTs in admixed populations. In the present study, extensive computer simulations were conducted to explore how population admixture affects the power of TDTs. It was found that (1) in newly admixed populations, the LD due to admixture makes no contribution to the power of TDTs, and it is the averaged background LD in the parental populations that determines the power of TDTs; but (2) after random mating between the admixed populations, the LD due to admixture becomes effective in increasing or decreasing the power of the tests, and (3) incomplete random mating can prolong the time for the LD due to admixture to become effective. This study clarifies the potential influence of population admixture on the performance of TDTs.

Exclusion mapping of chromosomes 1, 4, 6 and 14 with bone mineral density in 79 Caucasian pedigrees.

Low bone mineral density (BMD) is a major determinant of osteoporosis and is under strong genetic control. A large number of linkage and association studies for BMD variation have been conducted, with the results being largely inconsistent. Linkage exclusion analysis is a useful tool for gene mapping but has never been used on BMD. In the present study, we conducted a linkage exclusion mapping for BMD variation on chromosomes 1, 4, 6 and 17 in 79 Caucasian pedigrees. For hip BMD variation, several genomic regions were excluded for effect sizes of 10% or greater, including regions of 61-77 cM at 1p35-p34, 167-196 cM at 1q21-q23 and 261-291 cM at 1q42-q44; 85-112 cM at 4q21-q25 and 146-150 cM at 4q31; and 77-85 cM at 6p12-q13. For spine BMD, we were able to exclude the regions of 168-189 cM at 1q21-q23, 92-94 cM at 4q21 and 106-107 cM at 4q24 and 56-103 cM at 17q12-q25, as having effect sizes of 10% or greater. These results suggest that a number of candidate genes located in the excluded regions, such as interleukin 6 receptor (IL6R) gene, type I collagen alpha 1 (COL1A1) gene and bone morphogenetic protein-3 (BMP3) gene are unlikely to have a substantial effect on BMD variation in this Caucasian population. Along with previous studies searching for genes underlying BMD variation, the current study has further delineated the genetic basis of BMD variation and provided valuable information for future genetic studies.

FAPI: Fast and accurate P-value Imputation for genome-wide association study.

Imputing individual-level genotypes (or genotype imputation) is now a standard procedure in genome-wide association studies (GWAS) to examine disease associations at untyped common genetic variants. Meta-analysis of publicly available GWAS summary statistics can allow more disease-associated loci to be discovered, but these data are usually provided for various variant sets. Thus imputing these summary statistics of different variant sets into a common reference panel for meta-analyses is impossible using traditional genotype imputation methods. Here we develop a fast and accurate P-value imputation (FAPI) method that utilizes summary statistics of common variants only. Its computational cost is linear with the number of untyped variants and has similar accuracy compared with IMPUTE2 with prephasing, one of the leading methods in genotype imputation. In addition, based on the FAPI idea, we develop a metric to detect abnormal association at a variant and showed that it had a significantly greater power compared with LD-PAC, a method that quantifies the evidence of spurious associations based on likelihood ratio. Our method is implemented in a user-friendly software tool, which is available at http://statgenpro.psychiatry.hku.hk/fapi.

No major effect of the insulin-like growth factor I gene on bone mineral density in premenopausal Chinese women.

Osteoporosis is a major public health problem, mainly characterized by low bone mineral density (BMD). BMD is a complex trait that is determined by multiple genes. Insulin-like growth factor I (IGF-I) is an important growth factor of bone and thus IGF-I gene has been considered as an attractive candidate gene for osteoporosis. A few studies on the relationship between variants of the IGF-I gene and BMD variation, via traditional association and/or linkage methods, have yielded conflicting results. In this study, we simultaneously tested association and/or linkage of a cytosine-adenine (CA) repeat polymorphism at 1 kb upstream of the transcription initiation site of the IGF-I gene with BMD variation in a large cohort of premenopausal Chinese women. A total of 1263 subjects from 402 Chinese nuclear families were examined. Each family consists of both parents and at least one daughter aged between 20 and 45 years. BMDs (g/cm(2)) at the lumbar spine and hip were measured using dual-energy X-ray absorptiometry (DXA). Applying the QTDT (quantitative transmission disequilibrium tests) progam, we did not find significant evidence of association or linkage between the CA repeat polymorphism of the IGF-I gene and BMD variation at any skeletal site. Our data do not support the IGF-I gene having major effect on BMD variation in premenopausal Chinese women.

Potential effect of inter-genic action on peak bone mass (PBM) in Chinese females.

Peak bone mass (PBM) is a complex trait, determined by both genetic and environmental factors and also their interactions. Vitamin D receptor (VDR) estrogen receptor alpha (ERalpha), interleukin 6 (IL6), parathyroid hormone (PTH), collagen type I alpha 2 (COL1A2), bone Gla protein (BGP), alpha2-HS glycoprotein (AHSG) are among the important candidate genes of bone metabolism. The study aims to detect significant effect of potential inter-genic action underlying PBM in Chinese females. 361 unrelated healthy premenopausal Chinese females (aged 20 -44 years) with Han ethnicity were recruited from the Shanghai city in China. Bone mineral density (BMD) at the hip and the lumbar spine (L1-4) was measured using a Hologic QDR 2000 + dual-energy X-ray absorptiometry (DXA) scanner. Eight polymorphisms among the seven genes were genotyped, i. e. Apa I in VDR, Pvu II and Xba I in ERa (ERX and ERP, respectively), BsrB I in IL6, BstB I in PTH, Msp I in COL1A2, Hind III in BGP, and Sac I in AHSG, using PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) methods. Two-way analysis of variance (ANOVA) showed significant effects of IL6 x ERP interaction on PBM at the total hip (P = 0.019), intertrochanter (P = 0.016), and femoral neck (P =0. 019). The BMD difference between GGPp carriers and GGpp subjects (at these two loci) amounted to 18.0%, 19.5%, and 14.8% at the hip,intertrochanter,and femoral neck,respectively. The potential interaction effect of AHSG x IL6 was observed on femoral neck PBM (P = 0.046). GGSS individuals (at these two loci) had, on average, 18.8% higher femoral neck BMD than those subjects with GGSs genotype. The population-level statistical analysis indicates that IL6 x ERP and AHSG x IL6 have significant inter-genic effect on the genetic determination of PBM in Chinese females.

Alpha2-HS glycoprotein gene is associated with bone size at the hip in Chinese.

Bone size is an important risk factor, independent of bone mineral density (BMD), for osteoporotic fracture. Bone size has a high heritability. A better understanding of genetic factors regulating bone size will have important clinical implications. In this study, we explored the relationship between the alpha2-HS glycoprotein (AHSG) gene and bone size variation at the spine and hip in a Chinese population. The study sample comprised 1 260 subjects from 401 Chinese nuclear families (each including both parents and at least one female child). The Sac / polymorphism inside the exon 7 of the AHSG gene was genotyped and analyzed. This variant represents a nucleotide substitution of C to G at amino acid position 238 resulting in a translation polymorphism of threonine to serine and thus making a potential impact on gene function. We assessed population stratification but did not find significant evidence at any skeletal sites. We found significant association between the AHSG Sac / polymorphism and bone size at the intertrochanteric region (P = 0.019) and the total hip (P = 0.035). The polymorphisms explained 3.74% and 3.16% variations in bone size at the intertrochanteric region and total hip respectively. No significant evidence of linkage was detected, largely due to the limited number of sibpairs in this data set and less informative marker (AHSG Sac / polymorphism) (compared with microsatellite markers) for linkage analysis. Our results suggested that the AHSG gene may contribute to bone size variation at the hip in this Chinese population.

SNPTracker: A Swift Tool for Comprehensive Tracking and Unifying dbSNP rs IDs and Genomic Coordinates of Massive Sequence Variants.

The reference single nucleotide polymorphism (rs) ID in dbSNP (http://www.ncbi.nlm.nih.gov/SNP/) is a key resource identifier, which is widely used in human genetics and genomics studies. However, its application is often complicated by the varied IDs of different versions. Here, we developed a user-friendly tool, SNPTracker, for comprehensively tracking and unifying the rs IDs and genomic coordinates of massive sequence variants at a time. It worked perfectly, and had much higher accuracy and capacity than two alternative utilities in our proof-of-principle examples. SNPTracker will greatly facilitate genetic data exchange and integration in the postgenome-wide association study era.

Estrogen receptor alpha gene polymorphisms and peak bone density in Chinese nuclear families.

PBD is an important determinant of osteoporotic fractures. Few studies were performed to search for genes underlying PBD variation in Chinese populations. We tested linkage and/or association of the estrogen receptor alpha gene polymorphism with PBD in 401 Chinese nuclear families. This study suggests the ER-alpha gene may have some minor effects on PBM variation in the Chinese population. Low peak bone density (PBD) in adulthood is an important determinant of osteoporotic fractures in the elderly. PBD variation is mainly regulated by genetic factors. Extensive molecular genetics studies have been performed to search for genes underlying PBD variation, largely in whites. Few studies were performed in Chinese populations. In this study, we simultaneously test linkage and/or association of the estrogen receptor alpha (ER-alpha) gene polymorphism with PBD in 401 Chinese nuclear families (both parents plus their female children) of 1260 subjects, with the 458 children generally between 20 and 40 years of age. All the subjects were genotyped by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) at polymorphic PvuII and XbaI sites inside the ER-alpha gene. Bone mineral density was measured at the lumbar spine (L1-L4) and hip (femoral neck, trochanter, and intertrochanteric region). Raw bone mineral density values were adjusted by age, height, and weight as covariates. We detected marginally significant results for within-family association (transmission disequilibrium; p = 0.054) between the spine bone mineral density variation and the ER-alpha XbaI genotypes. For the hip bone mineral density variation, significant (p < 0.05) linkage results were generally found for the two intragenic markers. Analyses of the haplotypes defined by the two markers confer further evidence for linkage of the ER-alpha with the hip PBD variation. In conclusion, this study suggests that the ER-alpha gene may have minor effects on PBD variation in our Chinese population.