Yuying Xie

Genetic analysis of complex traits in the emerging Collaborative Cross

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

Mitochondrial aldehyde dehydrogenase-2 (ALDH2) Glu504Lys polymorphism contributes to the variation in efficacy of sublingual nitroglycerin

Glyceryl trinitrate (GTN), also known as nitroglycerin, has been used to treat angina and heart failure for more than 130 years. Recently, it was shown that mitochondrial aldehyde dehydrogenase-2 (ALDH2) is responsible for formation of NO, the metabolite needed for GTN efficacy. In the present study, we show that the common G-to-A polymorphism in exon 12 of ALDH2--resulting in a Glu504Lys replacement that virtually eliminates ALDH2 activity in both heterozygotes and homozygotes--is associated with a lack of efficacy of sublingual GTN in Chinese subjects. We also show that the catalytic efficiency (Vmax/Km) of GTN metabolism of the Glu504 protein is approximately 10-fold higher than that of the Lys504 enzyme. We conclude that the presence of the Lys504 allele contributes in large part to the lack of an efficacious clinical response to nitroglycerin; we recommend that this genetic factor be considered when administering nitroglycerin to patients, especially Asians, 30-50% of whom possess the inactive ALDH2*2 mutant allele.

Analyses of allele-specific gene expression in highly divergent mouse crosses identifies pervasive allelic imbalance

Complex human traits are influenced by variation in regulatory DNA through mechanisms that are not fully understood. Because regulatory elements are conserved between humans and mice, a thorough annotation of cis regulatory variants in mice could aid in further characterizing these mechanisms. Here we provide a detailed portrait of mouse gene expression across multiple tissues in a three-way diallel. Greater than 80% of mouse genes have cis regulatory variation. Effects from these variants influence complex traits and usually extend to the human ortholog. Further, we estimate that at least one in every thousand SNPs creates a cis regulatory effect. We also observe two types of parent-of-origin effects, including classical imprinting and a new global allelic imbalance in expression favoring the paternal allele. We conclude that, as with humans, pervasive regulatory variation influences complex genetic traits in mice and provide a new resource toward understanding the genetic control of transcription in mammals.

Effect of μ-Opioid Receptor Gene Polymorphisms on Heroin-Induced Subjective Responses in a Chinese Population

BACKGROUND Genetic factors that influence subjective responses to drug use (such as euphoria) contribute to the risk of addiction. mu-opioid receptor is the molecular target of heroin mediating its effects in both pain relief and euphoria. METHODS To evaluate the association of mu-opioid receptor gene (OPRM1) variants with heroin-induced positive responses on first use, we studied 336 Chinese Han heroin addicts recruited in Shanghai and divided heroin addicts into two groups (positive vs. negative) according to the self-reporting feeling on first use. Association analyses with the genotypes and alleles in nine tagging single nucleotide polymorphisms (tSNPs) in OPRM1 with subjective responses were performed. Similar analysis with haplotypes of these tSNPs was also performed. RESULTS Allele frequencies of three tSNPs were significantly different between the positive and negative groups. They were rs696522 (odds ratio [OR] = 3.06, p = .0013), rs1381376 (OR = 3.16, p = .0008), and rs3778151 (OR = 3.12, p = .0004). Such association remains after adjustment for demographic covariates and for multiple testing. The subjects with heroin-induced positive responses on first use consumed more drugs than the negative group (Mann-Whitney U = 224.0, Wilcoxon W = 16334.0, p <or= .0001). CONCLUSIONS Self-reported positive responses on first use of heroin were found to be associated with OPRM1. The findings suggest that heroin-induced positive responses are likely associated with more heroin consumption.

The Effect of Dopamine D2, D5 Receptor and Transporter (SLC6A3) Polymorphisms on the Cue-Elicited Heroin Craving in Chinese

Heroin dependence is resulted from the interaction between multiple genetic and environmental factors. Subjective craving is considered to be a central phenomenon, which contributes to the continuation of drug use in active abuser and the occurrence of relapse in detoxified abusers. Dopamine pathway has been implicated in the cue‐elicited craving for a variety of addictive substances. The objective of this study was to test the hypothesis that heroin addicts carrying specific variants in dopamine‐related genes would have higher levels of craving following exposure to a heroin‐related cue. Craving induced by a series of exposure to heroin‐related cue was assessed in a cohort of Chinese heroin abuser (n = 420) recruited from natural abstinence center at Shanghai. Significantly stronger cue‐elicited heroin craving was found in individuals carrying D2 dopamine receptor gene (DRD2) TaqI RFLP A1 allele than the non‐carriers (P < 0.001). Furthermore, we did not observed significant association of cue‐elicited craving with the nine‐repeat allelic variants in dopamine transporter gene (DAT) SLC6A3 and with the dinucleotide repeat polymorphism (DRP) 148bp allele in D5 dopamine receptor gene (DRD5). The results of our study suggest that human dopamine pathway be involved in cue‐induced heroin craving, and indicate a potential genetic risk factor for persistent heroin behavior and relapse.

COE: A general approach for efficient genome-wide two-locus epistasis test in disease association study

The availability of high-density single nucleotide polymorphisms (SNPs) data has made genome-wide association study computationally challenging. Two-locus epistasis (gene-gene interaction) detection has attracted great research interest as a promising method for genetic analysis of complex diseases. In this article, we propose a general approach, COE, for efficient large scale gene-gene interaction analysis, which supports a wide range of tests. In particular, we show that many commonly used statistics are convex functions. From the observed values of the events in two-locus association test, we can develop an upper bound of the test value. Such an upper bound only depends on single-locus test and the genotype of the SNP-pair. We thus group and index SNP-pairs by their genotypes. This indexing structure can benefit the computation of all convex statistics. Utilizing the upper bound and the indexing structure, we can prune most of the SNP-pairs without compromising the optimality of the result. Our approach is especially efficient for large permutation test. Extensive experiments demonstrate that our approach provides orders of magnitude performance improvement over the brute force approach.

