Catherine Chu

Association of common genetic variation in the insulin/IGF1 signaling pathway with human longevity.

The insulin/IGF1 signaling pathways affect lifespan in several model organisms, including worms, flies and mice. To investigate whether common genetic variation in this pathway influences lifespan in humans, we genotyped 291 common variants in 30 genes encoding proteins in the insulin/IGF1 signaling pathway in a cohort of elderly Caucasian women selected from the Study of Osteoporotic Fractures (SOF). The cohort included 293 long‐lived cases (lifespan ≥ 92 years (y), mean ± standard deviation (SD) = 95.3 ± 2.2y) and 603 average‐lifespan controls (lifespan ≤ 79y, mean = 75.7 ± 2.6y). Variants were selected for genotyping using a haplotype‐tagging approach. We found a modest excess of variants nominally associated with longevity. Nominally significant variants were then replicated in two additional Caucasian cohorts including both males and females: the Cardiovascular Health Study and Ashkenazi Jewish Centenarians. An intronic single nucleotide polymorphism in AKT1, rs3803304, was significantly associated with lifespan in a meta‐analysis across the three cohorts (OR = 0.78 95%CI = 0.68–0.89, adjusted P = 0.043); two intronic single nucleotide polymorphisms in FOXO3A demonstrated a significant lifespan association among women only (rs1935949, OR = 1.35, 95%CI = 1.15–1.57, adjusted P = 0.0093). These results demonstrate that common variants in several genes in the insulin/IGF1 pathway are associated with human lifespan.

Temporal Dissection of Tumorigenesis in Primary Cancers

Timely intervention for cancer requires knowledge of its earliest genetic aberrations. Sequencing of tumors and their metastases reveals numerous abnormalities occurring late in progression. A means to temporally order aberrations in a single cancer, rather than inferring them from serially acquired samples, would define changes preceding even clinically evident disease. We integrate DNA sequence and copy number information to reconstruct the order of abnormalities as individual tumors evolve for 2 separate cancer types. We detect vast, unreported expansion of simple mutations sharply demarcated by recombinative loss of the second copy of TP53 in cutaneous squamous cell carcinomas (cSCC) and serous ovarian adenocarcinomas, in the former surpassing 50 mutations per megabase. In cSCCs, we also report diverse secondary mutations in known and novel oncogenic pathways, illustrating how such expanded mutagenesis directly promotes malignant progression. These results reframe paradigms in which TP53 mutation is required later, to bypass senescence induced by driver oncogenes.

A hybrid approach for de novo human genome sequence assembly and phasing

Despite tremendous progress in genome sequencing, the basic goal of producing a phased (haplotype-resolved) genome sequence with end-to-end contiguity for each chromosome at reasonable cost and effort is still unrealized. In this study, we describe an approach to performing de novo genome assembly and experimental phasing by integrating the data from Illumina short-read sequencing, 10X Genomics linked-read sequencing, and BioNano Genomics genome mapping to yield a high-quality, phased, de novo assembled human genome.

Sequencing of TNFAIP3 and association of variants with multiple autoimmune diseases

The TNFAIP3 locus at 6q23, encoding A20, has been associated with multiple autoimmune diseases (AIDs). In this study, we sequence the coding portions of the gene to identify contributing causal polymorphisms that may explain some of the observed associations. A collection of 123 individuals from the Multiple Autoimmune Disease Genetics Consortium (MADGC) collection, each with multiple AIDs (mean=2.2 confirmed diagnoses), and 397 unrelated healthy controls were used for initial sequencing. A total of 32 polymorphisms were identified in the sequencing experiments, including 16 novel and 11 coding variants. Association testing in the entire MADGC collection (1,008 Caucasians with one or more AIDs and 770 unaffected family controls) revealed association of a novel intronic insertion–deletion polymorphism with rheumatoid arthritis (RA) (odds ratio (OR)=2.48, P=0.041). Genotyping of the most common coding polymorphism, rs2230926, in the MADGC collection and additional control individuals revealed a significant association with Sjögrens syndrome (OR=3.38, P=0.038), Crohns disease (OR=2.25, P=0.041), psoriasis (OR=0.037, P=0.036) and RA (OR=1.9, P=0.025). Finally, haplotype and additional testing of polymorphisms revealed that cases were enriched for 5′ and 3′ untranslated region variants (one-sided P-value=0.04), but not specifically for common (>2% minor allele frequency), rare, exonic, intronic, non-synonymous or synonymous variants.

