Search

Whole-genome sequencing (WGS) allows for a comprehensive view of the sequence of the human genome. We present and apply integrated methodologic steps for interrogating WGS data to characterize the genetic architecture of 10 heart- and blood-related traits in a sample of 1,860 African Americans. In order to evaluate the contribution of regulatory and non-protein coding regions of the genome, we conducted aggregate tests of rare variation across the entire genomic landscape using a sliding window, complemented by an annotation-based assessment of the genome using predefined regulatory elements and within the first intron of all genes. These tests were performed treating all variants equally as well as with individual variants weighted by a measure of predicted functional consequence. Significant findings were assessed in 1,705 individuals of European ancestry. After these steps, we identified and replicated components of the genomic landscape significantly associated with heart- and blood-related traits. For two traits, lipoprotein(a) levels and neutrophil count, aggregate tests of low-frequency and rare variation were significantly associated across multiple motifs. For a third trait, cardiac troponin T, investigation of regulatory domains identified a locus on chromosome 9. These practical approaches for WGS analysis led to the identification of informative genomic regions and also showed that defined non-coding regions, such as first introns of genes and regulatory domains, are associated with important risk factor phenotypes. This study illustrates the tractable nature of WGS data and outlines an approach for characterizing the genetic architecture of complex traits.

Characterizing meiotic recombination rates across the genomes of nonhuman primates is important for understanding the genetics of primate populations, performing genetic analyses of phenotypic variation and reconstructing the evolution of human recombination. Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primates in biomedical research. We constructed a high-resolution genetic map of the rhesus genome based on whole genome sequence data from Indian-origin rhesus macaques. The genetic markers used were approximately 18 million SNPs, with marker density 6.93 per kb across the autosomes. We report that the genome-wide recombination rate in rhesus macaques is significantly lower than rates observed in apes or humans, while the distribution of recombination across the macaque genome is more uniform. These observations provide new comparative information regarding the evolution of recombination in primates.

Background The decreasing costs of sequencing are driving the need for cost effective and real time variant calling of whole genome sequencing data. The scale of these projects are far beyond the capacity of typical computing resources available with most research labs. Other infrastructures like the cloud AWS environment and supercomputers also have limitations due to which large scale joint variant calling becomes infeasible, and infrastructure specific variant calling strategies either fail to scale up to large datasets or abandon joint calling strategies. Results We present a high throughput framework including multiple variant callers for single nucleotide variant (SNV) calling, which leverages hybrid computing infrastructure consisting of cloud AWS, supercomputers and local high performance computing infrastructures. We present a novel binning approach for large scale joint variant calling and imputation which can scale up to over 10,000 samples while producing SNV callsets with high sensitivity and specificity. As a proof of principle, we present results of analysis on Cohorts for Heart And Aging Research in Genomic Epidemiology (CHARGE) WGS freeze 3 dataset in which joint calling, imputation and phasing of over 5300 whole genome samples was produced in under 6 weeks using four state-of-the-art callers. The callers used were SNPTools, GATK-HaplotypeCaller, GATK-UnifiedGenotyper and GotCloud. We used Amazon AWS, a 4000-core in-house cluster at Baylor College of Medicine, IBM power PC Blue BioU at Rice and Rhea at Oak Ridge National Laboratory (ORNL) for the computation. AWS was used for joint calling of 180 TB of BAM files, and ORNL and Rice supercomputers were used for the imputation and phasing step. All other steps were carried out on the local compute cluster. The entire operation used 5.2 million core hours and only transferred a total of 6 TB of data across the platforms. Conclusions Even with increasing sizes of whole genome datasets, ensemble joint calling of SNVs for low coverage data can be accomplished in a scalable, cost effective and fast manner by using heterogeneous computing platforms without compromising on the quality of variants. Electronic supplementary material The online version of this article (doi:10.1186/s12859-016-1211-6) contains supplementary material, which is available to authorized users.

