Search

Congenital heart malformation (CHM) is the most common form of congenital human birth anomaly and is the leading cause of infant mortality. Although some causative genes have been identified, little progress has been made in identifying genes in which low-penetrance susceptibility variants occur in the majority of sporadic CHM cases. To identify common genetic variants associated with sporadic non-syndromic CHM in Han Chinese populations, we performed a multistage genome-wide association study (GWAS) in a total of 4,225 CHM cases and 5,112 non-CHM controls. The GWAS stage included 945 cases and 1,246 controls and was followed by 2-stage validation with 2,160 cases and 3,866 controls. The combined analyses identified significant associations (P < 5.0 × 10⁻⁸) at 1p12 (rs2474937 near TBX15; odds ratio (OR) = 1.40; P = 8.44 × 10⁻¹⁰) and 4q31.1 (rs1531070 in MAML3; OR = 1.40; P = 4.99 × 10⁻¹²). These results extend current knowledge of genetic contributions to CHM in Han Chinese populations.

To identify genetic risk factors underlying non-Hodgkin lymphomas (NHLs) from the B cell lineage, we conducted a genome-wide association study (GWAS) of 253 Chinese individuals with B cell NHL (cases) and 1,438 controls and further validation in 1,175 cases and 5,492 controls. We identified a new susceptibility locus, rs6773854, located between BCL6 (encoding B cell lymphoma protein 6) and LPP (encoding lipoma preferred partner) on oncogene-rich chromosome 3q27 that was significantly associated with increased risk of B cell NHL (meta-analysis P = 3.36 × 10⁻¹³, per-allele odds ratio (OR) = 1.44) and with diffuse large B cell lymphoma (DLBCL) in particular (meta-analysis P = 1.14 × 10⁻¹¹, OR = 1.47). We found no evidence of association of rs6773854 with non-B cell NHLs (T cell and natural killer (NK) lineages) (P = 0.17, OR = 1.12) and observed significant heterogeneity between B cell and non-B cell subtypes (Phet = 0.01, I² = 84%). Our results provide insight that germline variation in the intergenic region between BCL6 and LPP has a role in risk of B cell lymphomagenesis.

To extend understanding of the genetic architecture and molecular basis of type 2 diabetes (T2D), we conducted a meta-analysis of genetic variants on the Metabochip, including 34,840 cases and 114,981 controls, overwhelmingly of European descent. We identified ten previously unreported T2D susceptibility loci, including two showing sex-differentiated association. Genome-wide analyses of these data are consistent with a long tail of additional common variant loci explaining much of the variation in susceptibility to T2D. Exploration of the enlarged set of susceptibility loci implicates several processes, including CREBBP-related transcription, adipocytokine signaling and cell cycle regulation, in diabetes pathogenesis.

We sequenced and annotated the genomes of four P. vivax strains collected from disparate geographic locations, tripling the number of genome sequences available for this understudied parasite and providing the first genome-wide perspective of global variability in this species. We observe approximately twice as much SNP diversity among these isolates as we do among a comparable collection of isolates of P. falciparum, a malaria-causing parasite that results in higher mortality. This indicates a distinct history of global colonization and/or a more stable demographic history for P. vivax relative to P. falciparum, which is thought to have undergone a recent population bottleneck. The SNP diversity, as well as additional microsatellite and gene family variability, suggests a capacity for greater functional variation in the global population of P. vivax. These findings warrant a deeper survey of variation in P. vivax to equip disease interventions targeting the distinctive biology of this neglected but major pathogen.

Following a previous genome-wide association study (GWAS 1) including 744 cases and 895 controls, we analyzed genome-wide association data from a new cohort of Han Chinese (GWAS 2) with 1,510 polycystic ovary syndrome (PCOS) cases and 2,016 controls. We followed up significantly associated signals identified in the combined results of GWAS 1 and 2 in a total of 8,226 cases and 7,578 controls. In addition to confirming the three loci we previously reported, we identify eight new PCOS association signals at P < 5 × 10(-8): 9q22.32, 11q22.1, 12q13.2, 12q14.3, 16q12.1, 19p13.3, 20q13.2 and a second independent signal at 2p16.3 (the FSHR gene). These PCOS association signals show evidence of enrichment for candidate genes related to insulin signaling, sexual hormone function and type 2 diabetes (T2D). Other candidate genes were related to calcium signaling and endocytosis. Our findings provide new insight and direction for discovering the biological mechanisms of PCOS.