Genome-Wide Structural Variation Detection by Genome Mapping on Nanochannel Arrays.

Comprehensive whole-genome structural variation detection is challenging with current approaches. With diploid cells as DNA source and the presence of numerous repetitive elements, short-read DNA sequencing cannot be used to detect structural variation efficiently. In this report, we show that genome mapping with long, fluorescently labeled DNA molecules imaged on nanochannel arrays can be used for whole-genome structural variation detection without sequencing. While whole-genome haplotyping is not achieved, local phasing (across >150-kb regions) is routine, as molecules from the parental chromosomes are examined separately. In one experiment, we generated genome maps from a trio from the 1000 Genomes Project, compared the maps against that derived from the reference human genome, and identified structural variations that are >5 kb in size. We find that these individuals have many more structural variants than those published, including some with the potential of disrupting gene function or regulation.

Pungent General Anesthetics Activate Transient Receptor Potential-A1 to Produce Hyperalgesia and Neurogenic Bronchoconstriction

Background:Volatile anesthetics such as isoflurane and halothane have been in clinical use for many years and represent the group of drugs most commonly used to maintain general anesthesia. However, despite their widespread use, the molecular mechanisms by which these drugs exert their effects are not completely understood. Recently, a seemingly paradoxical effect of general anesthetics has been identified: the activation of peripheral nociceptors by irritant anesthetics. This mechanism may explain the hyperalgesic actions of inhaled anesthetics and their adverse effects in the airways. Methods:To test the hypothesis that irritant inhaled anesthetics activate the excitatory ion-channel transient receptor potential (TRP)-A1 and thereby contribute to hyperalgesia and irritant airway effects, we used the measurement of intracellular calcium concentration in isolated cells in culture. For our functional experiments, we used models of isolated guinea pig bronchi to measure bronchoconstriction and withdrawal threshold to mechanical stimulation with von Frey filaments in mice. Results:Irritant inhaled anesthetics activate TRPA1 expressed in human embryonic kidney cells and in nociceptive neurons. Isoflurane induces mechanical hyperalgesia in mice by a TRPA1-dependent mechanism. Isoflurane also induces TRPA1-dependent constriction of isolated bronchi. Nonirritant anesthetics do not activate TRPA1 and fail to produce hyperalgesia and bronchial constriction. Conclusions:General anesthetics induce a reversible loss of consciousness and render the patient unresponsive to painful stimuli. However, they also produce excitatory effects such as airway irritation and they contribute to postoperative pain. Activation of TRPA1 may contribute to these adverse effects, a hypothesis that remains to be tested in the clinical setting.

Manitoba-oculo-tricho-anal (MOTA) syndrome is caused by mutations in FREM1

Background Manitoba-oculo-tricho-anal (MOTA) syndrome is a rare condition defined by eyelid colobomas, cryptophthalmos and anophthalmia/microphthalmia, an aberrant hairline, a bifid or broad nasal tip, and gastrointestinal anomalies such as omphalocele and anal stenosis. Autosomal recessive inheritance had been assumed because of consanguinity in the Oji-Cre population of Manitoba and reports of affected siblings, but no locus or cytogenetic aberration had previously been described. Methods and results This study shows that MOTA syndrome is caused by mutations in FREM1, a gene previously mutated in bifid nose, renal agenesis, and anorectal malformations (BNAR) syndrome. MOTA syndrome and BNAR syndrome can therefore be considered as part of a phenotypic spectrum that is similar to, but distinct from and less severe than, Fraser syndrome. Re-examination of Frem1bat/bat mutant mice found new evidence that Frem1 is involved in anal and craniofacial development, with anal prolapse, eyelid colobomas, telecanthus, a shortened snout and reduced philtral height present in the mutant mice, similar to the human phenotype in MOTA syndrome. Conclusions The milder phenotypes associated with FREM1 deficiency in humans (MOTA syndrome and BNAR syndrome) compared to that resulting from FRAS1 and FREM2 loss of function (Fraser syndrome) are also consistent with the less severe phenotypes resulting from Frem1 loss of function in mice. Together, Fraser, BNAR and MOTA syndromes constitute a clinically overlapping group of FRAS–FREM complex diseases.

