A Flexible and Accurate Genotype Imputation Method for the Next Generation of Genome-Wide Association Studies

Genotype imputation methods are now being widely used in the analysis of genome-wide association studies. Most imputation analyses to date have used the HapMap as a reference dataset, but new reference panels (such as controls genotyped on multiple SNP chips and densely typed samples from the 1,000 Genomes Project) will soon allow a broader range of SNPs to be imputed with higher accuracy, thereby increasing power. We describe a genotype imputation method (IMPUTE version 2) that is designed to address the challenges presented by these new datasets. The main innovation of our approach is a flexible modelling framework that increases accuracy and combines information across multiple reference panels while remaining computationally feasible. We find that IMPUTE v2 attains higher accuracy than other methods when the HapMap provides the sole reference panel, but that the size of the panel constrains the improvements that can be made. We also find that imputation accuracy can be greatly enhanced by expanding the reference panel to contain thousands of chromosomes and that IMPUTE v2 outperforms other methods in this setting at both rare and common SNPs, with overall error rates that are 15%–20% lower than those of the closest competing method. One particularly challenging aspect of next-generation association studies is to integrate information across multiple reference panels genotyped on different sets of SNPs; we show that our approach to this problem has practical advantages over other suggested solutions

A combined genome-wide linkage and association approach to find susceptibility loci for platelet function phenotypes in European American and African American families with coronary artery disease

<p>Abstract</p> <p>Background</p> <p>The inability of aspirin (ASA) to adequately suppress platelet aggregation is associated with future risk of coronary artery disease (CAD). Heritability studies of agonist-induced platelet function phenotypes suggest that genetic variation may be responsible for ASA responsiveness. In this study, we leverage independent information from genome-wide linkage and association data to determine loci controlling platelet phenotypes before and after treatment with ASA.</p> <p>Methods</p> <p>Clinical data on 37 agonist-induced platelet function phenotypes were evaluated before and after a 2-week trial of ASA (81 mg/day) in 1231 European American and 846 African American healthy subjects with a family history of premature CAD. Principal component analysis was performed to minimize the number of independent factors underlying the covariance of these various phenotypes. Multi-point sib-pair based linkage analysis was performed using a microsatellite marker set, and single-SNP association tests were performed using markers from the Illumina 1 M genotyping chip from deCODE Genetics, Inc. All analyses were performed separately within each ethnic group.</p> <p>Results</p> <p>Several genomic regions appear to be linked to ASA response factors: a 10 cM region in African Americans on chromosome 5q11.2 had several STRs with suggestive (p-value < 7 × 10^-4) and significant (p-value < 2 × 10^-5) linkage to post aspirin platelet response to ADP, and ten additional factors had suggestive evidence for linkage (p-value < 7 × 10^-4) to thirteen genomic regions. All but one of these factors were aspirin <it>response </it>variables. While the strength of genome-wide SNP association signals for factors showing evidence for linkage is limited, especially at the strict thresholds of genome-wide criteria (N = 9 SNPs for 11 factors), more signals were considered significant when the association signal was weighted by evidence for linkage (N = 30 SNPs).</p> <p>Conclusions</p> <p>Our study supports the hypothesis that platelet phenotypes in response to ASA likely have genetic control and the combined approach of linkage and association offers an alternative approach to prioritizing regions of interest for subsequent follow-up.</p

Ancestral bias in the Hras1 gene and distal Chromosome 7 among inbred mice

Inbred strains of mice vary in their frequency of liver tumors initiated by a mutation in the Hras1 (H-ras) proto-oncogene. We sequenced 4.5 kb of the Hras1 gene on distal Chr 7 in a diverse set of 12 commonly used laboratory inbred strains of mice and detected no sequence variation to account for strain-specific differences in Hras1 mutation prevalence. Furthermore, the Hras1 sequence is essentially monoallelic for an ancestral gene derived from the M. m. domesticus species. To determine if the monoallelism and associated low rate of polymorphism are unique to Hras1 or representative of the general chromosomal locale, we extended the sequence analysis to 12 genes in the final 8 Mb of distal Chr 7. A region of at least 2.5 Mb that encompasses several genes, including Hras1 and the H19/Igf2 loci, demonstrates virtually no sequence variation. The 12 inbred strains share one dominant haplotype derived from the M. m. domesticus allele. Chromosomal regions flanking the monoallelic segment exhibit a significantly higher rate of variation and multiple haplotypes, a majority of which are attributed to M. m. domesticus or M. m. musculus ancestry

