9p21 DNA variants associated with coronary artery disease impair interferon-γ signalling response

Regulatory variation in a TBX5 enhancer leads to isolated congenital heart disease

Recent studies have identified the genetic underpinnings of a growing number of diseases through targeted exome sequencing. However, this strategy ignores the large component of the genome that does not code for proteins, but is nonetheless biologically functional. To address the possible involvement of regulatory variation in congenital heart diseases (CHDs), we searched for regulatory mutations impacting the activity of TBX5, a dosage-dependent transcription factor with well-defined roles in the heart and limb development that has been associated with the Holt-Oram syndrome (heart-hand syndrome), a condition that affects 1/100 000 newborns. Using a combination of genomics, bioinformatics and mouse genetic engineering, we scanned ∼700 kb of the TBX5 locus in search of cis-regulatory elements. We uncovered three enhancers that collectively recapitulate the endogenous expression pattern of TBX5 in the developing heart. We re-sequenced these enhancer elements in a cohort of non-syndromic patients with isolated atrial and/or ventricular septal defects, the predominant cardiac defects of the Holt-Oram syndrome, and identified a patient with a homozygous mutation in an enhancer ∼90 kb downstream of TBX5. Notably, we demonstrate that this single-base-pair mutation abrogates the ability of the enhancer to drive expression within the heart in vivo using both mouse and zebrafish transgenic models. Given the population-wide frequency of this variant, we estimate that 1/100 000 individuals would be homozygous for this variant, highlighting that a significant number of CHD associated with TBX5 dysfunction might arise from non-coding mutations in TBX5 heart enhancers, effectively decoupling the heart and hand phenotypes of the Holt-Oram syndrome.

Understanding cardiovascular disease: a journey through the genome (and what we found there)

Cardiovascular disease (CVD) is a major cause of mortality and hospitalization worldwide. Several risk factors have been identified that are strongly associated with the development of CVD. However, these explain only a fraction of cases, and the focus of research into the causes underlying the unexplained risk has shifted first to genetics and more recently to genomics. A genetic contribution to CVD has long been recognized; however, with the exception of certain conditions that show Mendelian inheritance, it has proved more challenging than anticipated to identify the precise genomic components responsible for the development of CVD. Genome-wide association studies (GWAS) have provided information about specific genetic variations associated with disease, but these are only now beginning to reveal the underlying molecular mechanisms. To fully understand the biological implications of these associations, we need to relate them to the exquisite, multilayered regulation of protein expression, which includes chromatin remodeling, regulatory elements, microRNAs and alternative splicing. Understanding how the information contained in the DNA relates to the operation of these regulatory layers will allow us not only to better predict the development of CVD but also to develop more effective therapies.

11q13 is a susceptibility locus for hormone receptor positive breast cancer

A recent two-stage genome-wide association study (GWAS) identified five novel breast cancer susceptibility loci on chromosomes 9, 10, and 11. To provide more reliable estimates of the relative risk associated with these loci and investigate possible heterogeneity by subtype of breast cancer, we genotyped the variants rs2380205, rs1011970, rs704010, rs614367, and rs10995190 in 39 studies from the Breast Cancer Association Consortium (BCAC), involving 49,608 cases and 48,772 controls of predominantly European ancestry. Four of the variants showed clear evidence of association (P ≤ 3 × 10(-9) ) and weak evidence was observed for rs2380205 (P = 0.06). The strongest evidence was obtained for rs614367, located on 11q13 (per-allele odds ratio 1.21, P = 4 × 10(-39) ). The association for rs614367 was specific to estrogen receptor (ER)-positive disease and strongest for ER plus progesterone receptor (PR)-positive breast cancer, whereas the associations for the other three loci did not differ by tumor subtype.

Genome-wide pharmacogenomic study of citalopram-induced side effects in STAR*D

Affecting about 1 in 12 Americans annually, depression is a leading cause of the global disease burden. While a range of effective antidepressants are now available, failure and relapse rates remain substantial, with intolerable side effect burden the most commonly cited reason for discontinuation. Thus, understanding individual differences in susceptibility to antidepressant therapy side effects will be essential to optimize depression treatment. Here we perform genome-wide association studies (GWAS) to identify genetic variation influencing susceptibility to citalopram-induced side effects. The analysis sample consisted of 1762 depression patients, successfully genotyped for 421K single-nucleotide polymorphisms (SNPs), from the Sequenced Treatment Alternatives to Relieve Depression (STAR(*)D) study. Outcomes included five indicators of citalopram side effects: general side effect burden, overall tolerability, sexual side effects, dizziness and vision/hearing side effects. Two SNPs met our genome-wide significance criterion (q<0.1), ensuring that, on average, only 10% of significant findings are false discoveries. In total, 12 additional SNPs demonstrated suggestive associations (q<0.5). The top finding was rs17135437, an intronic SNP within EMID2, mediating the effects of citalopram on vision/hearing side effects (P=3.27 × 10(-8), q=0.026). The second genome-wide significant finding, representing a haplotype spanning ∼30 kb and eight genotyped SNPs in a gene desert on chromosome 13, was associated with general side effect burden (P=3.22 × 10(-7), q=0.096). Suggestive findings were also found for SNPs at LAMA1, AOX2P, EGFLAM, FHIT and RTP2. Although our findings require replication and functional validation, this study demonstrates the potential of GWAS to discover genes and pathways that potentially mediate adverse effects of antidepressant medications.

The chromosome 9p21.3 coronary heart disease risk allele is associated with altered gene expression in normal heart and vascular tissues

