Defining the role of common variation in the genomic and biological architecture of adult human height

Using genome-wide data from 253,288 individuals, we identified 697 variants at genome-wide significance that together explained one-fifth of the heritability for adult height. By testing different numbers of variants in independent studies, we show that the most strongly associated 1/42,000, 1/43,700 and 1/49,500 SNPs explained 1/421%, 1/424% and 1/429% of phenotypic variance. Furthermore, all common variants together captured 60% of heritability. The 697 variants clustered in 423 loci were enriched for genes, pathways and tissue types known to be involved in growth and together implicated genes and pathways not highlighted in earlier efforts, such as signaling by fibroblast growth factors, WNT/I 2-catenin and chondroitin sulfate-related genes. We identified several genes and pathways not previously connected with human skeletal growth, including mTOR, osteoglycin and binding of hyaluronic acid. Our results indicate a genetic architecture for human height that is characterized by a very large but finite number (thousands) of causal variants.

New genetic loci link adipose and insulin biology to body fat distribution

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

A Candidate Gene Association Study Identifies DAPL1 as a Female-Specific Susceptibility Locus for Age-Related Macular Degeneration (AMD)

Age-related macular degeneration (AMD) is the leading cause of blindness among white caucasians over the age of 50 years with a prevalence rate expected to increase markedly with an anticipated increase in the life span of the world population. To further expand our knowledge of the genetic architecture of the disease, we pursued a candidate gene approach assessing 25 genes and a total of 109 variants. Of these, synonymous single nucleotide polymorphism (SNP) rs17810398 located in death-associated protein-like 1 (DAPL1) was found to be associated with AMD in a joint analysis of 3,229 cases and 2,835 controls from five studies [combined P ADJ = 1.15 × 10(-6), OR 1.332 (1.187-1.496)]. This association was characterized by a highly significant sex difference (P diff = 0.0032) in that it was clearly confined to females with genome-wide significance [P ADJ = 2.62 × 10(-8), OR 1.541 (1.324-1.796); males: P ADJ = 0.382, OR 1.084 (0.905-1.298)]. By targeted resequencing of risk and non-risk associated haplotypes in the DAPL1 locus, we identified additional potentially functional risk variants, namely a common 897-bp deletion and a SNP predicted to affect a putative binding site of an exonic splicing enhancer. We show that the risk haplotype correlates with a reduced retinal transcript level of two, less frequent, non-canonical DAPL1 isoforms. DAPL1 plays a role in epithelial differentiation and may be involved in apoptotic processes thereby suggesting a possible novel pathway in AMSaveD pathogenesis.

Genomic lesions associated with a different clinical outcome in diffuse large B-Cell lymphoma treated with R-CHOP-21

Despite recent therapeutic improvements, the clinical course of diffuse large B-cell lymphoma (DLBCL) still differs considerably among patients. We conducted this retrospective multi-centre study to evaluate the impact of genomic aberrations detected using a high-density genome wide-single nucleotide polymorphism-based array on clinical outcome in a population of DLBCL patients treated with R-CHOP-21 (rituximab, cyclophosphamide, doxorubicine, vincristine and prednisone repeated every 21_d). 166 DNA samples were analysed using the GeneChip Human Mapping 250K NspI. Genomic anomalies were analysed regarding their impact on the clinical course of 124 patients treated with R-CHOP-21. Unsupervised clustering was performed to identify genetically related subgroups of patients with different clinical outcomes. Twenty recurrent genetic lesions showed an impact on the clinical course. Loss of genomic material at 8p23.1 showed the strongest statistical significance and was associated with additional aberrations, such as 17p- and 15q-. Unsupervised clustering identified five DLBCL clusters with distinct genetic profiles, clinical characteristics and outcomes. Genetic features and clusters, associated with a different outcome in patients treated with R-CHOP, have been identified by arrayCGH.

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.

A population study of binocular function

As part of a genome-wide association study (GWAS) of perceptual traits in healthy adults, we measured stereo acuity, the duration of alternative percepts in binocular rivalry and the extent of dichoptic masking in 1060 participants. We present the distributions of the measures, the correlations between measures, and their relationships to other psychophysical traits. We report sex differences, and correlations with age, interpupillary distance, eye dominance, phorias, visual acuity and personality. The GWAS, using data from 988 participants, yielded one genetic association that passed a permutation test for significance: The variant rs1022907 in the gene VTI1A was associated with self-reported ability to see autostereograms. We list a number of other suggestive genetic associations (p<10^-5).