Systemic and Mucosal Immune Responses to Sublingual or Intramuscular Human Papilloma Virus Antigens in Healthy Female Volunteers

The sublingual route has been proposed as a needle-free option to induce systemic and mucosal immune protection against viral infections. In a translational study of systemic and mucosal humoral immune responses to sublingual or systemically administered viral antigens, eighteen healthy female volunteers aged 19–31 years received three immunizations with a quadravalent Human Papilloma Virus vaccine at 0, 4 and 16 weeks as sublingual drops (SL, n = 12) or intramuscular injection (IM, n = 6). IM antigen delivery induced or boosted HPV-specific serum IgG and pseudovirus-neutralizing antibodies, HPV-specific cervical and vaginal IgG, and elicited circulating IgG and IgA antibody secreting cells. SL antigens induced ∼38-fold lower serum and ∼2-fold lower cervical/vaginal IgG than IM delivery, and induced or boosted serum virus neutralizing antibody in only 3/12 subjects. Neither route reproducibly induced HPV-specific mucosal IgA. Alternative delivery systems and adjuvants will be required to enhance and evaluate immune responses following sublingual immunization in humans

The supernatural characters and powers of sacred trees in the Holy Land

This article surveys the beliefs concerning the supernatural characteristics and powers of sacred trees in Israel; it is based on a field study as well as a survey of the literature and includes 118 interviews with Muslims and Druze. Both the Muslims and Druze in this study attribute supernatural dimensions to sacred trees which are directly related to ancient, deep-rooted pagan traditions. The Muslims attribute similar divine powers to sacred trees as they do to the graves of their saints; the graves and the trees are both considered to be the abode of the soul of a saint which is the source of their miraculous powers. Any violation of a sacred tree would be strictly punished while leaving the opportunity for atonement and forgiveness. The Druze, who believe in the transmigration of souls, have similar traditions concerning sacred trees but with a different religious background. In polytheistic religions the sacred grove/forest is a centre of the community's official worship; any violation of the trees is regarded as a threat to the well being of the community. Punishments may thus be collective. In the monotheistic world (including Christianity, Islam and Druze) the pagan worship of trees was converted into the worship/adoration of saints/prophets; it is not a part of the official religion but rather a personal act and the punishments are exerted only on the violating individual

Comparison of Strategies to Detect Epistasis from eQTL Data

Genome-wide association studies have been instrumental in identifying genetic variants associated with complex traits such as human disease or gene expression phenotypes. It has been proposed that extending existing analysis methods by considering interactions between pairs of loci may uncover additional genetic effects. However, the large number of possible two-marker tests presents significant computational and statistical challenges. Although several strategies to detect epistasis effects have been proposed and tested for specific phenotypes, so far there has been no systematic attempt to compare their performance using real data. We made use of thousands of gene expression traits from linkage and eQTL studies, to compare the performance of different strategies. We found that using information from marginal associations between markers and phenotypes to detect epistatic effects yielded a lower false discovery rate (FDR) than a strategy solely using biological annotation in yeast, whereas results from human data were inconclusive. For future studies whose aim is to discover epistatic effects, we recommend incorporating information about marginal associations between SNPs and phenotypes instead of relying solely on biological annotation. Improved methods to discover epistatic effects will result in a more complete understanding of complex genetic effects

Genome-Wide Association Data Reveal a Global Map of Genetic Interactions among Protein Complexes

This work demonstrates how gene association studies can be analyzed to map a global landscape of genetic interactions among protein complexes and pathways. Despite the immense potential of gene association studies, they have been challenging to analyze because most traits are complex, involving the combined effect of mutations at many different genes. Due to lack of statistical power, only the strongest single markers are typically identified. Here, we present an integrative approach that greatly increases power through marker clustering and projection of marker interactions within and across protein complexes. Applied to a recent gene association study in yeast, this approach identifies 2,023 genetic interactions which map to 208 functional interactions among protein complexes. We show that such interactions are analogous to interactions derived through reverse genetic screens and that they provide coverage in areas not yet tested by reverse genetic analysis. This work has the potential to transform gene association studies, by elevating the analysis from the level of individual markers to global maps of genetic interactions. As proof of principle, we use synthetic genetic screens to confirm numerous novel genetic interactions for the INO80 chromatin remodeling complex