Genome-wide association studies have identified a coronary artery disease (CAD) risk locus in a non-coding region at 9p21.3, the nearest genes being CDKN2A and CDKN2B. To understand the pathways by which this locus might influence CAD susceptibility, we investigated associations between the 9p21.3 risk genotype and global gene expression in heart tissue from donors with no diagnosed heart disease (n = 108, predominant cause of death, cerebral vascular accident) and in carotid plaque (n = 106), aorta (n = 104) and mammary artery (n = 88) tissues from heart valve and carotid endarterectomy patients. Genotyping was performed with Taqman assays and Illumina arrays, and gene expression profiles generated with Affymetrix microarrays. Associations were analyzed with an additive genetic model. In heart tissue, 46 genes were putatively altered in association with the 9p21.3 risk allele (70% down-regulated, fold-change >1.1 per allele, p<0.05 adjusted for age, gender, ethnicity, cause of death). These genes were enriched for biomarkers of myocardial infarction (p = 1.53×10⁻⁹), response to wounding (p = 2.65×10⁻¹⁰) and inflammatory processes (p<1.97×10⁻⁷). Among the top 10 most down-regulated genes, 7 genes shared a set of transcription factor binding sites within conserved promoter regions (p<1.14×10⁻⁵), suggesting they may be co-regulated. Canonical pathway modelling of the most differentially expressed transcripts across all tissues (154 genes, 60% down-regulated, fold-change >1.1 per allele, p<0.01) showed that 75% of the genes could be transcriptionally regulated through the cell cycle G1 phase progression pathway (p<1.08×10⁻²⁵⁸), in which CDKN2A and CDKN2B play a regulatory role. These data suggest that the cell cycle G1 phase progression pathway is activated in individuals with the 9p21.3 risk allele. This may contribute to a proliferative phenotype that promotes adverse cardiac hypertrophy and vascular remodeling, leading to an increased CAD risk.

Genetic and functional analyses implicate the NUDT11, HNF1B, and SLC22A3 genes in prostate cancer pathogenesis

One of the central goals of human genetics is to discover the genes and pathways driving human traits. To date, most of the common risk alleles discovered through genome-wide association studies (GWAS) map to nonprotein-coding regions. Because of our relatively poorer understanding of this part of the genome, the functional consequences of trait-associated variants pose a considerable challenge. To identify the genes through which risk loci act, we hypothesized that the risk variants are regulatory elements. For each of 12 known risk polymorphisms, we evaluated the correlation between risk allele status and transcript abundance for all annotated protein-coding transcripts within a 1-Mb interval. A total of 103 transcripts were evaluated in 662 prostate tissue samples [normal (n = 407) and tumor (n = 255)] from 483 individuals [European Americans (n = 233), Japanese (n = 127), and African Americans (n = 123)]. In a pooled analysis, 4 of the 12 risk variants were strongly associated with five transcripts (NUDT11, MSMB, NCOA4, SLC22A3, and HNF1B) in histologically normal tissue (P ≤ 0.001). Although associations were also observed in tumor tissue, they tended to be more attenuated. Previously, we showed that MSMB and NCOA4 participate in prostate cancer pathogenesis. Suppressing the expression of NUDT11, SLC22A3, and HNF1B influences cellular phenotypes associated with tumor-related properties in prostate cancer cells. Taken together, the data suggest that these transcripts contribute to prostate cancer pathogenesis.

Genomic approaches towards finding cis-regulatory modules in animals

Differential gene expression is the fundamental mechanism underlying animal development and cell differentiation. However, it is a challenge to identify comprehensively and accurately the DNA sequences that are required to regulate gene expression: namely, cis-regulatory modules (CRMs). Three major features, either singly or in combination, are used to predict CRMs: clusters of transcription factor binding site motifs, non-coding DNA that is under evolutionary constraint and biochemical marks associated with CRMs, such as histone modifications and protein occupancy. The validation rates for predictions indicate that identifying diagnostic biochemical marks is the most reliable method, and understanding is enhanced by the analysis of motifs and conservation patterns within those predicted CRMs.

Variants at chromosome 20 (ASIP locus) and melanoma risk

Agouti signaling protein (ASIP) locus on chromosome 20q11 is implicated, as shown by genome-wide association studies, in phenotype variation and melanoma risk. We genotyped 837 melanoma cases and 1,154 controls for 21 single nucleotide polymorphisms (SNPs) informative for 495 polymorphisms at the locus. Our data showed an increased risk of melanoma (odds ratio [OR] 1.27, 95% confidence interval [95% CI] 1.03-1.57) in carriers of the rs4911414 variant, located 120 kb upstream of ASIP. The main effect of rs4911414, as reported previously, was in tandem with a 10 kb adjacent polymorphism rs1015362; two constituted risk-associated haplotype/diplotype. Except for rs1015363, none of the 12 tagging SNPs, genotyped to cover 239.9 kb region with polymorphisms linked to rs4911414 and rs1015362, were associated with melanoma. Our data confirmed a previous association of melanoma risk (OR 1.82, 95% CI 1.37-2.41) with rs4911442, located in intron 5 of the nuclear receptor coactivator 6 (NCOA6) gene. The rs910871, one of the six variants, genotyped to cover NCOA6, showed an association with melanoma risk (OR 1.33, 95% CI 1.04-1.70). Both, rs4911442 and rs910871 were in moderate linkage with a, previously reported, risk-associated rs910873 polymorphism. A haplotype from the variants within NCOA6 showed an association with risk of melanoma (OR 1.49, 95% CI 1.17-1.88). Interaction between risk-associated polymorphisms and previously genotyped melanocortin receptor 1 (MC1R) variants, in our study, was not statistically significant. Nevertheless, the carriers of the variant alleles over the background of MC1R variants were at a higher risk than the carriers not enriched for MC1R variants. Our data confirmed the association of different variants at chromosome 20q11 with melanoma risk.

Integrative functional genomics identifies an enhancer looping to the SOX9 gene disrupted by the 17q24.3 prostate cancer risk locus

