Search

Rearrangements involving the DUX4 gene (DUX4-r) define a subtype of paediatric and adult acute lymphoblastic leukaemia (ALL) with a favourable outcome. Currently, there is no 'standard of care' diagnostic method for their confident identification. Here, we present an open-source software tool designed to detect DUX4-r from short-read, whole-genome sequencing (WGS) data. Evaluation on a cohort of 210 paediatric ALL cases showed that our method detects all known, as well as previously unidentified, cases of IGH::DUX4 and rearrangements with other partner genes. These findings demonstrate the possibility of robustly detecting DUX4-r using WGS in the routine clinical setting., Competing Interests: Declarations. Ethical approval and consent to participate: Informed consent was obtained from all the participants and/or legal guardians for the study. Patients donated samples to the VIVO biobank. The study was approved by the National Research Ethics Service Committee South West - Central Bristol (11/SW/0153). The study adhered to the Declaration of Helsinki. Consent for publication: Not applicable. Competing interests: PG, SB, KJC, CF, IA, HN, DJM, PJC, JB, DRB, MM, and MTR are employees of Illumina, a public company that develops and markets systems for genetic analysis., (© 2024. The Author(s).)

Childhood B-cell acute lymphoblastic leukaemia (B-ALL) is characterised by recurrent genetic abnormalities that drive risk-directed treatment strategies. Using current techniques, accurate detection of such aberrations can be challenging, due to the rapidly expanding list of key genetic abnormalities. Whole genome sequencing (WGS) has the potential to improve genetic testing, but requires comprehensive validation. We performed WGS on 210 childhood B-ALL samples annotated with clinical and genetic data. We devised a molecular classification system to subtype these patients based on identification of key genetic changes in tumour-normal and tumour-only analyses. This approach detected 294 subtype-defining genetic abnormalities in 96% (202/210) patients. Novel genetic variants, including fusions involving genes in the MAP kinase pathway, were identified. WGS results were concordant with standard-of-care methods and whole transcriptome sequencing (WTS). We expanded the catalogue of genetic profiles that reliably classify PAX5alt and ETV6::RUNX1-like subtypes. Our novel bioinformatic pipeline improved detection of DUX4 rearrangements (DUX4-r): a good-risk B-ALL subtype with high survival rates. Overall, we have validated that WGS provides a standalone, reliable genetic test to detect all subtype-defining genetic abnormalities in B-ALL, accurately classifying patients for the risk-directed treatment stratification, while simultaneously performing as a research tool to identify novel disease biomarkers., (© 2023. The Author(s).)

Incorporating genetics into risk-stratification for treatment of childhood B-progenitor acute lymphoblastic leukaemia (B-ALL) has contributed significantly to improved survival. In about 30% B-ALL (B-other-ALL) without well-established chromosomal changes, new genetic subtypes have recently emerged, yet their true prognostic relevance largely remains unclear. We integrated next generation sequencing (NGS): whole genome sequencing (WGS) (n = 157) and bespoke targeted NGS (t-NGS) (n = 175) (overlap n = 36), with existing genetic annotation in a representative cohort of 351 B-other-ALL patients from the childhood ALL trail, UKALL2003. PAX5alt was most frequently observed (n = 91), whereas PAX5 P80R mutations (n = 11) defined a distinct PAX5 subtype. DUX4-r subtype (n = 80) was defined by DUX4 rearrangements and/or ERG deletions. These patients had a low relapse rate and excellent survival. ETV6::RUNX1-like subtype (n = 21) was characterised by multiple abnormalities of ETV6 and IKZF1, with no reported relapses or deaths, indicating their excellent prognosis in this trial. An inferior outcome for patients with ABL-class fusions (n = 25) was confirmed. Integration of NGS into genomic profiling of B-other-ALL within a single childhood ALL trial, UKALL2003, has shown the added clinical value of NGS-based approaches, through improved accuracy in detection and classification into the range of risk stratifying genetic subtypes, while validating their prognostic significance., (© 2022. The Author(s).)

The genomic landscape of breast cancer is complex, and inter- and intra-tumour heterogeneity are important challenges in treating the disease. In this study, we sequence 173 genes in 2,433 primary breast tumours that have copy number aberration (CNA), gene expression and long-term clinical follow-up data. We identify 40 mutation-driver (Mut-driver) genes, and determine associations between mutations, driver CNA profiles, clinical-pathological parameters and survival. We assess the clonal states of Mut-driver mutations, and estimate levels of intra-tumour heterogeneity using mutant-allele fractions. Associations between PIK3CA mutations and reduced survival are identified in three subgroups of ER-positive cancer (defined by amplification of 17q23, 11q13-14 or 8q24). High levels of intra-tumour heterogeneity are in general associated with a worse outcome, but highly aggressive tumours with 11q13-14 amplification have low levels of intra-tumour heterogeneity. These results emphasize the importance of genome-based stratification of breast cancer, and have important implications for designing therapeutic strategies.

