Association of putative enteroaggregative Escherichia coli virulence genes and biofilm production in isolates from travelers to developing countries.

Enteroaggregative Escherichia coli (EAEC) is an emerging enteric pathogen that causes acute and chronic diarrhea among children, human immunodeficiency virus-infected patients, and travelers to developing regions of the world. The pathogenesis of EAEC strains involves the production of biofilm. In this study, we determined the association between presence of putative EAEC virulence genes and biofilm formation in 57 EAEC isolates (as defined by HEp-2 adherence) from travelers with diarrhea and in 18 EAEC isolates from travelers without diarrhea. Twelve nondiarrheagenic E. coli isolates from healthy travelers were used as controls. Biofilm formation was measured by using a microtiter plate assay with the crystal violet staining method, and the presence of the putative EAEC virulence genes aap, aatA, aggR, astA, irp2, pet, set1A, and shf was determined by PCR. EAEC isolates were more likely to produce biofilm than nondiarrheagenic E. coli isolates (P = 0.027), and the production of biofilm was associated with the virulence genes aggR, set1A, aatA, and irp2, which were found in 16 (40%), 17 (43%), 10 (25%), and 27 (68%) of the biofilm producers versus only 4 (11%), 6 (6%), 2 (6%), and 15 (43%) in non-biofilm producers (P = 0.008 for aggR, P = 0.0004 for set1A, P = 0.029 for aatA, and P = 0.04 for irp2). Although the proportion of EAEC isolates producing biofilm in patients with diarrhea (51%) was similar to that in patients without diarrhea (61%), biofilm production was related to the carriage of aggR (P = 0.015), set1A (P = 0.001), and aatA (P = 0.025). Since aggR is a master regulator of EAEC, the presence of aap (P = 0.004), astA (P = 0.001), irp2 (P = 0.0006), pet (P = 0.002), and set1A (P = 0.014) in an aggR versus an aggR-lacking background was investigated and was also found to be associated with biofilm production. This study suggests that biofilm formation is a common phenomenon among EAEC isolates derived from travelers with or without diarrhea and that multiple genes associated with biofilm formation are regulated by aggR.

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Single nucleotide polymorphisms (SNPs) associated with TGF-beta pathway and their significance in systemic sclerosis - A multilevel analysis

Systemic sclerosis (SSc) or Scleroderma is a complex disease and its etiopathogenesis remains unelucidated. Fibrosis in multiple organs is a key feature of SSc and studies have shown that transforming growth factor-β (TGF-β) pathway has a crucial role in fibrotic responses. For a complex disease such as SSc, expression quantitative trait loci (eQTL) analysis is a powerful tool for identifying genetic variations that affect expression of genes involved in this disease. In this study, a multilevel model is described to perform a multivariate eQTL for identifying genetic variation (SNPs) specifically associated with the expression of three members of TGF-β pathway, CTGF, SPARC and COL3A1. The uniqueness of this model is that all three genes were included in one model, rather than one gene being examined at a time. A protein might contribute to multiple pathways and this approach allows the identification of important genetic variations linked to multiple genes belonging to the same pathway. In this study, 29 SNPs were identified and 16 of them located in known genes. Exploring the roles of these genes in TGF-β regulation will help elucidate the etiology of SSc, which will in turn help to better manage this complex disease. ^

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

MULTIPLE POSTTRANSCRIPTIONAL REGULATORY FEATURES CONTROL EXPRESSION OF ETHANOLAMINE UTILIZATION GENES IN ENTEROCOCCUS FAECALIS

Enterococcus faecalis is a Gram-positive bacterium that lives as a commensal organism in the mammalian gastrointestinal tract, but can behave as an opportunistic pathogen. Our lab discovered that mutation of the eutK gene attenuates virulence of E. faecalis in the C. elegans model host. eutK is part of the ethanolamine metabolic pathway which was previously unknown in E. faecalis. I discovered the presence of two unique posttranscriptional regulatory features that control expression of eut locus genes. The first feature I found is an AdoCBL riboswitch, a cis-acting RNA regulatory element that acts as a positive regulator of gene expression. The second feature I discovered is a unique two-component system, EutVW. The EutV response regulator contains an ANTAR family domain, which binds RNA to trigger transcriptional antitermination. I determined that induction of expression of several genes in the eut locus is dependent on ethanolamine, AdoCBL and the two-component system. AdoCBL and ethanolamine are both required for induction of eut locus gene expression. Additionally, I discovered eutG is regulated by a unique mechanism of antitermination. Both the AdoCBL riboswitch and EutV response regulator control the expression of the downstream gene eutG. EutV potentially acts through a novel antitermination mechanism in which a dimer of EutV binds to a pair of mRNA stem loops forming an antitermination complex. My data show a unique mechanism by which two environmental signals are integrated by two different posttranscriptional regulators to regulate a single locus.

