Single nucleotide polymorphism (SNP) selection for genotype-phenotype association studies

With hundreds of single nucleotide polymorphisms (SNPs) in a candidate gene and millions of SNPs across the genome, selecting an informative subset of SNPs to maximize the ability to detect genotype-phenotype association is of great interest and importance. In addition, with a large number of SNPs, analytic methods are needed that allow investigators to control the false positive rate resulting from large numbers of SNP genotype-phenotype analyses. This dissertation uses simulated data to explore methods for selecting SNPs for genotype-phenotype association studies. I examined the pattern of linkage disequilibrium (LD) across a candidate gene region and used this pattern to aid in localizing a disease-influencing mutation. The results indicate that the r2 measure of linkage disequilibrium is preferred over the common D′ measure for use in genotype-phenotype association studies. Using step-wise linear regression, the best predictor of the quantitative trait was not usually the single functional mutation. Rather it was a SNP that was in high linkage disequilibrium with the functional mutation. Next, I compared three strategies for selecting SNPs for application to phenotype association studies: based on measures of linkage disequilibrium, based on a measure of haplotype diversity, and random selection. The results demonstrate that SNPs selected based on maximum haplotype diversity are more informative and yield higher power than randomly selected SNPs or SNPs selected based on low pair-wise LD. The data also indicate that for genes with small contribution to the phenotype, it is more prudent for investigators to increase their sample size than to continuously increase the number of SNPs in order to improve statistical power. When typing large numbers of SNPs, researchers are faced with the challenge of utilizing an appropriate statistical method that controls the type I error rate while maintaining adequate power. We show that an empirical genotype based multi-locus global test that uses permutation testing to investigate the null distribution of the maximum test statistic maintains a desired overall type I error rate while not overly sacrificing statistical power. The results also show that when the penetrance model is simple the multi-locus global test does as well or better than the haplotype analysis. However, for more complex models, haplotype analyses offer advantages. The results of this dissertation will be of utility to human geneticists designing large-scale multi-locus genotype-phenotype association studies. ^

An Sp1/Sp3 site polymorphism associated with hypermethylation of the candidate tumor suppressor gene RIL in cancer

Gene silencing due to promoter methylation is an alternative to mutations and deletions, which inactivate tumor suppressor genes (TSG) in cancer. We identified RIL by Methylated CpG Island Amplification technique as a novel aberrantly methylated gene. RIL is expressed in normal tissues and maps to the 5q31 region, frequently deleted in leukemias. We found methylation of RIL in 55/80 (69%) cancer cell lines, with highest methylation in leukemia and colon. We also observed methylation in 46/80 (58%) primary tumors, whereas normal tissues showed substantially lower degrees of methylation. RIL expression was lost in 13/16 cancer cell lines and was restored by demethylating agent. Screening of 38 cell lines and 13 primary cancers by SSCP revealed no mutations in RIL, suggesting that methylation and LOH are the primary inactivation mechanisms. Stable transfection of RIL into colorectal cancer cells resulted in reduction in cell growth, clonogenicity, and increased apoptosis upon UVC treatment, suggesting that RIL is a good candidate TSG. ^ In searching for a cause of RIL hypermethylation, we identified a 12-bp polymorphic sequence around the transcription start site of the gene that creates a long allele containing 3CTC repeat. Evolutionary studies suggested that the long allele appeared late in evolution due to insertion. Using bisulfite sequencing, in cancers heterozygous for RIL, we found that the short allele is 4.4-fold more methylated than the long allele (P = 0.003). EMSA results suggested binding of factor(s) to the inserted region of the long allele, but not to the short. EMSA mutagenesis and competition studies, as well as supershifts using nuclear extracts or recombinant Sp1 strongly indicated that those DNA binding proteins are Sp1 and Sp3. Transient transfections of RIL allele-specific expression constructs showed less than 2-fold differences in luciferase activity, suggesting no major effects of the additional Sp1 site on transcription. However, stable transfection resulted in 3-fold lower levels of transcription from the short allele 60 days post-transfection, consistent with the concept that the polymorphic Sp1 site protects against time-dependent silencing. Thus, an insertional polymorphism in the RIL promoter creates an additional Sp1/Sp3 site, which appears to protect it from silencing and methylation in cancer. ^

