Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Fine-mapping analysis including over 254,000 East Asian and European descendants identifies 136 putative colorectal cancer susceptibility genes.

Genome-wide association studies (GWAS) have identified more than 200 common genetic variants independently associated with colorectal cancer (CRC) risk, but the causal variants and target genes are mostly unknown. We sought to fine-map all known CRC risk loci using GWAS data from 100,204 cases and 154,587 controls of East Asian and European ancestry. Our stepwise conditional analyses revealed 238 independent association signals of CRC risk, each with a set of credible causal variants (CCVs), of which 28 signals had a single CCV. Our cis-eQTL/mQTL and colocalization analyses using colorectal tissue-specific transcriptome and methylome data separately from 1299 and 321 individuals, along with functional genomic investigation, uncovered 136 putative CRC susceptibility genes, including 56 genes not previously reported. Analyses of single-cell RNA-seq data from colorectal tissues revealed 17 putative CRC susceptibility genes with distinct expression patterns in specific cell types. Analyses of whole exome sequencing data provided additional support for several target genes identified in this study as CRC susceptibility genes. Enrichment analyses of the 136 genes uncover pathways not previously linked to CRC risk. Our study substantially expanded association signals for CRC and provided additional insight into the biological mechanisms underlying CRC development.

Simulation-predicted and -explained inheritance model of pathogenicity confirmed by transgenic mice models.

Variants in the gap junction beta-2 (GJB2) gene are the most common cause of hereditary hearing impairment. However, how GJB2 variants lead to local physicochemical and structural changes in the hexameric ion channels of connexin 26 (Cx26), resulting in hearing impairment, remains elusive. In this study, using molecular dynamics (MD) simulations, we showed that detached inner-wall N-terminal "plugs" aggregated to reduce the channel ion flow in a highly prevalent V37I variant in humans. To examine the predictive ability of the computational platform, an artificial mutant, V37M, of which the effect was previously unknown in hearing loss, was created. Microsecond simulations showed that homo-hexameric V37M Cx26 hemichannels had an abnormal affinity between the inner edge and N-termini to block the narrower side of the cone-shaped Cx26, while the most stable hetero-hexameric channels did not. From the perspective of the conformational energetics of WT and variant Cx26 hexamers, we propose that unaffected carriers could result from a conformational predominance of the WT and pore-shrinkage-incapable hetero-hexamers, while mice with homozygous variants can only harbor an unstable and dysfunctional N-termini-blocking V37M homo-hexamer. Consistent with these predictions, homozygous V37M transgenic mice exhibited apparent hearing loss, but not their heterozygous counterparts, indicating a recessive inheritance mode. Reduced channel conductivity was found in Gjb2V37M/V37M outer sulcus and Claudius cells but not in Gjb2WT/WT cells. We view that the current computational platform could serve as an assessment tool for the pathogenesis and inheritance of GJB2-related hearing impairments and other diseases caused by connexin dysfunction.

The power of TOPMed imputation for the discovery of Latino-enriched rare variants associated with type 2 diabetes.

