Gene-centric genomewide association study via entropy

Genes are the functional units in most organisms. Compared to genetic variants located outside genes, genic variants are more likely to affect disease risk. The development of the human HapMap project provides an unprecedented opportunity for genetic association studies at the genomewide level for elucidating disease etiology. Currently, most association studies at the single-nucleotide polymorphism (SNP) or the haplotype level rely on the linkage information between SNP markers and disease variants, with which association findings are difficult to replicate. Moreover, variants in genes might not be sufficiently covered by currently available methods. In this article, we present a gene-centric approach via entropy statistics for a genomewide association study to identify disease genes. The new entropy-based approach considers genic variants within one gene simultaneously and is developed on the basis of a joint genotype distribution among genetic variants for an association test. A grouping algorithm based on a penalized entropy measure is proposed to reduce the dimension of the test statistic. Type I error rates and power of the entropy test are evaluated through extensive simulation studies. The results indicate that the entropy test has stable power under different disease models with a reasonable sample size. Compared to single SNP-based analysis, the gene-centric approach has greater power, especially when there is more than one disease variant in a gene. As the genomewide genic SNPs become available, our entropy-based gene-centric approach would provide a robust and computationally efficient way for gene-based genomewide association study.     

Mouse models and the interpretation of human GWAS in type 2 diabetes and obesity

Within the last 3 years, genome-wide association studies (GWAS) have had unprecedented success in identifying loci that are involved in common diseases. For example, more than 35 susceptibility loci have been identified for type 2 diabetes and 32 for obesity thus far. However, the causal gene and variant at a specific linkage disequilibrium block is often unclear. Using a combination of different mouse alleles, we can greatly facilitate the understanding of which candidate gene at a particular disease locus is associated with the disease in humans, and also provide functional analysis of variants through an allelic series, including analysis of hypomorph and hypermorph point mutations, and knockout and overexpression alleles. The phenotyping of these alleles for specific traits of interest, in combination with the functional analysis of the genetic variants, may reveal the molecular and cellular mechanism of action of these disease variants, and ultimately lead to the identification of novel therapeutic strategies for common human diseases. In this Commentary, we discuss the progress of GWAS in identifying common disease loci for metabolic disease, and the use of the mouse as a model to confirm candidate genes and provide mechanistic insights.     

Sequence diversity in 36 candidate genes for cardiovascular disorders.

Two strategies involving whole-genome association studies have been proposed for the identification of genes involved in complex diseases. The first one seeks to characterize all common variants of human genes and to test their association with disease. The second one seeks to develop dense maps of single-nucleotide polymorphisms (SNPs) and to detect susceptibility genes through linkage disequilibrium. We performed a molecular screening of the coding and/or flanking regions of 36 candidate genes for cardiovascular diseases. All polymorphisms identified by this screening were further genotyped in 750 subjects of European descent. In the whole set of genes, the lengths explored spanned 53.8 kb in the 5' regions, 68.4 kb in exonic regions, and 13 kb in the 3' regions. The strength of linkage disequilibrium within candidate regions suggests that genomewide maps of SNPs might be efficient ways to identify new disease-susceptibility genes, provided that the maps are sufficiently dense. However, the relatively large number of polymorphisms within coding and regulatory regions of candidate genes raises the possibility that several of them might be functional and that the pattern of genotype-phenotype association might be more complex than initially envisaged, as actually has been observed in some well-characterized genes. These results argue in favor of both genomewide association studies and detailed studies of the overall sequence variation of candidate genes, as complementary approaches.     

Molecular basis of human hypertension: role of angiotensinogen.

Essential hypertension is a common human disease believed to result from the interplay of multiple genetic and environmental determinants. In genetic studies of two large panels of hypertensive sibships from widely separated geographical areas, we obtained evidence of genetic linkage between the angiotensinogen gene (AGT) and hypertension, demonstrated association of AGT molecular variants with the disease, and found significant differences in plasma concentrations of angiotensinogen among hypertensive subjects with different AGT genotypes. The corroboration and replication afforded by these results support the interpretation that molecular variants of AGT constitute inherited predispositions to essential hypertension in humans.     

Genetics of sputum gene expression in chronic obstructive pulmonary disease

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.     

Linkage analysis in the next-generation sequencing era