Strong signatures of positive selection at newly arising genetic variants are well documented in humans(1-8), but this form of selection may not be widespread in recent human evolution(9). Because many human traits are highly polygenic and partly determined by common, ancient genetic variation, an alternative model for rapid genetic adaptation has been proposed: weak selection acting on many pre-existing (standing) genetic variants, or polygenic adaptation(10-12). By studying height, a classic polygenic trait, we demonstrate the first human signature of widespread selection on standing variation. We show that frequencies of alleles associated with increased height, both at known loci and genome wide, are systematically elevated in Northern Europeans compared with Southern Europeans (P < 4.3 × 10(-4)). This pattern mirrors intra-European height differences and is not confounded by ancestry or other ascertainment biases. The systematic frequency differences are consistent with the presence of widespread weak selection (selection coefficients ∼10(-3)-10(-5) per allele) rather than genetic drift alone (P < 10(-15)).

Thyrotoxic periodic paralysis (TPP) is a potentially life-threatening complication of thyrotoxicosis. We conducted a genome-wide association study (GWAS) and a replication study with a total of 123 southern Chinese with TPP (cases) and 1,170 healthy controls and identified a susceptibility locus on chromosome 17q24.3 near KCNJ2 (rs312691: odds ratio (OR) = 3.3; P(meta-analysis) = 1.8 × 10(-14)). All subjects with TPP also had Graves' disease, and subsequent TPP versus Graves' disease comparison confirmed that the association at 17q24.3 was specific to TPP. The area under the curve (AUC) of rs312691 genotype for risk prediction of TPP in subjects with Graves' disease was 0.73. Expression quantitative trait locus (eQTL) analysis identified SNPs in the region flanking rs312691 (±10 kb) that could potentially affect KCNJ2 expression (P = 0.0001). Our study has identified a susceptibility locus associated with TPP and provides insight into the causes of TPP.

Genome-wide association studies (GWAS) are a standard approach for studying the genetics of natural variation. A major concern in GWAS is the need to account for the complicated dependence structure of the data, both between loci as well as between individuals. Mixed models have emerged as a general and flexible approach for correcting for population structure in GWAS. Here, we extend this linear mixed-model approach to carry out GWAS of correlated phenotypes, deriving a fully parameterized multi-trait mixed model (MTMM) that considers both the within-trait and between-trait variance components simultaneously for multiple traits. We apply this to data from a human cohort for correlated blood lipid traits from the Northern Finland Birth Cohort 1966 and show greatly increased power to detect pleiotropic loci that affect more than one blood lipid trait. We also apply this approach to an Arabidopsis thaliana data set for flowering measurements in two different locations, identifying loci whose effect depends on the environment.

In addition to its activity in nicotinamide adenine dinucleotide (NAD(+)) synthesis, the nuclear nicotinamide mononucleotide adenyltransferase NMNAT1 acts as a chaperone that protects against neuronal activity-induced degeneration. Here we report that compound heterozygous and homozygous NMNAT1 mutations cause severe neonatal neurodegeneration of the central retina and early-onset optic atrophy in 22 unrelated individuals. Their clinical presentation is consistent with Leber congenital amaurosis and suggests that the mutations affect neuroprotection of photoreceptor cells.