AIMS/HYPOTHESIS: The Latino population has been systematically underrepresented in large-scale genetic analyses, and previous studies have relied on the imputation of ungenotyped variants based on the 1000 Genomes (1000G) imputation panel, which results in suboptimal capture of low-frequency or Latino-enriched variants. The National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) released the largest multi-ancestry genotype reference panel representing a unique opportunity to analyse rare genetic variations in the Latino population. We hypothesise that a more comprehensive analysis of low/rare variation using the TOPMed panel would improve our knowledge of the genetics of type 2 diabetes in the Latino population. METHODS: We evaluated the TOPMed imputation performance using genotyping array and whole-exome sequence data in six Latino cohorts. To evaluate the ability of TOPMed imputation to increase the number of identified loci, we performed a Latino type 2 diabetes genome-wide association study (GWAS) meta-analysis in 8150 individuals with type 2 diabetes and 10,735 control individuals and replicated the results in six additional cohorts including whole-genome sequence data from the All of Us cohort. RESULTS: Compared with imputation with 1000G, the TOPMed panel improved the identification of rare and low-frequency variants. We identified 26 genome-wide significant signals including a novel variant (minor allele frequency 1.7%; OR 1.37, p=3.4 × 10-9). A Latino-tailored polygenic score constructed from our data and GWAS data from East Asian and European populations improved the prediction accuracy in a Latino target dataset, explaining up to 7.6% of the type 2 diabetes risk variance. CONCLUSIONS/INTERPRETATION: Our results demonstrate the utility of TOPMed imputation for identifying low-frequency variants in understudied populations, leading to the discovery of novel disease associations and the improvement of polygenic scores. DATA AVAILABILITY: Full summary statistics are available through the Common Metabolic Diseases Knowledge Portal ( https://t2d.hugeamp.org/downloads.html ) and through the GWAS catalog ( https://www.ebi.ac.uk/gwas/ , accession ID: GCST90255648). Polygenic score (PS) weights for each ancestry are available via the PGS catalog ( https://www.pgscatalog.org , publication ID: PGP000445, scores IDs: PGS003443, PGS003444 and PGS003445).

Metabolomic profiling of glucose homeostasis in African Americans: the Insulin Resistance Atherosclerosis Family Study (IRAS-FS).

INTRODUCTION: African Americans are at increased risk for type 2 diabetes. OBJECTIVES: This work aimed to examine metabolomic signature of glucose homeostasis in African Americans. METHODS: We used an untargeted liquid chromatography-mass spectrometry metabolomic approach to comprehensively profile 727 plasma metabolites among 571 African Americans from the Insulin Resistance Atherosclerosis Family Study (IRAS-FS) and investigate the associations between these metabolites and both the dynamic (SI, insulin sensitivity; AIR, acute insulin response; DI, disposition index; and SG, glucose effectiveness) and basal (HOMA-IR and HOMA-B) measures of glucose homeostasis using univariate and regularized regression models. We also compared the results with our previous findings in the IRAS-FS Mexican Americans. RESULTS: We confirmed increased plasma metabolite levels of branched-chain amino acids and their metabolic derivatives, 2-aminoadipate, 2-hydroxybutyrate, glutamate, arginine and its metabolic derivatives, carbohydrate metabolites, and medium- and long-chain fatty acids were associated with insulin resistance, while increased plasma metabolite levels in the glycine, serine and threonine metabolic pathway were associated with insulin sensitivity. We also observed a differential ancestral effect of glutamate on glucose homeostasis with significantly stronger effects observed in African Americans than those previously observed in Mexican Americans. CONCLUSION: We extended the observations that metabolites are useful biomarkers in the identification of prediabetes in individuals at risk of type 2 diabetes in African Americans. We revealed, for the first time, differential ancestral effect of certain metabolites (i.e., glutamate) on glucose homeostasis traits. Our study highlights the need for additional comprehensive metabolomic studies in well-characterized multiethnic cohorts.

Generation of induced pluripotent stem cells (IBMSi027-A) from a patient with hearing loss carrying WFS1 c.2051C > T (p.Ala684Val) variant.

Pathogenic variants of the WFS1 gene can cause recessive-inherited Wolfram syndrome or dominant-inherited Wolfram-like syndrome with optic atrophy and hearing impairment. Using the Sendai virus delivery system, we generated induced pluripotent stem cells from the peripheral blood mononuclear cells of a female patient with the WFS1 pathogenic variant c.2051C > T (p.Ala684Val). The resulting induced pluripotent stem cells exhibited a normal karyotype and pluripotency, as confirmed using immunofluorescence staining, and differentiated into three germ layers in vivo. This cellular model provides a useful platform for investigating the pathogenic mechanisms of both blindness and deafness related to WFS1 variants.

Deciphering colorectal cancer genetics through multi-omic analysis of 100,204 cases and 154,587 controls of European and east Asian ancestries.

Colorectal cancer (CRC) is a leading cause of mortality worldwide. We conducted a genome-wide association study meta-analysis of 100,204 CRC cases and 154,587 controls of European and east Asian ancestry, identifying 205 independent risk associations, of which 50 were unreported. We performed integrative genomic, transcriptomic and methylomic analyses across large bowel mucosa and other tissues. Transcriptome- and methylome-wide association studies revealed an additional 53 risk associations. We identified 155 high-confidence effector genes functionally linked to CRC risk, many of which had no previously established role in CRC. These have multiple different functions and specifically indicate that variation in normal colorectal homeostasis, proliferation, cell adhesion, migration, immunity and microbial interactions determines CRC risk. Crosstissue analyses indicated that over a third of effector genes most probably act outside the colonic mucosa. Our findings provide insights into colorectal oncogenesis and highlight potential targets across tissues for new CRC treatment and chemoprevention strategies.

Gene therapy with a synthetic adeno-associated viral vector improves audiovestibular phenotypes in Pjvk-mutant mice.

Recessive PJVK mutations that cause a deficiency of pejvakin, a protein expressed in both sensory hair cells and first-order neurons of the inner ear, are an important cause of hereditary hearing impairment. Patients with PJVK mutations garner limited benefits from cochlear implantation; thus, alternative biological therapies may be required to address this clinical difficulty. The synthetic adeno-associated viral vector Anc80L65, with its wide tropism and high transduction efficiency in various inner ear cells, may provide a solution. We delivered the PJVK transgene to the inner ear of Pjvk mutant mice using the synthetic Anc80L65 vector. We observed robust exogenous pejvakin expression in the hair cells and neurons of the cochlea and vestibular organs. Subsequent morphologic and audiologic studies demonstrated significant restoration of spiral ganglion neuron density and hair cells in the cochlea, along with partial recovery of sensorineural hearing impairment. In addition, we observed a recovery of vestibular ganglion neurons and balance function to WT levels. Our study demonstrates the utility of Anc80L65-mediated gene delivery in Pjvk mutant mice and provides insights into the potential of gene therapy for PJVK-related inner ear deficits.

Generation of induced pluripotent stem cells from a patient with hearing loss carrying OPA1 c.1468T>C (p.Cys490Arg) variant.

Pathogenic variants of OPA1 have been associated with autosomal dominant optic atrophy (DOA), leading to optic, auditory, and other sensorineural neuropathies and myopathies. Using the Sendai virus delivery system, we generated induced pluripotent stem cells from the peripheral blood mononuclear cells of a female patient with the OPA1 pathogenic variant c.1468T>C (p.Cys490Arg). The resulting induced pluripotent stem cells exhibited a normal karyotype and pluripotency, as confirmed using immunofluorescence staining, and differentiated into three germ layers in vivo. This cellular model is a useful platform for investigating the pathogenic mechanisms of both blindness and deafness related to OPA1 variants.

Large-scale genome-wide association study of coronary artery disease in genetically diverse populations.

We report a genome-wide association study (GWAS) of coronary artery disease (CAD) incorporating nearly a quarter of a million cases, in which existing studies are integrated with data from cohorts of white, Black and Hispanic individuals from the Million Veteran Program. We document near equivalent heritability of CAD across multiple ancestral groups, identify 95 novel loci, including nine on the X chromosome, detect eight loci of genome-wide significance in Black and Hispanic individuals, and demonstrate that two common haplotypes at the 9p21 locus are responsible for risk stratification in all populations except those of African origin, in which these haplotypes are virtually absent. Moreover, in the largest GWAS for angiographically derived coronary atherosclerosis performed to date, we find 15 loci of genome-wide significance that robustly overlap with established loci for clinical CAD. Phenome-wide association analyses of novel loci and polygenic risk scores (PRSs) augment signals related to insulin resistance, extend pleiotropic associations of these loci to include smoking and family history, and precisely document the markedly reduced transferability of existing PRSs to Black individuals. Downstream integrative analyses reinforce the critical roles of vascular endothelial, fibroblast, and smooth muscle cells in CAD susceptibility, but also point to a shared biology between atherosclerosis and oncogenesis. This study highlights the value of diverse populations in further characterizing the genetic architecture of CAD.

Genetic discovery and risk characterization in type 2 diabetes across diverse populations.

Genomic discovery and characterization of risk loci for type 2 diabetes (T2D) have been conducted primarily in individuals of European ancestry. We conducted a multiethnic genome-wide association study of T2D among 53,102 cases and 193,679 control subjects from African, Hispanic, Asian, Native Hawaiian, and European population groups in the Population Architecture Genomics and Epidemiology (PAGE) and Diabetes Genetics Replication and Meta-analysis (DIAGRAM) Consortia. In individuals of African ancestry, we discovered a risk variant in the TGFB1 gene (rs11466334, risk allele frequency (RAF) = 6.8%, odds ratio [OR] = 1.27, p = 2.06 × 10-8), which replicated in independent studies of African ancestry (p = 6.26 × 10-23). We identified a multiethnic risk variant in the BACE2 gene (rs13052926, RAF = 14.1%, OR = 1.08, p = 5.75 × 10-9), which also replicated in independent studies (p = 3.45 × 10-4). We also observed a significant difference in the performance of a multiethnic genetic risk score (GRS) across population groups (pheterogeneity = 3.85 × 10-20). Comparing individuals in the top GRS risk category (40%-60%), the OR was highest in Asians (OR = 3.08) and European (OR = 2.94) ancestry populations, followed by Hispanic (OR = 2.39), Native Hawaiian (OR = 2.02), and African ancestry (OR = 1.57) populations. These findings underscore the importance of genetic discovery and risk characterization in diverse populations and the urgent need to further increase representation of non-European ancestry individuals in genetics research to improve genetic-based risk prediction across populations.

Insights into phenotypic differences between humans and mice with p.T721M and other C-terminal variants of the SLC26A4 gene.

Recessive variants of the SLC26A4 gene are an important cause of hereditary hearing impairment. Several transgenic mice with different Slc26a4 variants have been generated. However, none have recapitulated the auditory phenotypes in humans. Of the SLC26A4 variants identified thus far, the p.T721M variant is of interest, as it appears to confer a more severe pathogenicity than most of the other missense variants, but milder pathogenicity than non-sense and frameshift variants. Using a genotype-driven approach, we established a knock-in mouse model homozygous for p.T721M. To verify the pathogenicity of p.T721M, we generated mice with compound heterozygous variants by intercrossing Slc26a4+/T721M mice with Slc26a4919-2A>G/919-2A>G mice, which segregated the c.919-2A > G variant with abolished Slc26a4 function. We then performed serial audiological assessments, vestibular evaluations, and inner ear morphological studies. Surprisingly, both Slc26a4T721M/T721M and Slc26a4919-2A>G/T721M showed normal audiovestibular functions and inner ear morphology, indicating that p.T721M is non-pathogenic in mice and a single p.T721M allele is sufficient to maintain normal inner ear physiology. The evidence together with previous reports on mouse models with Slc26a4 p.C565Y and p.H723R variants, support our speculation that the absence of audiovestibular phenotypes in these mouse models could be attributed to different protein structures at the C-terminus of human and mouse pendrin.

Multiethnic Genome-Wide Association Study of Subclinical Atherosclerosis in Individuals With Type 2 Diabetes.

BACKGROUND: Coronary artery calcification (CAC) and carotid artery intima-media thickness (cIMT) are measures of subclinical atherosclerosis in asymptomatic individuals and strong risk factors for cardiovascular disease. Type 2 diabetes (T2D) is an independent cardiovascular disease risk factor that accelerates atherosclerosis. METHODS: We performed meta-analyses of genome-wide association studies in up to 2500 T2D individuals of European ancestry (EA) and 1590 T2D individuals of African ancestry with or without exclusion of prevalent cardiovascular disease, for CAC measured by cardiac computed tomography, and 3608 individuals of EA and 838 individuals of African ancestry with T2D for cIMT measured by ultrasonography within the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology) Consortium. RESULTS: We replicated 2 loci (rs9369640 and rs9349379 near PHACTR1 and rs10757278 near CDKN2B) for CAC and one locus for cIMT (rs7412 and rs445925 near APOE-APOC1) that were previously reported in the general EA populations. We identified one novel CAC locus (rs8000449 near CSNK1A1L/LINC00547/POSTN at 13q13.3) at P=2.0×10-8 in EA. No additional loci were identified with the meta-analyses of EA and African ancestry. The expression quantitative trait loci analysis with nearby expressed genes derived from arterial wall and metabolic tissues from the Genotype-Tissue Expression project pinpoints POSTN, encoding a matricellular protein involved in bone formation and bone matrix organization, as the potential candidate gene at this locus. In addition, we found significant associations (P<3.1×10-4) for 3 previously reported coronary artery disease loci for these subclinical atherosclerotic phenotypes (rs2891168 near CDKN2B-AS1 and rs11170820 near FLJ12825 for CAC, and rs7412 near APOE for cIMT). CONCLUSIONS: Our results provide potential biological mechanisms that could link CAC and cIMT to increased cardiovascular disease risk in individuals with T2D.

