Software

S-LDXR

(October 2019) The S-LDXR software can be downloaded here. S-LDXR is a method for stratifying squared trans-ethnic genetic correlation across genomic annotations, as described in our manuscript “Population-specific causal disease effect sizes in functionally important regions impacted by selection” (Shi et al. 2019 biorxiv).

PolyFun and PolyLoc

(October 2019) The PolyFun software can be downloaded here. PolyFun is a method that leverages genome-wide functional annotations to improve fine-mapping power, as described in our manuscript “Functionally-informed fine-mapping and polygenic localization of complex trait heritability” (Weissbrod et al. 2019 biorxiv).

(October 2019) The PolyLoc software can be downloaded here. PolyLoc is a method that leverages the results of PolyFun to perform polygenic localization of complex trait heritability, as described in our manuscript “Functionally-informed fine-mapping and polygenic localization of complex trait heritability (Weissbrod et al. 2019 biorxiv).

LT-FH

(July 2019) The LT-FH software can be downloaded here. LT-FH is a method for leveraging family history information to improve association power, as described in our manuscript “Combining case-controls status and family history of disease increases association power” (Hujoel et al. biorxiv)

LDpred-funct

(January 2019) The LDpred-funct software can be downloaded here. LDpred-funct is a method for leveraging functional enrichment to improve polygenic prediction accuracy, as described in our manuscript “Modeling functional enrichment improves polygenic prediction accuracy in UK Biobank and 23andMe data sets” (Marquez-Luna et al. biorxiv).

LD SCORE

(updated October 2018) The ldsc software can be downloaded here.  LD Score regression (Bulik-Sullivan et al. 2015a) is a method for distinguishing confounding from polygenicity in genome-wide association studies.  Stratified LD Score regression (Finucane et al. 2015; functional annotations here) is a method for partitioning heritability by functional category using GWAS summary statistics.  Cross-trait LD Score regression (Bulik-Sullivan et al. 2015b) is a method for estimating genetic correlations using GWAS summary statistics. We have extended stratified LD score regression to gene expression phenotypes (Liu et al. 2017). We have also extended stratified LD score regression to continuous annotations (Gazal et al. 2017). We have also developed an approach that uses stratified LD score regression to identify disease-relevant tissues and cell types with heritability enrichment in specifically expressed genes (Finucane et al. 2018; functional annotations here). We have also extended stratified LD score regression to low-frequency variants (Gazal et al. 2018; functional annotations here). We have also identified conditionally independent signals of disease heritability enrichment for ancient enhancers, enhancers that are conserved across species, ancient promoters, and promoters of loss-of-function intolerant ExAC genes (Hujoel et al. 2019; functional annotations here). We have also identified enriched pathway, network, and pathway+network annotations, concluding that genes with high network connectivity are enriched for disease heritability (Kim et al. 2019; annotations here). We have also identified conditionally independent signals of disease heritability for transposable elements (Hormozdiari et al. 2019; annotations here).

S-LD4M

(September 2018) The S-LD4M software can be downloaded here.  This software implements our Stratified LD 4th moments regression (S-LD4M) method for estimating polygenicity across allele frequencies and functional categories, as described in our manuscript “Polygenicity of complex traits is explained by negative selection” (O’Connor et al. 2019).

FINDOR

(January 2018) The FINDOR software can be downloaded here.  This software implements our Functionally Informed Novel Discovery Of Risk loci (FINDOR) method, as described in our manuscript “Leveraging polygenic functional enrichment to improve GWAS power” (Kichaev et al. 2019).

LCV

(October 2017) The LCV software can be downloaded here.  This software implements our Latent Causal Variable (LCV) causal inference method, as described in our manuscript “Distinguishing genetic correlation from causation across 52 diseases and complex traits” (O’Connor et al. 2018).

Signed LD profile regression

(October 2017) Signed LD profile regression software can be downloaded here.  Signed LD profile regression is a method for identifying genome-wide directional effects of signed functional annotations on diseases and complex traits, as described in our manuscript “Detecting genome-wide directional effects of transcription factor binding on polygenic disease risk” (Reshef et al. 2018).

Molecular QTL annotations

(October 2017) Molecular QTL annotations described in our manuscript “Leveraging molecular QTL to understand the genetic architecture of diseases and complex traits” (Hormozdiari et al. 2018) can be downloaded here.  The molecular QTL annotations include MaxCPP annotations constructed from GTEx eQTL data and MaxCPP annotations constructed from BLUEPRINT eQTL, hQTL (H3K27ac and H3K4me1), sQTL and meQTL data.