Genome-wide association studies (GWAS) are identifying genetic predisposition to various diseases. The 17q24.3 locus harbors the single nucleotide polymorphism (SNP) rs1859962 that is statistically associated with prostate cancer (PCa). It defines a 130-kb linkage disequilibrium (LD) block that lies in an ∼2-Mb gene desert area. The functional biology driving the risk associated with this LD block is unknown. Here, we integrate genome-wide chromatin landscape data sets, namely, epigenomes and chromatin openness from diverse cell types. This identifies a PCa-specific enhancer within the rs1859962 risk LD block that establishes a 1-Mb chromatin loop with the SOX9 gene. The rs8072254 and rs1859961 SNPs mapping to this enhancer impose allele-specific gene expression. The variant allele of rs8072254 facilitates androgen receptor (AR) binding driving increased enhancer activity. The variant allele of rs1859961 decreases FOXA1 binding while increasing AP-1 binding. The latter is key to imposing allele-specific gene expression. The rs8072254 variant in strong LD with the rs1859962 risk SNP can account for the risk associated with this locus, while rs1859961 is a rare variant less likely to contribute to the risk associated with this LD block. Together, our results demonstrate that multiple genetic variants mapping to a unique enhancer looping to the SOX9 oncogene can account for the risk associated with the PCa 17q24.3 locus. Allele-specific recruitment of the transcription factors androgen receptor (AR) and activating protein-1 (AP-1) account for the increased enhancer activity ascribed to this PCa-risk LD block. This further supports the notion that an integrative genomics approach can identify the functional biology disrupted by genetic risk variants.

Impact of Next-Generation Sequencing (NGS) technology on cardiovascular disease research

In recent years, hundreds of gene loci associated with multiple cardiovascular pathologies and traits have been identified through high-throughput Next-Generation Sequencing (NGS) technology. Due to the increasing efficiency and decreasing cost of NGS, rapid progresses anticipated in the field of CVD research. This review summarizes the main strategies of CV research with NGS at the level of genomics, transcriptomics, epigenetics, and proteomics.

The genetics of coronary heart disease

BACKGROUND

Coronary heart disease (CHD) is a leading cause of death worldwide, yet many areas of its pathogenesis remain unknown or poorly understood, leaving potential for novel preventive and therapeutic interventions. Recent major advances in genomic science and technology have opened new avenues of investigation in the pathogenesis of CHD, some of which are leading to clinical translation.

SOURCES OF DATA

The published literature in CHD genetics has burgeoned in the last 5 years with the reporting of genome-wide association studies (GWASs) and many other findings.

AREAS OF AGREEMENT

Identification of many genetic variants with small effects on CHD risk has been a common finding. These have included several predicted loci, such as those involved in conventional CHD risk factors (e.g. plasma lipids) and many novel loci, where their mechanism of action is unclear. The need for large, collaborative approaches to research has also become clear and is now an accepted modus operandi.

AREAS OF CONTROVERSY

The clinical utility of novel GWAS findings remains uncertain. In particular, the relative contribution of common variants of modest effect and rare variants of larger effects to risk of CHD or response to drugs is unclear.

GROWING POINTS

As a greater number of larger GWASs are conducted in CHD and its related phenotypes, much effort is being made to find translational applications for their findings. Therapeutics, prediction and pathology are major areas of research endeavour.

Human genetics and genomics a decade after the release of the draft sequence of the human genome

Substantial progress has been made in human genetics and genomics research over the past ten years since the publication of the draft sequence of the human genome in 2001. Findings emanating directly from the Human Genome Project, together with those from follow-on studies, have had an enormous impact on our understanding of the architecture and function of the human genome. Major developments have been made in cataloguing genetic variation, the International HapMap Project, and with respect to advances in genotyping technologies. These developments are vital for the emergence of genome-wide association studies in the investigation of complex diseases and traits. In parallel, the advent of high-throughput sequencing technologies has ushered in the 'personal genome sequencing' era for both normal and cancer genomes, and made possible large-scale genome sequencing studies such as the 1000 Genomes Project and the International Cancer Genome Consortium. The high-throughput sequencing and sequence-capture technologies are also providing new opportunities to study Mendelian disorders through exome sequencing and whole-genome sequencing. This paper reviews these major developments in human genetics and genomics over the past decade.

Genetics of human cardiovascular disease

Cardiovascular disease encompasses a range of conditions extending from myocardial infarction to congenital heart disease, most of which are heritable. Enormous effort has been invested in understanding the genes and specific DNA sequence variants that are responsible for this heritability. Here, we review the lessons learned for monogenic and common, complex forms of cardiovascular disease. We also discuss key challenges that remain for gene discovery and for moving from genomic localization to mechanistic insights, with an emphasis on the impact of next-generation sequencing and the use of pluripotent human cells to understand the mechanism by which genetic variation contributes to disease.

Landscapes and archipelagos: spatial organization of gene regulation in vertebrates

Vertebrate genes controlling critical developmental processes are often regulated by complex sets of global enhancer sequences, located at a distance, within neighboring gene deserts. Recent technological advances have made it possible to investigate the spatial organization of these 'regulatory landscapes'. The integration of such datasets with information on chromatin status, transcriptional activity and nuclear localization of these loci, as well as the effects of genetic modifications thereof, may bring a more comprehensive understanding of tissue- and/or stage-specific gene regulation in both normal and pathological contexts. Here, we review the impact of recent technological advances on our understanding of large-scale gene regulation in vertebrates, by focusing on paradigmatic gene loci.

Uncovering the role of genomic "dark matter" in human disease

