The chromosome 3q25 locus associated with fetal adiposity is not associated with childhood adiposity

Increased newborn adiposity is associated with later adverse metabolic outcomes. Previous genome-wide association studies (GWAS) demonstrated strong association of a locus on chromosome 3 (3q25.31) with newborn sum of skinfolds, a measure of overall adiposity. Whether this locus is associated with childhood adiposity is unknown. Genotype and sum of skinfolds data were available for 293 children at birth and age 2, and for 350 children at birth and age 6 from a European cohort (Belfast, UK) who participated in the Hyperglycemia and Adverse Pregnancy Outcome GWAS. We examined single nucleotide polymorphisms (SNPs) at the 3q25.31 locus associated with newborn adiposity. Linear regression analyses under an additive genetic model adjusting for maternal body mass index were performed. In both cohorts, a positive association was observed between all SNPs and sum of skinfolds at birth (P=2.3 × 10(-4), β=0.026 and P=4.8 × 10(-4), β=0.025). At the age of 2 years, a non-significant negative association was observed with sum of skinfolds (P=0.06; β =-0.015). At the age of 6 years, there was no evidence of association (P=0.86; β=0.002). The 3q25.31 locus strongly associated with newborn adiposity had no significant association with childhood adiposity suggesting that its impact may largely be limited to fetal fat accretion.

Genetic risk factors affecting mitochondrial function are associated with kidney disease in people with Type 1 diabetes

AIM: To evaluate the association with diabetic kidney disease of single nucleotide polymorphisms (SNPs) that may contribute to mitochondrial dysfunction.

METHODS: The mitochondrial genome and 1039 nuclear genes that are integral to mitochondrial function were investigated using a case (n=823 individuals with diabetic kidney disease) vs. control (n=903 individuals with diabetes and no renal disease) approach. All people included in the analysis were of white European origin and were diagnosed with Type 1 diabetes before the age of 31 years. Replication was conducted in 5093 people with similar phenotypes to those of the discovery collection. Association analyses were performed using the plink genetic analysis toolset, with adjustment for relevant covariates.

RESULTS: A total of 25 SNPs were evaluated in the mitochondrial genome, but none were significantly associated with diabetic kidney disease or end-stage renal disease. A total of 38 SNPs in nuclear genes influencing mitochondrial function were nominally associated with diabetic kidney disease and 16 SNPS were associated with end-stage renal disease, secondary to diabetic kidney disease, with meta-analyses confirming the same direction of effect. Three independent signals (seven SNPs) were common to the replication data for both phenotypes with Type 1 diabetes and persistent proteinuria or end-stage renal disease.

CONCLUSIONS: Our results suggest that SNPs in nuclear genes that influence mitochondrial function are significantly associated with diabetic kidney disease in a white European population

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

Introduction: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers.

Methods: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals.

Results: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk.

Conclusions: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.

Abundant genetic overlap between blood lipids and immune-mediated diseases indicates shared molecular genetic mechanisms

