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.

Meta-analyses identify 13 loci associated with age at menopause and highlight DNA repair and immune pathways.

To newly identify loci for age at natural menopause, we carried out a meta-analysis of 22 genome-wide association studies (GWAS) in 38,968 women of European descent, with replication in up to 14,435 women. In addition to four known loci, we identified 13 loci newly associated with age at natural menopause (at P < 5 × 10(-8)). Candidate genes located at these newly associated loci include genes implicated in DNA repair (EXO1, HELQ, UIMC1, FAM175A, FANCI, TLK1, POLG and PRIM1) and immune function (IL11, NLRP11 and PRRC2A (also known as BAT2)). Gene-set enrichment pathway analyses using the full GWAS data set identified exoDNase, NF-κB signaling and mitochondrial dysfunction as biological processes related to timing of menopause.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice.

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

De novo T cell generation and cell cycle analysis of CD4 and CD8 T cells during HIV-disease

RESUME POUR UN LARGE PUBLIC Parmi les globules blancs, les lymphocytes T 004 jouent un rôle primordial dans la coordination de la réponse immunitaire contre les pathogènes et les lymphocytes T CD8 dans leur élimination. Lors d'une infection par le virus de l'immunodéficience humaine (VIH-1), non seulement les cellules T CD4 sont les principales cibles d'infections, mais aussi elles disparaissent progressivement tout au long de la maladie. Ce phénomène, appelé aussi épuisement des lymphocytes T CD4, est la principale cause provoquant le Syndrome d'Immunodéficience Acquise (SIDA). Malgré de grands efforts de recherche, nous ne sommes toujours pas en mesure de dire si ce phénomène est dû à un défaut dans la production de nouvelles cellules ou à une destruction massive de cellules en circulation. Dans cette étude, nous nous proposions, dans un premier temps, de comparer la production de nouvelles cellules T CD4 et CD8 chez des individus VIH-négatifs et positifs. Les cellules nouvellement produites portent un marqueur commun que l'on appelle TREC et qui est facilement mesurable. En considérant des paramètres cliniques, nous étions en mesure de déterminer le niveau de TRECs de cellules T CD4 et CD8 dans différentes phases de la maladie. De là, nous avons pu déterminer que le niveau de TREC est toujours plus bas dans les cellules T CD8 de patients VIH-positifs comparativement à notre groupe contrôle. Nous avons pu déterminer par une analyse ultérieure que cette différence est due à une forte prolifération de ces cellules chez les patients VIH-positifs, ce qui a pour effet de diluer ce marqueur. En revanche, la production de nouvelles cellules T CD4 chez des patients VIH-positifs est accentuée lors de la phase précoce de la maladie et largement réprimée lors de la phase tardive. Dans un second temps, nous avons effectué une analyse à grande échelle de l'expression de gènes associés à la division cellulaire sur des lymphocytes T CD4 et CD8 d'individus VIH-¬positifs et négatifs, avec comme contrôle des cellules proliférant in vitro. De cette étude, nous avons pu conclure que les cellules T CD8 de patients VIH-positifs étaient en état de prolifération, alors que les lymphocytes T CD4 présentaient des défauts majeurs conduisant à un arrêt de la division cellulaire. Nos résultats montrent que la capacité à produire de nouvelles cellules chez des patients VIH¬positifs reste active longtemps pendant la maladie, mais que l'incapacité des cellules T CD4 à proliférer peut enrayer la reconstitution immunitaire chez ces individus. ABSTRACT The hallmark of HIV-1 infection is the depletion of CD4 T cells. Despite extensive investigation, the mechanisms responsible for the loss of CD4 T cells have been elucidated only partially. In particular, it remains controversial whether CD4 T cell depletion results from a defect in T cell production or from a massive peripheral destruction. In this study, de novo T cell generation has been investigated by measuring T cell receptor rearrangement excision circles (TRECs) on large cohorts of HIV-negative (N=120) and HIV-1 infected (N=298) individuals. Analysis of TREC levels was performed in HIV-infected subjects stratified by the stage of HIV disease based on CD4 T cell counts (early: >500 CD4 T cells/µl; intermediate: <500>200; late: <200) and by age (20 to 60 years, n = 259). Our data show that TREC levels in CD8 T cells were significantly lower in HIV-infected subjects at any stage of disease compared to the control group. In contrast, TREC levels in CD4 T cells were significantly higher in HIV-infected subjects at early stages disease while no significant differences were observed at intermediate stages of the disease and were severely reduced only at late stages of disease. To investigate further the status of cell cycle in peripheral CD4 and CD8 T cells in HIV-1 infections, we determined the pattern of gene expression with the microarray technology. In particular, CD4 and CD8 T cells of HIV-1 infected and HIV-negative subjects were analysed by Cell Cycle cDNA expression array. The patterns of gene expression were compared to in vitro stimulated CD4 and CD8 T cells and this analysis showed that CD8 T cells of HIV-1 infected subjects had a pattern of gene expression very similar to that of in vitro stimulated CD8 T cells thus indicating ongoing cell cycling. In contrast, CD4 T cells of HIV-1 infected subjects displayed a complex pattern of gene expression. In fact, CD4 T cells expressed high levels of genes typically associated with cell activation, but low levels of cell cycle genes. Therefore, these results indicated that activated CD4 T cells of HIV-1 infected subjects were in cell cycle arrest. Taking together these results indicate that thymus function is preserved for long time during HIV- 1 infection and the increase observed in early stage disease may represent a compensatory mechanism to the depletion of CD4 T cells. However, we provide evidence for a cell cycle arrest of peripheral CD4 T cells that may prevent potentially the replenishment of CD4 T cells. RESUME Les mécanismes responsables de la perte des lymphocytes T CD4 lors de l'infection pas VIH n'ont été élucidés que partiellement. Nous ne savons toujours pas si l'épuisement des lymphocytes T CD4 résulte d'un défaut dans la production de cellules ou d'une destruction périphérique massive. Dans cette étude, la production de cellules T a été étudiée en mesurant les cercles d'excision générés lors du réarrangement du récepteur au cellules T (TRECs) chez des individus VIH-négatifs (N=120) et VIH-1 positifs (N=298). L'analyse des niveaux de TREC a été faite chez sujets HIV-infectés en considérant les phases de la maladie sur la base des comptes CD4 (phase précoce: > 500 cellules CD4/µl; intermédiaire: < 500>200; tardive: < 200) et par âge. Nos données démontrent que les niveaux de TRECs des cellules T CD8 étaient significativement plus bas chez les sujets VIH-1 infectés, à tous les stades de la maladie comparativement au groupe contrôle. En revanche, les niveaux de TRECs des cellules T CD4 étaient significativement plus élevés chez les sujets VIH-1 infectés durant la phase précoce de la maladie, tandis qu'aucune différence significative n'était observée durant la phase intermédiaire et étaient très réduits dans la phase tardive. Dans une deuxième partie, nous avons utilisé la technique des biopuces à d'ADN complémentaire pour analyser la régulation du cycle cellulaire chez les lymphocytes T CD4 et CD8 périphériques lors d'une infection au VIH-1. Des profils d'expression ont été déterminés et comparés à ceux de cellules T CD4 et CD8 stimulées in vitro, démontrant que les cellules T CD8 des sujets VIH-positifs avaient un profil d'expression très semblable à celui des cellules stimulées in vitro en prolifération. En revanche, les lymphocytes T CD4 des sujets VIH-1 positifs avaient un profil d'expression de gène plus complexe. En fait, leur profil montrait une sur- expression de gènes associés à une activation cellulaire, mais une sous-expression de ceux induisant une division. Ainsi, ces résultats indiquent que les lymphocytes T CD4 d'individus VIH-positifs présentent des dérégulations qui conduisent à un arrêt du cycle cellulaire. Ces résultats montrent que la fonction thymique est préservée longtemps pendant l'infection au VIH-1 et que l'augmentation de la quantité de TRECs dans la phase précoce de la maladie peut représenter un mécanisme compensatoire à l'épuisement des cellules T CD4. Cependant, nous démontrons aussi un clair dysfonctionnement du cycle cellulaire chez les cellules T CD4 d'individus infectés par VIH-1 ce qui peut enrayer la reconstitution du système immunitaire.