Genetics of Diabetic Nephropathy: a Long Road of Discovery

The global prevalence of diabetic nephropathy is rising in parallel with the increasing incidence of diabetes in most countries. Unfortunately, up to 40 % of persons diagnosed with diabetes may develop kidney complications. Diabetic nephropathy is associated with substantially increased risks of cardiovascular disease and premature mortality. An inherited susceptibility to diabetic nephropathy exists, and progress is being made unravelling the genetic basis for nephropathy thanks to international research collaborations, shared biological resources and new analytical approaches. Multiple epidemiological studies have highlighted the clinical heterogeneity of nephropathy and the need for better phenotyping to help define important subgroups for analysis and increase the power of genetic studies. Collaborative genome-wide association studies for nephropathy have reported unique genes, highlighted novel biological pathways and suggested new disease mechanisms, but progress towards clinically relevant risk prediction models for diabetic nephropathy has been slow. This review summarises the current status, recent developments and ongoing challenges elucidating the genetics of diabetic nephropathy.

Polymorphisms near TBX5 and GDF7 are associated with increased risk for Barrett's esophagus

BACKGROUND & AIMS: Barrett's esophagus (BE) increases the risk of esophageal adenocarcinoma (EAC). We found the risk to be BE has been associated with single nucleotide polymorphisms (SNPs) on chromosome 6p21 (within the HLA region) and on 16q23, where the closest protein-coding gene is FOXF1. Subsequently, the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) identified risk loci for BE and esophageal adenocarcinoma near CRTC1 and BARX1, and within 100 kb of FOXP1. We aimed to identify further SNPs that increased BE risk and to validate previously reported associations.

METHODS: We performed a genome-wide association study (GWAS) to identify variants associated with BE and further analyzed promising variants identified by BEACON by genotyping 10,158 patients with BE and 21,062 controls.

RESULTS: We identified 2 SNPs not previously associated with BE: rs3072 (2p24.1; odds ratio [OR] = 1.14; 95% CI: 1.09-1.18; P = 1.8 × 10(-11)) and rs2701108 (12q24.21; OR = 0.90; 95% CI: 0.86-0.93; P = 7.5 × 10(-9)). The closest protein-coding genes were respectively GDF7 (rs3072), which encodes a ligand in the bone morphogenetic protein pathway, and TBX5 (rs2701108), which encodes a transcription factor that regulates esophageal and cardiac development. Our data also supported in BE cases 3 risk SNPs identified by BEACON (rs2687201, rs11789015, and rs10423674). Meta-analysis of all data identified another SNP associated with BE and esophageal adenocarcinoma: rs3784262, within ALDH1A2 (OR = 0.90; 95% CI: 0.87-0.93; P = 3.72 × 10(-9)).

CONCLUSIONS: We identified 2 loci associated with risk of BE and provided data to support a further locus. The genes we found to be associated with risk for BE encode transcription factors involved in thoracic, diaphragmatic, and esophageal development or proteins involved in the inflammatory response.

Plasma Clusterin Levels and the rs11136000 Genotype in Individuals with Mild Cognitive Impairment and Alzheimer's Disease

Aim: Substantial evidence links atherosclerosis and Alzheimer's disease (AD). Apolipoproteins, such as apolipoprotein E, have a causal relationship with both diseases. The rs11136000 SNP within the CLU gene, which encodes clusterin (apolipoprotein J), is also associated with increased AD risk. The aim of this study was to investigate the relationship between plasma clusterin and the rs11136000 genotype in mild cognitive impairment (MCI) and AD.

Methods: Plasma and DNA samples were collected from control, MCI and AD subjects (n=142, 111, 154, respectively). Plasma clusterin was determined by ELISA and DNA samples were genotyped for rs11136000 by TaqMan assay.