Epidemiological studies suggest a relationship between blood lipids and immune-mediated diseases, but the nature of these associations is not well understood. We used genome-wide association studies (GWAS) to investigate shared single nucleotide polymorphisms (SNPs) between blood lipids and immune-mediated diseases. We analyzed data from GWAS (n~200,000 individuals), applying new False Discovery Rate (FDR) methods, to investigate genetic overlap between blood lipid levels [triglycerides (TG), low density lipoproteins (LDL), high density lipoproteins (HDL)] and a selection of archetypal immune-mediated diseases (Crohn's disease, ulcerative colitis, rheumatoid arthritis, type 1 diabetes, celiac disease, psoriasis and sarcoidosis). We found significant polygenic pleiotropy between the blood lipids and all the investigated immune-mediated diseases. We discovered several shared risk loci between the immune-mediated diseases and TG (n = 88), LDL (n = 87) and HDL (n = 52). Three-way analyses differentiated the pattern of pleiotropy among the immune-mediated diseases. The new pleiotropic loci increased the number of functional gene network nodes representing blood lipid loci by 40%. Pathway analyses implicated several novel shared mechanisms for immune pathogenesis and lipid biology, including glycosphingolipid synthesis (e.g. FUT2) and intestinal host-microbe interactions (e.g. ATG16L1). We demonstrate a shared genetic basis for blood lipids and immune-mediated diseases independent of environmental factors. Our findings provide novel mechanistic insights into dyslipidemia and immune-mediated diseases and may have implications for therapeutic trials involving lipid-lowering and anti-inflammatory agents.

Shared common variants in prostate cancer and blood lipids

BACKGROUND: Epidemiological and clinical studies suggest comorbidity between prostate cancer (PCA) and cardiovascular disease (CVD) risk factors. However, the relationship between these two phenotypes is still not well understood. Here we sought to identify shared genetic loci between PCA and CVD risk factors.

METHODS: We applied a genetic epidemiology method based on conjunction false discovery rate (FDR) that combines summary statistics from different genome-wide association studies (GWAS), and allows identification of genetic overlap between two phenotypes. We evaluated summary statistics from large, multi-centre GWA studies of PCA (n=50 000) and CVD risk factors (n=200 000) [triglycerides (TG), low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, body mass index, waist-hip ratio and type 2 diabetes (T2D)]. Enrichment of single nucleotide polymorphisms (SNPs) associated with PCA and CVD risk factors was assessed with conditional quantile-quantile plots and the Anderson-Darling test. Moreover, we pinpointed shared loci using conjunction FDR.

RESULTS: We found the strongest enrichment of P-values in PCA was conditional on LDL and conditional on TG. In contrast, we found only weak enrichment conditional on HDL or conditional on the other traits investigated. Conjunction FDR identified altogether 17 loci; 10 loci were associated with PCA and LDL, 3 loci were associated with PCA and TG and additionally 4 loci were associated with PCA, LDL and TG jointly (conjunction FDR <0.01). For T2D, we detected one locus adjacent to HNF1B.

CONCLUSIONS: We found polygenic overlap between PCA predisposition and blood lipids, in particular LDL and TG, and identified 17 pleiotropic gene loci between PCA and LDL, and PCA and TG, respectively. These findings provide novel pathobiological insights and may have implications for trials using targeting lipid-lowering agents in a prevention or cancer setting.

Genome-wide association study for sight-threatening diabetic retinopathy reveals association with genetic variation near the GRB2 gene

AIMS/HYPOTHESIS: Diabetic retinopathy is a serious complication of diabetes mellitus and can lead to blindness. A genetic component, in addition to traditional risk factors, has been well described although strong genetic factors have not yet been identified. Here, we aimed to identify novel genetic risk factors for sight-threatening diabetic retinopathy using a genome-wide association study.

METHODS: Retinopathy was assessed in white Australians with type 2 diabetes mellitus. Genome-wide association analysis was conducted for comparison of cases of sight-threatening diabetic retinopathy (n = 336) with diabetic controls with no retinopathy (n = 508). Top ranking single nucleotide polymorphisms were typed in a type 2 diabetes replication cohort, a type 1 diabetes cohort and an Indian type 2 cohort. A mouse model of proliferative retinopathy was used to assess differential expression of the nearby candidate gene GRB2 by immunohistochemistry and quantitative western blot.

RESULTS: The top ranked variant was rs3805931 with p = 2.66 × 10(-7), but no association was found in the replication cohort. Only rs9896052 (p = 6.55 × 10(-5)) was associated with sight-threatening diabetic retinopathy in both the type 2 (p = 0.035) and the type 1 (p = 0.041) replication cohorts, as well as in the Indian cohort (p = 0.016). The study-wide meta-analysis reached genome-wide significance (p = 4.15 × 10(-8)). The GRB2 gene is located downstream of this variant and a mouse model of retinopathy showed increased GRB2 expression in the retina.

CONCLUSIONS/INTERPRETATION: Genetic variation near GRB2 on chromosome 17q25.1 is associated with sight-threatening diabetic retinopathy. Several genes in this region are promising candidates and in particular GRB2 is upregulated during retinal stress and neovascularisation.

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled 'reservoir' host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers. © 2012 Biek et al.

Integrin genetic variants and stage-specific tumor recurrence in patients with stage II and III colon cancer

Integrins (ITGs) are key elements in cancer biology, regulating tumor growth, angiogenesis and lymphangiogenesis through interactions of the tumor cells with the microenvironment. Moving from the hypothesis that ITGs could have different effects in stage II and III colon cancer, we tested whether a comprehensive panel of germline single-nucleotide polymorphisms (SNPs) in ITG genes could predict stage-specific time to tumor recurrence (TTR). A total of 234 patients treated with 5-fluorouracil-based chemotherapy at the University of Southern California were included in this study. Whole-blood samples were analyzed for germline SNPs in ITG genes using PCR-restriction fragment length polymorphism or direct DNA sequencing. In the multivariable analysis, stage II colon cancer patients with at least one G allele for ITGB3 rs4642 had higher risk of recurrence (hazard ratio (HR)=4.027, 95% confidence interval (95% CI) 1.556-10.421, P=0.004). This association was also significant in the combined stage II-III cohort (HR=1.975, 95% CI 1.194-3.269, P=0.008). The predominant role of ITGB3 rs4642 in stage II diseases was confirmed using recursive partitioning, showing that ITGB3 rs4642 was the most important factor in stage II diseases. In contrast, in stage III diseases the combined analysis of ITGB1 rs2298141 and ITGA4 rs7562325 allowed to identify three distinct prognostic subgroups (P=0.009). The interaction between stage and the combined ITGB1 rs2298141 and ITGA4 rs7562325 on TTR was significant (P=0.025). This study identifies germline polymorphisms in ITG genes as independent stage-specific prognostic markers for stage II and III colon cancer. These data may help to select subgroups of patients who may benefit from ITG-targeted treatments.

The Cyclin D1 (CCND1) A870G polymorphism predicts clinical outcome to lapatinib and capecitabine in HER2-positive metastatic breast cancer

BACKGROUND: Lapatinib plus capecitabine emerged as an efficacious therapy in metastatic breast cancer (mBC). We aimed to identify germline single-nucleotide polymorphisms (SNPs) in genes involved in capecitabine catabolism and human epidermal receptor signaling that were associated with clinical outcome to assist in selecting patients likely to benefit from this combination.

PATIENTS AND METHODS: DNA was extracted from 240 of 399 patients enrolled in EGF100151 clinical trial (NCT00078572; clinicaltrials.gov) and SNPs were successfully evaluated in 234 patients. The associations between SNPs and clinical outcome were analyzed using Fisher's exact test, Kaplan-Meier curves, log-rank tests, likelihood ratio test within logistic or Cox regression model, as appropriate.

RESULTS: There were significant interactions between CCND1 A870G and clinical outcome. Patients carrying the A-allele were more likely to benefit from lapatinib plus capecitabine versus capecitabine when compared with patients harboring G/G (P = 0.022, 0.024 and 0.04, respectively). In patients with the A-allele, the response rate (RR) was significantly higher with lapatinib plus capecitabine (35%) compared with capecitabine (11%; P = 0.001) but not between treatments in patients with G/G (RR = 24% and 32%, respectively; P = 0.85). Time to tumor progression (TTP) was longer in patients with the A-allele treated with lapatinib plus capecitabine compared with capecitabine (median TTP = 7.9 and 3.4 months; P < 0.001), but not in patients with G/G (median TTP = 6.1 and 6.6 months; P = 0.92).

CONCLUSION: Our findings suggest that CCND1A870G may be useful in predicting clinical outcome in HER2-positive mBC patients treated with lapatinib plus capecitabine.

Adiposity as a cause of cardiovascular disease: a Mendelian randomization study

Background: Adiposity, as indicated by body mass index (BMI), has been associated with risk of cardiovascular diseases in epidemiological studies. We aimed to investigate if these associations are causal, using Mendelian randomization (MR) methods.

Methods: The associations of BMI with cardiovascular outcomes [coronary heart disease (CHD), heart failure and ischaemic stroke], and associations of a genetic score (32 BMI single nucleotide polymorphisms) with BMI and cardiovascular outcomes were examined in up to 22 193 individuals with 3062 incident cardiovascular events from nine prospective follow-up studies within the ENGAGE consortium. We used random-effects meta-analysis in an MR framework to provide causal estimates of the effect of adiposity on cardiovascular outcomes.

Results: There was a strong association between BMI and incident CHD (HR = 1.20 per SD-increase of BMI, 95% CI, 1.12–1.28, P = 1.9·10−7), heart failure (HR = 1.47, 95% CI, 1.35–1.60, P = 9·10−19) and ischaemic stroke (HR = 1.15, 95% CI, 1.06–1.24, P = 0.0008) in observational analyses. The genetic score was robustly associated with BMI (β = 0.030 SD-increase of BMI per additional allele, 95% CI, 0.028–0.033, P = 3·10−107). Analyses indicated a causal effect of adiposity on development of heart failure (HR = 1.93 per SD-increase of BMI, 95% CI, 1.12–3.30, P = 0.017) and ischaemic stroke (HR = 1.83, 95% CI, 1.05–3.20, P = 0.034). Additional cross-sectional analyses using both ENGAGE and CARDIoGRAMplusC4D data showed a causal effect of adiposity on CHD.

Conclusions: Using MR methods, we provide support for the hypothesis that adiposity causes CHD, heart failure and, previously not demonstrated, ischaemic stroke.

Genetic Sharing with Cardiovascular Disease Risk Factors and Diabetes Reveals Novel Bone Mineral Density Loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR < 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.

SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Gene-associated markers provide tools for tackling illegal fishing and false eco-certification

Illegal, Unreported and Unregulated fishing has had a major role in the overexploitation of global fish populations. In response, international regulations have been imposed and many fisheries have been 'eco-certified' by consumer organizations, but methods for independent control of catch certificates and eco-labels are urgently needed. Here we show that, by using gene-associated single nucleotide polymorphisms, individual marine fish can be assigned back to population of origin with unprecedented high levels of precision. By applying high differentiation single nucleotide polymorphism assays, in four commercial marine fish, on a pan-European scale, we find 93-100% of individuals could be correctly assigned to origin in policy-driven case studies. We show how case-targeted single nucleotide polymorphism assays can be created and forensically validated, using a centrally maintained and publicly available database. Our results demonstrate how application of gene-associated markers will likely revolutionize origin assignment and become highly valuable tools for fighting illegal fishing and mislabelling worldwide.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Application of SNPs for population genetics of nonmodel organisms:new opportunities and challenges

Recent improvements in the speed, cost and accuracy of next generation sequencing are revolutionizing the discovery of single nucleotide polymorphisms (SNPs). SNPs are increasingly being used as an addition to the molecular ecology toolkit in nonmodel organisms, but their efficient use remains challenging. Here, we discuss common issues when employing SNP markers, including the high numbers of markers typically employed, the effects of ascertainment bias and the inclusion of nonneutral loci in a marker panel. We provide a critique of considerations specifically associated with the application and population genetic analysis of SNPs in nonmodel taxa, focusing specifically on some of the most commonly applied methods.