Linkage analysis was developed to detect excess co-segregation of the putative alleles underlying a phenotype with the alleles at a marker locus in family data. Many different variations of this analysis and corresponding study design have been developed to detect this co-segregation. Linkage studies have been shown to have high power to detect loci that have alleles (or variants) with a large effect size, i.e. alleles that make large contributions to the risk of a disease or to the variation of a quantitative trait. However, alleles with a large effect size tend to be rare in the population. In contrast, association studies are designed to have high power to detect common alleles which tend to have a small effect size for most diseases or traits. Although genome-wide association studies have been successful in detecting many new loci with common alleles of small effect for many complex traits, these common variants often do not explain a large proportion of disease risk or variation of the trait. In the past, linkage studies were successful in detecting regions of the genome that were likely to harbor rare variants with large effect for many simple Mendelian diseases and for many complex traits. However, identifying the actual sequence variant(s) responsible for these linkage signals was challenging because of difficulties in sequencing the large regions implicated by each linkage peak. Current 'next-generation' DNA sequencing techniques have made it economically feasible to sequence all exons or the whole genomes of a reasonably large number of individuals. Studies have shown that rare variants are quite common in the general population, and it is now possible to combine these new DNA sequencing methods with linkage studies to identify rare causal variants with a large effect size. A brief review of linkage methods is presented here with examples of their relevance and usefulness for the interpretation of whole-exome and whole-genome sequence data.     

Derived SNP alleles are used more frequently than ancestral alleles as risk-associated variants in common human diseases

Evolutionary aspects of the genetic architecture of common human diseases remain enigmatic. The results of more than 200 genome-wide association studies published to date were compiled in a catalog (). We used cataloged data to determine whether derived (mutant) alleles are associated with higher risk of human disease more frequently than ancestral alleles. We placed all allelic variants into ten categories of population frequency (0%-100%) in 10% increments. We then analyzed the relationship between allelic frequency, evolutionary status of the polymorphic site (ancestral versus derived), and disease risk status (risk versus protection). Given the same population frequency, derived alleles are more likely to be risk associated than ancestral alleles, as are rarer alleles. The common interpretation of this association is that negative selection prevents fixation of the risk variants. However, disease stratification as early or late onset suggests that weak selection against risk-associated alleles is unlikely a major factor shaping genetic architecture of common diseases. Our results clearly suggest that the duration of existence of an allele in a population is more important. Alleles existing longer tend to show weaker linkage disequilibrium with neighboring alleles, including the causal alleles, and are less likely to tag a SNP-disease association.     

Extent and distribution of linkage disequilibrium in three genomic regions

The positional cloning of genes underlying common complex diseases relies on the identification of linkage disequilibrium (LD) between genetic markers and disease. We have examined 127 polymorphisms in three genomic regions in a sample of 575 chromosomes from unrelated individuals of British ancestry. To establish phase, 800 individuals were genotyped in 160 families. The fine structure of LD was found to be highly irregular. Forty-five percent of the variation in disequilibrium measures could be explained by physical distance. Additional factors, such as allele frequency, type of polymorphism, and genomic location, explained <5% of the variation. Nevertheless, disequilibrium was occasionally detectable at 500 kb and was present for over one-half of marker pairs separated by <50 kb. Although these findings are encouraging for the prospects of a genomewide LD map, they suggest caution in interpreting localization due to allelic association.     

A Genome-Wide Association Study of Hypertension and Blood Pressure in African Americans

The evidence for the existence of genetic susceptibility variants for the common form of hypertension (“essential hypertension”) remains weak and inconsistent. We sought genetic variants underlying blood pressure (BP) by conducting a genome-wide association study (GWAS) among African Americans, a population group in the United States that is disproportionately affected by hypertension and associated complications, including stroke and kidney diseases. Using a dense panel of over 800,000 SNPs in a discovery sample of 1,017 African Americans from the Washington, D.C., metropolitan region, we identified multiple SNPs reaching genome-wide significance for systolic BP in or near the genes: PMS1, SLC24A4, YWHA7, IPO7, and CACANA1H. Two of these genes, SLC24A4 (a sodium/potassium/calcium exchanger) and CACNA1H (a voltage-dependent calcium channel), are potential candidate genes for BP regulation and the latter is a drug target for a class of calcium channel blockers. No variant reached genome wide significance for association with diastolic BP (top scoring SNP rs1867226, p = 5.8×10⁻⁷) or with hypertension as a binary trait (top scoring SNP rs9791170, p = 5.1×10⁻⁷). We replicated some of the significant SNPs in a sample of West Africans. Pathway analysis revealed that genes harboring top-scoring variants cluster in pathways and networks of biologic relevance to hypertension and BP regulation. This is the first GWAS for hypertension and BP in an African American population. The findings suggests that, in addition to or in lieu of relying solely on replicated variants of moderate-to-large effect reaching genome-wide significance, pathway and network approaches may be useful in identifying and prioritizing candidate genes/loci for further experiments.     

Genomewide search in Canadian families with inflammatory bowel disease reveals two novel susceptibility loci