Identification of genetic variation influencing apolipoprotein E levels: Follow-up of genome-wide linkage scans in the GENOA study

Apolipoprotein E (ApoE) plays a major role in the metabolism of high density and low density lipoproteins (HDL and LDL). Its common protein isoforms (E2, E3, E4) are risk factors for coronary artery disease (CAD) and explain between 16 to 23% of the inter-individual variation in plasma apoE levels. Linkage analysis has been completed for plasma apoE levels in the GENOA study (Genetic Epidemiology Network of Atherosclerosis). After stratification of the population by lipoprotein levels and body mass index (BMI) to create more homogeneity with regard to biological context for apoE levels, Hispanic families showed significant linkage on chromosome 17q for two strata (LOD=2.93 at 104 cM for a low cholesterol group, LOD=3.04 at 111 cM for a low cholesterol, high HDLC group). Replication of 17q linkage was observed for apoB and apoE levels in the unstratified Hispanic and African-American populations, and for apoE levels in African-American families. Replication of this 17q linkage in different populations and strata provides strong support for the presence of gene(s) in this region with significant roles in the determination of inter-individual variation in plasma apoE levels. Through a positional and functional candidate gene approach, ten genes were identified in the 17q linked region, and 62 polymorphisms in these genes were genotyped in the GENOA families. Association analysis was performed with FBAT, GEE, and variance-component based tests followed by conditional linkage analysis. Association studies with partial coverage of TagSNPs in the gene coding for apolipoprotein H (APOH) were performed, and significant results were found for 2 SNPs (APOH_20951 and APOH_05407) in the Hispanic low cholesterol strata accounting for 3.49% of the inter-individual variation in plasma apoE levels. Among the other candidate genes, we identified a haplotype block in the ACE1 gene that contains two major haplotypes associated with apoE levels as well as total cholesterol, apoB and LDLC levels in the unstratified Hispanic population. Identifying genes responsible for the remaining 60% of inter-individual variation in plasma apoE level, will yield new insights into the understanding of genetic interactions involved in the lipid metabolism, and a more precise understanding of the risk factors leading to CAD. ^

G870A polymorphism of cyclin D1 gene influences cancer risk - From epidemiological study to mechanism analysis

Alternate splicing of the cyclin D1 gene gives rise to transcript a and b which encode two protein isoforms cyclin D1a and cyclin D1b. Through testing transcript a and transcript b in a series of human samples, we found that cyclin D1 transcript b is ubiquitously expressed as transcript a but in the lower abundance compared to transcript a. Epidemiological studies have reported that the cyclin D1 gene (CCND1) G870A polymorphism influences the risk for a variety of cancer. In this investigation, we examined the cyclin D1b levels in tumor samples with different genotypes and found that higher levels of cyclin D1b are expressed from the A allele than the G allele. Cyclin D1 is known as a cell cycle regulator facilitating the progression of the cell cycle from G1 to S phase in response to the mitogenic signals. It also interacts with several transcription factors and transcriptional coregulators to modulate their activities. It has been reported that cyclin D1a can substitute for estrogen to activate estrogen receptor α (ERα) mediated transcription and can induce the proliferation of estrogen responsive tissues. However the biological role of cyclin D1b in ERα transcriptional regulation has not been previously explored. In this study, we determined that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription. Cell proliferation assays provided the evidence that cyclin D1b negatively regulates estrogen responsive breast cancer cell growth. Taken together, our findings show that the CCND1 G870A polymorphism is correlated with increased levels of cyclin D1b and that cyclin D1b antagonizes the action of cyclin D1a on ERα mediated transcription providing evidence for the mechanism by which the CCND1 G870A polymorphism may be protective in certain types of breast cancer. ^

THE NEW ROLE OF PROPROTEIN CONVERTASE SUBTILISIN/KEXIN TYPE 9: A CONNECTION OF PROPROTEIN CONVERTASE SUBTILISIN/KEXIN TYPE 9, APOLIPOPROTEIN B, and AUTOPHAGY

