Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin.

BACKGROUND: Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown. METHODS: We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing. RESULTS: Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion. CONCLUSIONS: Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche.

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.

Genetic determinants of circulating VEGF levels in major depressive disorder and electroconvulsive therapy response.

Alterations in peripheral vascular endothelial growth factor (VEGF) levels were observed in major depressive disorder and relative treatments and were shown to be influenced by genetic variants. The study objective was to explore, at a genome-wide level, possible interplaying effects between the genetic background and major depressive disorder in regulating VEGF levels. Moreover, we aimed to investigate the association between these variants and response to electroconvulsive therapy. A genome-wide association study was carried out both on controls and patients with major depressive disorder (n = 145; n = 121) in correlation with serum VEGF levels determined by ELISA. Five SNPs not included in SNP arrays were additionally genotyped. Seventy-one patients with treatment-resistant depression underwent electroconvulsive therapy and were evaluated as responders/nonresponders. An association between VEGF levels and a locus in 6p21.1, downstream the VEGF gene, was evidenced both in controls (best SNP: FDR-corrected p = 2.4 × 10-5 ) and in patients with major depressive disorder (best SNP: FDR-corrected p = 2.6 × 10-3 ). The alleles associated with lower VEGF concentrations in patients were also associated with nonresponse to electroconvulsive therapy (p = .01). These results confirm a role of SNPs in 6p21.1 locus as major influencers of circulating VEGF levels also in patients affected by major depressive disorder and indicate a possible implication in response to electroconvulsive therapy.

Genome-Wide Meta-Analysis Unravels Interactions between Magnesium Homeostasis and Metabolic Phenotypes.

Magnesium (Mg2+) homeostasis is critical for metabolism. However, the genetic determinants of the renal handling of Mg2+, which is crucial for Mg2+ homeostasis, and the potential influence on metabolic traits in the general population are unknown. We obtained plasma and urine parameters from 9099 individuals from seven cohorts, and conducted a genome-wide meta-analysis of Mg2+ homeostasis. We identified two loci associated with urinary magnesium (uMg), rs3824347 (P=4.4×10-13) near TRPM6, which encodes an epithelial Mg2+ channel, and rs35929 (P=2.1×10-11), a variant of ARL15, which encodes a GTP-binding protein. Together, these loci account for 2.3% of the variation in 24-hour uMg excretion. In human kidney cells, ARL15 regulated TRPM6-mediated currents. In zebrafish, dietary Mg2+ regulated the expression of the highly conserved ARL15 ortholog arl15b, and arl15b knockdown resulted in renal Mg2+ wasting and metabolic disturbances. Finally, ARL15 rs35929 modified the association of uMg with fasting insulin and fat mass in a general population. In conclusion, this combined observational and experimental approach uncovered a gene-environment interaction linking Mg2+ deficiency to insulin resistance and obesity.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

Common variants in CLDN14 are associated with differential excretion of magnesium over calcium in urine.

The nature and importance of genetic factors regulating the differential handling of Ca2+ and Mg2+ by the renal tubule in the general population are poorly defined. We conducted a genome-wide meta-analysis of urinary magnesium-to-calcium ratio to identify associated common genetic variants. We included 9320 adults of European descent from four genetic isolates and three urban cohorts. Urinary magnesium and calcium concentrations were measured centrally in spot urine, and each study conducted linear regression analysis of urinary magnesium-to-calcium ratio on ~2.5 million single-nucleotide polymorphisms (SNPs) using an additive model. We investigated, in mouse, the renal expression profile of the top candidate gene and its variation upon changes in dietary magnesium. The genome-wide analysis evidenced a top locus (rs172639, p = 1.7 × 10-12), encompassing CLDN14, the gene coding for claudin-14, that was genome-wide significant when using urinary magnesium-to-calcium ratio, but not either one taken separately. In mouse, claudin-14 is expressed in the distal nephron segments specifically handling magnesium, and its expression is regulated by chronic changes in dietary magnesium content. A genome-wide approach identified common variants in the CLDN14 gene exerting a robust influence on the differential excretion of Mg2+ over Ca2+ in urine. These data highlight the power of urinary electrolyte ratios to unravel genetic determinants of renal tubular function. Coupled with mouse experiments, these results support a major role for claudin-14, a gene associated with kidney stones, in the differential paracellular handling of divalent cations by the renal tubule.

A general approach for haplotype phasing across the full spectrum of relatedness.

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations.

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤ 21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.

