ExPheWas: a platform for cis-Mendelian randomization and gene-based association scans.

Establishing the relationship between protein-coding genes and phenotypes has the potential to inform on the molecular etiology of diseases. Here, we describe ExPheWas (exphewas.ca), a gene-based phenome-wide association study browser and platform that enables the conduct of gene-based Mendelian randomization. The ExPheWas data repository includes sex-stratified and sex-combined gene-based association results from 26 616 genes with 1746 phenotypes measured in up to 413 133 individuals from the UK Biobank. Interactive visualizations are provided through a browser to facilitate data exploration supported by false discovery rate control, and it includes tools for enrichment analysis. The interactive Mendelian randomization module in ExPheWas allows the estimation of causal effects of a genetically predicted exposure on an outcome by using genetic variation in a single gene as the instrumental variable.

An inherited duplication at the gene p21 Protein-Activated Kinase 7 (PAK7) is a risk factor for psychosis.

Identifying rare, highly penetrant risk mutations may be an important step in dissecting the molecular etiology of schizophrenia. We conducted a gene-based analysis of large (>100 kb), rare copy-number variants (CNVs) in the Wellcome Trust Case Control Consortium 2 (WTCCC2) schizophrenia sample of 1564 cases and 1748 controls all from Ireland, and further extended the analysis to include an additional 5196 UK controls. We found association with duplications at chr20p12.2 (P = 0.007) and evidence of replication in large independent European schizophrenia (P = 0.052) and UK bipolar disorder case-control cohorts (P = 0.047). A combined analysis of Irish/UK subjects including additional psychosis cases (schizophrenia and bipolar disorder) identified 22 carriers in 11 707 cases and 10 carriers in 21 204 controls [meta-analysis Cochran-Mantel-Haenszel P-value = 2 × 10(-4); odds ratio (OR) = 11.3, 95% CI = 3.7, ∞]. Nineteen of the 22 cases and 8 of the 10 controls carried duplications starting at 9.68 Mb with similar breakpoints across samples. By haplotype analysis and sequencing, we identified a tandem ~149 kb duplication overlapping the gene p21 Protein-Activated Kinase 7 (PAK7, also called PAK5) which was in linkage disequilibrium with local haplotypes (P = 2.5 × 10(-21)), indicative of a single ancestral duplication event. We confirmed the breakpoints in 8/8 carriers tested and found co-segregation of the duplication with illness in two additional family members of one of the affected probands. We demonstrate that PAK7 is developmentally co-expressed with another known psychosis risk gene (DISC1) suggesting a potential molecular mechanism involving aberrant synapse development and plasticity.

Comparison of genotype clustering tools with rare variants.

BACKGROUND: Along with the improvement of high throughput sequencing technologies, the genetics community is showing marked interest for the rare variants/common diseases hypothesis. While sequencing can still be prohibitive for large studies, commercially available genotyping arrays targeting rare variants prove to be a reasonable alternative. A technical challenge of array based methods is the task of deriving genotype classes (homozygous or heterozygous) by clustering intensity data points. The performance of clustering tools for common polymorphisms is well established, while their performance when conducted with a large proportion of rare variants (where data points are sparse for genotypes containing the rare allele) is less known. We have compared the performance of four clustering tools (GenCall, GenoSNP, optiCall and zCall) for the genotyping of over 10,000 samples using the Illumina's HumanExome BeadChip, which includes 247,870 variants, 90% of which have a minor allele frequency below 5% in a population of European ancestry. Different reference parameters for GenCall and different initial parameters for GenoSNP were tested. Genotyping accuracy was assessed using data from the 1000 Genomes Project as a gold standard, and agreement between tools was measured. RESULTS: Concordance of GenoSNP's calls with the gold standard was below expectations and was increased by changing the tool's initial parameters. While the four tools provided concordance with the gold standard above 99% for common alleles, some of them performed poorly for rare alleles. The reproducibility of genotype calls for each tool was assessed using experimental duplicates which provided concordance rates above 99%. The inter-tool agreement of genotype calls was high for approximately 95% of variants. Most tools yielded similar error rates (approximately 0.02), except for zCall which performed better with a 0.00164 mean error rate. CONCLUSIONS: The GenoSNP clustering tool could not be run straight "out of the box" with the HumanExome BeadChip, as modification of hard coded parameters was necessary to achieve optimal performance. Overall, GenCall marginally outperformed the other tools for the HumanExome BeadChip. The use of experimental replicates provided a valuable quality control tool for genotyping projects with rare variants.