Robustness of Drosophila early embryo and wing imaginal disc development

Within a developing organism, cells require information on where they are in order to differentiate into the correct cell-type. Pattern formation is the process by which cells acquire and process positional cues and thus determine their fate. This can be achieved by the production and release of a diffusible signaling molecule, called a morphogen, which forms a concentration gradient: exposure to different morphogen levels leads to the activation of specific signaling pathways. Thus, in response to the morphogen gradient, cells start to express different sets of genes, forming domains characterized by a unique combination of differentially expressed genes. As a result, a pattern of cell fates and specification emerges.Though morphogens have been known for decades, it is not yet clear how these gradients form and are interpreted in order to yield highly robust patterns of gene expression. During my PhD thesis, I investigated the properties of Bicoid (Bcd) and Decapentaplegic (Dpp), two morphogens involved in the patterning of the anterior-posterior axis of Drosophila embryo and wing primordium, respectively. In particular, I have been interested in understanding how the pattern proportions are maintained across embryos of different sizes or within a growing tissue. This property is commonly referred to as scaling and is essential for yielding functional organs or organisms. In order to tackle these questions, I analysed fluorescence images showing the pattern of gene expression domains in the early embryo and wing imaginal disc. After characterizing the extent of these domains in a quantitative and systematic manner, I introduced and applied a new scaling measure in order to assess how well proportions are maintained. I found that scaling emerged as a universal property both in early embryos (at least far away from the Bcd source) and in wing imaginal discs (across different developmental stages). Since we were also interested in understanding the mechanisms underlying scaling and how it is transmitted from the morphogen to the target genes down in the signaling cascade, I also quantified scaling in mutant flies where this property could be disrupted. While scaling is largely conserved in embryos with altered bcd dosage, my modeling suggests that Bcd trapping by the nuclei as well as pre-steady state decoding of the morphogen gradient are essential to ensure precise and scaled patterning of the Bcd signaling cascade. In the wing imaginal disc, it appears that as the disc grows, the Dpp response expands and scales with the tissue size. Interestingly, scaling is not perfect at all positions in the field. The scaling of the target gene domains is best where they have a function; Spalt, for example, scales best at the position in the anterior compartment where it helps to form one of the anterior veins of the wing. Analysis of mutants for pentagone, a transcriptional target of Dpp that encodes a secreted feedback regulator of the pathway, indicates that Pentagone plays a key role in scaling the Dpp gradient activity.

