Faculty Bibliography

The identification of common variants that contribute to the genesis of human inherited disorders remains a significant challenge. Hirschsprung disease (HSCR) is a multifactorial, non-mendelian disorder in which rare high-penetrance coding sequence mutations in the receptor tyrosine kinase RET contribute to risk in combination with mutations at other genes. We have used family-based association studies to identify a disease interval, and integrated this with comparative and functional genomic analysis to prioritize conserved and functional elements within which mutations can be sought. We now show that a common non-coding RET variant within a conserved enhancer-like sequence in intron 1 is significantly associated with HSCR susceptibility and makes a 20-fold greater contribution to risk than rare alleles do. This mutation reduces in vitro enhancer activity markedly, has low penetrance, has different genetic effects in males and females, and explains several features of the complex inheritance pattern of HSCR. Thus, common low-penetrance variants, identified by association studies, can underlie both common and rare diseases.

Elevated levels of the inflammatory cytokine IL-6 are associated with the development of disability, frailty, and mortality in older adults. These outcomes are likely mediated through inflammatory activity that alters hormones, skeletal muscle, and the immune system. Polymorphic variants in the IL-6 gene influence IL-6 expression. We hypothesized that IL-6 alleles associate with increased serum of IL-6, decreased muscle strength, and frailty, and tested this in the Women's Health and Aging cohorts. We genotyped 463 participants age 70-79, and identified three common IL-6 haplotype blocks for the Caucasian (n=363) and African American (n=100) subsets. Using linear and logistic regression, and adjusting for age, BMI, race, and osteoarthritis, we identified no significant or clinically meaningful relationship between any single IL-6 single nucleotide polymorphism (SNP) or any IL-6 haplotype and serum IL-6 level, grip, knee, or hip strength, or frailty. Given that the promoter SNP (rs1800795) has been reported to influence IL-6 levels and health outcomes, we performed a similar association study in the In Chianti population (n=266) and confirmed lack of association. These results suggest that IL-6 gene variation may not be an important factor in the determination of elevated IL-6 levels and related phenotypes found in older women.

OBJECTIVE: Function of the renin-angiotensin system is important to human hypertension, but its genetic etiology remains elusive. We set out to examine a hypothesis that multiple genetic variants in the system act together in blood pressure regulation, via intermediate phenotypes such as blood pressure reactivity. METHODS: A sample of 531 hypertensive cases and 417 controls was selected from the HyperGEN study. Hypertension-related traits including blood pressure responses to challenges to math test, handgrip and postural change (mathBP, gripBP, and postBP), and body mass index (BMI) were analyzed for association with 10 single nucleotide polymorphisms (SNPs) in the angiotensinogen (AGT) gene. Single-marker and haplotype analyses were performed to examine the effects of both individual and multiple variants. Multiple-trait profiling was used to assess interaction of latent intermediate factors with susceptible haplotypes. RESULTS: In Blacks, two SNPs in exon 5 and 3'UTR showed significant association with gripBP, and two promoter SNPs were strongly associated with postBP. In Whites, only borderline association was found for 2 promoter SNPs with mathBP. Haplotype analyses in Blacks confirmed association with gripBP, and detected significant association of a haplotype to BMI (p=0.029). With the interactions modeled, haplotype associations found in Blacks remain significant, while significant associations to BMI (p=0.009) and gripSBP emerged in Whites. CONCLUSION: Genetic variants in regulatory regions of AGT showed strong association with blood pressure reactivity. Interaction of promoter and genic SNPs in AGT revealed collective action of multiple variants on blood pressure reactivity and BMI both in Blacks and in Whites, possibly following different pathways.

Genome-wide disease-association mapping has been heralded as the study design of the next generation, but the lack of analytical methods to use genotype data fully is a large stumbling block. Here we describe an algorithm and statistical method that efficiently and exhaustively exploits haplotype information by subjecting alleles (a marker or contiguous sets of markers) from sliding windows of all sizes to transmission disequilibrium tests. By applying our method to simulated data and to Hirschsprung disease, we show that it can detect both common and rare disease variants of small effect. These results show that the theoretical benefits of genome-wide association studies are at last realizable.

Determining linkage phase from population samples with statistical methods is accurate only within regions of high linkage disequilibrium (LD). Yet, affected individuals in a genetic mapping study, including those involving cases and controls, may share sequences identical-by-descent stretching on the order of 10s to 100s of kilobases, quite possibly over regions of low LD in the population. At the same time, inferring phase from nuclear families may be hampered by missing family members, missing genotypes, and the noninformativity of certain genotype patterns. In this study, we reformulate our previous haplotype reconstruction algorithm, and its associated computer program, to phase parents with information derived from population samples as well as from their offspring. In applications of our algorithm to 100-kb stretches, simulated in accordance to a Wright-Fisher model with typical levels of LD in humans, we find that phase reconstruction for 160 trios with 10% missing data is highly accurate (>90%) over the entire length. Furthermore, our algorithm can estimate allelic status for missing data at high accuracy (>95%). Finally, the input capacity of the program is vast, easily handling thousands of segregating sites in > or = 1000 chromosomes.

SUMMARY: Three recent publications have examined the quality and completeness of public database single nucleotide polymorphism (dbSNP) and have come to dramatically different conclusions regarding dbSNPs false positive rate and the proportion of dbSNPs that are expected to be common. These studies employed different genotyping technologies and different protocols in determining minimum acceptable genotyping quality thresholds. Because heterozygous sites typically have lower quality scores than homozygous sites, a higher minimum quality threshold reduces the number of false positive SNPs, but yields fewer heterozygotes and leads to fewer confirmed SNPs. To account for the different confirmation rates and distributions of minor allele frequencies, we propose that the three confirmation studies have different false positive and false negative rates. We developed a mathematical model to predict SNP confirmation rates and the apparent distribution of minor allele frequencies under user-specified false positive and false negative rates. We applied this model to the three published studies and to our own resequencing effort. We conclude that the dbSNP false positive rate is approximately 15-17% and that the reported confirmation studies have vastly different genotyping error rates and patterns

Somatic mitochondrial mutations are common in human cancers, and can be used as a tool for early detection of cancer. We have developed a mitochondrial Custom Reseq microarray as an array-based sequencing platform for rapid and high-throughput analysis of mitochondrial DNA. The MitoChip contains oligonucleotide probes synthesized using standard photolithography and solid-phase synthesis, and is able to sequence >29 kb of double-stranded DNA in a single assay. Both strands of the entire human mitochondrial coding sequence (15,451 bp) are arrayed on the MitoChip; both strands of an additional 12,935 bp (84% of coding DNA) are arrayed in duplicate. We used 300 ng of genomic DNA to amplify the mitochondrial coding sequence in three overlapping long PCR fragments. We then sequenced >2 million base pairs of mitochondrial DNA, and successfully assigned base calls at 96.0% of nucleotide positions. Replicate experiments demonstrated >99.99% reproducibility. In matched fluid samples (urine and pancreatic juice, respectively) obtained from five patients with bladder cancer and four with pancreatic cancer, the MitoChip detected at least one cancer-associated mitochondrial mutation in six (66%) of nine samples. The MitoChip is a high-throughput sequencing tool for the reliable identification of mitochondrial DNA mutations from primary tumors in clinical samples.