BACKGROUND: DNA methylation is an important heritable epigenetic mark that plays a crucial role in transcriptional regulation and the pathogenesis of various human disorders. The commonly used DNA methylation measurement approaches, e.g., Illumina Infinium HumanMethylation-27 and -450 BeadChip arrays (27 K and 450 K arrays) and reduced representation bisulfite sequencing (RRBS), only cover a small proportion of the total CpG sites in the human genome, which considerably limited the scope of the DNA methylation analysis in those studies. RESULTS: We proposed a new computational strategy to impute the methylation value at the unmeasured CpG sites using the mixture of regression model (MRM) of radial basis functions, integrating information of neighboring CpGs and the similarities in local methylation patterns across subjects and across multiple genomic regions. Our method achieved a better imputation accuracy over a set of competing methods on both simulated and empirical data, particularly when the missing rate is high. By applying MRM to an RRBS dataset from subjects with low versus high bone mineral density (BMD), we recovered methylation values of ~ 300 K CpGs in the promoter regions of chromosome 17 and identified some novel differentially methylated CpGs that are significantly associated with BMD. CONCLUSIONS: Our method is well applicable to the numerous methylation studies. By expanding the coverage of the methylation dataset to unmeasured sites, it can significantly enhance the discovery of novel differential methylation signals and thus reveal the mechanisms underlying various human disorders/traits.

Genetic risk score (GRS, also known as polygenic risk score) analysis is an increasingly popular method for exploring genetic architectures and relationships of complex diseases. However, complex diseases are usually measured by multiple correlated phenotypes. Analyzing each disease phenotype individually is likely to reduce statistical power due to multiple testing correction. In order to conquer the disadvantage, we proposed a principal component analysis (PCA)-based GRS analysis approach. Extensive simulation studies were conducted to compare the performance of PCA-based GRS analysis and traditional GRS analysis approach. Simulation results observed significantly improved performance of PCA-based GRS analysis compared to traditional GRS analysis under various scenarios. For the sake of verification, we also applied both PCA-based GRS analysis and traditional GRS analysis to a real Caucasian genome-wide association study (GWAS) data of bone geometry. Real data analysis results further confirmed the improved performance of PCA-based GRS analysis. Given that GWAS have flourished in the past decades, our approach may help researchers to explore the genetic architectures and relationships of complex diseases or traits.

Motivation: Extreme phenotype sampling (EPS) is a broadly-used design to identify candidate genetic factors contributing to the variation of quantitative traits. By enriching the signals in extreme phenotypic samples, EPS can boost the association power compared to random sampling. Most existing statistical methods for EPS examine the genetic factors individually, despite many quantitative traits have multiple genetic factors underlying their variation. It is desirable to model the joint effects of genetic factors, which may increase the power and identify novel quantitative trait loci under EPS. The joint analysis of genetic data in high-dimensional situations requires specialized techniques, e.g. the least absolute shrinkage and selection operator (LASSO). Although there are extensive research and application related to LASSO, the statistical inference and testing for the sparse model under EPS remain unknown. Results: We propose a novel sparse model (EPS-LASSO) with hypothesis test for high-dimensional regression under EPS based on a decorrelated score function. The comprehensive simulation shows EPS-LASSO outperforms existing methods with stable type I error and FDR control. EPS-LASSO can provide a consistent power for both low- and high-dimensional situations compared with the other methods dealing with high-dimensional situations. The power of EPS-LASSO is close to other low-dimensional methods when the causal effect sizes are small and is superior when the effects are large. Applying EPS-LASSO to a transcriptome-wide gene expression study for obesity reveals 10 significant body mass index associated genes. Our results indicate that EPS-LASSO is an effective method for EPS data analysis, which can account for correlated predictors. Availability and implementation: The source code is available at https://github.com/xu1912/EPSLASSO. Supplementary information: Supplementary data are available at Bioinformatics online.

MOTIVATION: Pathway association analysis has made great achievements in elucidating the genetic basis of human complex diseases. However, current pathway association analysis approaches fail to consider tissue-specificity. RESULTS: We developed a tissue-specific pathway interaction enrichment analysis algorithm (TPIEA). TPIEA was applied to two large Caucasian and Chinese genome-wide association study summary datasets of bone mineral density (BMD). TPIEA identified several significant pathways for BMD [false discovery rate (FDR) < 0.05], such as KEGG FOCAL ADHESION and KEGG AXON GUIDANCE, which had been demonstrated to be involved in the development of osteoporosis. We also compared the performance of TPIEA and classical pathway enrichment analysis, and TPIEA presented improved performance in recognizing disease relevant pathways. TPIEA may help to fill the gap of classic pathway association analysis approaches by considering tissue specificity. AVAILABILITY AND IMPLEMENTATION: The online web tool of TPIEA is available at https://sourceforge.net/projects/tpieav1/files CONTACT: [email protected] information: Supplementary data are available at Bioinformatics online.

MOTIVATION: Several new de novo assembly tools have been developed recently to assemble short sequencing reads generated by next-generation sequencing platforms. However, the performance of these tools under various conditions has not been fully investigated, and sufficient information is not currently available for informed decisions to be made regarding the tool that would be most likely to produce the best performance under a specific set of conditions. RESULTS: We studied and compared the performance of commonly used de novo assembly tools specifically designed for next-generation sequencing data, including SSAKE, VCAKE, Euler-sr, Edena, Velvet, ABySS and SOAPdenovo. Tools were compared using several performance criteria, including N50 length, sequence coverage and assembly accuracy. Various properties of read data, including single-end/paired-end, sequence GC content, depth of coverage and base calling error rates, were investigated for their effects on the performance of different assembly tools. We also compared the computation time and memory usage of these seven tools. Based on the results of our comparison, the relative performance of individual tools are summarized and tentative guidelines for optimal selection of different assembly tools, under different conditions, are provided.