Novel FOXF1 Mutations in Sporadic and Familial Cases of Alveolar Capillary Dysplasia with Misaligned Pulmonary Veins Imply a Role for its DNA Binding Domain.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACD/MPV) is a rare and lethal developmental disorder of the lung defined by a constellation of characteristic histopathological features. Nonpulmonary anomalies involving organs of gastrointestinal, cardiovascular, and genitourinary systems have been identified in approximately 80% of patients with ACD/MPV. We have collected DNA and pathological samples from more than 90 infants with ACD/MPV and their family members. Since the publication of our initial report of four point mutations and 10 deletions, we have identified an additional 38 novel nonsynonymous mutations of FOXF1 (nine nonsense, seven frameshift, one inframe deletion, 20 missense, and one no stop). This report represents an up to date list of all known FOXF1 mutations to the best of our knowledge. Majority of the cases are sporadic. We report four familial cases of which three show maternal inheritance, consistent with paternal imprinting of the gene. Twenty five mutations (60%) are located within the putative DNA-binding domain, indicating its plausible role in FOXF1 function. Five mutations map to the second exon. We identified two additional genic and eight genomic deletions upstream to FOXF1. These results corroborate and extend our previous observations and further establish involvement of FOXF1 in ACD/MPV and lung organogenesis.

Genome-wide prediction of matrix attachment regions that increase gene expression in mammalian cells.

Gene transfer in eukaryotic cells and organisms suffers from epigenetic effects that result in low or unstable transgene expression and high clonal variability. Use of epigenetic regulators such as matrix attachment regions (MARs) is a promising approach to alleviate such unwanted effects. Dissection of a known MAR allowed the identification of sequence motifs that mediate elevated transgene expression. Bioinformatics analysis implied that these motifs adopt a curved DNA structure that positions nucleosomes and binds specific transcription factors. From these observations, we computed putative MARs from the human genome. Cloning of several predicted MARs indicated that they are much more potent than the previously known element, boosting the expression of recombinant proteins from cultured cells as well as mediating high and sustained expression in mice. Thus we computationally identified potent epigenetic regulators, opening new strategies toward high and stable transgene expression for research, therapeutic production or gene-based therapies.

