Erin N. Smith

Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors

Natural genetic variation can cause significant differences in gene expression, but little is known about the polymorphisms that affect gene regulation. We analyzed regulatory variation in a cross between laboratory and wild strains of Saccharomyces cerevisiae. Clustering and linkage analysis defined groups of coregulated genes and the loci involved in their regulation. Most expression differences mapped to trans-acting loci. Positional cloning and functional assays showed that polymorphisms in GPA1 and AMN1 affect expression of genes involved in pheromone response and daughter cell separation, respectively. We also asked whether particular classes of genes were more likely to contain trans-regulatory polymorphisms. Notably, transcription factors showed no enrichment, and trans-regulatory variation seems to be broadly dispersed across classes of genes with different molecular functions.

Integrating Large-Scale Functional Genomic Data to Dissect the Complexity of Yeast Regulatory Networks

A key goal of biology is to construct networks that predict complex system behavior. We combine multiple types of molecular data, including genotypic, expression, transcription factor binding site (TFBS), and protein–protein interaction (PPI) data previously generated from a number of yeast experiments, in order to reconstruct causal gene networks. Networks based on different types of data are compared using metrics devised to assess the predictive power of a network. We show that a network reconstructed by integrating genotypic, TFBS and PPI data is the most predictive. This network is used to predict causal regulators responsible for hot spots of gene expression activity in a segregating yeast population. We also show that the network can elucidate the mechanisms by which causal regulators give rise to larger-scale changes in gene expression activity. We then prospectively validate predictions, providing direct experimental evidence that predictive networks can be constructed by integrating multiple, appropriate data types.

Gene–Environment Interaction in Yeast Gene Expression

The effects of genetic variants on phenotypic traits often depend on environmental and physiological conditions, but such gene–environment interactions are poorly understood. Recently developed approaches that treat transcript abundances of thousands of genes as quantitative traits offer the opportunity to broadly characterize the architecture of gene–environment interactions. We examined the genetic and molecular basis of variation in gene expression between two yeast strains (BY and RM) grown in two different conditions (glucose and ethanol as carbon sources). We observed that most transcripts vary by strain and condition, with 2,996, 3,448, and 2,037 transcripts showing significant strain, condition, and strain–condition interaction effects, respectively. We expression profiled over 100 segregants derived from a cross between BY and RM in both growth conditions, and identified 1,555 linkages for 1,382 transcripts that show significant gene–environment interaction. At the locus level, local linkages, which usually correspond to polymorphisms in cis-regulatory elements, tend to be more stable across conditions, such that they are more likely to show the same effect or the same direction of effect across conditions. Distant linkages, which usually correspond to polymorphisms influencing trans-acting factors, are more condition-dependent, and often show effects in different directions in the two conditions. We characterized a locus that influences expression of many growth-related transcripts, and showed that the majority of the variation is explained by polymorphism in the gene IRA2. The RM allele of IRA2 appears to inhibit Ras/PKA signaling more strongly than the BY allele, and has undergone a change in selective pressure. Our results provide a broad overview of the genetic architecture of gene–environment interactions, as well as a detailed molecular example, and lead to key insights into how the effects of different classes of regulatory variants are modulated by the environment. These observations will guide the design of studies aimed at understanding the genetic basis of complex traits.

Genome-wide association study of bipolar disorder in European American and African American individuals

To identify bipolar disorder (BD) genetic susceptibility factors, we conducted two genome-wide association (GWA) studies: one involving a sample of individuals of European ancestry (EA; n=1001 cases; n=1033 controls), and one involving a sample of individuals of African ancestry (AA; n=345 cases; n=670 controls). For the EA sample, single-nucleotide polymorphisms (SNPs) with the strongest statistical evidence for association included rs5907577 in an intergenic region at Xq27.1 (P=1.6 × 10−6) and rs10193871 in NAP5 at 2q21.2 (P=9.8 × 10−6). For the AA sample, SNPs with the strongest statistical evidence for association included rs2111504 in DPY19L3 at 19q13.11 (P=1.5 × 10−6) and rs2769605 in NTRK2 at 9q21.33 (P=4.5 × 10−5). We also investigated whether we could provide support for three regions previously associated with BD, and we showed that the ANK3 region replicates in our sample, along with some support for C15Orf53; other evidence implicates BD candidate genes such as SLITRK2. We also tested the hypothesis that BD susceptibility variants exhibit genetic background-dependent effects. SNPs with the strongest statistical evidence for genetic background effects included rs11208285 in ROR1 at 1p31.3 (P=1.4 × 10−6), rs4657247 in RGS5 at 1q23.3 (P=4.1 × 10−6), and rs7078071 in BTBD16 at 10q26.13 (P=4.5 × 10−6). This study is the first to conduct GWA of BD in individuals of AA and suggests that genetic variations that contribute to BD may vary as a function of ancestry.

