Hundreds of variants clustered in genomic loci and biological pathways affect human height.

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.

Carriers of the fragile X mental retardation 1 (FMR1) premutation allele present with increased levels of cytokine IL-10.

BACKGROUND: Fragile X-associated tremor/ataxia syndrome (FXTAS) is an inherited late-onset neurodegenerative disorder, characterized both by neurological and cognitive deficits. It is caused by the expansion of CGG repeats (55 to 200 repeats) in the noncoding region of the fragile X mental retardation 1 (FMR1) gene. Abnormal immunological patterns are often associated with neurodegenerative disorders and implicated in their etiology. We therefore investigated the immune status of FXTAS patients, which had not been assessed prior to this study. METHOD: Peripheral blood mononuclear cells (PBMCs) were collected from 15 asymptomatic FMR1 premutation carriers and 20 age-matched controls. Concentrations of three cytokines (IL-6, IL-8, IL-10) were measured in PBMC supernatants using ELISA assays. RESULTS: We found a significant increase in the concentration of the major anti-inflammatory cytokine IL-10 in supernatants of PBMCs derived from premutation carriers, when compared with controls (P = 0.019). This increase correlated significantly with the number of CGG repeats (P = 0.002). CONCLUSIONS: Elevated IL-10 levels were observed in all premutation carriers, before appearance of the classical neurological symptoms; therefore, IL-10 may be one of the early biomarkers of FXTAS.

Association of genetic variation with systolic and diastolic blood pressure among African Americans: the Candidate Gene Association Resource study.

The prevalence of hypertension in African Americans (AAs) is higher than in other US groups; yet, few have performed genome-wide association studies (GWASs) in AA. Among people of European descent, GWASs have identified genetic variants at 13 loci that are associated with blood pressure. It is unknown if these variants confer susceptibility in people of African ancestry. Here, we examined genome-wide and candidate gene associations with systolic blood pressure (SBP) and diastolic blood pressure (DBP) using the Candidate Gene Association Resource (CARe) consortium consisting of 8591 AAs. Genotypes included genome-wide single-nucleotide polymorphism (SNP) data utilizing the Affymetrix 6.0 array with imputation to 2.5 million HapMap SNPs and candidate gene SNP data utilizing a 50K cardiovascular gene-centric array (ITMAT-Broad-CARe [IBC] array). For Affymetrix data, the strongest signal for DBP was rs10474346 (P= 3.6 × 10(-8)) located near GPR98 and ARRDC3. For SBP, the strongest signal was rs2258119 in C21orf91 (P= 4.7 × 10(-8)). The top IBC association for SBP was rs2012318 (P= 6.4 × 10(-6)) near SLC25A42 and for DBP was rs2523586 (P= 1.3 × 10(-6)) near HLA-B. None of the top variants replicated in additional AA (n = 11 882) or European-American (n = 69 899) cohorts. We replicated previously reported European-American blood pressure SNPs in our AA samples (SH2B3, P= 0.009; TBX3-TBX5, P= 0.03; and CSK-ULK3, P= 0.0004). These genetic loci represent the best evidence of genetic influences on SBP and DBP in AAs to date. More broadly, this work supports that notion that blood pressure among AAs is a trait with genetic underpinnings but also with significant complexity.

The hourglass and the early conservation models--co-existing patterns of developmental constraints in vertebrates.

Developmental constraints have been postulated to limit the space of feasible phenotypes and thus shape animal evolution. These constraints have been suggested to be the strongest during either early or mid-embryogenesis, which corresponds to the early conservation model or the hourglass model, respectively. Conflicting results have been reported, but in recent studies of animal transcriptomes the hourglass model has been favored. Studies usually report descriptive statistics calculated for all genes over all developmental time points. This introduces dependencies between the sets of compared genes and may lead to biased results. Here we overcome this problem using an alternative modular analysis. We used the Iterative Signature Algorithm to identify distinct modules of genes co-expressed specifically in consecutive stages of zebrafish development. We then performed a detailed comparison of several gene properties between modules, allowing for a less biased and more powerful analysis. Notably, our analysis corroborated the hourglass pattern at the regulatory level, with sequences of regulatory regions being most conserved for genes expressed in mid-development but not at the level of gene sequence, age, or expression, in contrast to some previous studies. The early conservation model was supported with gene duplication and birth that were the most rare for genes expressed in early development. Finally, for all gene properties, we observed the least conservation for genes expressed in late development or adult, consistent with both models. Overall, with the modular approach, we showed that different levels of molecular evolution follow different patterns of developmental constraints. Thus both models are valid, but with respect to different genomic features.

Novel Approach Identifies SNPs in SLC2A10 and KCNK9 with Evidence for Parent-of-Origin Effect on Body Mass Index.

The phenotypic effect of some single nucleotide polymorphisms (SNPs) depends on their parental origin. We present a novel approach to detect parent-of-origin effects (POEs) in genome-wide genotype data of unrelated individuals. The method exploits increased phenotypic variance in the heterozygous genotype group relative to the homozygous groups. We applied the method to >56,000 unrelated individuals to search for POEs influencing body mass index (BMI). Six lead SNPs were carried forward for replication in five family-based studies (of ∼4,000 trios). Two SNPs replicated: the paternal rs2471083-C allele (located near the imprinted KCNK9 gene) and the paternal rs3091869-T allele (located near the SLC2A10 gene) increased BMI equally (beta = 0.11 (SD), P<0.0027) compared to the respective maternal alleles. Real-time PCR experiments of lymphoblastoid cell lines from the CEPH families showed that expression of both genes was dependent on parental origin of the SNPs alleles (P<0.01). Our scheme opens new opportunities to exploit GWAS data of unrelated individuals to identify POEs and demonstrates that they play an important role in adult obesity.