High-Throughput Sequencing of Three Lemnoideae (Duckweeds) Chloroplast Genomes from Total DNA

BACKGROUND: Chloroplast genomes provide a wealth of information for evolutionary and population genetic studies. Chloroplasts play a particularly important role in the adaption for aquatic plants because they float on water and their major surface is exposed continuously to sunlight. The subfamily of Lemnoideae represents such a collection of aquatic species that because of photosynthesis represents one of the fastest growing plant species on earth. METHODS: We sequenced the chloroplast genomes from three different genera of Lemnoideae, Spirodela polyrhiza, Wolffiella lingulata and Wolffia australiana by high-throughput DNA sequencing of genomic DNA using the SOLiD platform. Unfractionated total DNA contains high copies of plastid DNA so that sequences from the nucleus and mitochondria can easily be filtered computationally. Remaining sequence reads were assembled into contiguous sequences (contigs) using SOLiD software tools. Contigs were mapped to a reference genome of Lemna minor and gaps, selected by PCR, were sequenced on the ABI3730xl platform. CONCLUSIONS: This combinatorial approach yielded whole genomic contiguous sequences in a cost-effective manner. Over 1,000-time coverage of chloroplast from total DNA were reached by the SOLiD platform in a single spot on a quadrant slide without purification. Comparative analysis indicated that the chloroplast genome was conserved in gene number and organization with respect to the reference genome of L. minor. However, higher nucleotide substitution, abundant deletions and insertions occurred in non-coding regions of these genomes, indicating a greater genomic dynamics than expected from the comparison of other related species in the Pooideae. Noticeably, there was no transition bias over transversion in Lemnoideae. The data should have immediate applications in evolutionary biology and plant taxonomy with increased resolution and statistical power

Longitudinal Replication Studies of GWAS Risk SNPs Influencing Body Mass Index over the Course of Childhood and Adulthood

Genome-wide association studies (GWAS) have identified multiple common variants associated with body mass index (BMI). In this study, we tested 23 genotyped GWAS-significant SNPs (p-value<5*10-8) for longitudinal associations with BMI during childhood (3–17 years) and adulthood (18–45 years) for 658 subjects. We also proposed a heuristic forward search for the best joint effect model to explain the longitudinal BMI variation. After using false discovery rate (FDR) to adjust for multiple tests, childhood and adulthood BMI were found to be significantly associated with six SNPs each (q-value<0.05), with one SNP associated with both BMI measurements: KCTD15 rs29941 (q-value<7.6*10-4). These 12 SNPs are located at or near genes either expressed in the brain (BDNF, KCTD15, TMEM18, MTCH2, and FTO) or implicated in cell apoptosis and proliferation (FAIM2, MAP2K5, and TFAP2B). The longitudinal effects of FAIM2 rs7138803 on childhood BMI and MAP2K5 rs2241423 on adulthood BMI decreased as age increased (q-value<0.05). The FTO candidate SNPs, rs6499640 at the 5 ′-end and rs1121980 and rs8050136 downstream, were associated with childhood and adulthood BMI, respectively, and the risk effects of rs6499640 and rs1121980 increased as birth weight decreased. The best joint effect model for childhood and adulthood BMI contained 14 and 15 SNPs each, with 11 in common, and the percentage of explained variance increased from 0.17% and 9.0*10−6% to 2.22% and 2.71%, respectively. In summary, this study evidenced the presence of long-term major effects of genes on obesity development, implicated in pathways related to neural development and cell metabolism, and different sets of genes associated with childhood and adulthood BMI, respectively. The gene effects can vary with age and be modified by prenatal development. The best joint effect model indicated that multiple variants with effects that are weak or absent alone can nevertheless jointly exert a large longitudinal effect on BMI