BOLT-LMM and BOLT-REML

(September 2017) The BOLT-LMM v2.3 software package (Loh et al. 2018), which includes multi-threaded support for the BGEN v1.2 imputed file format used by UK Biobank, can be downloaded here. Summary association statistics from BOLT-LMM analyses of all N=459K European-ancestry samples in UK Biobank are available here. The BOLT-LMM algorithm (Loh et al. 2015a) rapidly computes statistics for association between phenotype and genotypes using a linear mixed model (LMM). The BOLT-REML algorithm partitions SNP-heritability and estimates genetic correlations using a Monte Carlo algorithm for fast multi-component, multi-trait modeling (Loh et al. 2015b). By default, BOLT-LMM association analysis assumes a Bayesian mixture-of-normals prior for the random effect attributed to SNPs other than the one being tested. This model generalizes the standard “infinitesimal” mixed model used by previous mixed model association methods, providing an opportunity for increased power to detect associations while controlling false positives. Both algorithms are implemented in the BOLT-LMM v2.3 software package; see link for update log.

EIGENSOFT

(June 2017): EIGENSOFT version 7.2.1 is now available for download. The EIGENSOFT package combines functionality from our population genetics methods (Patterson et al. 2006) and our EIGENSTRAT stratification correction method (Price et al. 2006). The EIGENSTRAT method uses principal components analysis to explicitly model ancestry differences between cases and controls along continuous axes of variation; the resulting correction is specific to a candidate marker’s variation in frequency across ancestral populations, minimizing spurious associations while maximizing power to detect true associations. The EIGENSOFT package has a built-in plotting script and supports multiple file formats and quantitative phenotypes.

The latest version of EIGENSOFT (7.2.1) can be downloaded here. Source code, documentation and executables for using EIGENSOFT 6.1.4 on a Linux platform can be downloaded here. New features of EIGENSOFT 6.x include fastmode option which implements a very fast pca approximation (Galinsky et al. 2016a, Galinsky et al. 2016b) and support for multi-threading. EIGENSOFT 6.1.4 includes bug fixes for pcaselection and a better out of memory diagnostic message. Our previous release, version 6.1.3, can be downloaded here.

The EIGENSOFT FAQ (Frequently Asked Questions) is available here.
For further questions about the EIGENSOFT software, please write to Samuela Pollack (spollack@hsph.harvard.edu).

SOFTWARE REGISTRATION: we encourage EIGENSOFT users to register here. Registration is voluntary, but will allow us to send information about software updates.

SNP loadings computed from samples with European ancestry from the GERA cohort can be downloaded here.

LTSOFT

(January 2017): LTSOFT version 4.0 can be downloaded here. Changes to version 4.0 include the addition of LT-Fam to the LTMLM software implementing a multivariate liability threshold mixed linear model for settings with related individuals (Hayeck et al. 2017). Changes to version 3.0 include the addition of LTMLM, a new piece of software implementing a multivariate liability threshold mixed linear model association statistics for additional increase in power in settings of case control diseases (Hayeck et al. 2015).  LTSOFT is a software suite designed to more powerfully leverage clinical-covariates such as age, bmi, smoking status, and gender as well as genetic-covariates such as known associated variants when conducting case-control association studies. Including these covariates in standard regression models is not only suboptimal, but can in many instances reduce power. LTSOFT employs a liability threshold model approach that takes advantage of known epidemiological results to better model the covariates’ relationship to the phenotype of interest (Zaitlen et al. 2012 PLoS Genet and Zaitlen et al. 2012 Bioinformatics).

Eagle

(May 2016) The Eagle v2.0 software (Loh et al. 2016b) estimates haplotype phase either using a phased reference panel or within a genotyped cohort. Eagle2 is now the default phasing method on the Sanger and Michigan imputation servers and uses a new, very fast HMM-based algorithm that improves speed and accuracy over existing methods via two key ideas. : a new data structure based on the positional Burrows-Wheeler transform and a rapid search algorithm that explores only the most relevant paths through the HMM. Compared to the Eagle1 algorithm (Loh et al. 2016a), Eagle2 has similar speed but much greater accuracy at sample sizes <50,000. The Eagle software can be downloaded here.

haploSNP