Brain energy consumption induced by electrical stimulation promotes systemic glucose uptake.

BACKGROUND: Controlled transcranial stimulation of the brain is part of clinical treatment strategies in neuropsychiatric diseases such as depression, stroke, or Parkinson's disease. Manipulating brain activity by transcranial stimulation, however, inevitably influences other control centers of various neuronal and neurohormonal feedback loops and therefore may concomitantly affect systemic metabolic regulation. Because hypothalamic adenosine triphosphate-sensitive potassium channels, which function as local energy sensors, are centrally involved in the regulation of glucose homeostasis, we tested whether transcranial direct current stimulation (tDCS) causes an excitation-induced transient neuronal energy depletion and thus influences systemic glucose homeostasis and related neuroendocrine mediators.METHODS: In a crossover design testing 15 healthy male volunteers, we increased neuronal excitation by anodal tDCS versus sham and examined cerebral energy consumption with (31)phosphorus magnetic resonance spectroscopy. Systemic glucose uptake was determined by euglycemic-hyperinsulinemic glucose clamp, and neurohormonal measurements comprised the parameters of the stress systems.RESULTS: We found that anodic tDCS-induced neuronal excitation causes an energetic depletion, as quantified by (31)phosphorus magnetic resonance spectroscopy. Moreover, tDCS-induced cerebral energy consumption promotes systemic glucose tolerance in a standardized euglycemic-hyperinsulinemic glucose clamp procedure and reduces neurohormonal stress axes activity.CONCLUSIONS: Our data demonstrate that transcranial brain stimulation not only evokes alterations in local neuronal processes but also clearly influences downstream metabolic systems regulated by the brain. The beneficial effects of tDCS on metabolic features may thus qualify brain stimulation as a promising nonpharmacologic therapy option for drug-induced or comorbid metabolic disturbances in various neuropsychiatric diseases.

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits.

Genome-wide association and functional follow-up reveals new loci for kidney function.

Chronic kidney disease (CKD) is an important public health problem with a genetic component. We performed genome-wide association studies in up to 130,600 European ancestry participants overall, and stratified for key CKD risk factors. We uncovered 6 new loci in association with estimated glomerular filtration rate (eGFR), the primary clinical measure of CKD, in or near MPPED2, DDX1, SLC47A1, CDK12, CASP9, and INO80. Morpholino knockdown of mpped2 and casp9 in zebrafish embryos revealed podocyte and tubular abnormalities with altered dextran clearance, suggesting a role for these genes in renal function. By providing new insights into genes that regulate renal function, these results could further our understanding of the pathogenesis of CKD.

Distinct levels in Pom1 gradients limit Cdr2 activity and localization to time and position division.

Where and when cells divide are fundamental questions. In rod-shaped fission yeast cells, the DYRK-family kinase Pom1 is organized in concentration gradients from cell poles and controls cell division timing and positioning. Pom1 gradients restrict to mid-cell the SAD-like kinase Cdr2, which recruits Mid1/Anillin for medial division. Pom1 also delays mitotic commitment through Cdr2, which inhibits Wee1. Here, we describe quantitatively the distributions of cortical Pom1 and Cdr2. These reveal low profile overlap contrasting with previous whole-cell measurements and Cdr2 levels increase with cell elongation, raising the possibility that Pom1 regulates mitotic commitment by controlling Cdr2 medial levels. However, we show that distinct thresholds of Pom1 activity define the timing and positioning of division. Three conditions-a separation-of-function Pom1 allele, partial downregulation of Pom1 activity, and haploinsufficiency in diploid cells-yield cells that divide early, similar to pom1 deletion, but medially, like wild-type cells. In these cells, Cdr2 is localized correctly at mid-cell. Further, Cdr2 overexpression promotes precocious mitosis only in absence of Pom1. Thus, Pom1 inhibits Cdr2 for mitotic commitment independently of regulating its localization or cortical levels. Indeed, we show Pom1 restricts Cdr2 activity through phosphorylation of a C-terminal self-inhibitory tail. In summary, our results demonstrate that distinct levels in Pom1 gradients delineate a medial Cdr2 domain, for cell division placement, and control its activity, for mitotic commitment.

Using transcription modules to identify expression clusters perturbed in Williams-Beuren syndrome.

The genetic dissection of the phenotypes associated with Williams-Beuren Syndrome (WBS) is advancing thanks to the study of individuals carrying typical or atypical structural rearrangements, as well as in vitro and animal studies. However, little is known about the global dysregulations caused by the WBS deletion. We profiled the transcriptomes of skin fibroblasts from WBS patients and compared them to matched controls. We identified 868 differentially expressed genes that were significantly enriched in extracellular matrix genes, major histocompatibility complex (MHC) genes, as well as genes in which the products localize to the postsynaptic membrane. We then used public expression datasets from human fibroblasts to establish transcription modules, sets of genes coexpressed in this cell type. We identified those sets in which the average gene expression was altered in WBS samples. Dysregulated modules are often interconnected and share multiple common genes, suggesting that intricate regulatory networks connected by a few central genes are disturbed in WBS. This modular approach increases the power to identify pathways dysregulated in WBS patients, thus providing a testable set of additional candidates for genes and their interactions that modulate the WBS phenotypes.