The chronic inflammatory bowel diseases (IBDs)-Crohn disease (CD) and ulcerative colitis (UC)-are idiopathic, inflammatory disorders of the gastrointestinal tract. These conditions have a peak incidence in early adulthood and a combined prevalence of approximately 100-200/100,000. Although the etiology of IBD is multifactorial, a significant genetic contribution to disease susceptibility is implied by epidemiological data revealing a sibling risk of approximately 35-fold for CD and approximately 15-fold for UC. To elucidate the genetic basis for these disorders, we undertook a genomewide scan in 158 Canadian sib-pair families and identified three regions of suggestive linkage (3p, 5q31-33, and 6p) and one region of significant linkage to 19p13 (LOD score 4.6). Higher-density mapping in the 5q31-q33 region revealed a locus of genomewide significance (LOD score 3.9) that contributes to CD susceptibility in families with early-onset disease. Both of these genomic regions contain numerous genes that are important to the immune and inflammatory systems and that provide good targets for future candidate-gene studies.     

Allelic and phenotypic spectrum of plasma triglycerides

The genetic underpinnings of both normal and pathological variation in plasma triglyceride (TG) concentration are relatively well understood compared to many other complex metabolic traits. For instance, genome-wide association studies (GWAS) have revealed 32 common variants that are associated with plasma TG concentrations in healthy epidemiologic populations. Furthermore, GWAS in clinically ascertained hypertriglyceridemia (HTG) patients have shown that almost all of the same TG-raising alleles from epidemiologic samples are also associated with HTG disease status, and that greater accumulation of these alleles reflects the severity of the HTG phenotype. Finally, comprehensive resequencing studies show a burden of rare variants in some of these same genes - namely in LPL, GCKR, APOB and APOA5 - in HTG patients compared to normolipidemic controls. A more complete understanding of the genes and genetic variants associated with plasma TG concentration will enrich our understanding of the molecular pathways that modulate plasma TG metabolism, which may translate into clinical benefit. This article is part of a Special Issue entitled Triglyceride Metabolism and Disease.     

Combined analysis from eleven linkage studies of bipolar disorder provides strong evidence of susceptibility loci on chromosomes 6q and 8q

Several independent studies and meta-analyses aimed at identifying genomic regions linked to bipolar disorder (BP) have failed to find clear and consistent evidence of linkage regions. Our hypothesis is that combining the original genotype data provides benefits of increased power and control over sources of heterogeneity that outweigh the difficulty and potential pitfalls of the implementation. We conducted a combined analysis using the original genotype data from 11 BP genomewide linkage scans comprising 5,179 individuals from 1,067 families. Heterogeneity among studies was minimized in our analyses by using uniform methods of analysis and a common, standardized marker map and was assessed using novel methods developed for meta-analysis of genome scans. To date, this collaboration is the largest and most comprehensive analysis of linkage samples involving a psychiatric disorder. We demonstrate that combining original genome-scan data is a powerful approach for the elucidation of linkage regions underlying complex disease. Our results establish genomewide significant linkage to BP on chromosomes 6q and 8q, which provides solid information to guide future gene-finding efforts that rely on fine-mapping and association approaches.     

Large-scale analysis of non-synonymous coding region single nucleotide polymorphisms

MOTIVATION: Single nucleotide polymorphisms (SNPs) are the most common form of genetic variant in humans. SNPs causing amino acid substitutions are of particular interest as candidates for loci affecting susceptibility to complex diseases, such as diabetes and hypertension. To efficiently screen SNPs for disease association, it is important to distinguish neutral variants from deleterious ones. RESULTS: We describe the use of Pfam protein motif models and the HMMER program to predict whether amino acid changes in conserved domains are likely to affect protein function. We find that the magnitude of the change in the HMMER E-value caused by an amino acid substitution is a good predictor of whether it is deleterious. We provide internet-accessible display tools for a genomewide collection of SNPs, including 7391 distinct non-synonymous coding region SNPs in 2683 genes. AVAILABILITY: http://lpgws.nci.nih.gov/cgi-bin/GeneViewer.cgi     

Linkage analysis identifies a novel locus for restless legs syndrome on chromosome 2q in a South Tyrolean population isolate