Genomic determinants of the efficiency of internal ribosomal entry sites of viral and cellular origin

Variation in cellular gene expression levels has been shown to be inherited. Expression is controlled at transcriptional and post-transcriptional levels. Internal ribosome entry sites (IRES) are used by viruses to bypass inhibition of cap-dependent translation, and by eukaryotic cells to control translation under conditions when protein synthesis is inhibited. We aimed at identifying genomic determinants of variability in IRES-mediated translation of viral [Encephalomyocarditis virus (EMCV)] and cellular IRES [X-linked inhibitor-of-apoptosis (XIAP) and c-myc]. Bicistronic lentiviral constructs expressing two fluorescent reporters were used to transduce laboratory and B lymphoblastoid cell lines [15 CEPH pedigrees (n = 205) and 50 unrelated individuals]. IRES efficiency varied according to cell type and among individuals. Control of IRES activity has a significant genetic component (h(2) of 0.47 and 0.36 for EMCV and XIAP, respectively). Quantitative linkage analysis identified a suggestive locus (LOD 2.35) on chromosome 18q21.2, and genome-wide association analysis revealed of a cluster of SNPs on chromosome 3, intronic to the FHIT gene, marginally associated (P = 5.9E-7) with XIAP IRES function. This study illustrates the in vitro generation of intermediate phenotypes by using cell lines for the evaluation of genetic determinants of control of elements such as IRES.

Global transcriptional programs in peripheral nerve endoneurium and DRG are resistant to the onset of type 1 diabetic neuropathy in Ins2 mice.

While the morphological and electrophysiological changes underlying diabetic peripheral neuropathy (DPN) are relatively well described, the involved molecular mechanisms remain poorly understood. In this study, we investigated whether phenotypic changes associated with early DPN are correlated with transcriptional alterations in the neuronal (dorsal root ganglia [DRG]) or the glial (endoneurium) compartments of the peripheral nerve. We used Ins2(Akita/+) mice to study transcriptional changes underlying the onset of DPN in type 1 diabetes mellitus (DM). Weight, blood glucose and motor nerve conduction velocity (MNCV) were measured in Ins2(Akita/+) and control mice during the first three months of life in order to determine the onset of DPN. Based on this phenotypic characterization, we performed gene expression profiling using sciatic nerve endoneurium and DRG isolated from pre-symptomatic and early symptomatic Ins2(Akita/+) mice and sex-matched littermate controls. Our phenotypic analysis of Ins2(Akita/+) mice revealed that DPN, as measured by reduced MNCV, is detectable in affected animals already one week after the onset of hyperglycemia. Surprisingly, the onset of DPN was not associated with any major persistent changes in gene expression profiles in either sciatic nerve endoneurium or DRG. Our data thus demonstrated that the transcriptional programs in both endoneurial and neuronal compartments of the peripheral nerve are relatively resistant to the onset of hyperglycemia and hypoinsulinemia suggesting that either minor transcriptional alterations or changes on the proteomic level are responsible for the functional deficits associated with the onset of DPN in type 1 DM.

KCTD13 is a major driver of mirrored neuroanatomical phenotypes of the 16p11.2 copy number variant.

Copy number variants (CNVs) are major contributors to genetic disorders. We have dissected a region of the 16p11.2 chromosome--which encompasses 29 genes--that confers susceptibility to neurocognitive defects when deleted or duplicated. Overexpression of each human transcript in zebrafish embryos identified KCTD13 as the sole message capable of inducing the microcephaly phenotype associated with the 16p11.2 duplication, whereas suppression of the same locus yielded the macrocephalic phenotype associated with the 16p11.2 deletion, capturing the mirror phenotypes of humans. Analyses of zebrafish and mouse embryos suggest that microcephaly is caused by decreased proliferation of neuronal progenitors with concomitant increase in apoptosis in the developing brain, whereas macrocephaly arises by increased proliferation and no changes in apoptosis. A role for KCTD13 dosage changes is consistent with autism in both a recently reported family with a reduced 16p11.2 deletion and a subject reported here with a complex 16p11.2 rearrangement involving de novo structural alteration of KCTD13. Our data suggest that KCTD13 is a major driver for the neurodevelopmental phenotypes associated with the 16p11.2 CNV, reinforce the idea that one or a small number of transcripts within a CNV can underpin clinical phenotypes, and offer an efficient route to identifying dosage-sensitive loci.

