Genes to blood pressure: The effects of variation in genes of the renin-angiotensin system on the hormonal and renal control of blood pressure

Objective. Essential hypertension affects 25% of the US adult population and is a leading contributor to morbidity and mortality. Because BP is a multifactorial phenotype that resists simple genetic analysis, intermediate phenotypes within the complex network of BP regulatory systems may be more accessible to genetic dissection. The Renin-Angiotensin System (RAS) is known to influence intermediate and long-term blood pressure regulation through alterations in vascular tone and renal sodium and fluid resorption. This dissertation examines associations between renin (REN), angiotensinogen (AGT), angiotensin-converting enzyme (ACE) and angiotensin II type 1 receptor (AT1) gene variation and interindividual differences in plasma hormone levels, renal hemodynamics, and BP homeostasis.^ Methods. A total of 150 unrelated men and 150 unrelated women, between 20.0 and 49.9 years of age and free of acute or chronic illness except for a history of hypertension (11 men and 7 women, all off medications), were studied after one week on a controlled sodium diet. RAS plasma hormone levels, renal hemodynamics and BP were determined prior to and during angiotensin II (Ang II) infusion. Individuals were genotyped by PCR for a variable number tandem repeat (VNTR) polymorphism in REN, and for the following restriction fragment length polymorphisms (RFLP): AGT M235T, ACE I/D, and AT1 A1166C. Associations between clinical measurements and allelic variation were examined using multiple linear regression statistical models.^ Results. Women homozygous for the AT1 1166C allele demonstrated higher intracellular levels of sodium (p = 0.044). Men homozygous for the AGT T235 allele demonstrated a blunted decrement in renal plasma flow in response to Ang II infusion (p = 0.0002). There were no significant associations between RAS gene variation and interindividual variation in RAS plasma hormone levels or BP.^ Conclusions. Rather than identifying new BP controlling genes or alleles, the study paradigm employed in this thesis (i.e., measured genes, controlled environments and interventions) may provide mechanistic insight into how candidate genes affect BP homeostasis. ^

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

T-786C polymorphism of the endothelial nitric oxide synthase gene and neuralgia-inducing cavitational osteonecrosis of the jaws.

OBJECTIVE: We hypothesized that, similar to idiopathic hip osteonecrosis, the T-786C mutation of the endothelial nitric oxide synthase (eNOS) gene affecting nitric oxide (NO) production was associated with neuralgia-inducing cavitational osteonecrosis of the jaws (NICO). DESIGN: In 22 NICO patients, not having taken bisphosphonates, mutations affecting NO production (eNOS T-786C, stromelysin 5A6A) were measured by polymerase chain reaction. Two healthy normal control subjects were matched per case by race and gender. RESULTS: Homozygosity for the mutant eNOS allele (TT) was present in 6 out of 22 patients (27%) with NICO compared with 0 out of 44 (0%) race and gender-matched control subjects; heterozygosity (TC) was present in 8 patients (36%) versus 15 control subjects (34%); and the wild-type normal genotype (CC) was present in 9 patients (36%) versus 29 controls (66%) (P = .0008). The mutant eNOS T-786C allele was more common in cases (20 out of 44 [45%]) than in control subjects (15 out of 88 [17%]) (P = .0005). The distribution of the stromelysin 5A6A genotype in cases did not differ from control subjects (P = .13). CONCLUSIONS: The eNOS T-786C polymorphism affecting NO production is associated with NICO, may contribute to the pathogenesis of NICO, and may open therapeutic medical approaches to treatment of NICO through provision of L-arginine, the amino-acid precursor of NO.

Origins and evolution of VNTR loci: The apolipoprotein B 3$\sp\prime$ VNTR

I studied the apolipoprotein (apo) B 3$\sp\prime$ variable number tandem repeat (VNTR) and did computer simulations of the stepwise mutation model to address four questions: (1) How did the apo B VNTR originate? (2) What is the mutational mechanism of repeat number change at the apo B VNTR? (3) To what extent are population and molecular level events responsible for the determination of the contemporary apo B allele frequency distribution? (4) Can VNTR allele frequency distributions be explained by a simple and conservative mutation-drift model? I used three general approaches to address these questions: (1) I characterized the apo B VNTR region in non-human primate species; (2) I constructed haplotypes of polymorphic markers flanking the apo B VNTR in a sample of individuals from Lorrain, France and studied the associations between the flanking-marker haplotypes and apo B VNTR size; (3) I did computer simulations of the one-step stepwise mutation model and compared the results to real data in terms of four allele frequency distribution characteristics.^ The results of this work have allowed me to conclude that the apo B VNTR originated after an initial duplication of a sequence which is still present as a single copy sequence in New World monkey species. I conclude that this locus did not originate by the transposition of an array of repeats from somewhere else in the genome. It is unlikely that recombination is the primary mutational mechanism. Furthermore, the clustered nature of these associations implicates a stepwise mutational mechanism. From the high frequencies of certain haplotype-allele size combinations, it is evident that population level events have also been important in the determination of the apo B VNTR allele frequency distribution. Results from computer simulations of the one-step stepwise mutation model have allowed me to conclude that bimodal and multimodal allele frequency distributions are not unexpected at loci evolving via stepwise mutation mechanisms. Short tandem repeat loci fit the stepwise mutation model best, followed by microsatellite loci. I therefore conclude that there are differences in the mutational mechanisms of VNTR loci as classed by repeat unit size. (Abstract shortened by UMI.) ^