Background: Glioblastoma (GBM) is the most common malignant brain cancer occurring in adults, and is associated with dismal outcome and few therapeutic options. GBM has been shown to predominantly disrupt three core pathways through somatic aberrations, rendering it ideal for precision medicine approaches., Methods: We describe a 35-year-old female patient with recurrent GBM following surgical removal of the primary tumour, adjuvant treatment with temozolomide and a 3-year disease-free period. Rapid whole-genome sequencing (WGS) of three separate tumour regions at recurrence was carried out and interpreted relative to WGS of two regions of the primary tumour., Results: We found extensive mutational and copy-number heterogeneity within the primary tumour. We identified a TP53 mutation and two focal amplifications involving PDGFRA, KIT and CDK4, on chromosomes 4 and 12. A clonal IDH1 R132H mutation in the primary, a known GBM driver event, was detectable at only very low frequency in the recurrent tumour. After sub-clonal diversification, evidence was found for a whole-genome doubling event and a translocation between the amplified regions of PDGFRA, KIT and CDK4, encoded within a double-minute chromosome also incorporating miR26a-2. The WGS analysis uncovered progressive evolution of the double-minute chromosome converging on the KIT/PDGFRA/PI3K/mTOR axis, superseding the IDH1 mutation in dominance in a mutually exclusive manner at recurrence, consequently the patient was treated with imatinib. Despite rapid sequencing and cancer genome-guided therapy against amplified oncogenes, the disease progressed, and the patient died shortly after., Conclusion: This case sheds light on the dynamic evolution of a GBM tumour, defining the origins of the lethal sub-clone, the macro-evolutionary genomic events dominating the disease at recurrence and the loss of a clonal driver. Even in the era of rapid WGS analysis, cases such as this illustrate the significant hurdles for precision medicine success., (© The Author 2015. Published by Oxford University Press on behalf of the European Society for Medical Oncology.)

Coronary artery disease (CAD) is the commonest cause of death. Here, we report an association analysis in 63,746 CAD cases and 130,681 controls identifying 15 loci reaching genome-wide significance, taking the number of susceptibility loci for CAD to 46, and a further 104 independent variants (r(2) < 0.2) strongly associated with CAD at a 5% false discovery rate (FDR). Together, these variants explain approximately 10.6% of CAD heritability. Of the 46 genome-wide significant lead SNPs, 12 show a significant association with a lipid trait, and 5 show a significant association with blood pressure, but none is significantly associated with diabetes. Network analysis with 233 candidate genes (loci at 10% FDR) generated 5 interaction networks comprising 85% of these putative genes involved in CAD. The four most significant pathways mapping to these networks are linked to lipid metabolism and inflammation, underscoring the causal role of these activities in the genetic etiology of CAD. Our study provides insights into the genetic basis of CAD and identifies key biological pathways.

We conducted a genome-wide association study to identify novel associations between genetic variants and circulating plasminogen activator inhibitor-1 (PAI-1) concentration, and examined functional implications of variants and genes that were discovered. A discovery meta-analysis was performed in 19 599 subjects, followed by replication analysis of genome-wide significant (P < 5 × 10(-8)) single nucleotide polymorphisms (SNPs) in 10 796 independent samples. We further examined associations with type 2 diabetes and coronary artery disease, assessed the functional significance of the SNPs for gene expression in human tissues, and conducted RNA-silencing experiments for one novel association. We confirmed the association of the 4G/5G proxy SNP rs2227631 in the promoter region of SERPINE1 (7q22.1) and discovered genome-wide significant associations at 3 additional loci: chromosome 7q22.1 close to SERPINE1 (rs6976053, discovery P = 3.4 × 10(-10)); chromosome 11p15.2 within ARNTL (rs6486122, discovery P = 3.0 × 10(-8)); and chromosome 3p25.2 within PPARG (rs11128603, discovery P = 2.9 × 10(-8)). Replication was achieved for the 7q22.1 and 11p15.2 loci. There was nominal association with type 2 diabetes and coronary artery disease at ARNTL (P < .05). Functional studies identified MUC3 as a candidate gene for the second association signal on 7q22.1. In summary, SNPs in SERPINE1 and ARNTL and an SNP associated with the expression of MUC3 were robustly associated with circulating levels of PAI-1.

Through genome-wide association meta-analyses of up to 133,010 individuals of European ancestry without diabetes, including individuals newly genotyped using the Metabochip, we have increased the number of confirmed loci influencing glycemic traits to 53, of which 33 also increase type 2 diabetes risk (q < 0.05). Loci influencing fasting insulin concentration showed association with lipid levels and fat distribution, suggesting impact on insulin resistance. Gene-based analyses identified further biologically plausible loci, suggesting that additional loci beyond those reaching genome-wide significance are likely to represent real associations. This conclusion is supported by an excess of directionally consistent and nominally significant signals between discovery and follow-up studies. Functional analysis of these newly discovered loci will further improve our understanding of glycemic control.