Plasma low-density lipoprotein (LDL) levels are positively correlated with the incidence of coronary artery disease. In the circulation, the plasma LDL clearance is mainly achieved by the uptake via LDL receptor (LDLR). Proprotein convertase subtilisin/kexin type 9 (PCSK9) is a newly discovered gene, playing an important role in LDL metabolism. Gain-of-function mutations of PCSK9 lead to hypercholesterolemia and loss-of-function mutations of PCSK9 are associated with decrease of LDL cholesterol. The effects of PCSK9 on cholesterol levels are the consequence of a strong interaction between the catalytic domain of PCSK9 and epidermal growth factor-like repeat A (EGF-A) domain of LDLR on the cell surface of hepatocytes. This PCSK9/LDLR complex enters the cell via endocytosis, where both PCSK9 and LDLR are removed via the lysosome pathway, resulting in decreased levels of LDLR and accumulation of LDL in the plasma. However, whether this is the exclusive function of PCSK9 on LDL metabolism was challenged by us; we observed PCSK9 interacted with apolipoprotein B (apoB) and increased apoB production, irrespective of the LDLR. ApoB is the primary structure protein of LDL particle and it also serves as the ligand for the LDL receptor. There is ample evidence showing that the levels of apoB are a better indicator for heart disease than either total cholesterol or LDL cholesterol levels. We used a second-generation adenoviral vector to overexpress PCSK9 (Ad-PCSK9) in wild-type C57BL/6 and LDLR deficient mice (Ldlr-/- and Ldlr-/-Apobec1-/-). Our study revealed that overexpression of PCSK9 promoted the production and secretion of apoB in the form of very-low density lipoprotein (VLDL), which is the precursor of LDL, in the 3 mouse models studied (C57BL/6J, Ldlr-/-, and Ldlr-/-Apobec1-/-). The increased apoB production in mice was regulated at post-transcriptional levels, since there was no difference in apoB mRNA levels between mice treated with Ad-PCSK9 and control vector Ad-Null. By using pulse-chase experiment on primary hepatocytes, we showed that overexpression of PCSK9 increased the secretion of apoB, independent of LDLR. In the circulation, we showed that PCSK9 was associated with LDL particles. By using 3 different protein–protein interaction assays of co-immunoprecipitation, mammalian two-hybrid system, and in situ proximity ligation assay, we demonstrated a direct protein–protein interaction between PCSK9 and apoB. The impact of this interaction inhibited the physiological removal process of apoB via autophagosome/lysosome pathway in an LDLR-independent fashion, resulting in increased production and secretion of apoB-containing lipoproteins. The significance of this process was shown in the Pcsk9 knockout mice in the background of Ldlr-/-Apobec1-/- mice (triple knockout mice); in the absence of Pcsk9 (triple knockout mice) the levels of cholesterol, triacylglycerol, and apoB decreased significantly in comparison to that of Ldlr-/-Apobec1-/- mice. Taken together, our study demonstrated a direct intracellular interaction of PCSK9 with apoB, resulting in the inhibition of apoB degradation via the autophagosome/lysosome pathway independent of LDLR. This discovery provides a new concept of the importance of PCSK9 and suggests new approaches for the therapeutic intervention of hyperlipidemia.

HKT2;2/1, a K+-permeable transporter identified in a salt tolerant rice cultivar through surveys of natural genetic polymorphism

We have investigated OsHKT2;1 natural variation in a collection of 49 cultivars with different levels of salt tolerance and geographical origins. The effect of identified polymorphism on OsHKT2;1 activity was analysed through heterologous expression of variants in Xenopus oocytes. OsHKT2;1 appeared to be a highly conserved protein with only five possible amino acid substitutions that have no substantial effect on functional properties. Our study, however, also identified a new HKT isoform, No-OsHKT2;2/1 in Nona Bokra, a highly salt-tolerant cultivar. No-OsHKT2;2/1 probably originated from a deletion in chromosome 6, producing a chimeric gene. Its 5¢ region corresponds to that of OsHKT2;2, whose full-length sequence is not present in Nipponbare but has been identified in Pokkali, a salt-tolerant rice cultivar. Its 3¢ region corresponds to that of OsHKT2;1. No-OsHKT2;2/1 is essentially expressed in roots and displays a significant level of expression at high Na+ concentrations, in contrast to OsHKT2;1. Expressed in Xenopus oocytes or in Saccharomyces cerevisiae, No-OsHKT2;2/1 exhibited a strong permeability to Na+ and K+, even at high external Na+ concentrations, like OsHKT2;2, and in contrast to OsHKT2;1. Our results suggest that No-OsHKT2;2/1 can contribute to Nona Bokra salt tolerance by enabling root K+ uptake under saline conditions.