(July 2015) The haploSNP algorithm (Bhatia et al. biorxiv) constructs a set of haplotype polymorphisms (haploSNPs) from phased genotype data. haploSNPs are haplotypes of adjacent SNPs excluding a subset of masked sites that arise from skipped mismatches. Mismatches are skipped only if they can be potentially explained as mutations on a shared background. This is tested using a 4-gamete test between the haploSNP being extended and the mismatch SNP. If all 4 possible allelic combinations are observed, the mismatch cannot be explained as a mutation on a shared background, and the haploSNP is terminated.
Individuals are considered to carry 0,1, or 2 copies of the haploSNP if none, one or both of their chromosomes matches the haplotype at all unmasked sites. As haploSNPs are biallelic, they can be used in downstream analyses such as heritability estimation and association.
The haploSNP software can be downloaded here.

Efficient PCGC Regression

(June 2015) PCGC regression (Golan et al. 2014 PNAS) is designed to avoid biases in REML estimation of heritability in the context of ascertained case-control studies. We have released an efficient implementation of the PCGC regression method. Subject to the restriction that all GRMs be computed over identical lists of individuals (all *.grm.id files must be identical), this implementation eliminates in-memory storage of N x N matrices by accumulating dot products among regressors on-the-fly (i.e., streaming the GRM inputs), speeds up jackknife computations (storing partition results on the fly), and eliminates storage of “cleaned” GRMs (i.e., with PCs projected out) by projecting PCs on-the-fly. The Efficient PCGC Regression software is used in Loh et al. 2015b and Bhatia et al. biorxiv, and can be downloaded here.

LDpred

(March 2015) The LDpred software can be downloaded here.  LDpred (Vilhjalmsson et al. 2015) is a method for computing polygenic risk scores from summary association statistics while accounting for LD between markers. The method infers the posterior mean causal effect size of each marker using a non-infinitesimal prior distribution on effect sizes and LD information from an external reference panel.

SNPWEIGHTS

(May 2014): SNPweights version 2.1 can be downloaded here.  SNPweights is a software package for inferring genome-wide genetic ancestry using SNP weights precomputed from large external reference panels (Chen et al. 2013 Bioinformatics).  Changes to version 2.0 include new SNP weights for Native American reference samples, a new format for SNP weights files, and new software for users to derive SNP weights using their own reference samples.  Version 2.1 incorporates a bug fix in the inferanc program, which now works with all snpwt files.  SNP weights for European and West African ancestral populations can be downloaded here. SNP weights for European, West African and East Asian ancestral populations can be downloaded here. SNP weights for European, West African, East Asian and Native American ancestral populations can be downloaded here.  SNP weights for NW, SE and AJ ancestral populations of European Americans can be downloaded here.

FUNCTIONAL ANNOTATIONS

(March 2014): Functional annotations of SNPs and regions from our functional heritability paper “Regulatory variants explain much more heritability than coding variants across 11 common diseases” (Gusev et al. 2014) can be downloaded here.

IMPG-SUMMARY

(July 2013):  ImpG-Summary version 1.0 can be downloaded here.  ImpG-Summary is a software package for Gaussian imputation from summary association statistics, as described in our paper “Fast and accurate imputation from summary association statistics” (Pasaniuc et al. 2014).

MIXSCORE

(July 2012): MIXSCORE version 1.3 can be downloaded here. MIXSCORE is a method for combining SNP association and admixture association statistics to increase power in GWAS in admixed populations. For details, see the MIXSCORE paper (Pasaniuc et al. 2011 PLoS Genet, “Enhanced statistical tests for GWAS in admixed populations: assessment using African Americans from CARe and a breast cancer consortium”).

TREESELECT

(April 2012): TreeSelect version 1.1 can be downloaded here. TreeSelect is a software package for inferring natural selection from unusual population differentiation between closely related populations. For details, see our Africa selection paper (Bhatia et al. 2011 AJHG, “Genome-wide comparison of African-ancestry populations from CARe and other cohorts reveals signals of natural selection.”)

HAPMIX

(March 2011): HAPMIX version 1.2 can be downloaded here. Improvements to version 1.2 include an explicit check for discordance between admixed and reference population allele frequencies, and a script to interpolate estimates of local ancestry to a superset of SNPs. HAPMIX is an application for accurately inferring chromosomal segments of distinct continental ancestry in admixed populations, using dense genetic data.For details, see the HAPMIX paper (Price et al. 2009).

GENE EXPRESSION HERITABILITY

(January 2011): Source code and gene-by-gene results from our gene expression heritability paper “Single-tissue and cross-tissue heritability of gene expression via identity-by-descent in related or unrelated individuals” (Price et al. 2011) can be downloaded here.