Evidence of inbreeding depression on human height.

Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ(2) = 83.89, df = 1; p = 5.2 × 10(-20)). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.

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.

Common variants in UMOD associate with urinary uromodulin levels: a meta-analysis.

Uromodulin is expressed exclusively in the thick ascending limb and is the most abundant protein excreted in normal urine. Variants in UMOD, which encodes uromodulin, are associated with renal function, and urinary uromodulin levels may be a biomarker for kidney disease. However, the genetic factors regulating uromodulin excretion are unknown. We conducted a meta-analysis of urinary uromodulin levels to identify associated common genetic variants in the general population. We included 10,884 individuals of European descent from three genetic isolates and three urban cohorts. Each study measured uromodulin indexed to creatinine and conducted linear regression analysis of approximately 2.5 million single nucleotide polymorphisms using an additive model. We also tested whether variants in genes expressed in the thick ascending limb associate with uromodulin levels. rs12917707, located near UMOD and previously associated with renal function and CKD, had the strongest association with urinary uromodulin levels (P<0.001). In all cohorts, carriers of a G allele of this variant had higher uromodulin levels than noncarriers did (geometric means 10.24, 14.05, and 17.67 µg/g creatinine for zero, one, or two copies of the G allele). rs12446492 in the adjacent gene PDILT (protein disulfide isomerase-like, testis expressed) also reached genome-wide significance (P<0.001). Regarding genes expressed in the thick ascending limb, variants in KCNJ1, SORL1, and CAB39 associated with urinary uromodulin levels. These data indicate that common variants in the UMOD promoter region may influence urinary uromodulin levels. They also provide insights into uromodulin biology and the association of UMOD variants with renal function.

Integration of genome-wide association studies with biological knowledge identifies six novel genes related to kidney function.

In conducting genome-wide association studies (GWAS), analytical approaches leveraging biological information may further understanding of the pathophysiology of clinical traits. To discover novel associations with estimated glomerular filtration rate (eGFR), a measure of kidney function, we developed a strategy for integrating prior biological knowledge into the existing GWAS data for eGFR from the CKDGen Consortium. Our strategy focuses on single nucleotide polymorphism (SNPs) in genes that are connected by functional evidence, determined by literature mining and gene ontology (GO) hierarchies, to genes near previously validated eGFR associations. It then requires association thresholds consistent with multiple testing, and finally evaluates novel candidates by independent replication. Among the samples of European ancestry, we identified a genome-wide significant SNP in FBXL20 (P = 5.6 × 10(-9)) in meta-analysis of all available data, and additional SNPs at the INHBC, LRP2, PLEKHA1, SLC3A2 and SLC7A6 genes meeting multiple-testing corrected significance for replication and overall P-values of 4.5 × 10(-4)-2.2 × 10(-7). Neither the novel PLEKHA1 nor FBXL20 associations, both further supported by association with eGFR among African Americans and with transcript abundance, would have been implicated by eGFR candidate gene approaches. LRP2, encoding the megalin receptor, was identified through connection with the previously known eGFR gene DAB2 and extends understanding of the megalin system in kidney function. These findings highlight integration of existing genome-wide association data with independent biological knowledge to uncover novel candidate eGFR associations, including candidates lacking known connections to kidney-specific pathways. The strategy may also be applicable to other clinical phenotypes, although more testing will be needed to assess its potential for discovery in general.

Common variants in Mendelian kidney disease genes and their association with renal function.

Many common genetic variants identified by genome-wide association studies for complex traits map to genes previously linked to rare inherited Mendelian disorders. A systematic analysis of common single-nucleotide polymorphisms (SNPs) in genes responsible for Mendelian diseases with kidney phenotypes has not been performed. We thus developed a comprehensive database of genes for Mendelian kidney conditions and evaluated the association between common genetic variants within these genes and kidney function in the general population. Using the Online Mendelian Inheritance in Man database, we identified 731 unique disease entries related to specific renal search terms and confirmed a kidney phenotype in 218 of these entries, corresponding to mutations in 258 genes. We interrogated common SNPs (minor allele frequency >5%) within these genes for association with the estimated GFR in 74,354 European-ancestry participants from the CKDGen Consortium. However, the top four candidate SNPs (rs6433115 at LRP2, rs1050700 at TSC1, rs249942 at PALB2, and rs9827843 at ROBO2) did not achieve significance in a stage 2 meta-analysis performed in 56,246 additional independent individuals, indicating that these common SNPs are not associated with estimated GFR. The effect of less common or rare variants in these genes on kidney function in the general population and disease-specific cohorts requires further research.

