Faculty Bibliography

Coffee consumption is a model for addictive behavior. We performed a meta-analysis of genome-wide association studies (GWASs) on coffee intake from 8 Caucasian cohorts (N=18 176) and sought replication of our top findings in a further 7929 individuals. We also performed a gene expression analysis treating different cell lines with caffeine. Genome-wide significant association was observed for two single-nucleotide polymorphisms (SNPs) in the 15q24 region. The two SNPs rs2470893 and rs2472297 (P-values=1.6 × 10(-11) and 2.7 × 10(-11)), which were also in strong linkage disequilibrium (r(2)=0.7) with each other, lie in the 23-kb long commonly shared 5' flanking region between CYP1A1 and CYP1A2 genes. CYP1A1 was found to be downregulated in lymphoblastoid cell lines treated with caffeine. CYP1A1 is known to metabolize polycyclic aromatic hydrocarbons, which are important constituents of coffee, whereas CYP1A2 is involved in the primary metabolism of caffeine. Significant evidence of association was also detected at rs382140 (P-value=3.9 × 10(-09)) near NRCAM-a gene implicated in vulnerability to addiction, and at another independent hit rs6495122 (P-value=7.1 × 10(-09))-an SNP associated with blood pressure-in the 15q24 region near the gene ULK3, in the meta-analysis of discovery and replication cohorts. Our results from GWASs and expression analysis also strongly implicate CAB39L in coffee drinking. Pathway analysis of differentially expressed genes revealed significantly enriched ubiquitin proteasome (P-value=2.2 × 10(-05)) and Parkinson's disease pathways (P-value=3.6 × 10(-05)).

Next-generation sequencing (NGS) technologies can be a boon to human mutation detection given their high throughput: consequently, many genes and samples may be simultaneously studied with high coverage for accurate detection of heterozygotes. In circumstances requiring the intensive study of a few genes, particularly in clinical applications, a rapid turn around is another desirable goal. To this end, we assessed the performance of the bench-top 454 GS Junior platform as an optimized solution for mutation detection by amplicon sequencing of three type 3 semaphorin genes SEMA3A, SEMA3C, and SEMA3D implicated in Hirschsprung disease (HSCR). We performed mutation detection on 39 PCR amplicons totaling 14,014 bp in 47 samples studied in pools of 12 samples. Each 10-hr run was able to generate approximately 75,000 reads and approximately 28 million high-quality bases at an average read length of 371 bp. The overall sequencing error was 0.26 changes per kb at a coverage depth of >/=20 reads. Altogether, 37 sequence variants were found in this study of which 10 were unique to HSCR patients. We identified five missense mutations in these three genes that may potentially be involved in the pathogenesis of HSCR and need to be studied in larger patient samples.

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex and heterogeneous disease with an incidence of 1 in 5000 live births. Despite the multifactorial determination of HSCR in the vast majority of cases, there is a monogenic subgroup for which private rare RET coding sequence mutations with high penetrance are found (45% of HSCR familial cases). An asymmetrical parental origin is observed for RET coding sequence mutations with a higher maternal inheritance. A parent-of-origin effect is usually assumed. Here we show that a differential reproductive rate for males and females also leads to an asymmetrical parental origin, which was never considered as a possible explanation till now. In the case of HSCR, we show a positive association between penetrance of the mutation and parental transmission asymmetry: no parental transmission asymmetry is observed in sporadic RET CDS mutation carrier cases for which penetrance of the mutation is low, whereas a parental transmission asymmetry is observed in affected sib-pairs for which penetrance of the mutation is higher. This allows us to conclude that the explanation for this parental asymmetry is that more severe mutations have resulted in a differential reproductive rate between male and female carriers.

Atrial fibrillation is a highly prevalent arrhythmia and a major risk factor for stroke, heart failure and death. We conducted a genome-wide association study (GWAS) in individuals of European ancestry, including 6,707 with and 52,426 without atrial fibrillation. Six new atrial fibrillation susceptibility loci were identified and replicated in an additional sample of individuals of European ancestry, including 5,381 subjects with and 10,030 subjects without atrial fibrillation (P < 5 x 10(-8)). Four of the loci identified in Europeans were further replicated in silico in a GWAS of Japanese individuals, including 843 individuals with and 3,350 individuals without atrial fibrillation. The identified loci implicate candidate genes that encode transcription factors related to cardiopulmonary development, cardiac-expressed ion channels and cell signaling molecules.

We describe methods for rapid sequencing of the entire human mitochondrial genome (mtgenome), which involve long-range PCR for specific amplification of the mtgenome, pyrosequencing, quantitative mapping of sequence reads to identify sequence variants and heteroplasmy, as well as de novo sequence assembly. These methods have been used to study 40 publicly available HapMap samples of European (CEU) and African (YRI) ancestry to demonstrate a sequencing error rate <5.63x10(-4), nucleotide diversity of 1.6x10(-3) for CEU and 3.7x10(-3) for YRI, patterns of sequence variation consistent with earlier studies, but a higher rate of heteroplasmy varying between 10% and 50%. These results demonstrate that next-generation sequencing technologies allow interrogation of the mitochondrial genome in greater depth than previously possible which may be of value in biology and medicine.