Objectives: The purpose of this study is investigate the effects of variants in the apolipoprotein(a) gene (LPA) on vascular diseases with different atherosclerotic and thrombotic components., Background: It is unclear whether the LPA variants rs10455872 and rs3798220, which correlate with lipoprotein(a) levels and coronary artery disease (CAD), confer susceptibility predominantly via atherosclerosis or thrombosis., Methods: The 2 LPA variants were combined and examined as LPA scores for the association with ischemic stroke (and TOAST [Trial of Org 10172 in Acute Stroke Treatment] subtypes) (effective sample size [n(e)] = 9,396); peripheral arterial disease (n(e) = 5,215); abdominal aortic aneurysm (n(e) = 4,572); venous thromboembolism (n(e) = 4,607); intracranial aneurysm (n(e) = 1,328); CAD (n(e) = 12,716), carotid intima-media thickness (n = 3,714), and angiographic CAD severity (n = 5,588)., Results: LPA score was associated with ischemic stroke subtype large artery atherosclerosis (odds ratio [OR]: 1.27; p = 6.7 × 10(-4)), peripheral artery disease (OR: 1.47; p = 2.9 × 10(-14)), and abdominal aortic aneurysm (OR: 1.23; p = 6.0 × 10(-5)), but not with the ischemic stroke subtypes cardioembolism (OR: 1.03; p = 0.69) or small vessel disease (OR: 1.06; p = 0.52). Although the LPA variants were not associated with carotid intima-media thickness, they were associated with the number of obstructed coronary vessels (p = 4.8 × 10(-12)). Furthermore, CAD cases carrying LPA risk variants had increased susceptibility to atherosclerotic manifestations outside of the coronary tree (OR: 1.26; p = 0.0010) and had earlier onset of CAD (-1.58 years/allele; p = 8.2 × 10(-8)) than CAD cases not carrying the risk variants. There was no association of LPA score with venous thromboembolism (OR: 0.97; p = 0.63) or intracranial aneurysm (OR: 0.85; p = 0.15)., Conclusions: LPA sequence variants were associated with atherosclerotic burden, but not with primarily thrombotic phenotypes., (Copyright © 2012 American College of Cardiology Foundation. Published by Elsevier Inc. All rights reserved.)

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and β-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies., (Copyright © 2011 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Coronary artery disease (CAD) is the leading cause of death worldwide. Affected individuals cluster in families in patterns that reflect the sharing of numerous susceptibility genes. Genome-wide and large-scale gene-centric genotyping studies that involve tens of thousands of cases and controls have now mapped common disease variants to 34 distinct loci. Some coronary disease common variants show allelic heterogeneity or copy number variation. Some of the loci include candidate genes that imply conventional or emerging risk factor-mediated mechanisms of disease pathogenesis. Quantitative trait loci associations with risk factors have been informative in Mendelian randomization studies as well as fine-mapping of causative variants. But, for most loci, plausible mechanistic links are uncertain or obscure at present but provide potentially novel directions for research into this disease's pathogenesis. The common variants explain ~4% of inter-individual variation in disease risk and no more than 13% of the total heritability of coronary disease. Although many CAD genes are presently undiscovered, it is likely that larger collaborative genome-wide association studies will map further common/low-penetrance variants and hoped that low-frequency or rare high-penetrance variants will also be identified in medical resequencing experiments.

Numerous genetic loci have been associated with systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N = 74,064) and follow-up studies (N = 48,607), we identified at genome-wide significance (P = 2.7 × 10(-8) to P = 2.3 × 10(-13)) four new PP loci (at 4q12 near CHIC2, 7q22.3 near PIK3CG, 8q24.12 in NOV and 11q24.3 near ADAMTS8), two new MAP loci (3p21.31 in MAP4 and 10q25.3 near ADRB1) and one locus associated with both of these traits (2q24.3 near FIGN) that has also recently been associated with SBP in east Asians. For three of the new PP loci, the estimated effect for SBP was opposite of that for DBP, in contrast to the majority of common SBP- and DBP-associated variants, which show concordant effects on both traits. These findings suggest new genetic pathways underlying blood pressure variation, some of which may differentially influence SBP and DBP.

Blood pressure is a heritable trait influenced by several biological pathways and responsive to environmental stimuli. Over one billion people worldwide have hypertension (≥140 mm Hg systolic blood pressure or  ≥90 mm Hg diastolic blood pressure). Even small increments in blood pressure are associated with an increased risk of cardiovascular events. This genome-wide association study of systolic and diastolic blood pressure, which used a multi-stage design in 200,000 individuals of European descent, identified sixteen novel loci: six of these loci contain genes previously known or suspected to regulate blood pressure (GUCY1A3-GUCY1B3, NPR3-C5orf23, ADM, FURIN-FES, GOSR2, GNAS-EDN3); the other ten provide new clues to blood pressure physiology. A genetic risk score based on 29 genome-wide significant variants was associated with hypertension, left ventricular wall thickness, stroke and coronary artery disease, but not kidney disease or kidney function. We also observed associations with blood pressure in East Asian, South Asian and African ancestry individuals. Our findings provide new insights into the genetics and biology of blood pressure, and suggest potential novel therapeutic pathways for cardiovascular disease prevention.