Dopamine receptor 4 promoter polymorphism modulates memory and neuronal responses to salience

Animal models and human functional imaging data implicate the dopamine system in mediating enhanced encoding of novel stimuli into human memory. A separate line of investigation suggests an association between a functional polymorphism in the promoter region for the human dopamine 4 receptor gene (DRD4) and sensitivity to novelty. We demonstrate, in two independent samples, that the -521Cmayor queT DRD4 promoter polymorphism determines the magnitude of human memory enhancement for contextually novel, perceptual oddball stimuli in an allele dose-dependent manner. The genotype-dependent memory enhancement conferred by the C allele is associated with increased neuronal responses during successful encoding of perceptual oddballs in the ventral striatum, an effect which is again allele dose-dependent. Furthermore, with repeated presentations of oddball stimuli, this memory advantage decreases, an effect mirrored by adaptation of activation in the hippocampus and substantia nigra/ventral tegmental area in C carriers only. Thus, a dynamic modulation of human memory enhancement for perceptually salient stimuli is associated with activation of a dopaminergic-hippocampal system, which is critically dependent on a functional polymorphism in the DRD4 promoter region.

MCT1 T1470A polymorphism influences adherence to strength training in overweight and obese men following a weight loss program

The monocarboxylate transporter (MCT) family member MCT1 transports lactate into and out of myocytes. Oxidative cells import lactate through MCT1 as a substrate, being the role of MCT1 in glycolysis-derived lactate efflux less clear. MCT1 T1470A polymorphism (rs1049434), which has been related with lactate metabolism and sports specialty 1, 2, could be an influencing factor for exercise adherence. Therefore the aim of this study was to relate the adherence to different training modalities with the T1470A MCT1 polymorphism in overweight and obese men following a weight loss program (WLP).

Low density lipoprotein receptor-negative mice expressing human apolipoprotein B-100 develop complex atherosclerotic lesions on a chow diet: No accentuation by apolipoprotein(a)

We have generated mice with markedly elevated plasma levels of human low density lipoprotein (LDL) and reduced plasma levels of high density lipoprotein. These mice have no functional LDL receptors [LDLR−/−] and express a human apolipoprotein B-100 (apoB) transgene [Tg(apoB+/+)] with or without an apo(a) transgene [Tg(apoa+/−)]. Twenty animals (10 males and 10 females) of each of the following four genotypes were maintained on a chow diet: (i) LDLR−/−, (ii) LDLR−/−;Tg(apoa+/−), (iii) LDLR−/−;Tg(apoB+/+), and (iv)LDLR−/−;Tg(apoB+/+);Tg(apo+/−). The mice were killed at 6 mo, and the percent area of the aortic intimal surface that stained positive for neutral lipid was quantified. Mean percent areas of lipid staining were not significantly different between the LDLR−/− and LDLR−/−;Tg(apoa+/−) mice (1.0 ± 0.2% vs. 1.4 ± 0.3%). However, the LDLR−/−;Tg(apoB+/+) mice had ≈15-fold greater mean lesion area than the LDLR−/− mice. No significant difference was found in percent lesion area in the LDLR−/−;Tg(apoB+/+) mice whether or not they expressed apo(a) [18.5 ± 2.5%, without lipoprotein(a), Lp(a), vs. 16.0 ± 1.7%, with Lp(a)]. Histochemical analyses of the sections from the proximal aorta of LDLR−/−;Tg(apoB+/+) mice revealed large, complex, lipid-laden atherosclerotic lesions that stained intensely with human apoB-100 antibodies. In mice expressing Lp(a), large amounts of apo(a) protein colocalized with apoB-100 in the lesions. We conclude that LDLR−/−; Tg(apoB+/+) mice exhibit accelerated atherosclerosis on a chow diet and thus provide an excellent animal model in which to study atherosclerosis. We found no evidence that apo(a) increased atherosclerosis in this animal model.