We have conducted the first meta-analyses for nonsyndromic cleft lip with or without cleft palate (NSCL/P) using data from the two largest genome-wide association studies published to date. We confirmed associations with all previously identified loci and identified six additional susceptibility regions (1p36, 2p21, 3p11.1, 8q21.3, 13q31.1 and 15q22). Analysis of phenotypic variability identified the first specific genetic risk factor for NSCLP (nonsyndromic cleft lip plus palate) (rs8001641; P(NSCLP) = 6.51 × 10(-11); homozygote relative risk = 2.41, 95% confidence interval (CI) 1.84-3.16).

We performed a genome-wide association study of esophageal squamous cell carcinoma (ESCC) by genotyping 1,077 individuals with ESCC and 1,733 control subjects of Chinese Han descent. We selected 18 promising SNPs for replication in an additional 7,673 cases of ESCC and 11,013 control subjects of Chinese Han descent and 303 cases of ESCC and 537 control subjects of Chinese Uygur-Kazakh descent. We identified two previously unknown susceptibility loci for ESCC: PLCE1 at 10q23 (P(Han combined for ESCC) = 7.46 x 10(-56), odds ratio (OR) = 1.43; P(Uygur-Kazakh for ESCC) = 5.70 x 10(-4), OR = 1.53) and C20orf54 at 20p13 (P(Han combined for ESCC) = 1.21 x 10(-11), OR = 0.86; P(Uygur-Kazakh for ESCC) = 7.88 x 10(-3), OR = 0.66). We also confirmed association in 2,766 cases of gastric cardia adenocarcinoma cases and the same 11,013 control subjects (PLCE1, P(Han for GCA) = 1.74 x 10(-39), OR = 1.55 and C20orf54, P(Han for GCA) = 3.02 x 10(-3), OR = 0.91). PLCE1 and C20orf54 have important biological implications for both ESCC and GCA. PLCE1 might regulate cell growth, differentiation, apoptosis and angiogenesis. C20orf54 is responsible for transporting riboflavin, and deficiency of riboflavin has been documented as a risk factor for ESCC and GCA.

Atrial fibrillation (AF) is the most common sustained arrhythmia. Previous studies have identified several genetic loci associated with typical AF. We sought to identify common genetic variants underlying lone AF. This condition affects a subset of individuals without overt heart disease and with an increased heritability of AF. We report a meta-analysis of genome-wide association studies conducted using 1,335 individuals with lone AF (cases) and 12,844 unaffected individuals (referents). Cases were obtained from the German AF Network, Heart and Vascular Health Study, the Atherosclerosis Risk in Communities Study, the Cleveland Clinic and Massachusetts General Hospital. We identified an association on chromosome 1q21 to lone AF (rs13376333, adjusted odds ratio = 1.56; P = 6.3 x 10(-12)), and we replicated this association in two independent cohorts with lone AF (overall combined odds ratio = 1.52, 95% CI 1.40-1.64; P = 1.83 x 10(-21)). rs13376333 is intronic to KCNN3, which encodes a potassium channel protein involved in atrial repolarization.

Pheochromocytomas, which are catecholamine-secreting tumors of neural crest origin, are frequently hereditary. However, the molecular basis of the majority of these tumors is unknown. We identified the transmembrane-encoding gene TMEM127 on chromosome 2q11 as a new pheochromocytoma susceptibility gene. In a cohort of 103 samples, we detected truncating germline TMEM127 mutations in approximately 30% of familial tumors and about 3% of sporadic-appearing pheochromocytomas without a known genetic cause. The wild-type allele was consistently deleted in tumor DNA, suggesting a classic mechanism of tumor suppressor gene inactivation. Pheochromocytomas with mutations in TMEM127 are transcriptionally related to tumors bearing NF1 mutations and, similarly, show hyperphosphorylation of mammalian target of rapamycin (mTOR) effector proteins. Accordingly, in vitro gain-of-function and loss-of-function analyses indicate that TMEM127 is a negative regulator of mTOR. TMEM127 dynamically associates with the endomembrane system and colocalizes with perinuclear (activated) mTOR, suggesting a subcompartmental-specific effect. Our studies identify TMEM127 as a tumor suppressor gene and validate the power of hereditary tumors to elucidate cancer pathogenesis.