The human genome encodes thousands of long noncoding RNAs (lncRNAs). Although most remain functionally uncharacterized biological "dark matter," lncRNAs have garnered considerable attention for their diverse roles in human biology, including developmental programs and tumor suppressor gene networks. As the number of lncRNAs associated with human disease grows, ongoing research efforts are focusing on their regulatory mechanisms. New technologies that enable enumeration of lncRNA interaction partners and determination of lncRNA structure are well positioned to drive deeper understanding of their functions and involvement in pathogenesis. In turn, lncRNAs may become targets for therapeutic intervention or new tools for biotechnology.

Common variants at 9p21 and 8q22 are associated with increased susceptibility to optic nerve degeneration in glaucoma

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

A genome-wide association study of female sexual dysfunction

BACKGROUND

Female sexual dysfunction (FSD) is an important but controversial problem with serious negative impact on women's quality of life. Data from twin studies have shown a genetic contribution to the development and maintenance of FSD.

METHODOLOGY/PRINCIPAL FINDINGS

We performed a genome-wide association study (GWAS) on 2.5 million single-nucleotide polymorphisms (SNPs) in 1,104 female twins (25-81 years of age) in a population-based register and phenotypic data on lifelong sexual functioning. Although none reached conventional genome-wide level of significance (10 × -8), we found strongly suggestive associations with the phenotypic dimension of arousal (rs13202860, P = 1.2 × 10(-7); rs1876525, P = 1.2 × 10(-7); and rs13209281 P = 8.3 × 10(-7)) on chromosome 6, around 500 kb upstream of the locus HTR1E (5-hydroxytryptamine receptor 1E) locus, related to the serotonin brain pathways. We could not replicate previously reported candidate SNPs associated with FSD in the DRD4, 5HT2A and IL-1B loci.

CONCLUSIONS/SIGNIFICANCE

We report the first GWAS of FSD symptoms in humans. This has pointed to several "risk alleles" and the implication of the serotonin and GABA pathways. Ultimately, understanding key mechanisms via this research may lead to new FSD treatments and inform clinical practice and developments in psychiatric nosology.

Primate genome gain and loss: a bone dysplasia, muscular dystrophy, and bone cancer syndrome resulting from mutated retroviral-derived MTAP transcripts

Diaphyseal medullary stenosis with malignant fibrous histiocytoma (DMS-MFH) is an autosomal-dominant syndrome characterized by bone dysplasia, myopathy, and bone cancer. We previously mapped the DMS-MFH tumor-suppressing-gene locus to chromosomal region 9p21-22 but failed to identify mutations in known genes in this region. We now demonstrate that DMS-MFH results from mutations in the most proximal of three previously uncharacterized terminal exons of the gene encoding methylthioadenosine phosphorylase, MTAP. Intriguingly, two of these MTAP exons arose from early and independent retroviral-integration events in primate genomes at least 40 million years ago, and since then, their genomic integration has gained a functional role. MTAP is a ubiquitously expressed homotrimeric-subunit enzyme critical to polyamine metabolism and adenine and methionine salvage pathways and was believed to be encoded as a single transcript from the eight previously described exons. Six distinct retroviral-sequence-containing MTAP isoforms, each of which can physically interact with archetype MTAP, have been identified. The disease-causing mutations occur within one of these retroviral-derived exons and result in exon skipping and dysregulated alternative splicing of all MTAP isoforms. Our results identify a gene involved in the development of bone sarcoma, provide evidence of the primate-specific evolution of certain parts of an existing gene, and demonstrate that mutations in parts of this gene can result in human disease despite its relatively recent origin.

Genetics of atherosclerosis

Genome-wide association studies (GWAS) from the past several years have provided the first unbiased evidence of the genes contributing to common cardiovascular disease traits in European and some Asian populations. The results not only confirmed the importance of prior knowledge, such as the central role of lipoproteins, but also revealed that there is still much to learn about the underlying mechanisms of this disease, as most of the associated genes do not appear to be involved in pathways previously connected to atherosclerosis. In this review, I focus on the common forms of the disease and look at both human and animal model studies. I summarize what was known before GWAS, highlight how the field has been changed by GWAS, and discuss future considerations, such as the limitations of GWAS and strategies that may lead to a more complete, mechanistic understanding of atherosclerosis.

Integrated approaches to functionally characterize novel factors in lipoprotein metabolism

PURPOSE OF REVIEW

To discuss if and how the combined analysis of large-scale datasets from multiple independent sources benefits the mapping of novel genetic elements with relevance to lipoprotein metabolism and allows for conclusions on underlying molecular mechanisms.

RECENT FINDINGS

Genome-wide association studies (GWAS) have identified numerous genomic loci associated with plasma lipid levels and cardiovascular disease. Yet, despite being highly successful in mapping novel loci the GWAS approach falls short to systematically extract functional information from genomic data. With the aim to complement GWAS for a better insight into disease mechanisms and identification of the most promising targets for drug development, a number of high-throughput functional genomics strategies have now been applied. These include computational approaches, consideration of gene-gene and gene-environment interactions, as well as unbiased gene-expression analyses in relevant tissues. For a limited number of loci, mechanistic insight has been gained through in-vitro and in-vivo studies by knockdown and overexpression of candidate genes.

SUMMARY

The integration of GWAS data with existing functional genomics strategies has contributed to ascertain the relevance of a number of novel factors for lipoprotein biology and disease. However, technologies are warranted that provide a more systematic insight into the molecular function and pathogenic relevance of promising candidate genes.

Widespread site-dependent buffering of human regulatory polymorphism