Statistical Methods for Differential Expressions of Genes Detected in Multiple-Condition Experiment of Microarray

Most studies of differential gene-expressions have been conducted between two given conditions. The two-condition experimental (TCE) approach is simple in that all genes detected display a common differential expression pattern responsive to a common two-condition difference. Therefore, the genes that are differentially expressed under the other conditions other than the given two conditions are undetectable with the TCE approach. In order to address the problem, we propose a new approach called multiple-condition experiment (MCE) without replication and develop corresponding statistical methods including inference of pairs of conditions for genes, new t-statistics, and a generalized multiple-testing method for any multiple-testing procedure via a control parameter C. We applied these statistical methods to analyze our real MCE data from breast cancer cell lines and found that 85 percent of gene-expression variations were caused by genotypic effects and genotype-ANAX1 overexpression interactions, which agrees well with our expected results. We also applied our methods to the adenoma dataset of Notterman et al. and identified 93 differentially expressed genes that could not be found in TCE. The MCE approach is a conceptual breakthrough in many aspects: (a) many conditions of interests can be conducted simultaneously; (b) study of association between differential expressions of genes and conditions becomes easy; (c) it can provide more precise information for molecular classification and diagnosis of tumors; (d) it can save lot of experimental resources and time for investigators.^

ΔNp63 REGULATES A COMPLEX NETWORK OF TARGET GENES IN LIMB AND EPIDERMAL DEVELOPMENT

The skin is composed of two major compartments, the dermis and epidermis. The epidermis forms a barrier to protect the body. The stratified epithelium has self-renewing capacity throughout life, and continuous turnover is mediated by stem cells in the basal layer. p63 is structurally and functionally related to p53. In spite of their structural similarities, p63 is critical for the development and maintenance of stratified epithelial tissues, unlike p53. p63 is highly expressed in the epidermis and previously has been shown to play a critical role in the development and maintenance of the epidermis. The study of p63 has been complicated due to the existence of multiple isoforms: those with a transactivation domain (TAp63) and those lacking this domain (ΔNp63). Mice lacking p63 cannot form skin, have craniofacial and skeletal defects and die within hours after birth. These defects are due to the ability of p63 to regulate multiple processes in skin development including epithelial stem cell proliferation, differentiation, and adherence programs. To determine the roles of these isoforms in skin development and maintenance, isoform specific p63 conditional knock out mice were generated by our lab. TAp63-/- mice age prematurely, develop blisters, and display wound-healing defects that result from hyperproliferation of dermal stem cells. That results in premature depletion of these cells, which are necessary for wound repair, that indicates TAp63 plays a role in dermal/epidermal maintenance. To study the role of ΔNp63, I generated a ΔNp63-/- mouse and analyzed the skin by performing immunofluorescence for markers of epithelial differentiation. The ΔNp63-/- mice developed a thin, disorganized epithelium but differentiation markers were expressed. Interestingly, the epidermis from ΔNp63-/- mice co-expressed K14 and K10 in the same cell suggesting defects in epidermal differentiation and stratification. This phenotype is reminiscent of the DGCR8fl/fl;K14Cre and Dicerfl/fl;K14Cre mice skin. Importantly, DGCR8-/- embryonic stem cells (ESCs) display a hyperproliferation defect by failure to silence pluripotency genes. Furthermore, I have observed that epidermal cells lacking ΔNp63 display a phenotype reminiscent of embryonic stem cells instead of keratinocytes. Thus, I hypothesize that genes involved in maintaining pluripotency, like Oct4, may be upregulated in the absence of ΔNp63. To test this, q-RT PCR was performed for Oct4 mRNA with wild type and ΔNp63-/- 18.5dpc embryo skin. I found that the level of Oct4 was dramatically increased in the absence of ΔNp63-/-. Based on these results, I hypothesized that ΔNp63 induces differentiation by silencing pluripotency regulators, Oct4, Sox2 and Nanog directly through the regulation of DGCR8. I found that DGCR8 restoration resulted in repression of Oct4, Sox2 and Nanog in ΔNp63-/- epidermal cells and rescue differentiation defects. Loss of ΔNp63 resulted in pluripotency that caused defect in proper differentiation and stem cell like phenotype. This led me to culture the ΔNp63-/- epidermal cells in neuronal cell culture media in order to address whether restoration of DGCR8 can transform epidermal cells to neuronal cells. I found that DGCR8 restoration resulted in a change in cell fate. I also found that miR470 and miR145 play a role in the induction of pluripotency by repressing Oct4, Sox2 and Nanog. This indicates that ΔNp63 induces terminal differentiation through the regulation of DGCR8.