We conducted genome-wide association analyses of mean leukocyte telomere length in 2,917 individuals, with follow-up replication in 9,492 individuals. We identified an association with telomere length on 3q26 (rs12696304, combined P = 3.72 x 10(-14)) at a locus that includes TERC, which encodes the telomerase RNA component. Each copy of the minor allele of rs12696304 was associated with an approximately 75-base-pair reduction in mean telomere length, equivalent to approximately 3.6 years of age-related telomere-length attrition.

Eosinophils are pleiotropic multifunctional leukocytes involved in initiation and propagation of inflammatory responses and thus have important roles in the pathogenesis of inflammatory diseases. Here we describe a genome-wide association scan for sequence variants affecting eosinophil counts in blood of 9,392 Icelanders. The most significant SNPs were studied further in 12,118 Europeans and 5,212 East Asians. SNPs at 2q12 (rs1420101), 2q13 (rs12619285), 3q21 (rs4857855), 5q31 (rs4143832) and 12q24 (rs3184504) reached genome-wide significance (P = 5.3 x 10(-14), 5.4 x 10(-10), 8.6 x 10(-17), 1.2 x 10(-10) and 6.5 x 10(-19), respectively). A SNP at IL1RL1 associated with asthma (P = 5.5 x 10(-12)) in a collection of ten different populations (7,996 cases and 44,890 controls). SNPs at WDR36, IL33 and MYB that showed suggestive association with eosinophil counts were also associated with atopic asthma (P = 4.2 x 10(-6), 2.2 x 10(-5) and 2.4 x 10(-4), respectively). We also found that a nonsynonymous SNP at 12q24, in SH2B3, associated significantly (P = 8.6 x 10(-8)) with myocardial infarction in six different populations (6,650 cases and 40,621 controls).

Myocardial infarction is a common disease and among the leading causes of death in the world. We previously reported association of variants in LGALS2, encoding galectin-2, with myocardial infarction susceptibility in a case-control association study in a Japanese population. Here we identify BRAP (BRCA1-associated protein) as a galectin-2-binding protein. We report an association of SNPs in BRAP with myocardial infarction risk in a large Japanese cohort (P = 3.0 x 10(-18), OR = 1.48, 2,475 cases and 2,778 controls), with replication in additional Japanese and Taiwanese cohorts (P = 4.4 x 10(-6), 862 cases and 1,113 controls and P = 4.7 x 10(-3), 349 cases and 994 controls, respectively). BRAP expression was observed in smooth muscle cells (SMCs) and macrophages in human atherosclerotic lesions. BRAP knockdown by siRNA using cultured coronary endothelial cells suppressed activation of NF-kappaB, a central mediator of inflammation.

We identified a marker in LINGO1 showing genome-wide significant association (P = 1.2 x 10(-9), odds ratio = 1.55) with essential tremor. LINGO1 has potent, negative regulatory influences on neuronal survival and is also important in regulating both central-nervous-system axon regeneration and oligodendrocyte maturation. Increased axon integrity observed in Lingo1 mouse [corrected] knockout models highlights the potential role of LINGO1 in the pathophysiology of ET [corrected]

The common sequence variants that have recently been associated with cancer risk are particular to a single cancer type or at most two. Following up on our genome-wide scan of basal cell carcinoma, we found that rs401681[C] on chromosome 5p15.33 satisfied our threshold for genome-wide significance (OR = 1.25, P = 3.7 x 10(-12)). We tested rs401681 for association with 16 additional cancer types in over 30,000 cancer cases and 45,000 controls and found association with lung cancer (OR = 1.15, P = 7.2 x 10(-8)) and urinary bladder, prostate and cervix cancer (ORs = 1.07-1.31, all P < 4 x 10(-4)). However, rs401681[C] seems to confer protection against cutaneous melanoma (OR = 0.88, P = 8.0 x 10(-4)). Notably, most of these cancer types have a strong environmental component to their risk. Investigation of the region led us to rs2736098[A], which showed stronger association with some cancer types. However, neither variant could fully account for the association of the other. rs2736098 corresponds to A305A in the telomerase reverse transcriptase (TERT) protein and rs401681 is in an intron of the CLPTM1L gene.