Genome-wide association identifies OBFC1 as a locus involved in human leukocyte telomere biology

Telomeres are engaged in a host of cellular functions, and their length is regulated by multiple genes. Telomere shortening, in the course of somatic cell replication, ultimately leads to replicative senescence. In humans, rare mutations in genes that regulate telomere length have been identified in monogenic diseases such as dyskeratosis congenita and idiopathic pulmonary fibrosis, which are associated with shortened leukocyte telomere length (LTL) and increased risk for aplastic anemia. Shortened LTL is observed in a host of aging-related complex genetic diseases and is associated with diminished survival in the elderly. We report results of a genome-wide association study of LTL in a consortium of four observational studies (n = 3,417 participants with LTL and genome-wide genotyping). SNPs in the regions of the oligonucleotide/oligosaccharide-binding folds containing one gene (OBFC1; rs4387287; P = 3.9 × 10−9) and chemokine (C-X-C motif) receptor 4 gene (CXCR4; rs4452212; P = 2.9 × 10−8) were associated with LTL at a genome-wide significance level (P < 5 × 10−8). We attempted replication of the top SNPs at these loci through de novo genotyping of 1,893 additional individuals and in silico lookup in another observational study (n = 2,876), and we confirmed the association findings for OBFC1 but not CXCR4. In addition, we confirmed the telomerase RNA component (TERC) as a gene associated with LTL (P = 1.1 × 10−5). The identification of OBFC1 through genome-wide association as a locus for interindividual variation in LTL in the general population advances the understanding of telomere biology in humans and may provide insights into aging-related disorders linked to altered LTL dynamics.

Singleton deletions throughout the genome increase risk of bipolar disorder.

An overall burden of rare structural genomic variants has not been reported in bipolar disorder (BD), although there have been reports of cases with microduplication and microdeletion. Here, we present a genome-wide copy number variant (CNV) survey of 1001 cases and 1034 controls using the Affymetrix single nucleotide polymorphism (SNP) 6.0 SNP and CNV platform. Singleton deletions (deletions that appear only once in the dataset) more than 100 kb in length are present in 16.2% of BD cases in contrast to 12.3% of controls (permutation P=0.007). This effect was more pronounced for age at onset of mania ⩽18 years old. Our results strongly suggest that BD can result from the effects of multiple rare structural variants.

Enrichment of cis-regulatory gene expression SNPs and methylation quantitative trait loci among bipolar disorder susceptibility variants

We conducted a systematic study of top susceptibility variants from a genome-wide association (GWA) study of bipolar disorder to gain insight into the functional consequences of genetic variation influencing disease risk. We report here the results of experiments to explore the effects of these susceptibility variants on DNA methylation and mRNA expression in human cerebellum samples. Among the top susceptibility variants, we identified an enrichment of cis regulatory loci on mRNA expression (eQTLs), and a significant excess of quantitative trait loci for DNA CpG methylation, hereafter referred to as methylation quantitative trait loci (mQTLs). Bipolar disorder susceptibility variants that cis regulate both cerebellar expression and methylation of the same gene are a very small proportion of bipolar disorder susceptibility variants. This finding suggests that mQTLs and eQTLs provide orthogonal ways of functionally annotating genetic variation within the context of studies of pathophysiology in brain. No lymphocyte mQTL enrichment was found, suggesting that mQTL enrichment was specific to the cerebellum, in contrast to eQTLs. Separately, we found that using mQTL information to restrict the number of single-nucleotide polymorphisms studied enhances our ability to detect a significant association. With this restriction a priori informed by the observed functional enrichment, we identified a significant association (rs12618769, Pbonferroni<0.05) from two other GWA studies (TGen+GAIN; 2191 cases and 1434 controls) of bipolar disorder, which we replicated in an independent GWA study (WTCCC). Collectively, our findings highlight the importance of integrating functional annotation of genetic variants for gene expression and DNA methylation to advance the biological understanding of bipolar disorder.