Cotransfer of antibiotic resistance genes and a hylEfm-containing virulence plasmid in Enterococcus faecium.

The hyl(Efm) gene (encoding a putative hyaluronidase) has been found almost exclusively in Enterococcus faecium clinical isolates, and recently, it was shown to be on a plasmid which increased the ability of E. faecium strains to colonize the gastrointestinal tract. In this work, the results of mating experiments between hyl(Efm)-containing strains of E. faecium belonging to clonal cluster 17 and isolated in the United States and Colombia indicated that the hyl(Efm) gene of these strains is also carried on large plasmids (>145 kb) which we showed transfer readily from clinical strains to E. faecium hosts. Cotransfer of resistance to vancomycin and high-level resistance (HLR) to aminoglycosides (gentamicin and streptomycin) and erythromycin was also observed. The vanA gene cluster and gentamicin resistance determinants were genetically linked to hyl(Efm), whereas erm(B) and ant(6)-I, conferring macrolide-lincosamide-streptogramin B resistance and HLR to streptomycin, respectively, were not. A hyl(Efm)-positive transconjugant resulting from a mating between a well-characterized endocarditis strain [TX0016 (DO)] and a derivative of a fecal strain of E. faecium from a healthy human volunteer (TX1330RF) exhibited increased virulence in a mouse peritonitis model. These results indicate that E. faecium strains use a strategy which involves the recruitment into the same genetic unit of antibiotic resistance genes and determinants that increase the ability to produce disease. Our findings indicate that the acquisition of the hyl(Efm) plasmids may explain, at least in part, the recent successful emergence of some E. faecium strains as nosocomial pathogens.

Regulation of Spec genes in aboral ectoderm cells of sea urchin {\it Strongylocentrotus purpuratus}

The Spec genes of the sea urchin Stronylocentrotus purpuratus serves as an excellent model for studying cell type-specific gene expression during early embryogenesis. The Spec1/Spec2 genes encode cytosolic calcium-binding proteins related to the calmodulin/troponin C/myosin light chain superfamily. Members of the Spec gene family are activated shortly after the sixth cleavage as the lineage-specific founder cells giving rise to aboral ectoderm are established, and the accumulation of the Spec mRNAs is limited exclusively to aboral ectoderm cell lineages. In this dissertation, the transcriptional regulation of the Spec genes was studied. Sequence comparisons of the Spec gene 5$\sp\prime$ flanking regions showed that a DNA block of approximately 800 bp from the 3$\sp\prime$ end of the first exon to the 5$\sp\prime$ end of a repetitive DNA element, termed RSR, was highly conserved. In Spec2a, the conserved region was a continuous stretch of DNA, but in Spec1 and Spec2c, DNA insertions interrupt the conserved sequence block and alter the relative placement of the RSR element and other 5$\sp\prime$ flanking DNA. Thus, drastic rearrangements have occurred within the putative control regions of the Spec genes. In vivo expression experiments using the sea urchin embryo gene-transfer system showed that while the 5$\sp\prime$ flanking regions of all three Spec genes conferred proper temporal activation to the reporter CAT gene, only the Spec2a 5$\sp\prime$ flanking region could restrict lacZ gene expression to aboral ectoderm cells. However, the Spec2a conserved region alone was not sufficient to confer proper spatial expression, suggesting that negative spatial elements are also associated with the proper activation of Spec2a. A major positive regulatory region, defined as the RSR enhancer, was identified between base pairs $-$631 and $-$443 on Spec2a. The RSR enhancer was essential for maximal activity and conferred preferential aboral ectoderm expression to a lacZ reporter gene. DNaseI footprinting and band-shift analysis of the RSR enhancer revealed multiple DNA-elements. One of the elements, an A/T-rich sequence called the A/T palindrome was studied in detail. This element binds a single 45-kDa nuclear protein, the A/T palindrome binding protein (A/TBP), whose DNA-binding specificity suggests a possible relationship with the bicoid-class homeodomain proteins. Mutated A/T palindromes are incapable of binding the 45-kDa protein and lower promoter activity by 8-fold. DNA-binding activity for A/TBP is low in unfertilized eggs, increases by the 16-cell stage and continues rising in blastulae. These data suggest that A/TBP plays a major role in the activation of the Spec2a gene in aboral ectoderm cells. ^