Toward the Pathogenicity of the SLC26A4 p.C565Y Variant Using a Genetically Driven Mouse Model.

Recessive variants of the SLC26A4 gene are globally a common cause of hearing impairment. In the past, cell lines and transgenic mice were widely used to investigate the pathogenicity associated with SLC26A4 variants. However, discrepancies in pathogenicity between humans and cell lines or transgenic mice were documented for some SLC26A4 variants. For instance, the p.C565Y variant, which was reported to be pathogenic in humans, did not exhibit functional pathogenic consequences in cell lines. To address the pathogenicity of p.C565Y, we used a genotype-based approach in which we generated knock-in mice that were heterozygous (Slc26a4+/C565Y), homozygous (Slc26a4C565Y/C565Y), and compound heterozygous (Slc26a4919-2A>G/C565Y) for this variant. Subsequent phenotypic characterization revealed that mice with these genotypes demonstrated normal auditory and vestibular functions, and normal inner-ear morphology and pendrin expression. These findings indicate that the p.C565Y variant is nonpathogenic for mice, and that a single p.C565Y allele is sufficient to maintain normal inner-ear physiology in mice. Our results highlight the differences in pathogenicity associated with certain SLC26A4 variants between transgenic mice and humans, which should be considered when interpreting the results of animal studies for SLC26A4-related deafness.