Height is a classic complex trait with common variants in a growing list of genes known to contribute to the phenotype. Using a genecentric genotyping array targeted toward cardiovascular-related loci, comprising 49,320 SNPs across approximately 2000 loci, we evaluated the association of common and uncommon SNPs with adult height in 114,223 individuals from 47 studies and six ethnicities. A total of 64 loci contained a SNP associated with height at array-wide significance (p < 2.4 × 10(-6)), with 42 loci surpassing the conventional genome-wide significance threshold (p < 5 × 10(-8)). Common variants with minor allele frequencies greater than 5% were observed to be associated with height in 37 previously reported loci. In individuals of European ancestry, uncommon SNPs in IL11 and SMAD3, which would not be genotyped with the use of standard genome-wide genotyping arrays, were strongly associated with height (p < 3 × 10(-11)). Conditional analysis within associated regions revealed five additional variants associated with height independent of lead SNPs within the locus, suggesting allelic heterogeneity. Although underpowered to replicate findings from individuals of European ancestry, the direction of effect of associated variants was largely consistent in African American, South Asian, and Hispanic populations. Overall, we show that dense coverage of genes for uncommon SNPs, coupled with large-scale meta-analysis, can successfully identify additional variants associated with a common complex trait.

Objective: Glycated hemoglobin (HbA₁(c)), used to monitor and diagnose diabetes, is influenced by average glycemia over a 2- to 3-month period. Genetic factors affecting expression, turnover, and abnormal glycation of hemoglobin could also be associated with increased levels of HbA₁(c). We aimed to identify such genetic factors and investigate the extent to which they influence diabetes classification based on HbA₁(c) levels., Research Design and Methods: We studied associations with HbA₁(c) in up to 46,368 nondiabetic adults of European descent from 23 genome-wide association studies (GWAS) and 8 cohorts with de novo genotyped single nucleotide polymorphisms (SNPs). We combined studies using inverse-variance meta-analysis and tested mediation by glycemia using conditional analyses. We estimated the global effect of HbA₁(c) loci using a multilocus risk score, and used net reclassification to estimate genetic effects on diabetes screening., Results: Ten loci reached genome-wide significant association with HbA(1c), including six new loci near FN3K (lead SNP/P value, rs1046896/P = 1.6 × 10⁻²⁶), HFE (rs1800562/P = 2.6 × 10⁻²⁰), TMPRSS6 (rs855791/P = 2.7 × 10⁻¹⁴), ANK1 (rs4737009/P = 6.1 × 10⁻¹²), SPTA1 (rs2779116/P = 2.8 × 10⁻⁹) and ATP11A/TUBGCP3 (rs7998202/P = 5.2 × 10⁻⁹), and four known HbA₁(c) loci: HK1 (rs16926246/P = 3.1 × 10⁻⁵⁴), MTNR1B (rs1387153/P = 4.0 × 10⁻¹¹), GCK (rs1799884/P = 1.5 × 10⁻²⁰) and G6PC2/ABCB11 (rs552976/P = 8.2 × 10⁻¹⁸). We show that associations with HbA₁(c) are partly a function of hyperglycemia associated with 3 of the 10 loci (GCK, G6PC2 and MTNR1B). The seven nonglycemic loci accounted for a 0.19 (% HbA₁(c)) difference between the extreme 10% tails of the risk score, and would reclassify ∼2% of a general white population screened for diabetes with HbA₁(c)., Conclusions: GWAS identified 10 genetic loci reproducibly associated with HbA₁(c). Six are novel and seven map to loci where rarer variants cause hereditary anemias and iron storage disorders. Common variants at these loci likely influence HbA₁(c) levels via erythrocyte biology, and confer a small but detectable reclassification of diabetes diagnosis by HbA₁(c).

Waist-hip ratio (WHR) is a measure of body fat distribution and a predictor of metabolic consequences independent of overall adiposity. WHR is heritable, but few genetic variants influencing this trait have been identified. We conducted a meta-analysis of 32 genome-wide association studies for WHR adjusted for body mass index (comprising up to 77,167 participants), following up 16 loci in an additional 29 studies (comprising up to 113,636 subjects). We identified 13 new loci in or near RSPO3, VEGFA, TBX15-WARS2, NFE2L3, GRB14, DNM3-PIGC, ITPR2-SSPN, LY86, HOXC13, ADAMTS9, ZNRF3-KREMEN1, NISCH-STAB1 and CPEB4 (P = 1.9 × 10⁻⁹ to P = 1.8 × 10⁻⁴⁰) and the known signal at LYPLAL1. Seven of these loci exhibited marked sexual dimorphism, all with a stronger effect on WHR in women than men (P for sex difference = 1.9 × 10⁻³ to P = 1.2 × 10⁻¹³). These findings provide evidence for multiple loci that modulate body fat distribution independent of overall adiposity and reveal strong gene-by-sex interactions.

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Demographic and family studies support the existence of a genetic contribution to the pathogenesis of IgA nephropathy, but results from genetic association studies of candidate genes are inconsistent. To systematically survey common genetic variation in this disease, we performed a genome-wide analysis in a cohort of patients with IgA nephropathy selected from the UK Glomerulonephritis DNA Bank. We used two groups of controls: parents of affected individuals and previously genotyped, unaffected, ancestry-matched individuals from the 1958 British Birth Cohort and the UK Blood Service. We genotyped 914 affected or family controls for 318,127 single nucleotide polymorphisms (SNPs). Filtering for low genotype call rates and inferred non-European ancestry left 533 genotyped individuals (187 affected children) for the family-based association analysis and 244 cases and 4980 controls for the case-control analysis. A total of 286,200 SNPs with call rates >95% were available for analysis. Genome-wide analysis showed a strong signal of association on chromosome 6p in the region of the MHC (P = 1 × 10(-9)). The two most strongly associated SNPs showed consistent association in both family-based and case-control analyses. HLA imputation analysis showed that the strongest association signal arose from a combination of DQ loci with some support for an independent HLA-B signal. These results suggest that the HLA region contains the strongest common susceptibility alleles that predispose to IgA nephropathy in the European population.