{\it Hox\/} genes and specification of regional identity in the axial skeleton: Analysis of the individual and combined mutant phenotypes of the paralogous group 4 murine {\it Hox\/} genes; {\it hoxa\/}-4, {\it hoxb\/}-4 and {\it hoxd\/}-4

The Hox gene products are transcription factors involved in specifying regional identity along the anteroposterior body axis. In Drosophila, where these genes are known as HOM-C (Homeotic-complex) genes and where they have been most extensively studied, they are expressed in restricted domains along the anteroposterior axis with different anterior limits. Genetic analysis of a large number of gain- and loss-of-function alleles of these genes has revealed that these genes are important in specifying segmental identity at their anterior limits of expression. Furthermore, there is a functional dominance of posterior genes over anterior genes, such that posterior genes can dominantly specify their developmental programs in spite of the expression of more anterior genes in the same segment. In the mouse, there are four clusters of HOM-C genes, called Hox genes. Thus, there may be up to four genes, called paralogs, that are more highly homologous to each other and to their Drosophila homolog than they are to the other mouse Hox genes. The single mutants for two paralogous genes, hoxa-4 and hoxd-4, presented in this dissertation, are similar to several other mouse Hox mutants in that they show partial, incompletely penetrant homeotic transformations of vertebrae at their anterior limit of expression. These mutants were then bred with hoxb-4 mutants (Ramirez-Solis, et al. 1993) to generate the three possible double mutant combinations as well as the triple mutant. The skeletal phenotypes of these group 4 Hox compound mutants displayed clear alterations in regional identity, such that a nearly complete transformation towards the morphology of the first cervical vertebra occurs. These results suggest a certain degree of functional redundancy among paralogous genes in specifying regional identity. Furthermore, there was a remarkable dose-dependent increase in the number of vertebrae transformed to a first cervical vertebra identity, including the second through the fifth cervical vertebrae in the triple mutant. Thus, these genes are required in a larger anteroposterior domain than is revealed by the single mutant phenotypes alone, such that multiple mutations in these genes result in transformations of vertebrae that are not at their anterior limit of expression. ^

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

Evaluation of candidate genes for susceptibility to spina bifida

Neural tube defects (NTDs) are the most common severely disabling birth defects in the United States, with a frequency of approximately 1–2 of every 1,000 births. This text includes the identification and evaluation of candidate susceptibility genes that confer risk for the development of neural tube defects (NTDs). The project focused on isolated meningomyelocele, also termed spina bifida (SB). ^ Spina bifida is a complex disease with multifactorial inheritance, therefore the subject population (consisting of North American Caucasians and Hispanics of Mexicali-American descent) was composed of 459 simplex SB families who were tested for genetic associations utilizing the transmission disequilibrium test (TDT), a nonparametric linkage technique. Three categories of candidate genes were studied, including (1) human equivalents of genes determined in mouse models to cause NTDs, (2) HOX and PAX genes, and (3) the MTHFR gene involved in the metabolic pathway of folate. ^ The C677T variant of the 5,10-methylenetetrahydrofolate reductase (MTHFR) gene was the first mutation in this gene to be implicated as a risk factor for NTDs. Our evaluation of the MTHFR gene provides evidence that maternal C677T homozygosity is a risk factor for upper level spina bifida defects in Hispanics [OR = 2.3, P = 0.02]. This observed risk factor is of great importance due to the high prevalence of this homozygous genotype in the Hispanic population. Additionally, maternal C677T/A1298C compound heterozygosity is a risk factor for upper level spina bifida defects in non-Hispanic whites [OR = 3.6, P = 0.03]. ^ For TDT analysis, our total population of 1128 subjects were genotyped for 54 markers from within and/or flanking the 20 candidate genes/gene regions of interest. Significant TDT findings were obtained for 3 of the 54 analyzed markers: d20s101 flanking the PAX1 gene (P = 0.019), d1s228 within the PAX7 gene (P = 0.011), and d2s110 within the PAX8 gene (P = 0.013). These results were followed-up by testing the genes directly for mutations utilizing single-strand conformational analysis (SSCA) and direct sequencing. Multiple variations were detected in each of these PAX genes; however, these variations were not passed from parent to child in phase with the positively transmitted alleles. Therefore, these variations do not contribute to the susceptibility of spina bifida, but rather are previously unreported single nucleotide polymorphisms. ^

