Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus.

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

Genome-wide association study in almost 195,000 individuals identifies 50 previously unidentified genetic loci for eye color.

Human eye color is highly heritable, but its genetic architecture is not yet fully understood. We report the results of the largest genome-wide association study for eye color to date, involving up to 192,986 European participants from 10 populations. We identify 124 independent associations arising from 61 discrete genomic regions, including 50 previously unidentified. We find evidence for genes involved in melanin pigmentation, but we also find associations with genes involved in iris morphology and structure. Further analyses in 1636 Asian participants from two populations suggest that iris pigmentation variation in Asians is genetically similar to Europeans, albeit with smaller effect sizes. Our findings collectively explain 53.2% (95% confidence interval, 45.4 to 61.0%) of eye color variation using common single-nucleotide polymorphisms. Overall, our study outcomes demonstrate that the genetic complexity of human eye color considerably exceeds previous knowledge and expectations, highlighting eye color as a genetically highly complex human trait.

F cell numbers are associated with an X-linked genetic polymorphism and correlate with haematological parameters in patients with sickle cell disease.

Patients with sickle cell disease (SCD) with high fetal haemoglobin (HbF) tend to have a lower incidence of complications and longer survival due to inhibition of deoxyhaemoglobin S (HbS) polymerisation by HbF. HbF-containing cells, namely F cells, are strongly influenced by genetic factors. We measured the percentage of F cells (Fcells%) in 222 patients with SCD to evaluate the association of (i) Fcells% with genetic HbF-modifier variants and (ii) Fcells% with haematological parameters. There was a different distribution of Fcells% in females compared to males. The association of the B-cell lymphoma/leukaemia 11A (BCL11A) locus with Fcells% (ß = 8·238; P < 0·001) and with HbF% (ß = 2·490; P < 0·001) was significant. All red cell parameters except for Hb and mean corpuscular Hb concentration levels in males and females were significantly different. The Fcells% was positively associated with mean cell Hb, mean cell volume and reticulocytes. To explain the significant gender difference in Fcells%, we tested for associations with single nucleotide polymorphisms on the X chromosomal region Xp22.2, where a genetic determinant of HbF had been previously hypothesised. We found in males a significant association with a SNP in FERM and PDZ domain-containing protein 4 (FRMPD4) and adjacent to male-specific lethal complex subunit 3 (MSL3). Thus, we have identified an X-linked locus that could account for a significant fraction of the Fcells% variation in patients with SCD.

A bird's-eye view of Italian genomic variation through whole-genome sequencing.

The genomic variation of the Italian peninsula populations is currently under characterised: the only Italian whole-genome reference is represented by the Tuscans from the 1000 Genome Project. To address this issue, we sequenced a total of 947 Italian samples from three different geographical areas. First, we defined a new Italian Genome Reference Panel (IGRP1.0) for imputation, which improved imputation accuracy, especially for rare variants, and we tested it by GWAS analysis on red blood traits. Furthermore, we extended the catalogue of genetic variation investigating the level of population structure, the pattern of natural selection, the distribution of deleterious variants and occurrence of human knockouts (HKOs). Overall the results demonstrate a high level of genomic differentiation between cohorts, different signatures of natural selection and a distinctive distribution of deleterious variants and HKOs, confirming the necessity of distinct genome references for the Italian population.

Associations of autozygosity with a broad range of human phenotypes.

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

Target genes, variants, tissues and transcriptional pathways influencing human serum urate levels.

Elevated serum urate levels cause gout and correlate with cardiometabolic diseases via poorly understood mechanisms. We performed a trans-ancestry genome-wide association study of serum urate in 457,690 individuals, identifying 183 loci (147 previously unknown) that improve the prediction of gout in an independent cohort of 334,880 individuals. Serum urate showed significant genetic correlations with many cardiometabolic traits, with genetic causality analyses supporting a substantial role for pleiotropy. Enrichment analysis, fine-mapping of urate-associated loci and colocalization with gene expression in 47 tissues implicated the kidney and liver as the main target organs and prioritized potentially causal genes and variants, including the transcriptional master regulators in the liver and kidney, HNF1A and HNF4A. Experimental validation showed that HNF4A transactivated the promoter of ABCG2, encoding a major urate transporter, in kidney cells, and that HNF4A p.Thr139Ile is a functional variant. Transcriptional coregulation within and across organs may be a general mechanism underlying the observed pleiotropy between urate and cardiometabolic traits.