The average individual is expected to harbor thousands of variants within non-coding genomic regions involved in gene regulation. However, it is currently not possible to interpret reliably the functional consequences of genetic variation within any given transcription factor recognition sequence. To address this, we comprehensively analyzed heritable genome-wide binding patterns of a major sequence-specific regulator (CTCF) in relation to genetic variability in binding site sequences across a multi-generational pedigree. We localized and quantified CTCF occupancy by ChIP-seq in 12 related and unrelated individuals spanning three generations, followed by comprehensive targeted resequencing of the entire CTCF–binding landscape across all individuals. We identified hundreds of variants with reproducible quantitative effects on CTCF occupancy (both positive and negative). While these effects paralleled protein–DNA recognition energetics when averaged, they were extensively buffered by striking local context dependencies. In the significant majority of cases buffering was complete, resulting in silent variants spanning every position within the DNA recognition interface irrespective of level of binding energy or evolutionary constraint. The prevalence of complex partial or complete buffering effects severely constrained the ability to predict reliably the impact of variation within any given binding site instance. Surprisingly, 40% of variants that increased CTCF occupancy occurred at positions of human–chimp divergence, challenging the expectation that the vast majority of functional regulatory variants should be deleterious. Our results suggest that, even in the presence of “perfect” genetic information afforded by resequencing and parallel studies in multiple related individuals, genomic site-specific prediction of the consequences of individual variation in regulatory DNA will require systematic coupling with empirical functional genomic measurements.

A comprehensive understanding of the contribution of individual genome sequences to disease and quantitative traits will require the general ability to predict consequences of genetic variation in non-protein-coding regions, particularly those involved in gene regulation. Here we tested the power to predict such consequences when presented with “complete” information encompassing the genomic DNA binding site patterns of a well-studied regulatory protein across multiple related individuals, coupled with all individual genome sequences at the binding positions. We find that, while there is reasonable ability to predict the average effects of variation within the consensus recognition sequence of a transcriptional regulator, it is not possible to determine reliably the consequences of variation at any given genomic instance. This suggests that the interpretation of individual genome sequences will require comprehensive complementation with functional genomic studies.

Identification of a chronic obstructive pulmonary disease genetic determinant that regulates HHIP

Multiple intergenic single-nucleotide polymorphisms (SNPs) near hedgehog interacting protein (HHIP) on chromosome 4q31 have been strongly associated with pulmonary function levels and moderate-to-severe chronic obstructive pulmonary disease (COPD). However, whether the effects of variants in this region are related to HHIP or another gene has not been proven. We confirmed genetic association of SNPs in the 4q31 COPD genome-wide association study (GWAS) region in a Polish cohort containing severe COPD cases and healthy smoking controls (P = 0.001 to 0.002). We found that HHIP expression at both mRNA and protein levels is reduced in COPD lung tissues. We identified a genomic region located ∼85 kb upstream of HHIP which contains a subset of associated SNPs, interacts with the HHIP promoter through a chromatin loop and functions as an HHIP enhancer. The COPD risk haplotype of two SNPs within this enhancer region (rs6537296A and rs1542725C) was associated with statistically significant reductions in HHIP promoter activity. Moreover, rs1542725 demonstrates differential binding to the transcription factor Sp3; the COPD-associated allele exhibits increased Sp3 binding, which is consistent with Sp3's usual function as a transcriptional repressor. Thus, increased Sp3 binding at a functional SNP within the chromosome 4q31 COPD GWAS locus leads to reduced HHIP expression and increased susceptibility to COPD through distal transcriptional regulation. Together, our findings reveal one mechanism through which SNPs upstream of the HHIP gene modulate the expression of HHIP and functionally implicate reduced HHIP gene expression in the pathogenesis of COPD.

Strategies beyond genome-wide association studies for atherosclerosis

Atherosclerotic diseases, including coronary artery disease (CAD) and myocardial infarction (MI), are the leading causes of death in the world. The genetic basis of CAD and MI, which are caused by multiple interacting endogenous and exogenous factors, has gained considerable interest in the last years as genome-wide association studies (GWASs) have identified many new susceptibility loci for CAD and MI, and the underlying genes provide new insights into the genetic architecture of these diseases. Here we summarize the recent findings from GWASs of atherosclerosis and discuss their functional and biological implications. We also discuss the different post-GWAS strategies that are currently used for refining the location of causal variants, understanding their role, and shedding light on molecular mechanisms explaining their association to CAD. We finally discuss potential clinical translations of GWAS findings for individual risk prediction, advanced clinical strategies, and personalized treatments.

Common variants in CDKN2B-AS1 associated with optic-nerve vulnerability of glaucoma identified by genome-wide association studies in Japanese

BACKGROUND

To date, only a small portion of the genetic variation for primary open-angle glaucoma (POAG), the major type of glaucoma, has been elucidated.

METHODS AND PRINCIPAL FINDINGS

We examined our two data sets of the genome-wide association studies (GWAS) derived from a total of 2,219 Japanese subjects. First, we performed a GWAS by analyzing 653,519 autosomal common single-nucleotide polymorphisms (SNPs) in 833 POAG patients and 686 controls. As a result, five variants that passed the Bonferroni correction were identified in CDKN2B-AS1 on chromosome 9p21.3, which was already reported to be a significant locus in the Caucasian population. Moreover, we combined the data set with our previous GWAS data set derived from 411 POAG patients and 289 controls by the Mantel-Haenszel test, and all of the combined variants showed stronger association with POAG (P<5.8 × 10(-10)). We then subdivided the case groups into two subtypes based on the value of intraocular pressure (IOP)--POAG with high IOP (high pressure glaucoma, HPG) and that with normal IOP (normal pressure glaucoma, NPG)--and performed the GWAS using the two data sets, as the prevalence of NPG in Japanese is much higher than in Caucasians. The results suggested that the variants from the same CDKN2B-AS1 locus were likely to be significant for NPG patients.

CONCLUSIONS AND SIGNIFICANCE

In this study, we successfully identified POAG-associated variants in the CDKN2B-AS1 locus using a Japanese population, i.e., variants originally reported as being associated with the Caucasian population. Although we cannot rule out that the significance could be due to the differences in sample size between HPG and NPG, the variants could be associated specifically with the vulnerability of the optic nerve to IOP, which is useful for investigating the etiology of glaucoma.