Genome-wide association studies on Northern Italy isolated populations provide further support concerning genetic susceptibility for major depressive disorder.

OBJECTIVES: Major depressive disorder (MDD) is a psychiatric disorder with pathogenesis influenced by both genetic and environmental factors. To date, the molecular-level understanding of its aetiology remains unclear. Thus, we aimed to identify genetic variants and susceptibility genes for MDD with a genome-wide association study (GWAS) approach. METHODS: We performed a meta-analysis of GWASs and a gene-based analysis on two Northern Italy isolated populations (cases/controls n = 166/472 and 33/320), followed by replication and polygenic risk score (PRS) analyses in Italian independent samples (cases n = 464, controls n = 339). RESULTS: We identified two novel MDD-associated genes, KCNQ5 (lead SNP rs867262, p = 3.82 × 10-9) and CTNNA2 (rs6729523, p = 1.25 × 10-8). The gene-based analysis revealed another six genes (p < 2.703 × 10-6): GRM7, CTNT4, SNRK, SRGAP3, TRAPPC9, and FHIT. No replication of the genome-wide significant SNPs was found in the independent cohort, even if 14 SNPs around CTNNA2 showed association with MDD and related phenotypes at the nominal level of p (<0.05). Furthermore, the PRS model developed in the discovery cohort discriminated cases and controls in the replication cohort. CONCLUSIONS: Our work suggests new possible genes associated with MDD, and the PRS analysis confirms the polygenic nature of this disorder. Future studies are required to better understand the role of these findings in MDD.

Hearing function and thresholds: a genome-wide association study in European isolated populations identifies new loci and pathways.

BACKGROUND: Hearing is a complex trait, but until now only a few genes are known to contribute to variability of this process. In order to discover genes and pathways that underlie auditory function, a genome-wide association study was carried out within the International Consortium G-EAR. METHODS: Meta-analysis of genome-wide association study's data from six isolated populations of European ancestry for an overall number of 3417 individuals. RESULTS: Eight suggestive significant loci (p<10(-7)) were detected with a series of genes expressed within the inner ear such as: DCLK1, PTPRD, GRM8, CMIP. Additional biological candidates marked by a single nucleotide polymorphism (SNP) with a suggestive association (p<10(-6)) were identified, as well as loci encompassing 'gene desert regions'-genes of unknown function or genes whose function has not be linked to hearing so far. Some of these new loci map to already known hereditary hearing loss loci whose genes still need to be identified. Data have also been used to construct a highly significant 'in silico' pathway for hearing function characterised by a network of 49 genes, 34 of which are certainly expressed in the ear. CONCLUSION: These results provide new insights into the molecular basis of hearing function and may suggest new targets for hearing impairment treatment and prevention.

Case Report: Two cases of apparent discordance between non-invasive prenatal testing (NIPT) and amniocentesis resulting in feto-placental mosaicism of trisomy 21. Issues in diagnosis, investigation and counselling.

The sequencing of cell-free fetal DNA in the maternal plasma through non-invasive prenatal testing (NIPT) is an accurate genetic screening test to detect the most common fetal aneuploidies during pregnancy. The extensive use of NIPT, as a screening method, has highlighted the limits of the technique, including false positive and negative results. Feto-placental mosaicism is a challenging biological issue and is the most frequent cause of false positive and negative results in NIPT screening, and of discrepancy between NIPT and invasive test results. We are reporting on two cases of feto-placental mosaicism of trisomy 21, both with a low-risk NIPT result, identified by ultrasound signs and a subsequent amniocentesis consistent with a trisomy 21. In both cases, after the pregnancy termination, cytogenetic and/or cytogenomic analyses were performed on the placenta and fetal tissues, showing in the first case a mosaicism of trisomy 21 in both the placenta and the fetus, but a mosaicism in the placenta and a complete trisomy 21 in the fetus in the second case. These cases emphasize the need for accurate and complete pre-test NIPT counselling, as well as to identify situations at risk for a possible false negative NIPT result, which may underestimate a potential pathological condition, such as feto-placental mosaicism or fetal trisomy. Post-mortem molecular autopsy may discriminate between placental, fetal and feto-placental mosaicism, and between complete or mosaic fetal chromosomal anomalies. A multidisciplinary approach in counselling, as well as in the interpretation of biological events, is essential for the clarification of complex cases, such as feto-placental mosaicisms.