Individual response to small-molecule drugs is variable; a drug that provides a cure for some may confer no therapeutic benefit or trigger an adverse reaction in others. To begin to understand such differences systematically, we treated 104 genotyped segregants from a cross between two yeast strains with a collection of 100 diverse small molecules. We used linkage analysis to identify 124 distinct linkages between genetic markers and response to 83 compounds. The linked markers clustered at eight genomic locations, or quantitative-trait locus 'hotspots', that contain one or more polymorphisms that affect response to multiple small molecules. We also experimentally verified that a deficiency in leucine biosynthesis caused by a deletion of LEU2 underlies sensitivity to niguldipine, which is structurally related to therapeutic calcium channel blockers, and that a natural coding-region polymorphism in the inorganic phosphate transporter PHO84 underlies sensitivity to two polychlorinated phenols that uncouple oxidative phosphorylation. Our results provide a step toward a systematic understanding of small-molecule drug action in genetically distinct individuals.

Resequencing genes provides the opportunity to assess the full spectrum of variants that influence complex traits. Here we report the first application of resequencing to a large population (n = 3,551) to examine the role of the adipokine ANGPTL4 in lipid metabolism. Nonsynonymous variants in ANGPTL4 were more prevalent in individuals with triglyceride levels in the lowest quartile than in individuals with levels in the highest quartile (P = 0.016). One variant (E40K), present in approximately 3% of European Americans, was associated with significantly lower plasma levels of triglyceride and higher levels of high-density lipoprotein cholesterol in European Americans from the Atherosclerosis Risk in Communities Study and in Danes from the Copenhagen City Heart Study. The ratio of nonsynonymous to synonymous variants was higher in European Americans than in African Americans (4:1 versus 1.3:1), suggesting population-specific relaxation of purifying selection. Thus, resequencing of ANGPTL4 in a multiethnic population allowed analysis of the phenotypic effects of both rare and common variants while taking advantage of genetic variation arising from ethnic differences in population history.

Toll-like receptors (TLRs) and members of their signaling pathway are important in the initiation of the innate immune response to a wide variety of pathogens. The adaptor protein Mal (also known as TIRAP), encoded by TIRAP (MIM 606252), mediates downstream signaling of TLR2 and TLR4 (refs. 4-6). We report a case-control study of 6,106 individuals from the UK, Vietnam and several African countries with invasive pneumococcal disease, bacteremia, malaria and tuberculosis. We genotyped 33 SNPs, including rs8177374, which encodes a leucine substitution at Ser180 of Mal. We found that heterozygous carriage of this variant associated independently with all four infectious diseases in the different study populations. Combining the study groups, we found substantial support for a protective effect of S180L heterozygosity against these infectious diseases (N = 6,106; overall P = 9.6 x 10(-8)). We found that the Mal S180L variant attenuated TLR2 signal transduction.

Osteoarthritis (MIM 165720), characterized by degeneration of articular cartilage, is the most common form of human arthritis and a major concern for aging societies worldwide. Epidemiological and genetic studies have shown that osteoarthritis is a polygenic disease. Here, we report that the gene encoding growth differentiation factor 5 (GDF5) is associated with osteoarthritis in Asian populations. A SNP in the 5' UTR of GDF5 (+104T/C; rs143383) showed significant association (P = 1.8 x 10(-13)) with hip osteoarthritis in two independent Japanese populations. This association was replicated for knee osteoarthritis in Japanese (P = 0.0021) and Han Chinese (P = 0.00028) populations. This SNP, located in the GDF5 core promoter, exerts allelic differences on transcriptional activity in chondrogenic cells, with the susceptibility allele showing reduced activity. Our findings implicate GDF5 as a susceptibility gene for osteoarthritis and suggest that decreased GDF5 expression is involved in the pathogenesis of osteoarthritis.