Restless legs syndrome (RLS) is a common neurological condition with three loci (12q, 14q, and 9p) described so far, although none of these genes has yet been identified. We report a genomewide linkage scan of patients with RLS (n=37) assessed in a population isolate (n=530) of South Tyrol (Italy). Using both nonparametric and parametric analyses, we initially obtained suggestive evidence of a novel locus on chromosome 2q, with nominal evidence of linkage on chromosomes 5p and 17p. Follow-up genotyping yielded significant evidence of linkage (nonparametric LOD score 5.5, P相似文献   

A genome-wide linkage and association scan reveals novel loci for hypertension and blood pressure traits

Hypertension is caused by the interaction of environmental and genetic factors. The condition which is very common, with about 18% of the adult Hong Kong Chinese population and over 50% of older individuals affected, is responsible for considerable morbidity and mortality. To identify genes influencing hypertension and blood pressure, we conducted a combined linkage and association study using over 500,000 single nucleotide polymorphisms (SNPs) genotyped in 328 individuals comprising 111 hypertensive probands and their siblings. Using a family-based association test, we found an association with SNPs on chromosome 5q31.1 (rs6596140; P<9 × 10(-8)) for hypertension. One candidate gene, PDC, was replicated, with rs3817586 on 1q31.1 attaining P = 2.5 × 10(-4) and 2.9 × 10(-5) in the within-family tests for DBP and MAP, respectively. We also identified regions of significant linkage for systolic and diastolic blood pressure on chromosomes 2q22 and 5p13, respectively. Further family-based association analysis of the linkage peak on chromosome 5 yielded a significant association (rs1605685, P<7 × 10(-5)) for DBP. This is the first combined linkage and association study of hypertension and its related quantitative traits with Chinese ancestry. The associations reported here account for the action of common variants whereas the discovery of linkage regions may point to novel targets for rare variant screening.     

Imbalance production between interleukin-1beta (IL-1beta) and IL-1 receptor antagonist (IL-1Ra) in bronchial asthma

Genes of the IL-1 family encode three different peptides, IL-1alpha, IL-1beta, and IL-1Ra, respectively. IL-1 operates through IL-1RI, and is involved in airway inflammation in asthmatic subjects, whereas IL-1Ra appears to be a specific competitive inhibitor of IL-1. All genes are on chromosome 2q12-21 where genomewide searches have identified linkage for asthma. To test whether variants of IL-1 relate to asthma, we conducted a genetic association study in a Japanese population. We show that the A2 allele of IL1RN (encoding IL-1Ra) associates with nonatopic asthma [OR = 5.71, 95% CI: 1.63-19. 8, Pc = 0.007]. Both atopic and nonatopic asthmatics with the A2 allele had significantly lower serum IL-1Ra levels in both types of asthmatics. Peripheral blood cells from asthmatics with A2 alleles, however, produced as much IL-1 as did those with A1 homozygotes. Since Th1 and Th2 cytokines differentially regulate the ratio between IL-1beta and IL-1Ra, these findings suggest that dysregulation of IL-1beta/IL-1Ra, probably due to interaction between epithelium and immuno-competent cells in the airway, is important in asthma inflammation.     

Genomewide screen and identification of gene-gene interactions for asthma-susceptibility loci in three U.S. populations: collaborative study on the genetics of asthma

The genomewide screen to search for asthma-susceptibility loci, in the Collaborative Study on the Genetics of Asthma (CSGA), has been conducted in two stages and includes 266 families (199 nuclear and 67 extended pedigrees) from three U.S. populations: African American, European American, and Hispanic. Evidence for linkage with the asthma phenotype was observed for multiple chromosomal regions, through use of several analytical approaches that facilitated the identification of multiple disease loci. Ethnicity-specific analyses, which allowed for different frequencies of asthma-susceptibility genes in each ethnic population, provided the strongest evidence for linkage at 6p21 in the European American population, at 11q21 in the African American population, and at 1p32 in the Hispanic population. Both the conditional analysis and the affected-sib-pair two-locus analysis provided further evidence for linkage, at 5q31, 8p23, 12q22, and 15q13. Several of these regions have been observed in other genomewide screens and linkage or association studies, for asthma and related phenotypes. These results were used to develop a conceptual model to delineate asthma-susceptibility loci and their genetic interactions, which provides a promising basis for initiation of fine-mapping studies and, ultimately, for gene identification.     

Genomewide Linkage Analysis of Celiac Disease in Finnish Families

Celiac disease (CD), or gluten-sensitive enteropathy, is a common multifactorial disorder resulting from intolerance to cereal prolamins. The only established genetic susceptibility factor is HLA-DQ, which appears to explain only part of the overall genetic risk. We performed a genomewide scan of CD in 60 Finnish families. In addition to strong evidence for linkage to the HLA region at 6p21.3 (Z(max)>5), suggestive evidence for linkage was found for six other chromosomal regions--1p36, 4p15, 5q31, 7q21, 9p21-23, and 16q12. We further analyzed the three most convincing regions--4p15, 5q31, and 7q21--by evaluation of dense marker arrays across each region and by analysis of an additional 38 families. Although multipoint analysis with dense markers provided supportive evidence (multipoint LOD scores 3.25 at 4p15, 1.49 at 5q31, and 1.04 at 7q21) for the initial findings, the additional 38 families did not strengthen evidence for linkage. The role that HLA-DQ plays was studied in more detail by analysis of DQB1 alleles in all 98 families. All but one patient carried one or two HLA-DQ risk alleles, and 65% of HLA-DQ2 carriers were affected. Our study indicates that the HLA region harbors a predominant CD-susceptibility locus in these Finnish families.