Background: Plasma levels of coagulation factors VII (FVII), VIII (FVIII), and von Willebrand factor (vWF) influence risk of hemorrhage and thrombosis. We conducted genome-wide association studies to identify new loci associated with plasma levels., Methods and Results: The setting of the study included 5 community-based studies for discovery comprising 23 608 European-ancestry participants: Atherosclerosis Risk In Communities Study, Cardiovascular Health Study, British 1958 Birth Cohort, Framingham Heart Study, and Rotterdam Study. All subjects had genome-wide single-nucleotide polymorphism (SNP) scans and at least 1 phenotype measured: FVII activity/antigen, FVIII activity, and vWF antigen. Each study used its genotype data to impute to HapMap SNPs and independently conducted association analyses of hemostasis measures using an additive genetic model. Study findings were combined by meta-analysis. Replication was conducted in 7604 participants not in the discovery cohort. For FVII, 305 SNPs exceeded the genome-wide significance threshold of 5.0x10(-8) and comprised 5 loci on 5 chromosomes: 2p23 (smallest P value 6.2x10(-24)), 4q25 (3.6x10(-12)), 11q12 (2.0x10(-10)), 13q34 (9.0x10(-259)), and 20q11.2 (5.7x10(-37)). Loci were within or near genes, including 4 new candidate genes and F7 (13q34). For vWF, 400 SNPs exceeded the threshold and marked 8 loci on 6 chromosomes: 6q24 (1.2x10(-22)), 8p21 (1.3x10(-16)), 9q34 (<5.0x10(-324)), 12p13 (1.7x10(-32)), 12q23 (7.3x10(-10)), 12q24.3 (3.8x10(-11)), 14q32 (2.3x10(-10)), and 19p13.2 (1.3x10(-9)). All loci were within genes, including 6 new candidate genes, as well as ABO (9q34) and VWF (12p13). For FVIII, 5 loci were identified and overlapped vWF findings. Nine of the 10 new findings were replicated., Conclusions: New genetic associations were discovered outside previously known biological pathways and may point to novel prevention and treatment targets of hemostasis disorders.