Identification of a New Mutation in RSK2, the Gene for Coffin-Lowry Syndrome (CLS), in Two Related Patients with Mild and Atypical Phenotypes.

BACKGROUND: Coffin-Lowry syndrome (CLS) is a syndromic form of X-linked intellectual disability, in which specific associated facial, hand, and skeletal abnormalities are diagnostic features. METHODS: In the present study, an unreported missense genetic variant of the ribosomal S6 kinase 2 (RSK2) gene has been identified, by next-generation sequencing, in two related males with two different phenotypes of intellectual disability (ID) and peculiar facial dysmorphisms. We performed functional studies on this variant and another one, already reported in the literature, involving the same amino acid residue but, to date, without an efficient characterization. RESULTS: Our study demonstrated that the two variants involving residue 189 significantly impaired its kinase activity. CONCLUSIONS: We detected a loss-of-function RSK2 mutation with loss in kinase activity in a three-generation family with an X-linked ID.

Metabonomics and population studies: age-related amino acids excretion and inferring networks through the study of urine samples in two Italian isolated populations.

The study of two different Italian isolated populations was combined with a metabonomic approach to better understand tubular handling of amino acids. Levels of amino acids and metabolites have been analyzed by Nucleic Magnetic Resonance and expressed as ratio vs urinary creatinine concentration (mmol/mol). For most of the amino acids there is an age-related U shape pattern of excretion, with the peaks during childhood and old age, and a significant reduction in the adult age. Hierarchical cluster analysis has clearly identified three groups clustering the same amino acids: His, Thr and Ala (group one); Gly and Phe (group two) and a third larger one. Results have been further confirmed by factor and regression analysis, and used to confirm and, in some cases, infer new amino acids networks. As a matter of facts, the identification of strong evidences for clustering of urine excretion of several neutral amino acids suggests the predominant impact of relevant and common transporters.

Modeling the effect of 3 missense AGXT mutations on dimerization of the AGT enzyme in primary hyperoxaluria type 1.

INTRODUCTION: Mutations of the AGXT gene encoding the alanine:glyoxylate aminotransferase liver enzyme (AGT) cause primary hyperoxaluria type 1 (PH1). Here we report a molecular modeling study of selected missense AGXT mutations: the common Gly170Arg and the recently described Gly47Arg and Ser81Leu variants, predicted to be pathogenic using standard criteria. METHODS: Taking advantage of the refined 3D structure of AGT, we computed the dimerization energy of the wild-type and mutated proteins. RESULTS: Molecular modeling predicted that Gly47Arg affects dimerization with a similar effect to that shown previously for Gly170Arg through classical biochemical approaches. In contrast, no effect on dimerization was predicted for Ser81Leu. Therefore, this probably demonstrates pathogenic properties via a different mechanism, similar to that described for the adjacent Gly82Glu mutation that affects pyridoxine binding. CONCLUSION: This study shows that the molecular modeling approach can contribute to evaluating the pathogenicity of some missense variants that affect dimerization. However, in silico studies--aimed to assess the relationship between structural change and biological effects--require the integrated use of more than 1 tool.

Multi-ancestry GWAS of the electrocardiographic PR interval identifies 202 loci underlying cardiac conduction.

The electrocardiographic PR interval reflects atrioventricular conduction, and is associated with conduction abnormalities, pacemaker implantation, atrial fibrillation (AF), and cardiovascular mortality. Here we report a multi-ancestry (N = 293,051) genome-wide association meta-analysis for the PR interval, discovering 202 loci of which 141 have not previously been reported. Variants at identified loci increase the percentage of heritability explained, from 33.5% to 62.6%. We observe enrichment for cardiac muscle developmental/contractile and cytoskeletal genes, highlighting key regulation processes for atrioventricular conduction. Additionally, 8 loci not previously reported harbor genes underlying inherited arrhythmic syndromes and/or cardiomyopathies suggesting a role for these genes in cardiovascular pathology in the general population. We show that polygenic predisposition to PR interval duration is an endophenotype for cardiovascular disease, including distal conduction disease, AF, and atrioventricular pre-excitation. These findings advance our understanding of the polygenic basis of cardiac conduction, and the genetic relationship between PR interval duration and cardiovascular disease.