Population genetics of VNTR loci and their applications in evolutionary studies

Variable number of tandem repeats (VNTR) are genetic loci at which short sequence motifs are found repeated different numbers of times among chromosomes. To explore the potential utility of VNTR loci in evolutionary studies, I have conducted a series of studies to address the following questions: (1) What are the population genetic properties of these loci? (2) What are the mutational mechanisms of repeat number change at these loci? (3) Can DNA profiles be used to measure the relatedness between a pair of individuals? (4) Can DNA fingerprint be used to measure the relatedness between populations in evolutionary studies? (5) Can microsatellite and short tandem repeat (STR) loci which mutate stepwisely be used in evolutionary analyses?^ A large number of VNTR loci typed in many populations were studied by means of statistical methods developed recently. The results of this work indicate that there is no significant departure from Hardy-Weinberg expectation (HWE) at VNTR loci in most of the human populations examined, and the departure from HWE in some VNTR loci are not solely caused by the presence of population sub-structure.^ A statistical procedure is developed to investigate the mutational mechanisms of VNTR loci by studying the allele frequency distributions of these loci. Comparisons of frequency distribution data on several hundreds VNTR loci with the predictions of two mutation models demonstrated that there are differences among VNTR loci grouped by repeat unit sizes.^ By extending the ITO method, I derived the distribution of the number of shared bands between individuals with any kinship relationship. A maximum likelihood estimation procedure is proposed to estimate the relatedness between individuals from the observed number of shared bands between them.^ It was believed that classical measures of genetic distance are not applicable to analysis of DNA fingerprints which reveal many minisatellite loci simultaneously in the genome, because the information regarding underlying alleles and loci is not available. I proposed a new measure of genetic distance based on band sharing between individuals that is applicable to DNA fingerprint data.^ To address the concern that microsatellite and STR loci may not be useful for evolutionary studies because of the convergent nature of their mutation mechanisms, by a theoretical study as well as by computer simulation, I conclude that the possible bias caused by the convergent mutations can be corrected, and a novel measure of genetic distance that makes the correction is suggested. In summary, I conclude that hypervariable VNTR loci are useful in evolutionary studies of closely related populations or species, especially in the study of human evolution and the history of geographic dispersal of Homo sapiens. (Abstract shortened by UMI.) ^

The angiotensin-converting enzyme insertion/deletion polymorphism and its associations with carotid artery wall thickness and coronary heart disease: The atherosclerosis risk in communities study

Coronary heart disease (CHD) is the leading cause of death in the United States. Recently, renin-angiotensin system (RAS) was found associated with atherosclerosis formation, with angiotensin II inducing vascular smooth muscle cell growth and migration, platelet activation and aggregation, and stimulation of plasminogen activator inhibitor-1. Angiotensin II is converted from angiotensin I by angiotensin I-converting enzyme (ACE) and this enzyme is mainly genetically determined. The ACE gene has been assigned to chromosome 17q23 and an insertion/deletion (I/D)polymorphism has been characterized by the presence/absence of a 287 bp fragment in intron 16 of the gene. The two alleles form three genotypes, namely, DD, ID and II and the DD genotype has been linked to higher plasma ACE levels and cell ACE activity.^ In this study, the association between the ACE I/D polymorphism and carotid artery wall thickness measured by B-mode ultrasound was investigated in a biracial sample, and the association between the gene and incident CHD was investigated in whites and if the gene-CHD association in whites, if any, was due to the gene effect on atherosclerosis. The study participants are from the prospective Atherosclerosis Risk in Communities (ARIC) Study, including adults aged 45 to 65 years. The present dissertation used a matched case-control design for studying the associations of the ACE gene with carotid artery atherosclerosis and an unmatched case-control design for the association of the gene with CHD. A significant recessive effect of the D allele on carotid artery thickness was found in blacks (OR = 3.06, 95% C.I: 1.11-8.47, DD vs. ID and II) adjusting for age, gender, cigarette smoking, LDL-cholesterol and diabetes. No similar associations were found in whites. The ACE I/D polymorphism is significantly associated with coronary heart disease in whites, and while stratifying data by carotid artery wall thickness, the significant associations were only observed in thin-walled subgroups. Assuming a recessive effect of the D allele, odds ratio was 2.84 (95% C.I:1.17-6.90, DD vs. ID and II) and it was 2.30 (95% C.I:1.22-4.35, DD vs. ID vs. II) assuming a codominant effect of the D allele. No significant associations were observed while comparing thick-walled CHD cases with thin-walled controls. Following conclusions could be drawn: (1) The ACE I/D polymorphism is unlikely to confer appreciable increase in the risk of carotid atherosclerosis in US whites, but may increases the risk of carotid atherosclerosis in blacks. (2) ACE I/D polymorphism is a genetic risk factor for incident CHD in US whites and this effect is separate from the chronic process of atherosclerosis development. Finally, the associations observed here are not causal, since the I/D polymorphism is in an intron, where no ACE proteins are encoded. ^