Low intensity ultrasound enhances cisplatin uptake in vitro by cochlear hair cells.

Drug delivery to the inner ear has been challenging due to the blood-labyrinth barrier. Intracochlear drug delivery is an invasive alternative with less pharmacokinetic variables. In this study, the effect of low intensity ultrasound on drug uptake by hair cells is investigated. Cochlear explants harvested from newborn mice were cultured in a medium containing cisplatin to emulate drug delivered to the endolymph. The results demonstrated the exposure to ultrasound stimulation effectively enhanced cisplatin uptake by hair cells. The uptake started from the apical side of the hair cells and progressed inward as the exposure time increased.

Causal associations between urinary sodium with body mass, shape and composition: a Mendelian randomization study.

Observational studies have found associations between urinary sodium (UNa) with obesity, body shape and composition; but the findings may be biased by residual confounding. The objective of this two-sample Mendelian randomization (MR) study was to analyze their causal associations in both sex-combined and sex-specific models. Genome-wide association studies of UNa, body mass index (BMI), BMI-adjusted waist-to-hip ratio (WHR), body fat (BF) percentage and estimated glomerular filtration rate (eGFR) were identified. We initially extracted fifty SNPs associated with UNa at significance level of 5 × 10-8, but further removed those SNPs with potential horizontal pleiotropy. Univariable and multivariable MR with adjustment for eGFR were performed. Inverse-variance weighted MR was performed as the primary analysis, with MR-Egger methods as sensitivity analysis. The potential bidirectional association between BMI and UNa was investigated. All exposure and outcomes were continuous, and the effect measure was regression coefficients (beta) and their 95% confidence intervals (95% CI). The total sample size was up to 322 154. UNa was causally associated with increased BMI in both men [eGFR-adjusted beta 0.443 (0.163-0.724)] and women [0.594 (0.333-0.855)]. UNa caused BF percentage increase in men [0.622 (0.268-0.976)] and women [0.334 (0.007-0.662)]. UNa significantly elevated BMI-adjusted WHR in men [0.321 (0.094-0.548)], but not in women [0.170 (- 0.052 to 0.391)]. Additionally, we found that BMI causally increased UNa [0.043 (0.023-0.063)]. UNa increased BMI and BF percentage. Salt intake affects male body shape by increasing BMI-adjusted WHR, but showed no effects on female body shape. The bidirectional association between BMI and UNa suggested that salt reduction measures and weight reduction measures should be implemented simultaneously to break the vicious cycle and gain more health benefits.