Context: Plasma levels of C-reactive protein (CRP) are independently associated with risk of coronary heart disease, but whether CRP is causally associated with coronary heart disease or merely a marker of underlying atherosclerosis is uncertain., Objective: To investigate association of genetic loci with CRP levels and risk of coronary heart disease., Design, Setting, and Participants: We first carried out a genome-wide association (n = 17,967) and replication study (n = 13,615) to identify genetic loci associated with plasma CRP concentrations. Data collection took place between 1989 and 2008 and genotyping between 2003 and 2008. We carried out a mendelian randomization study of the most closely associated single-nucleotide polymorphism (SNP) in the CRP locus and published data on other CRP variants involving a total of 28,112 cases and 100,823 controls, to investigate the association of CRP variants with coronary heart disease. We compared our finding with that predicted from meta-analysis of observational studies of CRP levels and risk of coronary heart disease. For the other loci associated with CRP levels, we selected the most closely associated SNP for testing against coronary heart disease among 14,365 cases and 32,069 controls., Main Outcome Measure: Risk of coronary heart disease., Results: Polymorphisms in 5 genetic loci were strongly associated with CRP levels (% difference per minor allele): SNP rs6700896 in LEPR (-14.8%; 95% confidence interval [CI], -17.6% to -12.0%; P = 6.2 x 10(-22)), rs4537545 in IL6R (-11.5%; 95% CI, -14.4% to -8.5%; P = 1.3 x 10(-12)), rs7553007 in the CRP locus (-20.7%; 95% CI, -23.4% to -17.9%; P = 1.3 x 10(-38)), rs1183910 in HNF1A (-13.8%; 95% CI, -16.6% to -10.9%; P = 1.9 x 10(-18)), and rs4420638 in APOE-CI-CII (-21.8%; 95% CI, -25.3% to -18.1%; P = 8.1 x 10(-26)). Association of SNP rs7553007 in the CRP locus with coronary heart disease gave an odds ratio (OR) of 0.98 (95% CI, 0.94 to 1.01) per 20% lower CRP level. Our mendelian randomization study of variants in the CRP locus showed no association with coronary heart disease: OR, 1.00; 95% CI, 0.97 to 1.02; per 20% lower CRP level, compared with OR, 0.94; 95% CI, 0.94 to 0.95; predicted from meta-analysis of the observational studies of CRP levels and coronary heart disease (z score, -3.45; P < .001). SNPs rs6700896 in LEPR (OR, 1.06; 95% CI, 1.02 to 1.09; per minor allele), rs4537545 in IL6R (OR, 0.94; 95% CI, 0.91 to 0.97), and rs4420638 in the APOE-CI-CII cluster (OR, 1.16; 95% CI, 1.12 to 1.21) were all associated with risk of coronary heart disease., Conclusion: The lack of concordance between the effect on coronary heart disease risk of CRP genotypes and CRP levels argues against a causal association of CRP with coronary heart disease.

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N = 38,580) informative for adult waist circumference (WC) and waist-hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P = 1.9x10(-11)) and MSRA (WC, P = 8.9x10(-9)). A third locus, near LYPLAL1, was associated with WHR in women only (P = 2.6x10(-8)). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity., Competing Interests: The authors from Decode Genetics Inc. own stocks and stock options in the company. VM, NL, KS and DMW are full-time employees of GlaxoSmithKline, a pharmaceutical company. The other authors have no competing interests to declare.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5 million genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N ≤ 71,225 European ancestry, N ≤ 12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N = 29,136). We identified association between systolic or diastolic blood pressure and common variants in eight regions near the CYP17A1 (P = 7 × 10(-24)), CYP1A2 (P = 1 × 10(-23)), FGF5 (P = 1 × 10(-21)), SH2B3 (P = 3 × 10(-18)), MTHFR (P = 2 × 10(-13)), c10orf107 (P = 1 × 10(-9)), ZNF652 (P = 5 × 10(-9)) and PLCD3 (P = 1 × 10(-8)) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Genome-wide association studies have identified a region on chromosome 9p that is associated with coronary artery disease (CAD). The region is also associated with type 2 diabetes (T2D), a risk factor for CAD, although different SNPs were reported to be associated to each disease in separate studies. We have undertaken a case-control study in 4251 CAD cases and 4443 controls in four European populations using previously reported ('literature') and tagging SNPs. We replicated the literature SNPs (P = 8x10(-13); OR = 1.29; 95% CI: 1.20-1.38) and showed that the strong consistent association detected by these SNPs is a consequence of a 'yin-yang' haplotype pattern spanning 53 kb. There was no evidence of additional CAD susceptibility alleles over the major risk haplotype. CAD patients without myocardial infarction (MI) showed a trend towards stronger association than MI patients. The CAD susceptibility conferred by this locus did not differ by sex, age, smoking, obesity, hypertension or diabetes. A simultaneous test of CAD and diabetes susceptibility with CAD and T2D-associated SNPs indicated that these associations were independent of each other. Moreover, this region was not associated with differences in plasma levels of low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, fibrinogen, albumin, uric acid, bilirubin or homocysteine, although the CAD-high-risk allele was paradoxically associated with lower triglyceride levels. A large antisense non-coding RNA gene (ANRIL) collocates with the high-risk haplotype, is expressed in tissues and cell types that are affected by atherosclerosis and is a prime candidate gene for the chromosome 9p CAD locus.

Coronary artery disease (CAD) is a leading cause of death world-wide, and most cases have a complex, multifactorial aetiology that includes a substantial heritable component. Identification of new genes involved in CAD may inform pathogenesis and provide new therapeutic targets. The PROCARDIS study recruited 2,658 affected sibling pairs (ASPs) with onset of CAD before age 66 y from four European countries to map susceptibility loci for CAD. ASPs were defined as having CAD phenotype if both had CAD, or myocardial infarction (MI) phenotype if both had a MI. In a first study, involving a genome-wide linkage screen, tentative loci were mapped to Chromosomes 3 and 11 with the CAD phenotype (1,464 ASPs), and to Chromosome 17 with the MI phenotype (739 ASPs). In a second study, these loci were examined with a dense panel of grid-tightening markers in an independent set of families (1,194 CAD and 344 MI ASPs). This replication study showed a significant result on Chromosome 17 (MI phenotype; p = 0.009 after adjustment for three independent replication tests). An exclusion analysis suggests that further genes of effect size lambda(sib) > 1.24 are unlikely to exist in these populations of European ancestry. To our knowledge, this is the first genome-wide linkage analysis to map, and replicate, a CAD locus. The region on Chromosome 17 provides a compelling target within which to identify novel genes underlying CAD. Understanding the genetic aetiology of CAD may lead to novel preventative and/or therapeutic strategies., Competing Interests: Competing interests. The authors have declared that no competing interests exist.

Among bacteria, many species have synonymous codon usage patterns that have been influenced by natural selection for those codons that are translated more accurately and/or efficiently. However, in other species selection appears to have been ineffective. Here, we introduce a population genetics-based model for quantifying the extent to which selection has been effective. The approach is applied to 80 phylogenetically diverse bacterial species for which whole genome sequences are available. The strength of selected codon usage bias, S, is found to vary substantially among species; in 30% of the genomes examined, there was no significant evidence that selection had been effective. Values of S are highly positively correlated with both the number of rRNA operons and the number of tRNA genes. These results are consistent with the hypothesis that species exposed to selection for rapid growth have more rRNA operons, more tRNA genes and more strongly selected codon usage bias. For example, Clostridium perfringens, the species with the highest value of S, can have a generation time as short as 7 min.