Genome-wide association meta-analysis identifies five novel loci for age-related hearing impairment.

Previous research has shown that genes play a substantial role in determining a person's susceptibility to age-related hearing impairment. The existing studies on this subject have different results, which may be caused by difficulties in determining the phenotype or the limited number of participants involved. Here, we have gathered the largest sample to date (discovery n = 9,675; replication n = 10,963; validation n = 356,141), and examined phenotypes that represented low/mid and high frequency hearing loss on the pure tone audiogram. We identified 7 loci that were either replicated and/or validated, of which 5 loci are novel in hearing. Especially the ILDR1 gene is a high profile candidate, as it contains our top SNP, is a known hearing loss gene, has been linked to age-related hearing impairment before, and in addition is preferentially expressed within hair cells of the inner ear. By verifying all previously published SNPs, we can present a paper that combines all new and existing findings to date, giving a complete overview of the genetic architecture of age-related hearing impairment. This is of importance as age-related hearing impairment is highly prevalent in our ageing society and represents a large socio-economic burden.

Publisher Correction: Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A catalog of genetic loci associated with kidney function from analyses of a million individuals.

Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through trans-ancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these, 147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.

Effects of Calcium, Magnesium, and Potassium Concentrations on Ventricular Repolarization in Unselected Individuals.

BACKGROUND: Subclinical changes on the electrocardiogram are risk factors for cardiovascular mortality. Recognition and knowledge of electrolyte associations in cardiac electrophysiology are based on only in vitro models and observations in patients with severe medical conditions. OBJECTIVES: This study sought to investigate associations between serum electrolyte concentrations and changes in cardiac electrophysiology in the general population. METHODS: Summary results collected from 153,014 individuals (54.4% women; mean age 55.1 ± 12.1 years) from 33 studies (of 5 ancestries) were meta-analyzed. Linear regression analyses examining associations between electrolyte concentrations (mmol/l of calcium, potassium, sodium, and magnesium), and electrocardiographic intervals (RR, QT, QRS, JT, and PR intervals) were performed. The study adjusted for potential confounders and also stratified by ancestry, sex, and use of antihypertensive drugs. RESULTS: Lower calcium was associated with longer QT intervals (-11.5 ms; 99.75% confidence interval [CI]: -13.7 to -9.3) and JT duration, with sex-specific effects. In contrast, higher magnesium was associated with longer QT intervals (7.2 ms; 99.75% CI: 1.3 to 13.1) and JT. Lower potassium was associated with longer QT intervals (-2.8 ms; 99.75% CI: -3.5 to -2.0), JT, QRS, and PR durations, but all potassium associations were driven by use of antihypertensive drugs. No physiologically relevant associations were observed for sodium or RR intervals. CONCLUSIONS: The study identified physiologically relevant associations between electrolytes and electrocardiographic intervals in a large-scale analysis combining cohorts from different settings. The results provide insights for further cardiac electrophysiology research and could potentially influence clinical practice, especially the association between calcium and QT duration, by which calcium levels at the bottom 2% of the population distribution led to clinically relevant QT prolongation by >5 ms.

Next Generation Sequencing and Animal Models Reveal SLC9A3R1 as a New Gene Involved in Human Age-Related Hearing Loss.

Age-related hearing loss (ARHL) is the most common sensory impairment in the elderly affecting millions of people worldwide. To shed light on the genetics of ARHL, a large cohort of 464 Italian patients has been deeply characterized at clinical and molecular level. In particular, 46 candidate genes, selected on the basis of genome-wide association studies (GWAS), animal models and literature updates, were analyzed by targeted re-sequencing. After filtering and prioritization steps, SLC9A3R1 has been identified as a strong candidate and then validated by "in vitro" and "in vivo" studies. Briefly, a rare (MAF: 2.886e-5) missense variant c.539G > A, p.(R180Q) was detected in two unrelated male patients affected by ARHL characterized by a severe to profound high-frequency hearing loss. The variant, predicted as damaging, was not present in healthy matched controls. Protein modeling confirmed the pathogenic effect of p.(R180Q) variant on protein's structure leading to a change in the total number of hydrogen bonds. In situ hybridization showed slc9a3r1 expression in zebrafish inner ear. A zebrafish knock-in model, generated by CRISPR-Cas9 technology, revealed a reduced auditory response at all frequencies in slc9a3r1 R180Q/R180Q mutants compared to slc9a3r1 +/+ and slc9a3r1 +/R180Q animals. Moreover, a significant reduction (5.8%) in the total volume of the saccular otolith (which is responsible for sound detection) was observed in slc9a3r1 R180Q/R180Q compared to slc9a3r1 +/+ (P = 0.0014), while the utricular otolith, necessary for balance, was not affected in agreement with the human phenotype. Overall, these data strongly support the role of SLC9A3R1 gene in the pathogenesis of ARHL opening new perspectives in terms of diagnosis, prevention and treatment.