Genetic Studies of Leptin Concentrations Implicate Leptin in the Regulation of Early Adiposity.

Leptin influences food intake by informing the brain about the status of body fat stores. Rare LEP mutations associated with congenital leptin deficiency cause severe early-onset obesity that can be mitigated by administering leptin. However, the role of genetic regulation of leptin in polygenic obesity remains poorly understood. We performed an exome-based analysis in up to 57,232 individuals of diverse ancestries to identify genetic variants that influence adiposity-adjusted leptin concentrations. We identify five novel variants, including four missense variants, in LEP, ZNF800, KLHL31, and ACTL9, and one intergenic variant near KLF14. The missense variant Val94Met (rs17151919) in LEP was common in individuals of African ancestry only, and its association with lower leptin concentrations was specific to this ancestry (P = 2 × 10-16, n = 3,901). Using in vitro analyses, we show that the Met94 allele decreases leptin secretion. We also show that the Met94 allele is associated with higher BMI in young African-ancestry children but not in adults, suggesting that leptin regulates early adiposity.

Generation and pathological characterization of a transgenic mouse model carrying a missense PJVK mutation.

Hearing loss is the most prevalent hereditary sensory disorder in children. Approximately 2 in 1000 infants are affected by genetic hearing loss. The PJVK gene, which encodes the pejvakin protein, has been linked to autosomal recessive non-syndromic hearing loss DFNB59. Previous clinical studies have revealed that PJVK mutations might be associated with a wide spectrum of auditory manifestations, ranging from hearing loss of pure cochlear origin to that involving the retrocochlear central auditory pathway. The phenotypic variety makes the pathogenesis of this disease difficult to determine. Similarly, mouse models carrying different Pjvk defects show phenotypic variability and inconsistency. In this study, we generated a knockin mouse model carrying the c.874G > A (p.G292R) variant to model and investigate the auditory and vestibular phenotypes of DFNB59.

Efficient in Utero Gene Transfer to the Mammalian Inner Ears by the Synthetic Adeno-Associated Viral Vector Anc80L65.

Sensorineural hearing loss is one of the most common sensory disorders worldwide. Recent advances in vector design have paved the way for investigations into the use of adeno-associated vectors (AAVs) for hearing disorder gene therapy. Numerous AAV serotypes have been discovered to be applicable to inner ears, constituting a key advance for gene therapy for sensorineural hearing loss, where transduction efficiency of AAV in inner ear cells is critical for success. One such viral vector, AAV2/Anc80L65, has been shown to yield high expression in the inner ears of mice treated as neonates or adults. Here, to evaluate the feasibility of prenatal gene therapy for deafness, we assessed the transduction efficiency of AAV2/Anc80L65-eGFP (enhanced green fluorescent protein) after microinjection into otocysts in utero. This embryonic delivery method achieved high transduction efficiency in both inner and outer hair cells of the cochlea. Additionally, the transduction efficiency was high in the hair cells of the vestibules and semicircular canals and in spiral ganglion neurons. Our results support the potential of Anc80L65 as a gene therapy vehicle for prenatal inner ear disorders.