A single epidermal stem cell strategy for safe ex vivo gene therapy.

There is a widespread agreement from patient and professional organisations alike that the safety of stem cell therapeutics is of paramount importance, particularly for ex vivo autologous gene therapy. Yet current technology makes it difficult to thoroughly evaluate the behaviour of genetically corrected stem cells before they are transplanted. To address this, we have developed a strategy that permits transplantation of a clonal population of genetically corrected autologous stem cells that meet stringent selection criteria and the principle of precaution. As a proof of concept, we have stably transduced epidermal stem cells (holoclones) obtained from a patient suffering from recessive dystrophic epidermolysis bullosa. Holoclones were infected with self-inactivating retroviruses bearing a COL7A1 cDNA and cloned before the progeny of individual stem cells were characterised using a number of criteria. Clonal analysis revealed a great deal of heterogeneity among transduced stem cells in their capacity to produce functional type VII collagen (COLVII). Selected transduced stem cells transplanted onto immunodeficient mice regenerated a non-blistering epidermis for months and produced a functional COLVII. Safety was assessed by determining the sites of proviral integration, rearrangements and hit genes and by whole-genome sequencing. The progeny of the selected stem cells also had a diploid karyotype, was not tumorigenic and did not disseminate after long-term transplantation onto immunodeficient mice. In conclusion, a clonal strategy is a powerful and efficient means of by-passing the heterogeneity of a transduced stem cell population. It guarantees a safe and homogenous medicinal product, fulfilling the principle of precaution and the requirements of regulatory affairs. Furthermore, a clonal strategy makes it possible to envision exciting gene-editing technologies like zinc finger nucleases, TALENs and homologous recombination for next-generation gene therapy.

Genome-wide profiling of the cardiac transcriptome after myocardial infarction identifies novel heart-specific long non-coding RNAs.

AIM: Heart disease is recognized as a consequence of dysregulation of cardiac gene regulatory networks. Previously, unappreciated components of such networks are the long non-coding RNAs (lncRNAs). Their roles in the heart remain to be elucidated. Thus, this study aimed to systematically characterize the cardiac long non-coding transcriptome post-myocardial infarction and to elucidate their potential roles in cardiac homoeostasis. METHODS AND RESULTS: We annotated the mouse transcriptome after myocardial infarction via RNA sequencing and ab initio transcript reconstruction, and integrated genome-wide approaches to associate specific lncRNAs with developmental processes and physiological parameters. Expression of specific lncRNAs strongly correlated with defined parameters of cardiac dimensions and function. Using chromatin maps to infer lncRNA function, we identified many with potential roles in cardiogenesis and pathological remodelling. The vast majority was associated with active cardiac-specific enhancers. Importantly, oligonucleotide-mediated knockdown implicated novel lncRNAs in controlling expression of key regulatory proteins involved in cardiogenesis. Finally, we identified hundreds of human orthologues and demonstrate that particular candidates were differentially modulated in human heart disease. CONCLUSION: These findings reveal hundreds of novel heart-specific lncRNAs with unique regulatory and functional characteristics relevant to maladaptive remodelling, cardiac function and possibly cardiac regeneration. This new class of molecules represents potential therapeutic targets for cardiac disease. Furthermore, their exquisite correlation with cardiac physiology renders them attractive candidate biomarkers to be used in the clinic.

A Potential Contributory Role for Ciliary Dysfunction in the 16p11.2 600 kb BP4-BP5 Pathology.