Genetic basis of atherosclerosis: insights from mice and humans

Atherosclerosis is a complex and heritable disease involving multiple cell types and the interactions of many different molecular pathways. The genetic and molecular mechanisms of atherosclerosis have, in part, been elucidated by mouse models; at least 100 different genes have been shown to influence atherosclerosis in mice. Importantly, unbiased genome-wide association studies have recently identified a number of novel loci robustly associated with atherosclerotic coronary artery disease. Here, we review the genetic data elucidated from mouse models of atherosclerosis, as well as significant associations for human coronary artery disease. Furthermore, we discuss in greater detail some of these novel human coronary artery disease loci. The combination of mouse and human genetics has the potential to identify and validate novel genes that influence atherosclerosis, some of which may be candidates for new therapeutic approaches.

Sex differences in disease risk from reported genome-wide association study findings

Men and women differ in susceptibility to many diseases and in responses to treatment. Recent advances in genome-wide association studies (GWAS) provide a wealth of data for associating genetic profiles with disease risk; however, in general, these data have not been systematically probed for sex differences in gene-disease associations. Incorporating sex into the analysis of GWAS results can elucidate new relationships between single nucleotide polymorphisms (SNPs) and human disease. In this study, we performed a sex-differentiated analysis on significant SNPs from GWAS data of the seven common diseases studied by the Wellcome Trust Case Control Consortium. We employed and compared three methods: logistic regression, Woolf's test of heterogeneity, and a novel statistical metric that we developed called permutation method to assess sex effects (PMASE). After correction for false discovery, PMASE finds SNPs that are significantly associated with disease in only one sex. These sexually dimorphic SNP-disease associations occur in Coronary Artery Disease and Crohn's Disease. GWAS analyses that fail to consider sex-specific effects may miss discovering sexual dimorphism in SNP-disease associations that give new insights into differences in disease mechanism between men and women.

Advancing cardiovascular research

Over the past 50 years, we have seen dramatic changes in cardiovascular science and clinical care, accompanied by marked declines in the morbidity and mortality. Nonetheless, cardiovascular disease remains the leading cause of death and disability in the world, and its nature is changing as Americans become older, fatter, and ethnically more diverse. Instead of young or middle-aged men with ST-segment elevation myocardial infarction, the "typical" cardiac patient now presents with acute coronary syndrome or with complications related to chronic hypertension or ischemic heart disease, including heart failure, sudden death, and atrial fibrillation. Analogously, structural heart disease is now dominated by degenerative valve or congenital disease, far more common than rheumatic disease. The changing clinical scene presents cardiovascular scientists with a number of opportunities and challenges, including taking advantage of high-throughput technologies to elucidate complex disease mechanisms, accelerating development and implementation of evidence-based strategies, assessing evolving technologies of unclear value, addressing a global epidemic of cardiovascular disease, and maintaining high levels of innovation in a time of budgetary constraint and economic turmoil.

Molecular genetic studies of complex phenotypes

The approach to molecular genetic studies of complex phenotypes evolved considerably during the recent years. The candidate gene approach, which is restricted to an analysis of a few single-nucleotide polymorphisms (SNPs) in a modest number of cases and controls, has been supplanted by the unbiased approach of genome-wide association studies (GWAS), wherein a large number of tagger SNPs are typed in many individuals. GWAS, which are designed on the common disease-common variant hypothesis (CD-CV), identified several SNPs and loci for complex phenotypes. However, the alleles identified through GWAS are typically not causative but rather in linkage disequilibrium (LD) with the true causal variants. The common alleles, which may not capture the uncommon and rare variants, account only for a fraction of heritability of the complex traits. Hence, the focus is being shifted to rare variants-common disease (RV-CD) hypothesis, surmising that rare variants exert large effect sizes on the phenotype. In conjunctional with this conceptual shift, technologic advances in DNA sequencing techniques have dramatically enhanced whole genome or whole exome sequencing capacity. The sequencing approach affords identification of not only the rare but also the common variants. The approach-whether used in complementation with GWAS or as a stand-alone approach-could define the genetic architecture of the complex phenotypes. Robust phenotyping and large-scale sequencing studies are essential to extract the information content of the vast number of DNA sequence variants (DSVs) in the genome. To garner meaningful clinical information and link the genotype to a phenotype, the identification and characterization of a large number of causal fields beyond the information content of DNA sequence variants would be necessary. This review provides an update on the current progress and limitations in identifying DSVs that are associated with phenotypic effects.

Association between chromosome 9p21 variants and the ankle-brachial index identified by a meta-analysis of 21 genome-wide association studies

BACKGROUND

Genetic determinants of peripheral arterial disease (PAD) remain largely unknown. To identify genetic variants associated with the ankle-brachial index (ABI), a noninvasive measure of PAD, we conducted a meta-analysis of genome-wide association study data from 21 population-based cohorts.

METHODS AND RESULTS

Continuous ABI and PAD (ABI ≤0.9) phenotypes adjusted for age and sex were examined. Each study conducted genotyping and imputed data to the ≈2.5 million single nucleotide polymorphisms (SNPs) in HapMap. Linear and logistic regression models were used to test each SNP for association with ABI and PAD using additive genetic models. Study-specific data were combined using fixed effects inverse variance weighted meta-analyses. There were a total of 41 692 participants of European ancestry (≈60% women, mean ABI 1.02 to 1.19), including 3409 participants with PAD and with genome-wide association study data available. In the discovery meta-analysis, rs10757269 on chromosome 9 near CDKN2B had the strongest association with ABI (β=-0.006, P=2.46×10(-8)). We sought replication of the 6 strongest SNP associations in 5 population-based studies and 3 clinical samples (n=16 717). The association for rs10757269 strengthened in the combined discovery and replication analysis (P=2.65×10(-9)). No other SNP associations for ABI or PAD achieved genome-wide significance. However, 2 previously reported candidate genes for PAD and 1 SNP associated with coronary artery disease were associated with ABI: DAB21P (rs13290547, P=3.6×10(-5)), CYBA (rs3794624, P=6.3×10(-5)), and rs1122608 (LDLR, P=0.0026).