Longitudinal genome-wide association of cardiovascular disease risk factors in the Bogalusa heart study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10−24), and rs445925 at APOE with LDL levels (combined P = 8.7×10−19). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.

Genome-wide meta-analysis points to CTC1 and ZNF676 as genes regulating telomere homeostasis in humans

Leukocyte telomere length (LTL) is associated with a number of common age-related diseases and is a heritable trait. Previous genome-wide association studies (GWASs) identified two loci on chromosomes 3q26.2 (TERC) and 10q24.33 (OBFC1) that are associated with the inter-individual LTL variation. We performed a meta-analysis of 9190 individuals from six independent GWAS and validated our findings in 2226 individuals from four additional studies. We confirmed previously reported associations with OBFC1 (rs9419958 P = 9.1 × 10−11) and with the telomerase RNA component TERC (rs1317082, P = 1.1 × 10−8). We also identified two novel genomic regions associated with LTL variation that map near a conserved telomere maintenance complex component 1 (CTC1; rs3027234, P = 3.6 × 10−8) on chromosome17p13.1 and zinc finger protein 676 (ZNF676; rs412658, P = 3.3 × 10−8) on 19p12. The minor allele of rs3027234 was associated with both shorter LTL and lower expression of CTC1. Our findings are consistent with the recent observations that point mutations in CTC1 cause short telomeres in both Arabidopsis and humans affected by a rare Mendelian syndrome. Overall, our results provide novel insights into the genetic architecture of inter-individual LTL variation in the general population.

Influence of genetic polymorphisms on the effect of high- and standard-dose clopidogrel after percutaneous coronary intervention: the GIFT (Genotype Information and Functional Testing) study.

OBJECTIVES This study sought to evaluate the influence of single nucleotide polymorphisms (SNPs) on the pharmacodynamic effect of high- or standard-dose clopidogrel after percutaneous coronary intervention (PCI). BACKGROUND There is a lack of prospective, multicenter data regarding the effect of different genetic variants on clopidogrel pharmacodynamics over time in patients undergoing PCI. METHODS The GRAVITAS (Gauging Responsiveness with A VerifyNow assay-Impact on Thrombosis And Safety) trial screened patients with platelet function testing after PCI and randomly assigned those with high on-treatment reactivity (OTR) to either high- or standard-dose clopidogrel; a cohort of patients without high OTR were also followed. DNA samples obtained from 1,028 patients were genotyped for 41 SNPs in 17 genes related to platelet reactivity. After adjusting for clinical characteristics, the associations between the SNPs and OTR using linear regression were evaluated. RESULTS CYP2C19*2 was significantly associated with OTR at 12 to 24 h (R(2) = 0.07, p = 2.2 × 10(-15)), 30 days (R(2) = 0.10, p = 1.3 × 10(-7)), and 6 months after PCI (R(2) = 0.07, p = 1.9 × 10(-11)), whereas PON1, ABCB1 3435 C→T, and other candidate SNPs were not. Carriers of 1 and 2 reduced-function CYP2C19 alleles were significantly more likely to display persistently high OTR at 30 days and 6 months, irrespective of treatment assignment. The portion of the risk of persistently high OTR at 30 days attributable to reduced-function CYP2C19 allele carriage was 5.2% in the patients randomly assigned to high-dose clopidogrel. CONCLUSIONS CYP2C19, but not PON1 or ABCB1, is a significant determinant of the pharmacodynamic effects of clopidogrel, both early and late after PCI. In patients with high OTR identified by platelet function testing, the CYP2C19 genotype provides limited incremental information regarding the risk of persistently high reactivity with clopidogrel 150-mg maintenance dosing. (Genotype Information and Functional Testing Study [GIFT]; NCT00992420).