Background: Horizontal gene transfer is central to evolution in most bacterial species. The detection of exchanged regions is often based upon analysis of compositional characteristics and their comparison to the organism as a whole. In this study we describe a new methodology combining aspects of established signature analysis with textual analysis approaches. This approach has been used to analyze the two available genome sequences of H. pylori., Results: This gene-by-gene analysis reveals a wide range of genes related to both virulence behaviour and the strain differences that have been relatively recently acquired from other sequence backgrounds. These frequently involve single genes or small numbers of genes that are not associated with transposases or bacteriophage genes, nor with inverted repeats typically used as markers for horizontal transfer. In addition, clear examples of horizontal exchange in genes associated with 'core' metabolic functions were identified, supported by differences between the sequenced strains, including: ftsK, xerD and polA. In some cases it was possible to determine which strain represented the 'parent' and 'altered' states for insertion-deletion events. Different signature component lengths showed different sensitivities for the detection of some horizontally transferred genes, which may reflect different amelioration rates of sequence components., Conclusion: New implementations of signature analysis that can be applied on a gene-by-gene basis for the identification of horizontally acquired sequences are described. These findings highlight the central role of the availability of homologous substrates in evolution mediated by horizontal exchange, and suggest that some components of the supposedly stable 'core genome' may actually be favoured targets for integration of foreign sequences because of their degree of conservation.

Unlabelled: The high cost of genotyping single nucleotide polymorphisms (SNPs) generally prohibits the systematic mapping of entire genetic linkage regions in order to find the polymorphisms associated with increased risk of disease. In practice, SNPs are selected at approximately equal spacing across the linkage region to try to locate a SNP lying in the haplotype block of the disease SNP. The size of the haplotype block may not be known, however, and SNPs taken from public domain sources may not in fact be polymorphic. Our program will choose a subset of the SNPs in a linkage region so as to maximize the expected proportion of the sequence that lies within a given distance of a real SNP., Availability: The software is available, free of charge, for academic use on request from the authors., Supplementary Information: www.oxagen.co.uk

We have examined the phenotypic effects of 21 independent deletions from the fully sequenced and annotated 356 kb telomeric region of the short arm of chromosome 16 (16p13.3). Fifteen genes contained within this region have been highly conserved throughout evolution and encode proteins involved in important housekeeping functions, synthesis of haemoglobin, signalling pathways and critical developmental pathways. Although a priori many of these genes would be considered candidates for critical haploinsufficient genes, none of the deletions within the 356 kb interval cause any discernible phenotype other than alpha thalassaemia whether inherited via the maternal or paternal line. These findings contrast with previous observations on patients with larger (> 1 Mb) deletions from the 16p telomere and therefore address the mechanisms by which monosomy gives rise to human genetic disease.

We have sequenced 1949 kb from the terminal Giemsa light band of human chromosome 16p, enabling us to fully annotate the region extending from the telomeric repeats to the previously published tuberous sclerosis disease 2 (TSC2) and polycystic kidney disease 1 (PKD1) genes. This region can be subdivided into two GC-rich, Alu-rich domains and one GC-rich, Alu-poor domain. The entire region is extremely gene rich, containing 100 confirmed genes and 20 predicted genes. Many of the genes encode widely expressed proteins orchestrating basic cellular processes (e.g. DNA recombination, repair, transcription, RNA processing, signal transduction, intracellular signalling and mRNA translation). Others, such as the alpha globin genes (HBA1 and HBA2), PDIP and BAIAP3, are specialized tissue-restricted genes. Some of the genes have been previously implicated in the pathophysiology of important human genetic diseases (e.g. asthma, cataracts and the ATR-16 syndrome). Others are known disease genes for alpha thalassaemia, adult polycystic kidney disease and tuberous sclerosis. There is also linkage evidence for bipolar affective disorder, epilepsy and autism in this region. Sixty-three chromosomal deletions reported here and elsewhere allow us to interpret the results of removing progressively larger numbers of genes from this well defined human telomeric region.

Phase variation, mediated through variation in the length of simple sequence repeats, is recognized as an important mechanism for controlling the expression of factors involved in bacterial virulence. Phase variation is associated with most of the currently recognized virulence determinants of Neisseria meningitidis. Based upon the complete genome sequence of the N. meningitidis serogroup B strain MC58, we have identified tracts of potentially unstable simple sequence repeats and their potential functional significance determined on the basis of sequence context. Of the 65 potentially phase variable genes identified, only 13 were previously recognized. Comparison with the sequences from the other two pathogenic Neisseria sequencing projects shows differences in the length of the repeats in 36 of the 65 genes identified, including 25 of those not previously known to be phase variable. Six genes that did not have differences in the length of the repeat instead had polymorphisms such that the gene would not be expected to be phase variable in at least one of the other strains. A further 12 candidates did not have homologues in either of the other two genome sequences. The large proportion of these genes that are associated with frameshifts and with differences in repeat length between the neisserial genome sequences is further corroborative evidence that they are phase variable. The number of potentially phase variable genes is substantially greater than for any other species studied to date, and would allow N. meningitidis to generate a very large repertoire of phenotypes through expression of these genes in different combinations. Novel phase variable candidates identified in the strain MC58 genome sequence include a spectrum of genes encoding glycosyltransferases, toxin related products, and metabolic activities as well as several restriction/modification and bacteriocin-related genes and a number of open reading frames (ORFs) for which the function is currently unknown. This suggests that the potential role of phase variation in mediating bacterium-host interactions is much greater than has been appreciated to date. Analysis of the distribution of homopolymeric tract lengths indicates that this species has sequence-specific mutational biases that favour the instability of sequences associated with phase variation.