Global patterns of DNA polymorphism in noncoding regions in human populations

DNA sequence variation is currently a major source of data for studying human origins, evolution, and demographic history, and for detecting linkage association of complex diseases. In this dissertation, I investigated DNA variation in worldwide populations from two ∼10 kb autosomal regions on 22q11.2 (noncoding) and 1q24 (introns). A total of 75 variant sites were found among 128 human sequences in the 22q11.2 region, yielding an estimate of 0.088% for nucleotide diversity (π), and a total of 52 variant sites were found among 122 human sequences in the 1q24 region with an estimated π value of 0.057%. The data from these two regions and a 10 kb noncoding region on Xq13.3 all show a strong excess of low-frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The effective population sizes estimated from the three regions were 11,000, 12,700, and 8,600, respectively, which are close to the commonly used value of 10,000. In each of the two autosomal regions, the age of the most recent common ancestor (MRCA) was estimated to be older than 1 million years among all the sequences and ∼600,000 years among non-African sequences, providing first evidence from autosomal noncoding or intronic regions for a genetic history of humans much more ancient than the emergence of modern humans. The ancient genetic history of humans indicates no severe bottleneck during the evolution of humans in the last half million years; otherwise, much of the ancient genetic history would have been lost during a severe bottleneck. This study strongly suggests that both the “out of Africa” and the multiregional models are too simple for explaining the evolution of modern humans. A compilation of genome-wide data revealed that nucleotide diversity is highest in autosomal regions, intermediate in X-linked regions, and lowest in Y-linked regions. The data suggest the existence of background selection or selective sweep on Y-linked loci. In general, the nucleotide diversity in humans is low compared to that in chimpanzee and Drosophila populations. ^

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. ^

Are functional and genetic components of platelets linked to coronary artery disease: A case-control study

Coronary artery disease (CAD) is a multifactorial disease process involving behavioral, inflammatory, clinical, thrombotic, and genetic components. Previous epidemiologic studies focused on identifying behavioral and demographic risk factors of CAD, but none focused on platelets. Current platelet literature lacks the known effects of platelet function and platelet receptor polymorphisms on CAD. This case-control analysis addressed these issues by analyzing data collected for a previous study. Cases were individuals who had undergone CABG and thus had been diagnosed with CAD, while the controls were volunteers presumed to be CAD free. The platelet function variables analyzed included fibrinogen Von Willebrand Factor activity (VWF), shear-induced platelet aggregation (SIPA), sCD40L, and mean platelet volume; and the platelet polymorphisms studied included PIA, α2 807, Ko, Kozak, and VNTR. Univariate analysis found fibrinogen, VWF, SIPA, and PIA to be independent risk factors of CAD. Logistic regression was used to build a predictive model for CAD using the platelet function and platelet polymorphism data adjusted for age, sex, race, and current smoking status. A model containing only platelet polymorphisms and their respective receptor densities, found polymorphisms within GPIbα to be associated with CAD, yielding an 86% (95% C.I. 0.97–3.55) increased risk with the presence of at least 1 polymorphism in Ko, Kozak, or VNTR. Another model included both platelet function and platelet polymorphism data. Fibrinogen, the receptor density of GPIbα, and the polymorphism in GPIa-IIa (α2 807) were all associated with CAD with odds ratios of 1.10, 1.04, and 2.30 for fibrinogen (10mg/dl increase), GPIbα receptors (1 MFI increase), and GPIa-IIa, respectively. In addition, risk estimates and 99% confidence intervals adjusted for race were calculated to determine if the presence of a platelet receptor polymorphism was associated with CAD. The results were as follows: PIA (1.64, 0.74–3.65); α2 807 (1.35, 0.77–2.37); Ko (1.71, 0.70–4.16); Kozak (1.17, 0.54–2.52); and VNTR (1.24, 0.52–2.91). Although not statistically significant, all platelet polymorphisms were associated with an increased risk for CAD. These exploratory findings indicate that platelets do appear to have a role in atherosclerosis and that anti-platelet drugs targeting GPI-IIa and GPIbα may be better treatment candidates for individuals with CAD. ^

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. ^

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. ^