The goal of the International HapMap Project is to determine the common patterns of DNA sequence variation in the human genome and to make this information freely available in the public domain. An international consortium is developing a map of these patterns across the genome by determining the genotypes of one million or more sequence variants, their frequencies and the degree of association between them, in DNA samples from populations with ancestry from parts of Africa, Asia and Europe. The HapMap will allow the discovery of sequence variants that affect common disease, will facilitate development of diagnostic tools, and will enhance our ability to choose targets for therapeutic intervention. © 2003 Nature Publishing Group.

Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primate in biomedical research, have the largest natural geographic distribution of any nonhuman primate, and have been the focus of much evolutionary and behavioral investigation. Consequently, rhesus macaques are one of the most thoroughly studied nonhuman primate species. However, little is known about genome-wide genetic variation in this species. A detailed understanding of extant genomic variation among rhesus macaques has implications for the use of this species as a model for studies of human health and disease, as well as for evolutionary population genomics. Whole-genome sequencing analysis of 133 rhesus macaques revealed more than 43.7 million single-nucleotide variants, including thousands predicted to alter protein sequences, transcript splicing, and transcription factor binding sites. Rhesus macaques exhibit 2.5-fold higher overall nucleotide diversity and slightly elevated putative functional variation compared with humans. This functional variation in macaques provides opportunities for analyses of coding and noncoding variation, and its cellular consequences. Despite modestly higher levels of nonsynonymous variation in the macaques, the estimated distribution of fitness effects and the ratio of nonsynonymous to synonymous variants suggest that purifying selection has had stronger effects in rhesus macaques than in humans. Demographic reconstructions indicate this species has experienced a consistently large but fluctuating population size. Overall, the results presented here provide new insights into the population genomics of nonhuman primates and expand genomic information directly relevant to primate models of human disease.