Mitochondrial DNA Sequence Variation and Risk of Pancreatic Cancer

Although the mitochondrial genome exhibits high mutation rates, common mitochondrial DNA (mtDNA) variation has not been consistently associated with pancreatic cancer. Here, we comprehensively examined mitochondrial genomic variation by sequencing the mtDNA of participants (cases = 286, controls = 283) in a San Francisco Bay Area pancreatic cancer case-control study. Five common variants were associated with pancreatic cancer at nominal statistical significance (P < 0.05) with the strongest finding for mt5460g in the ND2 gene [OR = 3.9; 95% confidence interval (CI), 1.5-10; P = 0.004] which encodes an A331T substitution. Haplogroup K was nominally associated with reduced pancreatic cancer risk (OR = 0.32; 95% CI, 0.13-0.76; P = 0.01) when compared with the most common haplogroup, H. A total of 19 haplogroup-specific rare variants yielded nominal statistically significant associations (P < 0.05) with pancreatic cancer risk, with the majority observed in genes involved in oxidative phosphorylation. Weighted-sum statistics were used to identify an aggregate effect of variants in the 22 mitochondrial tRNAs on pancreatic cancer risk (P = 0.02). While the burden of singleton variants in the HV2 and 12S RNA regions was three times higher among European haplogroup N cases than controls, the prevalence of singleton variants in ND4 and ND5 was two to three times higher among African haplogroup L cases than in controls. Together, the results of this study provide evidence that aggregated common and rare variants and the accumulation of singleton variants are important contributors to pancreatic cancer risk.

Direct determination of haplotypes from single DNA molecules

Determining the long-range haplotypes in a diploid individual is a major technical challenge. Here we report a method of molecular haplotyping by directly imaging multiple polymorphic sites on individual human DNA molecules simultaneously. We demonstrate the utility of this technology by accurately determining the haplotypes consisting of up to 16 single-nucleotide polymorphisms in genomic regions up to 50 kilobases.

OATP1B1-related drug–drug and drug–gene interactions as potential risk factors for cerivastatin-induced rhabdomyolysis

Objective Genetic variation in drug metabolizing enzymes and membrane transporters as well as concomitant drug therapy can modulate the beneficial and the deleterious effects of drugs. We investigated whether patients exhibiting rhabdomyolysis who were taking cerivastatin possess functional genetic variants in SLCO1B1 and whether they were on concomitant medications that inhibit OATP1B1, resulting in accumulation of cerivastatin. Methods This study had three components: (a) resequencing the SLCO1B1 gene in 122 patients who developed rhabdomyolysis while on cerivastatin; (b) functional evaluation of the identified SLCO1B1 nonsynonymous variants and haplotypes in in-vitro HEK293/FRT cells stably transfected with pcDNA5/FRT empty vector, SLCO1B1 reference, variants, and haplotypes; and (c) in-vitro screening of 15 drugs commonly used among the rhabdomyolysis cases for inhibition of OATP1B1-mediated uptake of cerivastatin in HEK293/FRT cells stably transfected with reference SLCO1B1. Results The resequencing of the SLCO1B1 gene identified 54 variants. In-vitro functional analysis of SLCO1B1 nonsynonymous variants and haplotypes showed that the V174A, R57Q, and P155T variants, a novel frameshift insertion, OATP1B1*14 and OATP1B1*15 haplotype were associated with a significant reduction (P<0.001) in cerivastatin uptake (32, 18, 72, 3.4, 2.1 and 5.7% of reference, respectively). Furthermore, clopidogrel and seven other drugs were shown to inhibit OATP1B1-mediated uptake of cerivastatin. Conclusion Reduced function of OATP1B1 related to genetic variation and drug–drug interactions likely contributed to cerivastatin-induced rhabdomyolysis. Although cerivastatin is no longer in clinical use, these findings may translate to related statins and other substrates of OATP1B1.