Computerized polymorphic marker identification: Experimental validation and a predicted human polymorphism catalog

A computational system for the prediction of polymorphic loci directly and efficiently from human genomic sequence was developed and verified. A suite of programs, collectively called pompous (polymorphic marker prediction of ubiquitous simple sequences) detects tandem repeats ranging from dinucleotides up to 250 mers, scores them according to predicted level of polymorphism, and designs appropriate flanking primers for PCR amplification. This approach was validated on an approximately 750-kilobase region of human chromosome 3p21.3, involved in lung and breast carcinoma homozygous deletions. Target DNA from 36 paired B lymphoblastoid and lung cancer lines was amplified and allelotyped for 33 loci predicted by pompous to be variable in repeat size. We found that among those 36 predominately Caucasian individuals 22 of the 33 (67%) predicted loci were polymorphic with an average heterozygosity of 0.42. Allele loss in this region was found in 27/36 (75%) of the tumor lines using these markers. pompous provides the genetic researcher with an additional tool for the rapid and efficient identification of polymorphic markers, and through a World Wide Web site, investigators can use pompous to identify polymorphic markers for their research. A catalog of 13,261 potential polymorphic markers and associated primer sets has been created from the analysis of 141,779,504 base pairs of human genomic sequence in GenBank. This data is available on our Web site (pompous.swmed.edu) and will be updated periodically as GenBank is expanded and algorithm accuracy is improved.

Low-dose expression of a human apolipoprotein E transgene in macrophages restores cholesterol efflux capacity of apolipoprotein E-deficient mouse plasma

Apolipoprotein E- (apoE) deficient (E−/−) mice develop severe hyperlipidemia and diffuse atherosclerosis. Low-dose expression of a human apoE3 transgene in macrophages of apoE-deficient mice (E−/−hTgE+/0), which results in about 5% of wild-type apoE plasma levels, did not correct hyperlipidemia but significantly reduced the extent of atherosclerotic lesions. To investigate the contribution of apoE to reverse cholesterol transport, we compared plasmas of wild-type (E+/+), E−/−, and E−/−hTgE+/0 mice for the appearance of apoE-containing lipoproteins by electrophoresis and their capacity to take up and esterify 3H-labeled cholesterol from radiolabeled fibroblasts or J774 macrophages. Wild-type plasma displayed lipoproteins containing apoE that were the size of high density lipoprotein and that had either electrophoretic α or γ mobilities. Similar particles were also present in E−/−hTgE+/0 plasma. Depending on incubation time, E−/− plasma released 48–74% less 3H-labeled cholesterol from fibroblasts than E+/+ plasma, whereas cholesterol efflux into E−/−hTgE+/0 plasma was only 11–25% lower than into E+/+ plasma. E−/−hTgE+/0 plasma also released 10% more 3H-labeled cholesterol from radiolabeled J774 macrophages than E−/− plasma. E+/+ and E−/−hTgE+/0 plasma each esterified significantly more cell-derived 3H-labeled cholesterol than E−/− plasma. Moreover, E−/− plasma accumulated much smaller proportions of fibroblast-derived 3H-labeled cholesterol in fractions with electrophoretic γ and α mobility than E+/+ and E−/−hTgE+/0 plasma. Thus, low-dose expression of apoE in macrophages nearly restored the cholesterol efflux capacity of apoE-deficient plasma through the formation of apoE-containing particles, which efficiently take up cell-derived cholesterol, and through the increase of cholesterol esterification activity. Thus, macrophage-derived apoE may protect against atherosclerosis by increasing cholesterol efflux from arterial wall cells.

Epitopes close to the apolipoprotein B low density lipoprotein receptor-binding site are modified by advanced glycation end products