Pharmacogenomic markers of metoprolol and α-OH-metoprolol concentrations: a genome-wide association study.

Aim: Few genome-wide association studies (GWASs) have been conducted to identify predictors of drug concentrations. The authors therefore sought to discover the pharmacogenomic markers involved in metoprolol pharmacokinetics. Patients & methods: The authors performed a GWAS of a cross-sectional study of 993 patients from the Montreal Heart Institute Biobank taking metoprolol. Results: A total of 391 and 444 SNPs reached the significance threshold of 5 × 10-8 for metoprolol and α-OH-metoprolol concentrations, respectively. All were located on chromosome 22 at or near the CYP2D6 gene, encoding CYP450 2D6, metoprolol's main metabolizing enzyme. Conclusion: The results reinforce previous findings of the importance of the CYP2D6 locus for metoprolol concentrations and confirm that large biobanks can be used to identify genetic determinants of drug pharmacokinetics at a GWAS significance level.

Do variants in the coding regions of FOXP2, a gene implicated in speech disorder, confer a risk for congenital amusia?

Congenital amusia is a lifelong disorder that compromises the normal development of musical abilities in 1.5-4% of the general population. There is a substantial genetic contribution to congenital amusia, and it bears similarities to neurodevelopmental disorders of language. Here, we examine the extent to which variants in the forkhead box P2 gene (FOXP2)-the first gene to be identified as causal in developmental speech deficits-are associated with the amusic trait. Using a cohort of 49 individuals with amusia, of which 27 were unrelated, the role of FOXP2 variants in amusia was evaluated. Fourteen variants were examined in the cohort. None segregated with the amusic trait among participants for whom family information was available; nor were they predicted to be deleterious to protein function. Thus, variants in FOXP2 are not likely to cause amusia. Implications for ongoing debates about the distinction between musicality and language are discussed.

Partitioning of copy-number genotypes in pedigrees.

BACKGROUND: Copy number variations (CNVs) and polymorphisms (CNPs) have only recently gained the genetic community's attention. Conservative estimates have shown that CNVs and CNPs might affect more than 10% of the genome and that they may be at least as important as single nucleotide polymorphisms in assessing human variability. Widely used tools for CNP analysis have been implemented in Birdsuite and PLINK for the purpose of conducting genetic association studies based on the unpartitioned total number of CNP copies provided by the intensities from Affymetrix's Genome-Wide Human SNP Array. Here, we are interested in partitioning copy number variations and polymorphisms in extended pedigrees for the purpose of linkage analysis on familial data. RESULTS: We have developed CNGen, a new software for the partitioning of copy number polymorphism using the integrated genotypes from Birdsuite with the Affymetrix platform. The algorithm applied to familial trios or extended pedigrees can produce partitioned copy number genotypes with distinct parental alleles. We have validated the algorithm using simulations on a complex pedigree structure using frequencies calculated from a real dataset of 300 genotyped samples from 42 pedigrees segregating a congenital heart defect phenotype. CONCLUSIONS: CNGen is the first published software for the partitioning of copy number genotypes in pedigrees, making possible the use CNPs and CNVs for linkage analysis. It was implemented with the Python interpreter version 2.5.2. It was successfully tested on current Linux, Windows and Mac OS workstations.

Rare and low-frequency coding variants in CXCR2 and other genes are associated with hematological traits.

Hematological traits are important clinical parameters. To test the effects of rare and low-frequency coding variants on hematological traits, we analyzed hemoglobin concentration, hematocrit levels, white blood cell (WBC) counts and platelet counts in 31,340 individuals genotyped on an exome array. We identified several missense variants in CXCR2 associated with reduced WBC count (gene-based P = 2.6 × 10(-13)). In a separate family-based resequencing study, we identified a CXCR2 frameshift mutation in a pedigree with congenital neutropenia that abolished ligand-induced CXCR2 signal transduction and chemotaxis. We also identified missense or splice-site variants in key hematopoiesis regulators (EPO, TFR2, HBB, TUBB1 and SH2B3) associated with blood cell traits. Finally, we were able to detect associations between a rare somatic JAK2 mutation (encoding p.Val617Phe) and platelet count (P = 3.9 × 10(-22)) as well as hemoglobin concentration (P = 0.002), hematocrit levels (P = 9.5 × 10(-7)) and WBC count (P = 3.1 × 10(-5)). In conclusion, exome arrays complement genome-wide association studies in identifying new variants that contribute to complex human traits.