Familial aggregation, candidate genes and genome scans: Analyzing the role of genetics in stroke

The genetic etiology of stroke likely reflects the influence of multiple loci with small effects, each modulating different pathophysiological processes. This research project utilized three analytical strategies to address the paucity of information related to the identification and characterization of genetic variation associated with stroke in the general population. ^ First, the general contribution of familial factors to stroke susceptibility was evaluated in a population-based sample of unrelated individuals. Increased risk of subclinical cerebral infarction was observed among individuals with a positive parental history of stroke. This association did not appear to be mediated by established stroke risk factors, specifically blood pressure levels or hypertension status. ^ The need to identify specific gene variation associated with stroke in the general population was addressed by evaluating seven candidate gene polymorphisms in a population-based sample of unrelated individuals. Three polymorphisms were significantly associated with increased subclinical cerebral infarction or incident clinical ischemic stroke risk. These relationships include the G-protein β3 subunit 825C/T polymorphism and clinical stroke in Whites, the lipoprotein lipase S/X447 polymorphism and subclinical and clinical stroke in men, and the angiotensin I-converting enzyme Ins/Del polymorphism and subclinical stroke in White men. These associations did not appear to be obfuscated by the stroke risk factors adjusted for in the analysis models specifically blood pressure levels or anti-hypertensive medication use. ^ The final research strategy considered, on a genome-wide scale, the idea that genetic variation may contribute to the occurrence of hypertension or stroke through a common etiologic pathway. Genomic regions were identified for which significant evidence of heterogeneity was observed among hypertensive sibpairs stratified by family history of stroke information. Regions identified on chromosome 15 in African Americans, and chromosome 13 in Whites and African Americans, suggest the presence of genes influencing hypertension and stroke susceptibility. ^ Insight into the role of genetics in stroke is useful for the potential early identification of individuals at increased risk for stroke and improved understanding of the etiology of the disease. The ultimate goal of these endeavors is to guide the development of therapeutic intervention and informed prevention to provide a lasting and positive impact on public health. ^

Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Diabetes genes and risk of prostate cancer in the Atherosclerosis Risk in Communities study

Prostate cancer (PrCa) is a leading cause of morbidity and mortality, yet the etiology remains uncertain. Meta-analyses show that PrCa risk is reduced by 16% in men with type 2 diabetes (T2D), but the mechanism is unknown. Recent genome-wide association studies and meta-analyses have found single nucleotide polymorphisms (SNPs) that consistently predict T2D risk. We evaluated associations of incident PrCa with 14 T2D SNPs in the Atherosclerosis Risk in Communities (ARIC) study. From 1987-2000, there were 397 incident PrCa cases ascertained from state or local cancer registries among 6,642 men (1,560 blacks and 5,082 whites) aged 45-64 years at baseline. Genotypes were determined by TaqMan assay. Cox proportional hazards models were used to assess the association between PrCa and increasing number of T2D risk-raising alleles for individual SNPs and for genetic risk scores (GRS) comprised of the number of T2D risk-raising alleles across SNPs. Two-way gene-gene interactions were evaluated with likelihood ratio tests. Using additive genetic models, the T2D risk-raising allele was associated with significantly reduced risk of PrCa for IGF2BP2 rs4402960 (hazard ratio [HR]=0.79; P=0.07 among blacks only), SLC2A2 rs5400 (race-adjusted HR=0.85; P=0.05) and UCP2 rs660339 (race-adjusted HR=0.84; P=0.02), but significantly increased risk of PrCa for CAPN10 rs3792267 (race-adjusted HR=1.20; P=0.05). No other SNPs were associated with PrCa using an additive genetic model. However, at least one copy of the T2D risk-raising allele for TCF7L2 rs7903146 was associated with reduced PrCa risk using a dominant genetic model (race-adjusted HR=0.79; P=0.03). These results imply that the T2D-PrCa association may be partly due to shared genetic variation, but these results should be verified since multiple tests were performed. When the combined, additive effects of these SNPs were tested using a GRS, there was nearly a 10% reduction in risk of PrCa per T2D risk-raising allele (race-adjusted HR=0.92; P=0.02). SNPs in IGF2BP2, KCNJ11 and SLC2A2 were also involved in multiple synergistic gene-gene interactions on a multiplicative scale. In conclusion, it appears that the T2D-PrCa association may be due, in part, to common genetic variation. Further knowledge of T2D gene-PrCa mechanisms may improve understanding of PrCa etiology and may inform PrCa prevention and treatment.^