CONCLUSIONS

Genome-wide association studies in more than 40 000 individuals identified 1 genome wide significant association on chromosome 9p21 with ABI. Two candidate genes for PAD and 1 SNP for coronary artery disease are associated with ABI.

Genomic research to identify novel pathways in the development of abdominal aortic aneurysm

Abdominal aortic aneurysm (AAA) is a common disease with a large heritable component. There is a need to improve our understanding of AAA pathogenesis in order to develop novel treatment paradigms. Genomewide association studies have revolutionized research into the genetic variants that underpin the development of many complex diseases including AAA. This article reviews the progress that has been made to date in this regard, including mechanisms by which loci identified by GWAS may contribute to the development of AAA. It also highlights potential post-GWAS analytical strategies to improve our understanding of the disease further.

Cdkn2a is an atherosclerosis modifier locus that regulates monocyte/macrophage proliferation

OBJECTIVE

Common genetic variants in a 58-kb region of chromosome 9p21, near the CDKN2A/CDKN2B tumor suppressor locus, are strongly associated with coronary artery disease. However, the underlying mechanism of action remains unknown.

METHODS AND RESULTS

We previously reported a congenic mouse model harboring an atherosclerosis susceptibility locus and the region of homology with the human 9p21 locus. Microarray and transcript-specific expression analyses showed markedly decreased Cdkn2a expression, including both p16(INK4a) and p19(ARF), but not Cdkn2b (p15(INK4b)), in macrophages derived from congenic mice compared with controls. Atherosclerosis studies in subcongenic strains revealed genetic complexity and narrowed 1 locus to a small interval including Cdkn2a/b. Bone marrow (BM) transplantation studies implicated myeloid lineage cells as the culprit cell type, rather than resident vascular cells. To directly test the role of BM-derived Cdkn2a transcripts in atherogenesis and inflammatory cell proliferation, we performed a transplantation study using Cdkn2a(-/-) cells in the Ldlr(-/-) mouse model. Cdkn2a-deficient BM recipients exhibited accelerated atherosclerosis, increased Ly6C proinflammatory monocytes, and increased monocyte/macrophage proliferation compared with controls.

CONCLUSION

These data provide a plausible mechanism for accelerated atherogenesis in susceptible congenic mice, involving decreased expression of Cdkn2a and increased proliferation of monocyte/macrophages, with possible relevance to the 9p21 human locus.

Transcriptional enhancers in development and disease

Distal transcription enhancers are cis-regulatory elements that promote gene expression, enabling spatiotemporal control of genetic programs such as those required in metazoan developmental processes. Because of their importance, their disruption can lead to disease.

Long-range DNA looping and gene expression analyses identify DEXI as an autoimmune disease candidate gene

The chromosome 16p13 region has been associated with several autoimmune diseases, including type 1 diabetes (T1D) and multiple sclerosis (MS). CLEC16A has been reported as the most likely candidate gene in the region, since it contains the most disease-associated single-nucleotide polymorphisms (SNPs), as well as an imunoreceptor tyrosine-based activation motif. However, here we report that intron 19 of CLEC16A, containing the most autoimmune disease-associated SNPs, appears to behave as a regulatory sequence, affecting the expression of a neighbouring gene, DEXI. The CLEC16A alleles that are protective from T1D and MS are associated with increased expression of DEXI, and no other genes in the region, in two independent monocyte gene expression data sets. Critically, using chromosome conformation capture (3C), we identified physical proximity between the DEXI promoter region and intron 19 of CLEC16A, separated by a loop of >150 kb. In reciprocal experiments, a 20 kb fragment of intron 19 of CLEC16A, containing SNPs associated with T1D and MS, as well as with DEXI expression, interacted with the promotor region of DEXI but not with candidate DNA fragments containing other potential causal genes in the region, including CLEC16A. Intron 19 of CLEC16A is highly enriched for transcription-factor-binding events and markers associated with enhancer activity. Taken together, these data indicate that although the causal variants in the 16p13 region lie within CLEC16A, DEXI is an unappreciated autoimmune disease candidate gene, and illustrate the power of the 3C approach in progressing from genome-wide association studies results to candidate causal genes.

HERC2 rs12913832 modulates human pigmentation by attenuating chromatin-loop formation between a long-range enhancer and the OCA2 promoter

Pigmentation of skin, eye, and hair reflects some of the most evident common phenotypes in humans. Several candidate genes for human pigmentation are identified. The SNP rs12913832 has strong statistical association with human pigmentation. It is located within an intron of the nonpigment gene HERC2, 21 kb upstream of the pigment gene OCA2, and the region surrounding rs12913832 is highly conserved among animal species. However, the exact functional role of HERC2 rs12913832 in human pigmentation is unknown. Here we demonstrate that the HERC2 rs12913832 region functions as an enhancer regulating OCA2 transcription. In darkly pigmented human melanocytes carrying the rs12913832 T-allele, we detected binding of the transcription factors HLTF, LEF1, and MITF to the HERC2 rs12913832 enhancer, and a long-range chromatin loop between this enhancer and the OCA2 promoter that leads to elevated OCA2 expression. In contrast, in lightly pigmented melanocytes carrying the rs12913832 C-allele, chromatin-loop formation, transcription factor recruitment, and OCA2 expression are all reduced. Hence, we demonstrate that allelic variation of a common noncoding SNP located in a distal regulatory element not only disrupts the regulatory potential of this element but also affects its interaction with the relevant promoter. We provide the key mechanistic insight that allele-dependent differences in chromatin-loop formation (i.e., structural differences in the folding of gene loci) result in differences in allelic gene expression that affects common phenotypic traits. This concept is highly relevant for future studies aiming to unveil the functional basis of genetically determined phenotypes, including diseases.