The 16p11.2 600 kb copy-number variants (CNVs) are associated with mirror phenotypes on BMI, head circumference, and brain volume and represent frequent genetic lesions in autism spectrum disorders (ASDs) and schizophrenia. Here we interrogated the transcriptome of individuals carrying reciprocal 16p11.2 CNVs. Transcript perturbations correlated with clinical endophenotypes and were enriched for genes associated with ASDs, abnormalities of head size, and ciliopathies. Ciliary gene expression was also perturbed in orthologous mouse models, raising the possibility that ciliary dysfunction contributes to 16p11.2 pathologies. In support of this hypothesis, we found structural ciliary defects in the CA1 hippocampal region of 16p11.2 duplication mice. Moreover, by using an established zebrafish model, we show genetic interaction between KCTD13, a key driver of the mirrored neuroanatomical phenotypes of the 16p11.2 CNV, and ciliopathy-associated genes. Overexpression of BBS7 rescues head size and neuroanatomical defects of kctd13 morphants, whereas suppression or overexpression of CEP290 rescues phenotypes induced by KCTD13 under- or overexpression, respectively. Our data suggest that dysregulation of ciliopathy genes contributes to the clinical phenotypes of these CNVs.

Copy number variations and cognitive phenotypes in unselected populations

IMPORTANCE: The association of copy number variations (CNVs), differing numbers of copies of genetic sequence at locations in the genome, with phenotypes such as intellectual disability has been almost exclusively evaluated using clinically ascertained cohorts. The contribution of these genetic variants to cognitive phenotypes in the general population remains unclear. OBJECTIVE: To investigate the clinical features conferred by CNVs associated with known syndromes in adult carriers without clinical preselection and to assess the genome-wide consequences of rare CNVs (frequency ≤0.05%; size ≥250 kilobase pairs [kb]) on carriers' educational attainment and intellectual disability prevalence in the general population. DESIGN, SETTING, AND PARTICIPANTS: The population biobank of Estonia contains 52,000 participants enrolled from 2002 through 2010. General practitioners examined participants and filled out a questionnaire of health- and lifestyle-related questions, as well as reported diagnoses. Copy number variant analysis was conducted on a random sample of 7877 individuals and genotype-phenotype associations with education and disease traits were evaluated. Our results were replicated on a high-functioning group of 993 Estonians and 3 geographically distinct populations in the United Kingdom, the United States, and Italy. MAIN OUTCOMES AND MEASURES: Phenotypes of genomic disorders in the general population, prevalence of autosomal CNVs, and association of these variants with educational attainment (from less than primary school through scientific degree) and prevalence of intellectual disability. RESULTS: Of the 7877 in the Estonian cohort, we identified 56 carriers of CNVs associated with known syndromes. Their phenotypes, including cognitive and psychiatric problems, epilepsy, neuropathies, obesity, and congenital malformations are similar to those described for carriers of identical rearrangements ascertained in clinical cohorts. A genome-wide evaluation of rare autosomal CNVs (frequency, ≤0.05%; ≥250 kb) identified 831 carriers (10.5%) of the screened general population. Eleven of 216 (5.1%) carriers of a deletion of at least 250 kb (odds ratio [OR], 3.16; 95% CI, 1.51-5.98; P = 1.5e-03) and 6 of 102 (5.9%) carriers of a duplication of at least 1 Mb (OR, 3.67; 95% CI, 1.29-8.54; P = .008) had an intellectual disability compared with 114 of 6819 (1.7%) in the Estonian cohort. The mean education attainment was 3.81 (P = 1.06e-04) among 248 (≥250 kb) deletion carriers and 3.69 (P = 5.024e-05) among 115 duplication carriers (≥1 Mb). Of the deletion carriers, 33.5% did not graduate from high school (OR, 1.48; 95% CI, 1.12-1.95; P = .005) and 39.1% of duplication carriers did not graduate high school (OR, 1.89; 95% CI, 1.27-2.8; P = 1.6e-03). Evidence for an association between rare CNVs and lower educational attainment was supported by analyses of cohorts of adults from Italy and the United States and adolescents from the United Kingdom. CONCLUSIONS AND RELEVANCE: Known pathogenic CNVs in unselected, but assumed to be healthy, adult populations may be associated with unrecognized clinical sequelae. Additionally, individually rare but collectively common intermediate-size CNVs may be negatively associated with educational attainment. Replication of these findings in additional population groups is warranted given the potential implications of this observation for genomics research, clinical care, and public health.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHRadjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHRadjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHRadjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHRadjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHRadjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHRadjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.