Structural variants are implicated in numerous diseases and make up the majority of varying nucleotides among human genomes. Here we describe an integrated set of eight structural variant classes comprising both balanced and unbalanced variants, which we constructed using short-read DNA sequencing data and statistically phased onto haplotype blocks in 26 human populations. Analysing this set, we identify numerous gene-intersecting structural variants exhibiting population stratification and describe naturally occurring homozygous gene knockouts that suggest the dispensability of a variety of human genes. We demonstrate that structural variants are enriched on haplotypes identified by genome-wide association studies and exhibit enrichment for expression quantitative trait loci. Additionally, we uncover appreciable levels of structural variant complexity at different scales, including genic loci subject to clusters of repeated rearrangement and complex structural variants with multiple breakpoints likely to have formed through individual mutational events. Our catalogue will enhance future studies into structural variant demography, functional impact and disease association.

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies., Wellcome Trust (London, England) (Core Award 090532/Z/09/Z), Wellcome Trust (London, England) (Senior Investigator Award 095552/Z/11/Z ), Wellcome Trust (London, England) (WT095908), Wellcome Trust (London, England) (WT109497), Wellcome Trust (London, England) (WT098051), Wellcome Trust (London, England) (WT086084/Z/08/Z), Wellcome Trust (London, England) (WT100956/Z/13/Z ), Wellcome Trust (London, England) (WT097307), Wellcome Trust (London, England) (WT0855322/Z/08/Z ), Wellcome Trust (London, England) (WT090770/Z/09/Z ), Wellcome Trust (London, England) (Major Overseas program in Vietnam grant 089276/Z.09/Z), Medical Research Council (Great Britain) (grant G0801823), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I02593X/1), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I021213/1), Zhongguo ke xue ji shu qing bao yan jiu suo. Office of 863 Programme of China (2012AA02A201), National Basic Research Program of China (2011CB809201), National Basic Research Program of China (2011CB809202), National Basic Research Program of China (2011CB809203), National Natural Science Foundation of China (31161130357), Shenzhen Municipal Government of China (grant ZYC201105170397A), Canadian Institutes of Health Research (grant 136855), Quebec Ministry of Economic Development, Innovation, and Exports (PSR-SIIRI-195), Germany. Bundesministerium für Bildung und Forschung (0315428A), Germany. Bundesministerium für Bildung und Forschung (01GS08201), Germany. Bundesministerium für Bildung und Forschung (BMBF-EPITREAT grant 0316190A), Deutsche Forschungsgemeinschaft (Emmy Noether Grant KO4037/1-1), Beatriu de Pinos Program (2006 BP-A 10144), Beatriu de Pinos Program (2009 BP-B 00274), Spanish National Institute for Health (grant PRB2 IPT13/0001-ISCIII-SGEFI/FEDER), Japan Society for the Promotion of Science (fellowship number PE13075), Marie Curie Actions Career Integration (grant 303772), Fonds National Suisse del la Recherche, SNSF, Scientifique (31003A_130342), National Center for Biotechnology Information (U.S.) (U54HG3067), National Center for Biotechnology Information (U.S.) (U54HG3273), National Center for Biotechnology Information (U.S.) (U01HG5211), National Center for Biotechnology Information (U.S.) (U54HG3079), National Center for Biotechnology Information (U.S.) (R01HG2898), National Center for Biotechnology Information (U.S.) (R01HG2385), National Center for Biotechnology Information (U.S.) (RC2HG5552), National Center for Biotechnology Information (U.S.) (U01HG6513), National Center for Biotechnology Information (U.S.) (U01HG5214), National Center for Biotechnology Information (U.S.) (U01HG5715), National Center for Biotechnology Information (U.S.) (U01HG5718), National Center for Biotechnology Information (U.S.) (U01HG5728), National Center for Biotechnology Information (U.S.) (U41HG7635), National Center for Biotechnology Information (U.S.) (U41HG7497), National Center for Biotechnology Information (U.S.) (R01HG4960), National Center for Biotechnology Information (U.S.) (R01HG5701), National Center for Biotechnology Information (U.S.) (R01HG5214), National Center for Biotechnology Information (U.S.) (R01HG6855), National Center for Biotechnology Information (U.S.) (R01HG7068), National Center for Biotechnology Information (U.S.) (R01HG7644), National Center for Biotechnology Information (U.S.) (DP2OD6514), National Center for Biotechnology Information (U.S.) (DP5OD9154), National Center for Biotechnology Information (U.S.) (R01CA166661), National Center for Biotechnology Information (U.S.) (R01CA172652), National Center for Biotechnology Information (U.S.) (P01GM99568), National Center for Biotechnology Information (U.S.) (R01GM59290), National Center for Biotechnology Information (U.S.) (R01GM104390), National Center for Biotechnology Information (U.S.) (T32GM7790), National Center for Biotechnology Information (U.S.) (R01HL87699), National Center for Biotechnology Information (U.S.) (R01HL104608), National Center for Biotechnology Information (U.S.) (T32HL94284), National Center for Biotechnology Information (U.S.) (HHSN268201100040C), National Center for Biotechnology Information (U.S.) (HHSN272201000025C), Lundbeck Foundation (grant R170-2014-1039, Simons Foundation (SFARI award SF51), National Science Foundation (U.S.) (Research Fellowship DGE-1147470)

We deployed the genetic variant pipeline SNPTools in the cloud utilizing the Amazon Web Service (AWS). With the cloud SNPTools pipeline, we performed the SNP calling and genotype imputation on the 1000 Genomes Project Phase 3 data and assessed the quality of SNPs. We also explored different strategies of exploiting Amazon Elastic Cloud Compute instances and the Amazon Simple Storage Service in order to optimize the performance and cost of cloud computing. Our analysis shows that cloud computing will be indispensable to the Next Generation Sequencing data processing.

The synthesis, secretion and clearance of von Willebrand factor (VWF) are regulated by genetic variations in coding and promoter regions of the VWF gene. We have previously identified 19 single nucleotide polymorphisms (SNPs), primarily in introns that are associated with VWF antigen levels in subjects of European descent. In this study, we conducted race by gender analyses to compare the association of VWF SNPs with VWF antigen among 10,434 healthy Americans of European (EA) or African (AA) descent from the Atherosclerosis Risk in Communities (ARIC) study. Among 75 SNPs analyzed, 13 and 10 SNPs were associated with VWF antigen levels in EA male and EA female subjects, respectively. However, only one SNP (RS1063857) was significantly associated with VWF antigen in AA females and none was in AA males. Haplotype analysis of the ARIC samples and studying racial diversities in the VWF gene from the 1000 genomes database suggest a greater degree of variations in the VWF gene in AA subjects as compared to EA subjects. Together, these data suggest potential race and gender divergence in regulating VWF expression by genetic variations.

As DNA sequencing becomes more affordable and data ‘tsunami’ is practically here, it is clear that the inforamtics and analysis are the rate-limiting step currently for scientific discoveries as well as medical actions in the genomics enterprise. In this chapter, soon after the brief historical overview on the genomics field, we describe details on alignment and variation analysis algorithms and software packages. Remaining challenges and potential directions are discussed at the end of the chapter.

We describe initial steps for interrogating whole-genome sequence data to characterize the genetic architecture of a complex trait, levels of high-density lipoprotein cholesterol (HDL-C). We report whole-genome sequencing and analysis of 962 individuals from the Cohorts for Heart and Aging Research in Genetic Epidemiology (CHARGE) studies. From this analysis, we estimate that common variation contributes more to heritability of HDL-C levels than rare variation, and screening for mendelian variants for dyslipidemia identified individuals with extreme HDL-C levels. Whole-genome sequencing analyses highlight the value of regulatory and non-protein-coding regions of the genome in addition to protein-coding regions.

VWF is extensively glycosylated with biantennary core fucosylated glycans. Most N-linked and O-linked glycans on VWF are sialylated. FVIII is also glycosylated, with a glycan structure similar to that of VWF. ST3GAL sialyltransferases catalyze the transfer of sialic acids in the α2,3 linkage to termini of N- and O-glycans. This sialic acid modification is critical for VWF synthesis and activity. We analyzed genetic and phenotypic data from the Atherosclerosis Risk in Communities (ARIC) study for the association of single nucleotide polymorphisms (SNPs) in the ST3GAL4 gene with plasma VWF levels and FVIII activity in 12,117 subjects. We also analyzed ST3GAL4 SNPs found in 2,535 subjects of 26 ethnicities from the 1000 Genomes (1000G) project for ethnic diversity, SNP imputation, and ST3GAL4 haplotypes. We identified 14 and 1,714 ST3GAL4 variants in the ARIC GWAS and 1000G databases respectively, with 46% being ethnically diverse in their allele frequencies. Among the 14 ST3GAL4 SNPs found in ARIC GWAS, the intronic rs2186717, rs7928391, and rs11220465 were associated with VWF levels and with FVIII activity after adjustment for age, BMI, hypertension, diabetes, ever-smoking status, and ABO. This study illustrates the power of next-generation sequencing in the discovery of new genetic variants and a significant ethnic diversity in the ST3GAL4 gene. We discuss potential mechanisms through which these intronic SNPs regulate ST3GAL4 biosynthesis and the activity that affects VWF and FVIII.

Background Next-generation DNA sequencing is opening new avenues for genetic association studies in common diseases that, like deep vein thrombosis (DVT), have a strong genetic predisposition still largely unexplained by currently identified risk variants. In order to develop sequencing and analytical pipelines for the application of next-generation sequencing to complex diseases, we conducted a pilot study sequencing the coding area of 186 hemostatic/proinflammatory genes in 10 Italian cases of idiopathic DVT and 12 healthy controls. Results A molecular-barcoding strategy was used to multiplex DNA target capture and sequencing, while retaining individual sequence information. Genomic libraries with barcode sequence-tags were pooled (in pools of 8 or 16 samples) and enriched for target DNA sequences. Sequencing was performed on ABI SOLiD-4 platforms. We produced > 12 gigabases of raw sequence data to sequence at high coverage (average: 42X) the 700-kilobase target area in 22 individuals. A total of 1876 high-quality genetic variants were identified (1778 single nucleotide substitutions and 98 insertions/deletions). Annotation on databases of genetic variation and human disease mutations revealed several novel, potentially deleterious mutations. We tested 576 common variants in a case-control association analysis, carrying the top-5 associations over to replication in up to 719 DVT cases and 719 controls. We also conducted an analysis of the burden of nonsynonymous variants in coagulation factor and anticoagulant genes. We found an excess of rare missense mutations in anticoagulant genes in DVT cases compared to controls and an association for a missense polymorphism of FGA (rs6050; p = 1.9 × 10-5, OR 1.45; 95% CI, 1.22-1.72; after replication in > 1400 individuals). Conclusions We implemented a barcode-based strategy to efficiently multiplex sequencing of hundreds of candidate genes in several individuals. In the relatively small dataset of our pilot study we were able to identify bona fide associations with DVT. Our study illustrates the potential of next-generation sequencing for the discovery of genetic variation predisposing to complex diseases.

Background Whole exome capture sequencing allows researchers to cost-effectively sequence the coding regions of the genome. Although the exome capture sequencing methods have become routine and well established, there is currently a lack of tools specialized for variant calling in this type of data. Results Using statistical models trained on validated whole-exome capture sequencing data, the Atlas2 Suite is an integrative variant analysis pipeline optimized for variant discovery on all three of the widely used next generation sequencing platforms (SOLiD, Illumina, and Roche 454). The suite employs logistic regression models in conjunction with user-adjustable cutoffs to accurately separate true SNPs and INDELs from sequencing and mapping errors with high sensitivity (96.7%). Conclusion We have implemented the Atlas2 Suite and applied it to 92 whole exome samples from the 1000 Genomes Project. The Atlas2 Suite is available for download at http://sourceforge.net/projects/atlas2/. In addition to a command line version, the suite has been integrated into the Genboree Workbench, allowing biomedical scientists with minimal informatics expertise to remotely call, view, and further analyze variants through a simple web interface. The existing genomic databases displayed via the Genboree browser also streamline the process from variant discovery to functional genomics analysis, resulting in an off-the-shelf toolkit for the broader community.

Until recently, sequencing has primarily been carried out in large genome centers who also invested heavily in developing the computational infrastructure to enable post sequencing analysis. The recent advancements in sequencing technologies have lead to a wide dissemination of sequencing and we are now seeing many sequencing projects being undertaken in small laboratories. However, the limited accessibility to the computational infrastructure and high quality bioinformatic tools needed to enable analysis remains a serious road-block. The cloud computing and Software-as-a-Service (SaaS) technologies can help address this barrier. We deploy the Atlas2 Cloud Pipeline for personal genome analysis via the Genboree Workbench using software-as-a-service model. We report on a successful case study of personal genome analysis using this pipeline.

Summary Genomic structural variants (SVs) are abundant in humans, differing from other variation classes in extent, origin, and functional impact. Despite progress in SV characterization, the nucleotide resolution architecture of most SVs remains unknown. We constructed a map of unbalanced SVs (i.e., copy number variants) based on whole genome DNA sequencing data from 185 human genomes, integrating evidence from complementary SV discovery approaches with extensive experimental validations. Our map encompassed 22,025 deletions and 6,000 additional SVs, including insertions and tandem duplications. Most SVs (53%) were mapped to nucleotide resolution, which facilitated analyzing their origin and functional impact. We examined numerous whole and partial gene deletions with a genotyping approach and observed a depletion of gene disruptions amongst high frequency deletions. Furthermore, we observed differences in the size spectra of SVs originating from distinct formation mechanisms, and constructed a map constructed a map of SV hotspots formed by common mechanisms. Our analytical framework and SV map serves as a resource for sequencing-based association studies.

Background Generation of long (>5 Kb) DNA sequencing reads provides an approach for interrogation of complex regions in the human genome. Currently, large-insert whole genome sequencing (WGS) technologies from Pacific Biosciences (PacBio) enable analysis of chromosomal structural variations (SVs), but the cost to achieve the required sequence coverage across the entire human genome is high. Results We developed a method (termed PacBio-LITS) that combines oligonucleotide-based DNA target-capture enrichment technologies with PacBio large-insert library preparation to facilitate SV studies at specific chromosomal regions. PacBio-LITS provides deep sequence coverage at the specified sites at substantially reduced cost compared with PacBio WGS. The efficacy of PacBio-LITS is illustrated by delineating the breakpoint junctions of low copy repeat (LCR)-associated complex structural rearrangements on chr17p11.2 in patients diagnosed with Potocki–Lupski syndrome (PTLS; MIM#610883). We successfully identified previously determined breakpoint junctions in three PTLS cases, and also were able to discover novel junctions in repetitive sequences, including LCR-mediated breakpoints. The new information has enabled us to propose mechanisms for formation of these structural variants. Conclusions The new method leverages the cost efficiency of targeted capture-sequencing as well as the mappability and scaffolding capabilities of long sequencing reads generated by the PacBio platform. It is therefore suitable for studying complex SVs, especially those involving LCRs, inversions, and the generation of chimeric Alu elements at the breakpoints. Other genomic research applications, such as haplotype phasing and small insertion and deletion validation could also benefit from this technology. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1370-2) contains supplementary material, which is available to authorized users.

Background Massively parallel DNA sequencing generates staggering amounts of data. Decreasing cost, increasing throughput, and improved annotation have expanded the diversity of genomics applications in research and clinical practice. This expanding scale creates analytical challenges: accommodating peak compute demand, coordinating secure access for multiple analysts, and sharing validated tools and results. Results To address these challenges, we have developed the Mercury analysis pipeline and deployed it in local hardware and the Amazon Web Services cloud via the DNAnexus platform. Mercury is an automated, flexible, and extensible analysis workflow that provides accurate and reproducible genomic results at scales ranging from individuals to large cohorts. Conclusions By taking advantage of cloud computing and with Mercury implemented on the DNAnexus platform, we have demonstrated a powerful combination of a robust and fully validated software pipeline and a scalable computational resource that, to date, we have applied to more than 10,000 whole genome and whole exome samples.