Transcriptional and epigenetic control of T helper cell specification: molecular mechanisms underlying commitment and plasticity

T helper cell differentiation occurs in the context of the extracellular cytokine milieu evoked by diverse microbes and other pathogenic stimuli along with T cell receptor stimulation. The culmination of these signals results in specification of T helper lineages, which occurs through the combinatorial action of multiple transcription factors that establish distinctive transcriptomes. In this manner, inducible, but constitutively active, master regulators work in conjunction with factors such as the signal transducer and activator of transcriptions (STATs) that sense the extracellular environment. The acquisition of a distinctive transcriptome also depends on chromatin modifications that impact key cis elements as well as the changes in global genomic organization. Thus, signal transduction and epigenetics are linked in these processes of differentiation. In this review, recent advances in understanding T helper lineage specification and deciphering the action of transcription factors are summarized with emphasis on comprehensive views of the dynamic T cell epigenome.

9p21 and the Genetic Revolution for Coronary Artery Disease

BACKGROUND

It has long been recognized that 50% of the susceptibility for coronary artery disease (CAD) is due to predisposing genetic factors. Comprehensive prevention is likely to require knowledge of these genetic factors.

CONTENT

Using a genomewide association study (GWAS), the Ottawa Heart Genomic Study and the deCODE group simultaneously identified the first genetic risk variant, at chromosome 9p21. The 9p21 variant became the first risk factor to be identified since 1964. 9p21 occurs in 75% of the population except for African Americans and is associated with a 25% increased risk for CAD with 1 copy and a 50% increased risk with 2 copies. Perhaps the most remarkable finding is that 9p21 is independent of all known risk factors, indicating there are factors contributing to the pathogenesis of CAD that are yet unknown. 9p21 in individuals with premature CAD is associated with a 2-fold increase in risk, similar to that of smoking and cholesterol. Routine genetic testing will probably remain controversial until a specific treatment is developed. Over a period of 5 years, however, GWASs have identified 30 genetic variants for CAD risk, of which only 6 act through the known risk factors.

SUMMARY

The 9p21 variant has now been established as an independent risk factor for CAD and, along with the additional 29 risk genetic variants recently identified, is likely to provide the thrust for genetic testing and personalized medicine in the near future.

Genomewide Association Studies in Cardiovascular Disease—An Update 2011

BACKGROUND

Genomewide association studies have led to an enormous boost in the identification of susceptibility genes for cardiovascular diseases. This review aims to summarize the most important findings of recent years.

CONTENT

We have carefully reviewed the current literature (PubMed search terms: “genome wide association studies,” “genetic polymorphism,” “genetic risk factors,” “association study” in connection with the respective diseases, “risk score,” “transcriptome”).

SUMMARY

Multiple novel genetic loci for such important cardiovascular diseases as myocardial infarction, hypertension, heart failure, stroke, and hyperlipidemia have been identified. Given that many novel genetic risk factors lie within hitherto-unsuspected genes or influence gene expression, these findings have inspired discoveries of biological function. Despite these successes, however, only a fraction of the heritability for most cardiovascular diseases has been explained thus far. Forthcoming techniques such as whole-genome sequencing will be important to close the gap of missing heritability.

Variability of platelet indices and function: acquired and genetic factors

Each individual has an inherent variable risk of bleeding linked to genetic or acquired abnormal platelet number or platelet dysfunction. In contrast, it is less obvious that the variability of platelet phenotypes (number, mean platelet volume, function) may contribute to the variable individual risk of thrombosis. Interindividual variability of platelet indices or function may be either due to acquired factors, such as age, sex, metabolic variables, smoke, dietary habits, and ongoing inflammation, or due to genetic factors. Acquired variables explain a small portion of the heterogeneity of platelet parameters. Genetic factors, instead, appear to play a major role, although a consistent portion of such a genetic variance has not yet been attributed to any specific genetic factor, possibly due to the high number of DNA loci potentially involved and to the limited effect size of each individual SNP. A portion of variance remains thus unexplained, also due to variability of test performance. A major contradiction in present platelet knowledge is, indeed, the difficulty to reconcile the universally accepted importance of platelet indices or function and the lack of reliable platelet parameters in cardiovascular risk prediction models. Trials on antiplatelet drugs were generally designed to select a homogeneous sample, whose results could be applied to an "average subject," tending to exclude the deviation/extreme values. As the current indications for antiplatelet treatment in primary or secondary prevention of ischemic vascular disease still derive from the results of such clinical trials where platelet function and its variability was not investigated, we cannot at present rely upon any current platelet test to either initiate, or monitor, or modify or stop treatment with any antiplatelet drug. Evidence is, however, increasing that traditional platelet aggregometry and other more recently developed platelet function assays could be useful to optimize antiplatelet therapy and to predict major adverse cardiac events.The observation of interindividual differences in platelet response to antiplatelet drugs has enlarged the spectrum and the possible clinical relevance of the variability of platelet indices or function. The development of "personalized medicine" will benefit from the concepts discussed in this chapter.

Genomics: is it ready for primetime?

The next decade will focus on identifying the missing heritability of coronary artery disease (CAD). This process will involve a more comprehensive interrogation of common single nucleotide polymorphisms (SNPs) that impart modest biologic effect and an interrogation of rare SNPs that impart profound biologic effect. In parallel, an investigation of the underlying biology of the described association will likely yield novel pathways that provide therapeutic targets. Once we obtain a more complete inventory of sequence variation that predisposes to CAD, a more realistic assessment of the role of genetic risk scoring allied with standard risk algorithms will be possible.