The 2,272,351-base pair genome of Neisseria meningitidis strain MC58 (serogroup B), a causative agent of meningitis and septicemia, contains 2158 predicted coding regions, 1158 (53.7%) of which were assigned a biological role. Three major islands of horizontal DNA transfer were identified; two of these contain genes encoding proteins involved in pathogenicity, and the third island contains coding sequences only for hypothetical proteins. Insights into the commensal and virulence behavior of N. meningitidis can be gleaned from the genome, in which sequences for structural proteins of the pilus are clustered and several coding regions unique to serogroup B capsular polysaccharide synthesis can be identified. Finally, N. meningitidis contains more genes that undergo phase variation than any pathogen studied to date, a mechanism that controls their expression and contributes to the evasion of the host immune system.

Uptake sequences are abundant sequence motifs, often located downstream of ORFs, that are used to facilitate the within-species horizontal transfer of DNA. A frequent word analysis of the complete genome sequence of Helicobacter pylori strain 26685 was performed to search for and determine the identity of an uptake sequence in this species. The results demonstrated that Hel. pylori does not possess an uptake sequence. This is the first naturally transformable Gram-negative species shown to lack such a transformation-targeting system.

The transposon Tn916 was evaluated as a tool for generalized mutagenesis of the genome of Haemophilus influenzae. This was achieved in silico by searching the genome sequence of H. influenzae Rd for the published Tn916 target site consensus sequence 5' TT/ATTTT(N)6AAAAAA/TA. This search identified 16 putative target sites. In subsequent experiments, integration of Tn916 did not occur at any of these sites. Using the nucleotide sequences of these observed integration sites, a new consensus sequence, 5' TTTTT(N)xAAAAA (4 < or = x < or = 7), was derived. This sequence reflects the curve-twist-curve DNA topology which is a feature common to all Tn916 integration sites. A search of the H. influenzae Rd genome using the new consensus sequence identified 167 potential target sites, representing approximately 1% of the total genome. Only 80 of these sites were located within ORFs. The presence of such a limited number of target sites places severe constraints on the use of Tn916 as a tool for generalized mutagenesis of the genome of H. influenzae.

We describe an integrated system for the analysis of DNA sequence motifs within complete bacterial genome sequences. This system is based around ACeDB, a genome database with an integrated graphical user interface; we identify and display motifs in the context of genetic, sequence and bibliographic data. Tomb et aL (1997) previously reported the identification of contingency genes in Helicobacter pylori through their association with homopolymeric tracts and dinucleotide repeats. With this as a starting point, we validated the system by a search for this type of repeat and used the contextual information to assess the likelihood that they mediate phase variation in the associated open reading frames (ORFs). We found all of the repeats previously described, and identified 27 putative phase-variable genes (including 17 previously described). These could be divided into three groups: lipopolysaccharide (LPS) biosynthesis, cell-surface-associated proteins and DNA restriction/modification systems. Five of the putative genes did not have obvious homologues in any of the public domain sequence databases. The reading frame of some ORFs was disrupted by the presence of the repeats, including the alpha(1-2) fucosyltransferase gene, necessary for the synthesis of the Lewis Y epitope. An additional benefit of this approach is that the results of each search can be analysed further and compared with those from other genomes. This revealed that H. pylori has an unusually high frequency of homopurine:homopyrimidine repeats suggesting mechanistic biases that favour their presence and instability.

Silent sites (positions that can undergo synonymous substitutions) in protein-coding genes can illuminate two evolutionary processes. First, despite being silent, they may be subject to natural selection. Among eukaryotes this is exemplified by yeast, where synonymous codon usage patterns are shaped by selection for particular codons that are more efficiently and/or accurately translated by the most abundant tRNAs; codon usage across the genome, and the abundance of different tRNA species, are highly co-adapted. Second, in the absence of selection, silent sites reveal underlying mutational patterns. Codon usage varies enormously among human genes, and yet silent sites do not appear to be influenced by natural selection, suggesting that mutation patterns vary among regions of the genome. At first, the yeast and human genomes were thought to reflect a dichotomy between unicellular and multicellular organisms. However, it now appears that natural selection shapes codon usage in some multicellular species (e.g. Drosophila and Caenorhabditis), and that regional variations in mutation biases occur in yeast. Silent sites (in serine codons) also provide evidence for mutational events changing adjacent nucleotides simultaneously.