Results: Plasma clusterin levels were higher in MCI and AD subjects vs. controls (222.3 +/- 61.3 and 193.6 +/- 58.2 vs. 178.6 +/- 52.3 mu g/ml, respectively; p

Conclusion: This study examined control, MCI and AD subjects, identifying for the first time that plasma clusterin levels were influenced, not only by the presence of AD, but also the transitional stage of MCI, while rs11136000 genotype only influenced plasma clusterin levels in the control group. The increase in plasma clusterin in MCI and AD subjects may occur in response to the disease process and would be predicted to increase binding capacity for amyloid-beta peptides in plasma, enhancing their removal from the brain.

Exposure to aflatoxin B1 in utero is associated with DNA methylation in white blood cells of infants in The Gambia

BACKGROUND: Exposure to environmental toxins during embryonic development may lead to epigenetic changes that influence disease risk in later life. Aflatoxin is a contaminant of staple foods in sub-Saharan Africa, is a known human liver carcinogen and has been associated with stunting in infants.

METHODS: We have measured aflatoxin exposure in 115 pregnant women in The Gambia and examined the DNA methylation status of white blood cells from their infants at 2-8 months old (mean 3.6 ± 0.9). Aflatoxin exposure in women was assessed using an ELISA method to measure aflatoxin albumin (AF-alb) adducts in plasma taken at 1-16 weeks of pregnancy. Genome-wide DNA methylation of infant white blood cells was measured using the Illumina Infinium HumanMethylation450beadchip.

RESULTS: AF-alb levels ranged from 3.9 to 458.4 pg/mg albumin. We found that aflatoxin exposure in the mothers was associated to DNA methylation in their infants for 71 CpG sites (false discovery rate < 0.05), with an average effect size of 1.7% change in methylation. Aflatoxin-associated differential methylation was observed in growth factor genes such as FGF12 and IGF1, and immune-related genes such as CCL28, TLR2 and TGFBI. Moreover, one aflatoxin-associated methylation region (corresponding to the miR-4520b locus) was identified.

CONCLUSIONS: This study shows that maternal exposure to aflatoxin during the early stages of pregnancy is associated with differential DNA methylation patterns of infants, including in genes related to growth and immune function. This reinforces the need for interventions to reduce aflatoxin exposure, especially during critical periods of fetal and infant development.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data

BACKGROUND: Urothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.

METHODS: In our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.

RESULTS: We found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.

CONCLUSION: To our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

Background: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

Results: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

Conclusions: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.

Convergent genetic and expression data implicate immunity in Alzheimer's disease

Background: Late-onset Alzheimer's disease (AD) is heritable with 20 genes showing genome-wide association in the International Genomics of Alzheimer's Project (IGAP). To identify the biology underlying the disease, we extended these genetic data in a pathway analysis.

Methods: The ALIGATOR and GSEA algorithms were used in the IGAP data to identify associated functional pathways and correlated gene expression networks in human brain.

Results: ALIGATOR identified an excess of curated biological pathways showing enrichment of association. Enriched areas of biology included the immune response (P = 3.27 X 10(-12) after multiple testing correction for pathways), regulation of endocytosis (P = 1.31 X 10(-11)), cholesterol transport (P = 2.96 X 10(-9)), and proteasome-ubiquitin activity (P = 1.34 X 10(-6)). Correlated gene expression analysis identified four significant network modules, all related to the immune response (corrected P = .002-.05).

Conclusions: The immime response, regulation of endocytosis, cholesterol transport, and protein ubiquitination represent prime targets for AD therapeutics. (C) 2015 Published by Elsevier Inc. on behalf of The Alzheimer's Association.

Integrative analysis of chemo-transcriptomic profiles for drug-feature specific gene expression signatures

One of the major challenges in systems biology is to understand the complex responses of a biological system to external perturbations or internal signalling depending on its biological conditions. Genome-wide transcriptomic profiling of cellular systems under various chemical perturbations allows the manifestation of certain features of the chemicals through their transcriptomic expression profiles. The insights obtained may help to establish the connections between human diseases, associated genes and therapeutic drugs. The main objective of this study was to systematically analyse cellular gene expression data under various drug treatments to elucidate drug-feature specific transcriptomic signatures. We first extracted drug-related information (drug features) from the collected textual description of DrugBank entries using text-mining techniques. A novel statistical method employing orthogonal least square learning was proposed to obtain drug-feature-specific signatures by integrating gene expression with DrugBank data. To obtain robust signatures from noisy input datasets, a stringent ensemble approach was applied with the combination of three techniques: resampling, leave-one-out cross validation, and aggregation. The validation experiments showed that the proposed method has the capacity of extracting biologically meaningful drug-feature-specific gene expression signatures. It was also shown that most of signature genes are connected with common hub genes by regulatory network analysis. The common hub genes were further shown to be related to general drug metabolism by Gene Ontology analysis. Each set of genes has relatively few interactions with other sets, indicating the modular nature of each signature and its drug-feature-specificity. Based on Gene Ontology analysis, we also found that each set of drug feature (DF)-specific genes were indeed enriched in biological processes related to the drug feature. The results of these experiments demonstrated the pot- ntial of the method for predicting certain features of new drugs using their transcriptomic profiles, providing a useful methodological framework and a valuable resource for drug development and characterization.