Using the emerging Collaborative Cross to probe the immune system

The Collaborative Cross (CC) is an emerging panel of recombinant inbred (RI) mouse strains. Each strain is genetically distinct but all descended from the same eight inbred founders. In 66 strains from incipient lines of the CC (pre-CC), as well as the 8 CC founders and some of their F1 offspring, we examined subsets of lymphocytes and antigen-presenting cells. We found significant variation among the founders, with even greater diversity in the pre-CC. Genome-wide association using inferred haplotypes detected highly significant loci controlling B-to-T cell ratio, CD8 T-cell numbers, CD11c and CD23 expression. Comparison of overall strain effects in the CC founders with strain effects at QTL in the pre-CC revealed sharp contrasts in the genetic architecture of two traits with significant loci: variation in CD23 can be explained largely by additive genetics at one locus, whereas variation in B-to-T ratio has a more complex etiology. For CD23, we found a strong QTL whose confidence interval contained the CD23 structural gene Fcer2a. Our data on the pre-CC demonstrate the utility of the CC for studying immunophenotypes and the value of integrating founder, CC and F1 data. The extreme immunophenotypes observed could have pleiotropic effects in other CC experiments.

Fine-mapping diabetes-related traits, including insulin resistance, in heterogeneous stock rats

Type 2 diabetes (T2D) is a disease of relative insulin deficiency resulting from both insulin resistance and beta cell failure. We have previously used heterogeneous stock (HS) rats to fine-map a locus for glucose tolerance. We show here that glucose intolerance in the founder strains of the HS colony is mediated by different mechanisms: insulin resistance in WKY and an insulin secretion defect in ACI, and we demonstrate a high degree of variability for measures of insulin resistance and insulin secretion in HS rats. As such, our goal was to use HS rats to fine-map several diabetes-related traits within a region on rat chromosome 1. We measured blood glucose and plasma insulin levels after a glucose tolerance test in 782 male HS rats. Using 97 SSLP markers, we genotyped a 68 Mb region on rat chromosome 1 previously implicated in glucose and insulin regulation. We used linkage disequilibrium mapping by mixed model regression with inferred descent to identify a region from 198.85 to 205.9 that contains one or more quantitative trait loci (QTL) for fasting insulin and a measure of insulin resistance, the quantitative insulin sensitivity check index. This region also encompasses loci identified for fasting glucose and Insulin_AUC (area under the curve). A separate <3 Mb QTL was identified for body weight. Using a novel penalized regression method we then estimated effects of alternative haplotype pairings under each locus. These studies highlight the utility of HS rats for fine-mapping genetic loci involved in the underlying causes of T2D.

Candidate risk factors and mechanisms for tolvaptan-induced liver injury are identified using a collaborative cross approach

Clinical trials of tolvaptan showed it to be a promising candidate for the treatment of Autosomal Dominant Polycystic Kidney Disease (ADPKD) but also revealed potential for idiosyncratic drug-induced liver injury (DILI) in this patient population. To identify risk factors and mechanisms underlying tolvaptan DILI, 8 mice in each of 45 strains of the genetically diverse Collaborative Cross (CC) mouse population were treated with a single oral dose of either tolvaptan or vehicle. Significant elevations in plasma alanine aminotransferase (ALT) were observed in tolvaptan-treated animals in 3 of the 45 strains. Genetic mapping coupled with transcriptomic analysis in the liver was used to identify several candidate susceptibility genes including epoxide hydrolase 2, interferon regulatory factor 3, and mitochondrial fission factor. Gene pathway analysis revealed that oxidative stress and immune response pathways were activated in response to tolvaptan treatment across all strains, but genes involved in regulation of bile acid homeostasis were most associated with tolvaptan-induced elevations in ALT. Secretory leukocyte peptidase inhibitor (Slpi) mRNA was also induced in the susceptible strains and was associated with increased plasma levels of Slpi protein, suggesting a potential serum marker for DILI susceptibility. In summary, tolvaptan induced signs of oxidative stress, mitochondrial dysfunction, and innate immune response in all strains, but variation in bile acid homeostasis was most associated with susceptibility to the liver response. This CC study has indicated potential mechanisms underlying tolvaptan DILI and biomarkers of susceptibility that may be useful in managing the risk of DILI in ADPKD patients.

Telltale tumor infiltrating lymphocytes (TIL) in oral, head & neck cancer

Evidence gleaned from recent studies on the role of tumor-infiltrating lymphocytes (TILs) suggests that cancer is not only a genetic disease but also an immunologic disease. Head and Neck Squamous Cell Carcinoma (HNSCC) has been a significant model to study cancer cell-immune cell interactions. First, immune cell infiltration is an important feature of these tumors. Second, HNSCC frequently develops resistance to immunogenic cytotoxicity, which provides a window to decipher how tumors engage the immune system to establish immune tolerance. Finally, chemoradiation therapy, as a central modality for HNSCC treatment, has been shown to elicit immune activation. The presence of effector immune cells in the tumor microenvironment is often associated with superior clinical response to adjuvant therapy. On the other hand, an activated immune system, in addition to limiting tumor initiation and progression, could also exert selective pressure to promote the growth of less immunogenic tumors, as a pivotal immunoediting process. But it remains unclear how cancer cell signaling regulates tumor immunogenicity and how to mitigate HNSCC-potentiated TIL suppression. In this review, we will revisit the prognostic role of TILs in HNSCC, and collectively discuss how cancer cell machinery impacts upon the plasticity of TILs.