Cx26 partial loss causes accelerated presbycusis by redox imbalance and dysregulation of Nfr2 pathway.

Mutations in GJB2, the gene that encodes connexin 26 (Cx26), are the most common cause of sensorineural hearing impairment. The truncating variant 35delG, which determines a complete loss of Cx26 protein function, is the prevalent GJB2 mutation in several populations. Here, we generated and analyzed Gjb2+/- mice as a model of heterozygous human carriers of 35delG. Compared to control mice, auditory brainstem responses (ABRs) and distortion product otoacoustic emissions (DPOAEs) worsened over time more rapidly in Gjb2+/- mice, indicating they were affected by accelerated age-related hearing loss (ARHL), or presbycusis. We linked causally the auditory phenotype of Gjb2+/- mice to apoptosis and oxidative damage in the cochlear duct, reduced release of glutathione from connexin hemichannels, decreased nutrient delivery to the sensory epithelium via cochlear gap junctions and deregulated expression of genes that are under transcriptional control of the nuclear factor erythroid 2-related factor 2 (Nrf2), a pivotal regulator of tolerance to redox stress. Moreover, a statistically significant genome-wide association with two genes (PRKCE and TGFB1) related to the Nrf2 pathway (p-value < 4â€¯× 10-2) was detected in a very large cohort of 4091 individuals, originating from Europe, Caucasus and Central Asia, with hearing phenotype (including 1076 presbycusis patients and 1290 healthy matched controls). We conclude that (i) elements of the Nrf2 pathway are essential for hearing maintenance and (ii) their dysfunction may play an important role in the etiopathogenesis of human presbycusis.

Whole-genome sequencing reveals new insights into age-related hearing loss: cumulative effects, pleiotropy and the role of selection.

Age-related hearing loss (ARHL) is the most common sensory disorder in the elderly. Although not directly life threatening, it contributes to loss of autonomy and is associated with anxiety, depression and cognitive decline. To search for genetic risk factors underlying ARHL, a large whole-genome sequencing (WGS) approach has been carried out in a cohort of 212 cases and controls, both older than 50 years to select genes characterized by a burden of variants specific to cases or controls. Accordingly, the total variation load per gene was compared and two groups were detected: 375 genes more variable in cases and 371 more variable in controls. In both cases, Gene Ontology analysis showed that the largest enrichment for biological processes (fold > 5, p-value = 0.042) was the "sensory perception of sound", suggesting cumulative genetic effects were involved. Replication confirmed 141 genes, while additional analysis based on natural selection led to a prioritization of 21 genes. The majority of them (20 out of 21) showed positive expression in mouse cochlea cDNA and were associated with two functional pathways. Among them, two genes were previously associated with hearing (CSMD1 and PTRPD) and re-sequenced in a large Italian cohort of ARHL patients (N = 389). Results led to the identification of six coding variants not detected in cases so far, suggesting a possible protective role, which requires investigation. In conclusion, we show that this multistep strategy (WGS, selection, expression, pathway analysis and targeted re-sequencing) can provide major insights into the molecular characterization of complex diseases such as ARHL.

Genome-wide association meta-analysis of individuals of European ancestry identifies new loci explaining a substantial fraction of hair color variation and heritability.

Hair color is one of the most recognizable visual traits in European populations and is under strong genetic control. Here we report the results of a genome-wide association study meta-analysis of almost 300,000 participants of European descent. We identified 123 autosomal and one X-chromosome loci significantly associated with hair color; all but 13 are novel. Collectively, single-nucleotide polymorphisms associated with hair color within these loci explain 34.6% of red hair, 24.8% of blond hair, and 26.1% of black hair heritability in the study populations. These results confirm the polygenic nature of complex phenotypes and improve our understanding of melanin pigment metabolism in humans.

Genomic Studies in a Large Cohort of Hearing Impaired Italian Patients Revealed Several New Alleles, a Rare Case of Uniparental Disomy (UPD) and the Importance to Search for Copy Number Variations.

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.