Advanced glycation end products (AGEs) are thought to contribute to the abnormal lipoprotein profiles and increased risk of cardiovascular disease of patients with diabetes and renal failure, in part by preventing apolipoprotein B (apoB)-mediated cellular uptake of low density lipoproteins (LDL) by LDL receptors (LDLr). It has been proposed that AGE modification at one site in apoB, almost 1,800 residues from the putative apoB LDLr-binding domain, may be sufficient to induce an apoB conformational change that prevents binding to the LDLr. To further explore this hypothesis, we used 29 anti-human apoB mAbs to identify other potential sites on apoB that may be modified by in vitro advanced glycation of LDL. Glycation of LDL caused a time-dependent decrease in its ability to bind to the LDLr and in the immunoreactivity of six distinct apoB epitopes, including two that flank the apoB LDLr-binding domain. ApoB appears to be modified at multiple sites by these criteria, as the loss of glycation-sensitive epitopes was detected on both native glycated LDL and denatured, delipidated glycated apoB. Moreover, residues directly within the putative apoB LDLr-binding site are not apparently modified in glycated LDL. We propose that the inability of LDL modified by AGEs to bind to the LDLr is caused by modification of residues adjacent to the putative LDLr-binding site that were undetected by previous immunochemical studies. AGE modification either eliminates the direct participation of the residues in LDLr binding or indirectly alters the conformation of the apoB LDLr-binding site.

Fibrinogen deficiency reduces vascular accumulation of apolipoprotein(a) and development of atherosclerosis in apolipoprotein(a) transgenic mice

To test directly whether fibrin(ogen) is a key binding site for apolipoprotein(a) [apo(a)] in vessel walls, apo(a) transgenic mice and fibrinogen knockout mice were crossed to generate fibrin(ogen)-deficient apo(a) transgenic mice and control mice. In the vessel wall of apo(a) transgenic mice, fibrin(ogen) deposition was found to be essentially colocalized with focal apo(a) deposition and fatty-streak type atherosclerotic lesions. Fibrinogen deficiency in apo(a) transgenic mice decreased the average accumulation of apo(a) in vessel walls by 78% and the average lesion (fatty streak type) development by 81%. Fibrinogen deficiency in wild-type mice did not significantly reduce lesion development. Our results suggest that fibrin(ogen) provides one of the major sites to which apo(a) binds to the vessel wall and participates in the generation of atherosclerosis.

Both apolipoprotein E and A-I genes are present in a nonmammalian vertebrate and are highly expressed during embryonic development

Apolipoprotein E (apoE) is associated with several classes of plasma lipoproteins and mediates uptake of lipoproteins through its ability to interact with specific cell surface receptors. Besides its role in cardiovascular diseases, accumulating evidence has suggested that apoE could play a role in neurodegenerative diseases, such as Alzheimer disease. In vertebrates, apoA-I is the major protein of high-density lipoprotein. ApoA-I may play an important role in regulating the cholesterol content of peripheral tissues through the reverse cholesterol transport pathway. We have isolated cDNA clones that code for apoE and apoA-I from a zebrafish embryo library. Analysis of the deduced amino acid sequences showed the presence of a region enriched in basic amino acids in zebrafish apoE similar to the lipoprotein receptor-binding region of human apoE. We demonstrated by whole-mount in situ hybridization that apoE and apoA-I genes are highly expressed in the yolk syncytial layer, an extraembryonic structure implicated in embryonic and larval nutrition. ApoE transcripts were also observed in the deep cell layer during blastula stage, in numerous ectodermal derivatives after gastrulation, and after 3 days of development in a limited number of cells both in brain and in the eyes. Our data indicate that apoE can be found in a nonmammalian vertebrate and that the duplication events, from which apoE and apoA-I genes arose, occurred before the divergence of the tetrapod and teleost ancestors. Zebrafish can be used as a simple and useful model for studying the role of apolipoproteins in embryonic and larval nutrition and of apoE in brain morphogenesis and regeneration.

A natural polymorphism in β-lactamase is a global suppressor

A M182T substitution was discovered as a second-site suppressor of a missense mutation in TEM-1 β-lactamase. The combination of the M182T substitution with other substitutions in the enzyme indicates the M182T substitution is a global suppressor of missense mutations in β-lactamase. The M182T substitution also is found in natural variants of TEM-1 β-lactamase with altered substrate specificity that have evolved in response to antibiotic therapy. The M182T substitution may have been selected in natural isolates as a suppressor of folding or stability defects resulting from mutations associated with drug resistance. This pathway of protein evolution may occur in other targets of antimicrobial drugs such as the HIV protease.