Genome-wide association studies (GWAS) have identified >250 loci for body mass index (BMI), implicating pathways related to neuronal biology. Most GWAS loci represent clusters of common, noncoding variants from which pinpointing causal genes remains challenging. Here we combined data from 718,734 individuals to discover rare and low-frequency (minor allele frequency (MAF) < 5%) coding variants associated with BMI. We identified 14 coding variants in 13 genes, of which 8 variants were in genes (ZBTB7B, ACHE, RAPGEF3, RAB21, ZFHX3, ENTPD6, ZFR2 and ZNF169) newly implicated in human obesity, 2 variants were in genes (MC4R and KSR2) previously observed to be mutated in extreme obesity and 2 variants were in GIPR. The effect sizes of rare variants are ~10 times larger than those of common variants, with the largest effect observed in carriers of an MC4R mutation introducing a stop codon (p.Tyr35Ter, MAF = 0.01%), who weighed ~7 kg more than non-carriers. Pathway analyses based on the variants associated with BMI confirm enrichment of neuronal genes and provide new evidence for adipocyte and energy expenditure biology, widening the potential of genetically supported therapeutic targets in obesity.

Melinda Mills, Nicola Barban, Harold Snieder, Marcel den Hoed and colleagues perform a meta-analysis of data from over 300,000 individuals for age at first birth and number of children ever born. They identify 12 significant loci that associate with these traits, providing insights into the genetic basis of human reproductive behavior. The genetic architecture of human reproductive behavior—age at first birth (AFB) and number of children ever born (NEB)—has a strong relationship with fitness, human development, infertility and risk of neuropsychiatric disorders. However, very few genetic loci have been identified, and the underlying mechanisms of AFB and NEB are poorly understood. We report a large genome-wide association study of both sexes including 251,151 individuals for AFB and 343,072 individuals for NEB. We identified 12 independent loci that are significantly associated with AFB and/or NEB in a SNP-based genome-wide association study and 4 additional loci associated in a gene-based effort. These loci harbor genes that are likely to have a role, either directly or by affecting non-local gene expression, in human reproduction and infertility, thereby increasing understanding of these complex traits.

An analysis of 16 health-related quantitative traits in approximately 350,000 individuals reveals statistically significant associations between genome-wide homozygosity and four complex traits (height, lung function, cognitive ability and educational attainment); in each case increased homozygosity associates with a decreased trait value, but no evidence was seen of an influence on blood pressure, cholesterol, or ten other cardio-metabolic traits. This consortium meta-analysis of 102 cohorts and more than 350,000 individuals investigates the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity, ROH), which are inferred to be homozygous along their complete length. Focusing on 16 health-related quantitative traits, the authors find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in a second, general cognitive ability and educational attainment. In each case increased homozygosity associates with decreased trait value. No evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Homozygosity has long been associated with rare, often devastating, Mendelian disorders1, and Darwin was one of the first to recognize that inbreeding reduces evolutionary fitness2. However, the effect of the more distant parental relatedness that is common in modern human populations is less well understood. Genomic data now allow us to investigate the effects of homozygosity on traits of public health importance by observing contiguous homozygous segments (runs of homozygosity), which are inferred to be homozygous along their complete length. Given the low levels of genome-wide homozygosity prevalent in most human populations, information is required on very large numbers of people to provide sufficient power3,4. Here we use runs of homozygosity to study 16 health-related quantitative traits in 354,224 individuals from 102 cohorts, and find statistically significant associations between summed runs of homozygosity and four complex traits: height, forced expiratory lung volume in one second, general cognitive ability and educational attainment (P < 1 × 10−300, 2.1 × 10−6, 2.5 × 10−10 and 1.8 × 10−10, respectively). In each case, increased homozygosity was associated with decreased trait value, equivalent to the offspring of first cousins being 1.2 cm shorter and having 10 months’ less education. Similar effect sizes were found across four continental groups and populations with different degrees of genome-wide homozygosity, providing evidence that homozygosity, rather than confounding, directly contributes to phenotypic variance. Contrary to earlier reports in substantially smaller samples5,6, no evidence was seen of an influence of genome-wide homozygosity on blood pressure and low density lipoprotein cholesterol, or ten other cardio-metabolic traits. Since directional dominance is predicted for traits under directional evolutionary selection7, this study provides evidence that increased stature and cognitive function have been positively selected in human evolution, whereas many important risk factors for late-onset complex diseases may not have been.

Genome-wide association studies conducted on QRS duration, an electrocardiographic measurement associated with heart failure and sudden cardiac death, have led to novel biological insights into cardiac function. However, the variants identified fall predominantly in non-coding regions and their underlying mechanisms remain unclear. Here, we identify putative functional coding variation associated with changes in the QRS interval duration by combining Illumina HumanExome BeadChip genotype data from 77,898 participants of European ancestry and 7695 of African descent in our discovery cohort, followed by replication in 111,874 individuals of European ancestry from the UK Biobank and deCODE cohorts. We identify ten novel loci, seven within coding regions, including ADAMTS6, significantly associated with QRS duration in gene-based analyses. ADAMTS6 encodes a secreted metalloprotease of currently unknown function. In vitro validation analysis shows that the QRS-associated variants lead to impaired ADAMTS6 secretion and loss-of function analysis in mice demonstrates a previously unappreciated role for ADAMTS6 in connexin 43 gap junction expression, which is essential for myocardial conduction. Our approach identifies novel coding and non-coding variants underlying ventricular depolarization and provides a possible mechanism for the ADAMTS6-associated conduction changes.

Hereditary hearing loss (HHL) is a common disorder characterized by a huge genetic heterogeneity. The definition of a correct molecular diagnosis is essential for proper genetic counseling, recurrence risk estimation, and therapeutic options. From 20 to 40% of patients carry mutations in GJB2 gene, thus, in more than half of cases it is necessary to look for causative variants in the other genes so far identified (~100). In this light, the use of next-generation sequencing technologies has proved to be the best solution for mutational screening, even though it is not always conclusive. Here we describe a combined approach, based on targeted re-sequencing (TRS) of 96 HHL genes followed by high-density SNP arrays, aimed at the identification of the molecular causes of non-syndromic HHL (NSHL). This strategy has been applied to study 103 Italian unrelated cases, negative for mutations in GJB2, and led to the characterization of 31% of them (i.e., 37% of familial and 26.3% of sporadic cases). In particular, TRS revealed TECTA and ACTG1 genes as major players in the Italian population. Furthermore, two de novo missense variants in ACTG1 have been identified and investigated through protein modeling and molecular dynamics simulations, confirming their likely pathogenic effect. Among the selected patients analyzed by SNP arrays (negative to TRS, or with a single variant in a recessive gene) a molecular diagnosis was reached in ~36% of cases, highlighting the importance to look for large insertions/deletions. Moreover, copy number variants analysis led to the identification of the first case of uniparental disomy involving LOXHD1 gene. Overall, taking into account the contribution of GJB2, plus the results from TRS and SNP arrays, it was possible to reach a molecular diagnosis in ~51% of NSHL cases. These data proved the usefulness of a combined approach for the analysis of NSHL and for the definition of the epidemiological picture of HHL in the Italian population.