ATAV: a comprehensive platform for population-scale genomic analyses
BMC Bioinformatics volume 22, Article number: 149 (2021)
A common approach for sequencing studies is to do joint-calling and store variants of all samples in a single file. If new samples are continually added or controls are re-used for several studies, the cost and time required to perform joint-calling for each analysis can become prohibitive.
We present ATAV, an analysis platform for large-scale whole-exome and whole-genome sequencing projects. ATAV stores variant and per site coverage data for all samples in a centralized database, which is efficiently queried by ATAV to support diagnostic analyses for trios and singletons, as well as rare-variant collapsing analyses for finding disease associations in complex diseases. Runtime logs ensure full reproducibility and the modularized ATAV framework makes it extensible to continuous development. Besides helping with the identification of disease-causing variants for a range of diseases, ATAV has also enabled the discovery of disease-genes by rare-variant collapsing on datasets containing more than 20,000 samples. Analyses to date have been performed on data of more than 110,000 individuals demonstrating the scalability of the framework.
To allow users to easily access variant-level data directly from the database, we provide a web-based interface, the ATAV data browser (http://atavdb.org/). Through this browser, summary-level data for more than 40,000 samples can be queried by the general public representing a mix of cases and controls of diverse ancestries. Users have access to phenotype categories of variant carriers, as well as predicted ancestry, gender, and quality metrics. In contrast to many other platforms, the data browser is able to show data of newly-added samples in real-time and therefore evolves rapidly as more and more samples are sequenced.
Through ATAV, users have public access to one of the largest variant databases for patients sequenced at a tertiary care center and can look up any genes or variants of interest. Additionally, since the entire code is freely available on GitHub, ATAV can easily be deployed by other groups that wish to build their own platform, database, and user interface.
Diagnostic and cohort sequencing studies benefit from the analysis of a large number of samples combined with similarly processed controls. A common approach to reach the necessary scale for analysis is to use a joint-calling procedure and store all samples in a single VCF file [1, 2]. While effective in allowing a single analysis of all samples included in the single VCF file, this approach has significant limitations. Perhaps most importantly, this approach is not amenable to ongoing analyses as new samples become available. Moreover, when projects combine multiple cohorts that were not sequenced together and in which controls might be re-used for several studies, the cost and time required to perform joint-calling for each analysis can become prohibitive. In addition to these considerations, typical sequencing file formats (VCF, BAM) place a sizeable overhead in moving these data from physical storage to the compute nodes for dynamic and multi user analysis needs. Furthermore, standard diagnostic and case–control studies leverage a range of filtering parameters, including variant calling (genotype quality, read coverage), variant annotation (gene, effect), internal population frequencies (minor allele frequency, genotype frequency) and external dataset filters (gnomAD , RVIS ) to identify "qualifying variants" that meet a specific set of user-defined criteria. These sophisticated needs require systematic logging and version control for re-analyses and reproducibility. As the data size and number of simultaneous users increase, ad-hoc analyses become prohibitively inefficient in the conventional single joint-genotyped VCF framework.
To address these constraints and dynamic analyses needs, we have developed ATAV (Analysis Tool for Annotated Variants, see Fig. 1) to streamline genomic analysis needs ranging from the standard diagnostic case interpretation to large-scale cohort analyses for disease-associated gene discovery. The ATAV platform is built on an open-source relational database. The database (ATAVDB) is configured with a feature allowing data replication across a cluster of nodes. ATAVDB contains sample level variant data, read coverage data, variant annotation data, external annotation data, and metadata. A data pipeline toolkit provided with the code extracts variants, annotations and associated quality data from VCF files and the coverage and genotype quality from BAM files. The Institute for Genomic Medicine (IGM) at Columbia University currently has data for more than 100 K whole exomes, and the coding-regions for over 10 K whole genomes stored in ATAVDB. It contains over 24 billion variant calls from over 220 million distinct genomic co-ordinates and read coverage information for all samples.
Several other tools and systems have been developed that address some of the same needs. For example, GORpipe  offers a fast way to query variant data, but does not offer easy pipelines for performing diagnostic variant prioritization or association studies. TileDB-VCF (https://github.com/TileDB-Inc/TileDB-VCF) offers a solution for simple processing, storage, and querying of data derived from single sample VCF files, but additional analyses have to be implemented by the user. Other frameworks such as Glow (https://github.com/projectglow/glow) and Hail (Hail Team. Hail 0.2. https://github.com/hail-is/hail) recommend jointly-called variant files as input, but offer functions for more complex analyses such as association tests. ATAV is, to the best of our knowledge, the only one that takes single sample files as input and still offers tools that can easily be used to perform full diagnostic analyses and rare-variant association tests. Furthermore, ATAV offers a web-interface and API that offers users access to a rich catalogue of variants detected in patients with a variety of diseases.
We use Percona Server for MySQL and its high-performance storage engine Percona TokuDB to improve scalability and operational efficiency. In the database, we store a universal variant list across all samples, annotation data that is annotated through ClinEff , sample-level variant calls and associated quality metrics, as well as all sites’ coverage data for inferring reference alleles at non-call sites. In addition, ATAV has standardized code for incorporating any external data that is gene-based, site-based, or variant-based. ATAVDB currently stores external databases such as allele frequencies from gnomAD , ExAC , or, DiscovEHR ; scores such as TraP , LIMBR , MTR , RVIS , subRVIS , REVEL , PrimateAI , CCR ; and clinical annotations from ClinVar [16, 17], ClinGen , HGMD , and OMIM (see Fig. 2).
Variants, annotations and associated quality data are extracted from single-sample VCF files. A custom script checks whether a different representation of the same variant already exists in the database before adding a new variant to the database to ensure the same variant is represented identically in all samples.
For efficiently storing coverage information for every site and every sample, the ATAV data pipeline parses the BAM files to generate read coverage data and converts site coverage values into binned values: a [0–9]; b [10–19]; c [20–29]; d [30–49]; e [50–199]; f ≥ 200. A run-length encoding procedure is used to further compress data within fixed 1000 bp block regions (see Fig. 3). This method reduces the data size by about 1000 times making it possible to store the coverage information for more than 100 K samples. Years of applied use have helped us to identify the information that is most often required for the standard genetic analyses performed as part of both diagnostic genetic studies and gene discovery. For example, in diagnostic analyses for identifying de novo mutations in affected children, it is necessary to know that the parental samples have sufficient coverage at the relevant site, but not necessary to know the precise number of reads, leading to the binning strategy for coverage described above. For the vast majority of applications, we have found that the necessary information can be economically stored and retrieved as described.
The platform architecture is depicted in Fig. 4. In order to run ATAV jobs, users log in to the head node, which automatically allocates resources and submits jobs to the cluster. A standard setup with a 6 node Sun Grid Engine (SGE) cluster (2 × 10 Cores, 128 GB RAM) allows at least 100 jobs to be run simultaneously. Each job queries a replica database with minimum database connections thus optimizing speed and workload. Using a local customized bioinformatics pipeline, it is possible to continue loading new samples into the master database which will automatically replicate to all replica databases.
The ATAV command line tool is the programmatic user interface to ATAVDB. Written in Java, ATAV consists of three modules. (1) The command line parser and query engine translate user defined parameters and the input sample list (in PLINK’s PED format ) into an efficient SQL query for interrogating the relational database, (2) A runtime variant object creator parses SQL output into a collection of variant objects. Each variant object includes variant information (genomic coordinates, annotation), variant calls in the sample list, sample genotype calls at coordinates without a called variant and external annotation data. (3) A statistical analysis module iterates over the variant objection collection to perform downstream analyses. ATAV currently supports tests for diagnostic analyses such as identifying putative de novo and inherited genotypes of interest in trios, and a framework for performing region-based rare-variant collapsing analyses that identify genes or other genomic units that carry an excess of qualifying variants among cases in comparison to the background variation observed in internal controls in ATAVDB. Furthermore, ATAV provides commands to output other popular formats such as PLINK’s PED/MAP files or multi-sample VCF files that can be used as input for many tools outside the ATAV framework. The modularized ATAV framework enables the continuous development of new functions that operate on sequencing/variant datasets. Critical to data integrity, all ATAV analyses include an auditable log of the software and database versions, the filter parameters used, the input sample lists used in the specific run and all runtime logs that ensure full reproducibility.
The analysts and researchers at the IGM, have run about 33,000 ATAV jobs within the last year. From a runtime perspective, 22,000 jobs completed in minutes, 8000 jobs completed in hours, and the remaining 3000 jobs completed within two days.
To allow access to variant-level data directly from the full dataset (for authorized users) or the publicly available dataset (for anonymous users) in ATAVDB, we provide a web-based user interface, the ATAV data browser (http://atavdb.org/). It supports the search of variants by gene, region, and variant ID. The gene or region view displays a list of variants with allele count, allele frequency, number of samples, effect, gene etc. The variant view (see Fig. 5) displays a set of annotations (effect, gene, transcript, PolyPhen ) and details about variant carriers (gender, predicted ancestry, phenotype, and quality metrics). It includes links to other public data resources such as Ensembl, gnomAD , ClinVar [16, 17], and others, and directly integrates additional annotations via APIs such as. the Genoox Franklin API for clinical variant interpretation (see Fig. 5). The data browser has several advanced filters that restrict results to rare or ultra-rare variants (using a maximum allele frequency threshold), high quality variants, or variants from samples with a specific phenotype. This variant information can also be easily queried programmatically through a REST API. The public view currently contains more than 40,000 samples representing a mix of cases and healthy controls of diverse ancestries. Users can look up potential disease-causing variants and check whether the phenotype of variant carriers in ATAVDB matches their phenotype of interest. In contrast to many other platforms, the data browser is able to show data of newly added samples in real-time and is therefore evolving rapidly as more and more samples are sequenced. While the current version only supports hg19, future updates will also include a version for hg38 once enough data has been generated.
ATAV provides functions for all recommended steps of the rare-variant collapsing workflow recently summarized in Povysil et al. . For the sample pruning steps, ATAV creates the necessary input files by querying data from ATAVDB and automatically calls existing standard tools such as KING  or FlashPCA . Because the coverage information for every sample and site is already efficiently stored in ATAVDB, ATAV can efficiently compare coverage between cases and controls and provides two different tests to perform coverage harmonization: sites can be removed if cases and controls show differing proportions of individuals with enough coverage ; or if a binomial test shows that the case/control status and coverage are not independent . The outputs of the sample pruning and coverage harmonization steps can be used as inputs for dominant or recessive collapsing models. Within the collapsing model call, ATAV selects qualifying variants (QVs) that pass filters based on variant quality (Phred quality (QUAL), genotype Phred quality (GQ), quality by depth (QD), mapping quality (MQ), and variant quality score log-odds (VQSLOD)), variant annotation (effect, pathogenicity prediction scores, intolerance scores), as well as internal and external minor allele frequencies (MAFs). All QVs are used for building the collapsing matrix, a gene-by-individual indicator matrix with a value of 0 if there is no qualifying variant found in that gene in that individual, and a value of 1 if there is at least one. This collapsing matrix is used for examining associations between genes with QVs and the phenotype of interest by using a Fisher’s exact test or Firth-based logistic regression. Finally, quantile–quantile (QQ) plots are created and the genomic inflation factor lambda is estimated using a permutation-based expected distribution of p-values . A standard collapsing analysis usually consists of several different models that all capture specific types of QVs. While quality control (QC) filters are used for all models, other filters, such as the predicted variant effects or population allele frequencies, depend on the specific model in use. In order to speed up computation, ATAV provides the option of running a general collapsing model first, using the QC filters that are shared by all models and relaxed allele frequency thresholds. The output of this initial model can be used as input for a “collapsing-lite” function that makes it possible to run the individual collapsing models within minutes because additional filters can simply be applied to the previous output and the variant database does not have to be queried again. Example commands detailing the full workflow can be found on GitHub (https://github.com/nickzren/atav/wiki/Collapsing-Workflow).
ATAV also supports a diagnostic workflow, which highlights all variants in an individual that are candidates to satisfy criteria for “Pathogenic” or “Likely Pathogenic” according to the American College of Medical Genetics and Genomics (ACMG). These candidate variants can then undergo further examination by genetic counselors and a laboratory director. All annotations and filters mentioned previously, such as QC filters or internal or external MAFs, are also important for diagnostic analyses—especially for singletons where we cannot use additional family information. In addition, ATAV provides special functions for trios and families to reduce the number of potential disease-causing variants in the final output. ATAV leverages information about family structure and affectedness status that is provided by the sample file (PLINK-style PED file). Multiple families can be analyzed at once and related controls are automatically removed when calculating control frequencies. Furthermore, the affectedness status is used to decide whether to look for inherited or de novo variants. In the standard trio case of one affected offspring and two unaffected parents, ATAV uses a series of functions to extract all novel genotypes: de novo variants, newly compound-heterozygous, and newly homozygous variants. For distinguishing compound-heterozygosity from variants that are in-phase, ATAV checks that both parents carry one of the qualifying variants. ATAV not only considers the genotype of the individuals, but also their coverage. If the coverage at a variant site is below a minimum threshold of 10 for any of the individuals the variant is still included in the output, but flagged as possibly de novo, possibly newly compound-heterozygous or possibly newly homozygous. Furthermore, ATAV identifies putative parental-mosaic variant transmissions. For each parent–child pair, it extracts all variants that were transmitted from parent to child where the variant in the parent has a low proportion of alternate alleles indicating mosaicism.
ATAV also combines information from multiple variant and disease databases (e.g. ClinVar [16, 17], HGMD , OMIM, ClinGen ) into an external annotation dataset called KnownVar. The data is stored in ATAVDB and regularly updated. KnownVar annotations are not only included if the "exact" variant has been reported before, but also if a different variant at the same site has been linked to disease. Typical annotations include the associated disease, ClinVar clinical significance, HGMD Class and Pubmed IDs of relevant papers. Variants are also annotated with information extracted from HGMD and ClinVar for any disease-associated variants in close proximity. On a gene level, annotations include the total number of likely pathogenic or pathogenic variants of each category (copy number variation, small insertion/deletion, splice, nonsense, missense) in ClinVar, disease associations and inheritance from OMIM, and dosage sensitivity from ClinGen. All the information provided by KnownVar can be used as additional information in the diagnostic setting to evaluate whether a variant can be considered as diagnostic for a specific patient. Example commands for running diagnostic analyses with ATAV can be found on GitHub (https://github.com/nickzren/atav/wiki/Diagnostic-Workflow).
The collapsing framework of ATAV has enabled the confirmation of known and the discovery of novel genes in a wide range of diseases such as epilepsies [27, 28], sudden unexplained death in epilepsy , congenital kidney malformations , chronic kidney disease , amyotrophic lateral sclerosis [32, 33], Alzheimer’s disease , retinal dystrophy , idiopathic pulmonary fibrosis , and heart failure . Cirulli et al. 2015  used ATAV’s rare-variant collapsing framework described above to look for genes with an excess of rare, presumably deleterious variants in patients with amyotrophic lateral sclerosis compared to controls. Among the genes that reached study-wide significance were known ALS genes, such as SOD1, but also a novel one called TBK1. In a recent study , ATAV was used to detect a significant enrichment of rare protein-truncating variants in the TTN gene in patients with heart failure of mostly ischemic etiology compared with controls.
Furthermore the diagnostic framework has helped to identify both diagnostic genotypes in known genes and candidate genotypes in novel genes in a wide range of diseases including stillbirth , rare undiagnosed genetic disorders [37, 38], epilepsies [39,40,41], alternating hemiplegia of childhood , and chronic kidney disease . Zhu et al.  used ATAV’s trio diagnostic pipeline to analyze 119 patients with undiagnosed genetic diseases. By restricting the analysis to de novo variants, newly compound-heterozygous or newly homozygous variants and integrating data from ClinVar, HGMD, and OMIM, the authors were able to obtain a genetic diagnosis for 29 (24%) of patients. Furthermore, they identified an enrichment of damaging de novo mutations in intolerant genes highlighting a possible way of identifying novel disease genes and expanding phenotypes for known disease genes. A recent study on causal genetic variants in stillbirth  used ATAV’s non-trio workflow to prioritized variants thought to be enriched for pathogenicity by focusing on variants in Mendelian disease genes from OMIM that are rare in the general population. The authors identified a diagnosis in 15 of 246 cases of stillbirth (6.1%) involving both genes that had been previously implicated in stillbirth and genes that are potential candidates for phenotypic expansion.
Examining the distribution of allele frequencies (AFs) can help characterize the composition of variants of the more than 100 k samples in our version of ATAVDB. Figure 6 shows AF distributions for bi-allelic coding variants, restricted to coding or splice regions that are well covered in the majority of our samples and using data of roughly 101,000 unrelated individuals. We also applied basic quality control filters including the removal of variants marked as potential artifacts by any of gnomAD’s filters. As expected, the vast majority of variants are singletons in our dataset, reflecting a frequency of 4.95 × 10−6 (left-most bar in panel e) and absent in the gnomAD v2.1 exome data (left-most bar in panel f).
We present ATAV as an analysis platform for large-scale whole-exome and whole-genome sequencing projects. In order to encourage the widespread adoption of ATAV, we provide detailed guidelines on GitHub for how to set up the whole framework, including the database, command line tool, and data browser. This gives any user the ability to both create a local version of ATAVDB into which they can load their data and use all functionality of the command line tool and data browser in order to run analyses and queries provided by ATAV. The advantages of the ATAV framework are that (1) it allows continuous real-time analyses of all samples loaded into the database without the need for computationally demanding joint calling preceding each analysis and (2) it allows convenient logging of specific analyses performed. The newly added ATAV data browser provides easy access even to users with little computational experience by providing an intuitive web interface to query variant-level data directly from the database.
Our experience with this platform on a database carrying more than 100,000 samples indicates that a relational database can be optimized in a way that makes it possible to analyze large-scale genomic datasets. Our current data processing and storage framework is robust and flexible when combining data from multiple projects and analyzing exomes and genomes together. ATAV supports diagnostic analyses for trios and singletons, as well as rare-variant collapsing analyses for finding disease associations in complex diseases. Further optimizations are possible such as database sharding, a horizontal partition of data in a database or search engine. Other potential solutions include storing the data in Hadoop Distributed File System (HDFS) and utilizing Apache Spark to do distributed cluster computing. This would allow the processing of large amounts of variant data in parallel at once, speeding up computations and enabling an even further increase in sample sizes.
The goal of ATAV is to standardize and optimize storage and data processing for large scale sequencing data across multiple studies and to provide an easy to use interface for users with little computational experience while ensuring full reproducibility.
All code for building ATAV is publicly available, providing a convenient way for other groups to build up their own analysis platform, database, and user interface. Additionally, since we provide general access to part of our database via the ATAV browser, users can also query one of the largest variant databases available for patients sequenced at a tertiary care center. Currently, public users have access to summary-level data for more than 40,000 samples, but since the data of newly sequenced samples are added in real-time, this number grows steadily, increasing the value of the database even further as more and more samples are sequenced.
Availability and requirements
Project name: ATAV
Project home page: https://github.com/nickzren/atav
Operating system(s): Platform independent
Other requirements: Java 1.8 or higher, Percona Server 5.6 or higher
License: MIT License
Any restrictions to use by non-academics: No restrictions
Web links and URLs
ATAV data browser, http://www.atavdb.org/
Ensembl GRCh37, https://grch37.ensembl.org/
Genome Asia, https://browser.genomeasia100k.org/
GME Variome, http://igm.ucsd.edu/gme/
TOPMed hg19, https://bravo.sph.umich.edu/freeze3a/hg19/
UCSC Genome Browser, https://genome.ucsc.edu/index.html
Variant call format
Binary alignment map
The Genome aggregation database
Residual variation intolerance score
Institute for genomic medicine
Clinical Variant Annotations Software
The Exome Aggregation Consortium
Genomic evolutionary rate profiling
The transcript-inferred pathogenicity
The localized intolerance model using Bayesian regression
The missense tolerance ratio
The human gene mutation database
Online Mendelian inheritance in man
Structured query language
Genotype quality score
Quality by depth score
Mapping quality score
Variant quality score log-odds
Minor allele frequencies
Hadoop distributed file system
Hout CV Van, Tachmazidou I, Backman JD, Hoffman JX, Ye B, Pandey AK, et al. Whole exome sequencing and characterization of coding variation in 49,960 individuals in the UK Biobank. bioRxiv. 2019; 572347. https://doi.org/10.1101/572347.
Taliun D, Harris DN, Kessler MD, Carlson J, Szpiech ZA, Torres R, et al. Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program. bioRxiv. 2019; 563866. https://doi.org/10.1101/563866.
Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature. 2020;581:434–43. https://doi.org/10.1038/s41586-020-2308-7.
Petrovski S, Gussow AB, Wang Q, Halvorsen M, Han Y, Weir WH, et al. The intolerance of regulatory sequence to genetic variation predicts gene dosage sensitivity. PLoS Genet. 2015;11:e1005492. https://doi.org/10.1371/journal.pgen.1005492.
Guðbjartsson H, Georgsson GF, Guðjónsson SA, Valdimarsson RÞ, Sigurðsson JH, Stefánsson SK, et al. GORpipe: a query tool for working with sequence data based on a Genomic Ordered Relational (GOR) architecture. Bioinformatics. 2016;32:3081–8. https://doi.org/10.1093/bioinformatics/btw199.
Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms. SnpEff Fly (Austin). 2012;6:80–92. https://doi.org/10.4161/fly.19695.
Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536:285–91. https://doi.org/10.1038/nature19057.
Dewey FE, Murray MF, Overton JD, Habegger L, Leader JB, Fetterolf SN, et al. Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study. Science. 2016;354:aaf6814. https://doi.org/10.1126/science.aaf6814.
Gelfman S, Wang Q, McSweeney KM, Ren Z, La Carpia F, Halvorsen M, et al. Annotating pathogenic non-coding variants in genic regions. Nat Commun. 2017;8:236. https://doi.org/10.1038/s41467-017-00141-2.
Hayeck TJ, Stong N, Wolock CJ, Copeland B, Kamalakaran S, Goldstein DB, et al. Improved pathogenic variant localization via a hierarchical model of sub-regional intolerance. Am J Hum Genet. 2019;104:299–309. https://doi.org/10.1016/J.AJHG.2018.12.020.
Traynelis J, Silk M, Wang Q, Berkovic SF, Liu L, Ascher DB, et al. Optimizing genomic medicine in epilepsy through a gene-customized approach to missense variant interpretation. Genome Res. 2017;27:1715–29. https://doi.org/10.1101/gr.226589.117.
Gussow AB, Petrovski S, Wang Q, Allen AS, Goldstein DB. The intolerance to functional genetic variation of protein domains predicts the localization of pathogenic mutations within genes. Genome Biol. 2016;17:9. https://doi.org/10.1186/s13059-016-0869-4.
Ioannidis NM, Rothstein JH, Pejaver V, Middha S, McDonnell SK, Baheti S, et al. REVEL: an ensemble method for predicting the pathogenicity of rare missense variants. Am J Hum Genet. 2016;99:877–85. https://doi.org/10.1016/J.AJHG.2016.08.016.
Sundaram L, Gao H, Padigepati SR, McRae JF, Li Y, Kosmicki JA, et al. Predicting the clinical impact of human mutation with deep neural networks. Nat Genet. 2018;50:1161–70. https://doi.org/10.1038/s41588-018-0167-z.
Havrilla JM, Pedersen BS, Layer RM, Quinlan AR. A map of constrained coding regions in the human genome. Nat Genet. 2019;51:88–95. https://doi.org/10.1038/s41588-018-0294-6.
Landrum MJ, Lee JM, Riley GR, Jang W, Rubinstein WS, Church DM, et al. ClinVar: public archive of relationships among sequence variation and human phenotype. Nucleic Acids Res. 2014;42(Database issue):D980–5. https://doi.org/10.1093/nar/gkt1113.
Landrum MJ, Lee JM, Benson M, Brown GR, Chao C, Chitipiralla S, et al. ClinVar: improving access to variant interpretations and supporting evidence. Nucleic Acids Res. 2018;46:D1062–7. https://doi.org/10.1093/nar/gkx1153.
Rehm HL, Berg JS, Brooks LD, Bustamante CD, Evans JP, Landrum MJ, et al. ClinGen: the clinical genome resource. N Engl J Med. 2015;372:2235–42. https://doi.org/10.1056/NEJMsr1406261.
Stenson PD, Ball EV, Mort M, Phillips AD, Shiel JA, Thomas NST, et al. Human gene mutation database (HGMD®): 2003 update. Hum Mutat. 2003;21:577–81. https://doi.org/10.1002/humu.10212.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–9. https://doi.org/10.1038/nmeth0410-248.
Povysil G, Petrovski S, Hostyk J, Aggarwal V, Allen AS, Goldstein DB. Rare-variant collapsing analyses for complex traits: guidelines and applications. Nat Rev Genet. 2019;20:747–59.
Manichaikul A, Mychaleckyj JC, Rich SS, Daly K, Sale M, Chen W-M. Robust relationship inference in genome-wide association studies. Bioinformatics. 2010;26:2867–73. https://doi.org/10.1093/bioinformatics/btq559.
Abraham G, Qiu Y, Inouye M. FlashPCA2: principal component analysis of Biobank-scale genotype datasets. Bioinformatics. 2017;33:2776–8. https://doi.org/10.1093/bioinformatics/btx299.
Petrovski S, Todd JL, Durheim MT, Wang Q, Chien JW, Kelly FL, et al. An exome sequencing study to assess the role of rare genetic variation in pulmonary fibrosis. Am J Respir Crit Care Med. 2017;196:82–93. https://doi.org/10.1164/rccm.201610-2088OC.
Raghavan NS, Brickman AM, Andrews H, Manly JJ, Schupf N, Lantigua R, et al. Whole-exome sequencing in 20,197 persons for rare variants in Alzheimer’s disease. Ann Clin Transl Neurol. 2018;5:832–42. https://doi.org/10.1002/acn3.582.
Allen AS, Bellows ST, Berkovic SF, Bridgers J, Burgess R, Cavalleri G, et al. Ultra-rare genetic variation in common epilepsies: a case-control sequencing study. Lancet Neurol. 2017;16:135–43. https://doi.org/10.1016/S1474-4422(16)30359-3.
Zhu X, Padmanabhan R, Copeland B, Bridgers J, Ren Z, Kamalakaran S, et al. A case-control collapsing analysis identifies epilepsy genes implicated in trio sequencing studies focused on de novo mutations. PLoS Genet. 2017;13:e1007104. https://doi.org/10.1371/journal.pgen.1007104.
Bagnall RD, Crompton DE, Petrovski S, Lam L, Cutmore C, Garry SI, et al. Exome-based analysis of cardiac arrhythmia, respiratory control, and epilepsy genes in sudden unexpected death in epilepsy. Ann Neurol. 2016;79:522–34. https://doi.org/10.1002/ana.24596.
Sanna-Cherchi S, Khan K, Westland R, Krithivasan P, Fievet L, Rasouly HM, et al. Exome-wide association study identifies GREB1L mutations in congenital kidney malformations. Am J Hum Genet. 2017;101:789–802. https://doi.org/10.1016/J.AJHG.2017.09.018.
Cameron-Christie S, Wolock CJ, Groopman E, Petrovski S, Kamalakaran S, Povysil G, et al. Exome-based rare-variant analyses in CKD. J Am Soc Nephrol. 2019;30:1109–22. https://doi.org/10.1681/ASN.2018090909.
Cirulli ET, Lasseigne BN, Petrovski S, Sapp PC, Dion PA, Leblond CS, et al. Exome sequencing in amyotrophic lateral sclerosis identifies risk genes and pathways. Science. 2015;347:1436–41. https://doi.org/10.1126/science.aaa3650.
Gelfman S, Dugger S, de Araujo Martins-Moreno C, Ren Z, Wolock CJ, Shneider NA, et al. A new approach for rare variation collapsing on functional protein domains implicates specific genic regions in ALS. Genome Res. 2019;29:809–18. https://doi.org/10.1101/gr.243592.118.
Wolock CJ, Stong N, Ma CJ, Nagasaki T, Lee W, Tsang SH, et al. A case–control collapsing analysis identifies retinal dystrophy genes associated with ophthalmic disease in patients with no pathogenic ABCA4 variants. Genet Med. 2019;21:2336–44.
Povysil G, Chazara O, Carss KJ, Deevi SVV, Wang Q, Armisen J, et al. Assessing the role of rare genetic variation in patients with heart failure. JAMA Cardiol. 2020. https://doi.org/10.1001/jamacardio.2020.6500.
Stanley KE, Giordano J, Thorsten V, Buchovecky C, Thomas A, Ganapathi M, et al. Causal Genetic Variants in Stillbirth. N Engl J Med. 2020. https://doi.org/10.1056/nejmoa1908753.
Zhu X, Petrovski S, Xie P, Ruzzo EK, Lu Y-F, McSweeney KM, et al. Whole-exome sequencing in undiagnosed genetic diseases: interpreting 119 trios. Genet Med. 2015;17:774. https://doi.org/10.1038/gim.2014.191.
Petrovski S, Shashi V, Petrou S, Schoch K, McSweeney KM, Dhindsa RS, et al. Exome sequencing results in successful riboflavin treatment of a rapidly progressive neurological condition. Mol Case Stud. 2015;1:a000257.
Allen AS, Berkovic SF, Cossette P, Delanty N, Dlugos D, et al. De novo mutations in epileptic encephalopathies. Nature. 2013;501:217–21. https://doi.org/10.1038/nature12439.
Myers CT, Stong N, Mountier EI, Helbig KL, Freytag S, Sullivan JE, et al. De Novo mutations in PPP3CA cause severe neurodevelopmental disease with seizures. Am J Hum Genet. 2017;101:516–24.
Petrovski S, Küry S, Myers CT, Anyane-Yeboa K, Cogné B, Bialer M, et al. Germline de Novo mutations in GNB1 cause severe neurodevelopmental disability, hypotonia, and seizures. Am J Hum Genet. 2016;98:1001–10.
Heinzen EL, Swoboda KJ, Hitomi Y, Gurrieri F, De Vries B, Tiziano FD, et al. De novo mutations in ATP1A3 cause alternating hemiplegia of childhood. Nat Genet. 2012;44:1030–4.
Groopman EE, Marasa M, Cameron-Christie S, Petrovski S, Aggarwal VS, Milo-Rasouly H, et al. Diagnostic utility of exome sequencing for kidney disease. N Engl J Med. 2019;380:142–51.
We thank Slavé Petrovski for contributing to the original design of the analysis framework and Quanli Wang for contributing to the original platform development.
This project was funded by the Institute for Genomic Medicine, Columbia University Irving Medical Center. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Consent for publication
DBG is a founder of and holds equity in Praxis, holds equity in Q-State Biosciences, serves as a consultant to AstraZeneca, and has received research support from Janssen, Gilead, Biogen, AstraZeneca, and Union Chimique Belge (UCB). ZR, GP, JAH, HC, and NB declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ren, Z., Povysil, G., Hostyk, J.A. et al. ATAV: a comprehensive platform for population-scale genomic analyses. BMC Bioinformatics 22, 149 (2021). https://doi.org/10.1186/s12859-021-04071-1
- Association testing
- Genome analysis
- Gene discovery
- Web platform