Search

Abstract Asthma has striking disparities across ancestral groups, but the molecular underpinning of these differences is poorly understood and minimally studied. A goal of the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) is to understand multi-omic signatures of asthma focusing on populations of African ancestry. RNASeq and DNA methylation data are generated from nasal epithelium including cases (current asthma, N = 253) and controls (never-asthma, N = 283) from 7 different geographic sites to identify differentially expressed genes (DEGs) and gene networks. We identify 389 DEGs; the top DEG, FN1, was downregulated in cases (q = 3.26 × 10−9) and encodes fibronectin which plays a role in wound healing. The top three gene expression modules implicate networks related to immune response (CEACAM5; p = 9.62 × 10−16 and CPA3; p = 2.39 × 10−14) and wound healing (FN1; p = 7.63 × 10−9). Multi-omic analysis identifies FKBP5, a co-chaperone of glucocorticoid receptor signaling known to be involved in drug response in asthma, where the association between nasal epithelium gene expression is likely regulated by methylation and is associated with increased use of inhaled corticosteroids. This work reveals molecular dysregulation on three axes – increased Th2 inflammation, decreased capacity for wound healing, and impaired drug response – that may play a critical role in asthma within the African Diaspora.

Abstract Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation.

Abstract Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Abstract Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods focus only on participants with distinct primary continental ancestry without accommodating recently-admixed individuals with mosaic continental ancestry backgrounds for different segments of their genomes. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals. GAUDI explicitly models ancestry-differential effects while borrowing information across segments with shared ancestry in admixed genomes. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses for traits with associated variants exhibiting ancestral-differential effects. Leveraging data from the Women’s Health Initiative study, we show that GAUDI improves PRS prediction of white blood cell count and C-reactive protein in African Americans by > 64% compared to alternative methods, and even outperforms PRS-CSx with large European GWAS for some scenarios. We believe GAUDI will be a valuable tool to mitigate disparities in PRS performance in admixed individuals.

Abstract Telomere length (TL) attrition, epigenetic age acceleration, and mitochondrial DNA copy number (mtDNAcn) decline are established hallmarks of aging. Each has been individually associated with Alzheimer’s dementia, cognitive function, and pathologic Alzheimer’s disease (AD). Epigenetic age and mtDNAcn have been studied in brain tissue directly but prior work on TL in brain is limited to small sample sizes and most studies have examined leukocyte TL. Importantly, TL, epigenetic age clocks, and mtDNAcn have not been studied jointly in brain tissue from an AD cohort. We examined dorsolateral prefrontal cortex (DLPFC) tissue from N = 367 participants of the Religious Orders Study (ROS) or the Rush Memory and Aging Project (MAP). TL and mtDNAcn were estimated from whole genome sequencing (WGS) data and cortical clock age was computed on 347 CpG sites. We examined dementia, MCI, and level of and change in cognition, pathologic AD, and three quantitative AD traits, as well as measures of other neurodegenerative diseases and cerebrovascular diseases (CVD). We previously showed that mtDNAcn from DLPFC brain tissue was associated with clinical and pathologic features of AD. Here, we show that those associations are independent of TL. We found TL to be associated with β-amyloid levels (beta = − 0.15, p = 0.023), hippocampal sclerosis (OR = 0.56, p = 0.0015) and cerebral atherosclerosis (OR = 1.44, p = 0.0007). We found strong associations between mtDNAcn and clinical measures of AD. The strongest associations with pathologic measures of AD were with cortical clock and there were associations of mtDNAcn with global AD pathology and tau tangles. Of the other pathologic traits, mtDNAcn was associated with hippocampal sclerosis, macroscopic infarctions and CAA and cortical clock was associated with Lewy bodies. Multi-modal age acceleration, accelerated aging on both mtDNAcn and cortical clock, had greater effect size than a single measure alone. These findings highlight for the first time that age acceleration determined on multiple genomic measures, mtDNAcn and cortical clock may have a larger effect on AD/AD related disorders (ADRD) pathogenesis than single measures.

Most genetic association studies have been done on single nucleotide polymorphisms and small indels, while other types of variants have been less studied. Here, the authors use whole genome sequencing in a diverse population to identify and provide experimental evidence for associations between structural variants and blood-cell traits.

Although the common genetic variants contributing to blood lipid levels have been studied, the contribution of rare variants is less understood. Here, the authors perform a rare coding and noncoding variant association study of blood lipid levels using whole genome sequencing data.

Abstract Background Asthma is the most common chronic disease in children, occurring at higher frequencies and with more severe disease in children with African ancestry. Methods We tested for association with haplotypes at the most replicated and significant childhood-onset asthma locus at 17q12-q21 and asthma in European American and African American children. Following this, we used whole-genome sequencing data from 1060 African American and 100 European American individuals to identify novel variants on a high-risk African American–specific haplotype. We characterized these variants in silico using gene expression and ATAC-seq data from airway epithelial cells, functional annotations from ENCODE, and promoter capture (pc)Hi-C maps in airway epithelial cells. Candidate causal variants were then assessed for correlation with asthma-associated phenotypes in African American children and adults. Results Our studies revealed nine novel African-specific common variants, enriched on a high-risk asthma haplotype, which regulated the expression of GSDMA in airway epithelial cells and were associated with features of severe asthma. Using ENCODE annotations, ATAC-seq, and pcHi-C, we narrowed the associations to two candidate causal variants that are associated with features of T2 low severe asthma. Conclusions Previously unknown genetic variation at the 17q12-21 childhood-onset asthma locus contributes to asthma severity in individuals with African ancestries. We suggest that many other population-specific variants that have not been discovered in GWAS contribute to the genetic risk for asthma and other common diseases.

This study of 23,000 non-diabetic individuals highlights loci associated with fasting glucose and fasting insulin in diverse cohorts with whole genome sequence data.

Abstract Background While large genome-wide association studies have identified nearly one thousand loci associated with variation in blood pressure, rare variant identification is still a challenge. In family-based cohorts, genome-wide linkage scans have been successful in identifying rare genetic variants for blood pressure. This study aims to identify low frequency and rare genetic variants within previously reported linkage regions on chromosomes 1 and 19 in African American families from the Trans-Omics for Precision Medicine (TOPMed) program. Genetic association analyses weighted by linkage evidence were completed with whole genome sequencing data within and across TOPMed ancestral groups consisting of 60,388 individuals of European, African, East Asian, Hispanic, and Samoan ancestries. Results Associations of low frequency and rare variants in RCN3 and multiple other genes were observed for blood pressure traits in TOPMed samples. The association of low frequency and rare coding variants in RCN3 was further replicated in UK Biobank samples (N = 403,522), and reached genome-wide significance for diastolic blood pressure (p = 2.01 × 10− 7). Conclusions Low frequency and rare variants in RCN3 contributes blood pressure variation. This study demonstrates that focusing association analyses in linkage regions greatly reduces multiple-testing burden and improves power to identify novel rare variants associated with blood pressure traits.

Summary: Anthropometric traits, measuring body size and shape, are highly heritable and significant clinical risk factors for cardiometabolic disorders. These traits have been extensively studied in genome-wide association studies (GWASs), with hundreds of genome-wide significant loci identified. We performed a whole-exome sequence analysis of the genetics of height, body mass index (BMI) and waist/hip ratio (WHR). We meta-analyzed single-variant and gene-based associations of whole-exome sequence variation with height, BMI, and WHR in up to 22,004 individuals, and we assessed replication of our findings in up to 16,418 individuals from 10 independent cohorts from Trans-Omics for Precision Medicine (TOPMed). We identified four trait associations with single-nucleotide variants (SNVs; two for height and two for BMI) and replicated the LECT2 gene association with height. Our expression quantitative trait locus (eQTL) analysis within previously reported GWAS loci implicated CEP63 and RFT1 as potential functional genes for known height loci. We further assessed enrichment of SNVs, which were monogenic or syndromic variants within loci associated with our three traits. This led to the significant enrichment results for height, whereas we observed no Bonferroni-corrected significance for all SNVs. With a sample size of ∼20,000 whole-exome sequences in our discovery dataset, our findings demonstrate the importance of genomic sequencing in genetic association studies, yet they also illustrate the challenges in identifying effects of rare genetic variants.

Konigsberg et al. profile DNA methylation in blood samples from SARS-CoV-2 cases and controls. The authors use machine learning to classify infected vs. non-infected individuals and predict clinical outcomes related to disease severity.

RationalePrevious studies identified an interaction between HLA and oral peanut exposure. HLA-DQA1*01:02 had a protective role with the induction of Ara h 2 epitope-specific IgG4 associated with peanut consumption during the LEAP clinical trial for prevention of peanut allergy, while it was a risk allele for peanut allergy in the peanut avoidance group. We have now evaluated this gene-environment interaction in two subsequent peanut oral immunotherapy (OIT) trials - IMPACT and POISED - to better understand the potential for the HLA-DQA1*01:02 allele as an indicator of higher likelihood of desensitization, sustained unresponsiveness, and peanut allergy remission.MethodsWe determined HLA-DQA1*01:02 carrier status using genome sequencing from POISED (N=118, age: 7-55yr) and IMPACT (N=126, age: 12-

Yang et al. examine whether genetic ancestry is associated with genetic variation in fatty acid desaturases and plasma phospholipid levels of long chain polyunsaturated fatty acids (LC-PUFAs) in Hispanic Americans. They find strong associations between Amerind genetic ancestry and LC-PUFA levels; and report that the well-known FADS rs174537 variant was associated with several metabolic, inflammatory and anthropomorphic traits including circulating triglycerides.

Platelet aggregation is associated with myocardial infarction and stroke. Here, the authors have conducted a whole genome sequencing association study on platelet aggregation, discovering a locus in RGS18, where enhancer assays suggest an effect on activity of haematopoeitic lineage transcription factors.

The influence of X chromosome genetic variation on blood lipids and coronary heart disease (CHD) is not well understood. Here, the authors analyse X chromosome sequencing data across 65,322 multi-ancestry individuals, identifying associations of the Xq23 locus with lipid changes and reduced risk of CHD and diabetes mellitus.

Drugs targeting cardiovascular disease (CVD) can have negative consequences for liver function. Here, the authors combine genome wide analyses on 69,479 individuals to identify loss-of-function variants with beneficial effects on CVD-related traits without negative impacts on liver function.

Polygenic risk scores (PRSs) aggregate the effects of genetic variants across the genome and are used to predict risk of complex diseases, such as obesity. Current PRSs only include common variants (minor allele frequency (MAF) ≥1%), whereas the contribution of rare variants in PRSs to predict disease remains unknown. Here, we examine whether augmenting the standard common variant PRS (PRScommon) with a rare variant PRS (PRSrare) improves prediction of obesity. We used genome-wide genotyped and imputed data on 451,145 European-ancestry participants of the UK Biobank, as well as whole exome sequencing (WES) data on 184,385 participants. We performed single variant analyses (for both common and rare variants) and gene-based analyses (for rare variants) for association with BMI (kg/m2), obesity (BMI ≥ 30 kg/m2), and extreme obesity (BMI ≥ 40 kg/m2). We built PRSscommon and PRSsrare using a range of methods (Clumping+Thresholding [C+T], PRS-CS, lassosum, gene-burden test). We selected the best-performing PRSs and assessed their performance in 36,757 European-ancestry unrelated participants with whole genome sequencing (WGS) data from the Trans-Omics for Precision Medicine (TOPMed) program. The best-performing PRScommon explained 10.1% of variation in BMI, and 18.3% and 22.5% of the susceptibility to obesity and extreme obesity, respectively, whereas the best-performing PRSrare explained 1.49%, and 2.97% and 3.68%, respectively. The PRSrare was associated with an increased risk of obesity and extreme obesity (ORobesity = 1.37 per SDPRS, Pobesity = 1.7x10-85; ORextremeobesity = 1.55 per SDPRS, Pextremeobesity = 3.8x10-40), which was attenuated, after adjusting for PRScommon (ORobesity = 1.08 per SDPRS, Pobesity = 9.8x10-6; ORextremeobesity= 1.09 per SDPRS, Pextremeobesity = 0.02). When PRSrare and PRScommon are combined, the increase in explained variance attributed to PRSrare was small (incremental Nagelkerke R2 = 0.24% for obesity and 0.51% for extreme obesity). Consistently, combining PRSrare to PRScommon provided little improvement to the prediction of obesity (PRSrare AUC = 0.591; PRScommon AUC = 0.708; PRScombined AUC = 0.710). In summary, while rare variants show convincing association with BMI, obesity and extreme obesity, the PRSrare provides limited improvement over PRScommon in the prediction of obesity risk, based on these large populations.

Human diets in developed countries such as the US have changed dramatically over the past 75 years, leading to increased obesity, inflammation, and cardiometabolic dysfunction. Evidence over the past decade indicates that the interaction of genetic variation with changes in the intake of 18-carbon essential dietary omega-6 (n-6) and omega-3 (n-3) polyunsaturated fatty acids (PUFA), linoleic acid (LA) and α-linolenic acid (ALA), respectively, has impacted numerous molecular and clinical phenotypes. Interactions are particularly relevant with the FADS1 and FADS2 genes, which encode key fatty acid desaturases in the pathway that converts LA and ALA to their long chain (≥20 carbons), highly unsaturated fatty acid (HUFA) counterparts. These gene by nutrient interactions affect the levels and balance of n-6 and n-3 HUFA that in turn are converted to a wide array of lipids with signaling roles, including eicosanoids, docosanoids, other oxylipins and endocannabinoids. With few exceptions, n-6 HUFA are precursors of pro-inflammatory/pro-thrombotic signaling lipids, and n-3 HUFA are generally anti-inflammatory/anti-thrombotic. We and others have demonstrated that African ancestry populations have much higher frequencies (vs. European-, Asian- or indigenous Americas-ancestry populations) of a FADS “derived” haplotype that is associated with the efficient conversion of high levels of dietary n-6 PUFA to pro-inflammatory n-6 HUFA. By contrast, an “ancestral” haplotype, carrying alleles associated with a limited capacity to synthesize HUFA, which can lead to n-3 HUFA deficiency, is found at high frequency in certain Hispanic populations and is nearly fixed in several indigenous populations from the Americas. Based on these observations, a focused secondary subgroup analysis of the VITAL n-3 HUFA supplementation trial stratifying the data based on self-reported ancestry revealed that African Americans may benefit from n-3 HUFA supplementation, and both ancestry and FADS variability should be factored into future clinical trials design.

We investigated the interplay between genetics and oral peanut protein exposure in the determination of the immunological response to peanut using the targeted intervention in the LEAP clinical trial. We identified an association between peanut-specific IgG4 and HLA-DQA1*01:02 that was only observed in the presence of sustained oral peanut protein exposure. The association between IgG4 and HLA-DQA1*01:02 was driven by IgG4 specific for the Ara h 2 component. Once peanut consumption ceased, the association between IgG4-specific Ara h 2 and HLA-DQA1*01:02 was attenuated. The association was validated by observing expanded IgG4-specific epitopes in people who carried HLA-DQA1*01:02. Notably, we confirmed the previously reported associations with HLA-DQA1*01:02 and peanut allergy risk in the absence of oral peanut protein exposure. Interaction between HLA and presence or absence of exposure to peanut in an allergen- and epitope-specific manner implicates a mechanism of antigen recognition that is fundamental to driving immune responses related to allergy risk or protection.

Summary: Genetic studies on telomere length are important for understanding age-related diseases. Prior GWASs for leukocyte TL have been limited to European and Asian populations. Here, we report the first sequencing-based association study for TL across ancestrally diverse individuals (European, African, Asian, and Hispanic/Latino) from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. We used whole-genome sequencing (WGS) of whole blood for variant genotype calling and the bioinformatic estimation of telomere length in n = 109,122 individuals. We identified 59 sentinel variants (p < 5 × 10−9) in 36 loci associated with telomere length, including 20 newly associated loci (13 were replicated in external datasets). There was little evidence of effect size heterogeneity across populations. Fine-mapping at OBFC1 indicated that the independent signals colocalized with cell-type-specific eQTLs for OBFC1 (STN1). Using a multi-variant gene-based approach, we identified two genes newly implicated in telomere length, DCLRE1B (SNM1B) and PARN. In PheWAS, we demonstrated that our TL polygenic trait scores (PTSs) were associated with an increased risk of cancer-related phenotypes.

Telomere length (TL) is widely studied as a possible biomarker for stress-related cellular aging and decreased longevity. There have been conflicting findings about the relationship between family caregiving stress and TL. Several initial cross-sectional studies have found associations between longer duration of caregiving or perceived stressfulness of caregiving and shortened TL, suggesting that caregiving poses grave risks to health. Previous reviews have suggested the need for longitudinal methods to investigate this topic. This study examined the association between the transition to family caregiving and change in TL across ~9 years. Data was utilized from the Caregiving Transitions Study, an ancillary study to the Reasons for Geographic and Racial Differences in Stroke (REGARDS) study. TL was assayed using qPCR and analyzed as the telomere-to-single copy gene ratio for each participant at baseline and follow-up. General linear models examined the association between caregiving status and the change in TL for 208 incident caregivers and 205 controls, as well as associations between perceived stress and TL among caregivers. No association was found between TL change and caregiving (p = 0.494), and fully adjusted models controlling for health and socioeconomic factors did not change the null relationship (p = 0.305). Among caregivers, no association was found between perceived caregiving stress and change in TL (p = 0.336). In contrast to earlier cross-sectional studies, this longitudinal, population-based study did not detect a significant relationship between the transition into a family caregiving role and changes in TL over time. Given the widespread citation of previous findings suggesting that caregiving shortens telomeres and places caregivers at risk of early mortality, these results demonstrate the potential need of a more balanced narrative about caregiving.

Coronary artery disease (CAD) remains a major cause of mortality and morbidity worldwide. The aggregation of activated platelets on a ruptured atherosclerotic plaque is a critical step in most acute cardiovascular events like myocardial infarction. Platelet aggregation both at baseline and after aspirin is highly heritable. Genome-wide association studies (GWAS) have identified a common variant within the first intron of the platelet endothelial aggregation receptor1 (PEAR1), to be robustly associated with platelet aggregation. In this study, we used targeted deep sequencing to fine-map the prior GWAS peak and identify additional rare variants of PEAR1 that account for missing heritability in platelet aggregation within the GeneSTAR families. In this study, 1709 subjects (1043 European Americans, EA and 666 African Americans, AA) from families in the GeneSTAR study were included. In vitro platelet aggregation in response to collagen, ADP and epinephrine was measured at baseline and 14 days after aspirin therapy (81 mg/day). Targeted deep sequencing of PEAR1 in addition to 2kb of upstream and downstream of the gene was performed. Under an additive genetic model, the association of single variants of PEAR1 with platelet aggregation phenotypes were examined. Additionally, we examined the association between the burden of PEAR1 rare non-synonymous variants and platelet aggregation phenotypes. Of 532 variants identified through sequencing, the intron 1 variant, rs12041331, was significantly associated with all platelet aggregation phenotypes at baseline and after platelet inhibition with aspirin therapy. rs12566888, which is in linkage disequilibrium with rs12041331, was associated with platelet aggregation phenotypes but to a lesser extent. In the EA families, the burden of PEAR1 missense variants was associated with platelet aggregation after aspirin therapy when the platelets were stimulated with epinephrine (p = 0.0009) and collagen (p = 0.03). In AAs, the burden of PEAR1 missense variants was associated, to a lesser degree, with platelet aggregation in response to epinephrine (p = 0.02) and ADP (p = 0.04). Our study confirmed that the GWAS-identified variant, rs12041331, is the strongest variant associated with platelet aggregation both at baseline and after aspirin therapy in our GeneSTAR families in both races. We identified additional association of rare missense variants in PEAR1 with platelet aggregation following aspirin therapy. However, we observed a racial difference in the contribution of these rare variants to the platelet aggregation, most likely due to higher residual missing heritability of platelet aggregation after accounting for rs12041331 in the EAs compared to AAs.

Previous genome-wide association studies (GWAS) have identified several variants associated with platelet function phenotypes; however, the proportion of variance explained by the identified variants is mostly small. Rare coding variants, particularly those with high potential for impact on protein structure/function, may have substantial impact on phenotype but are difficult to detect by GWAS. The main purpose of this study was to identify low frequency or rare variants associated with platelet function using genotype data from the Illumina HumanExome Bead Chip. Three family-based cohorts of European ancestry, including ~4,000 total subjects, comprised the discovery cohort and two independent cohorts, one of European and one of African American ancestry, were used for replication. Optical aggregometry in platelet-rich plasma was performed in all the discovery cohorts in response to adenosine diphosphate (ADP), epinephrine, and collagen. Meta-analyses were performed using both gene-based and single nucleotide variant association methods. The gene-based meta-analysis identified a significant association (P = 7.13 × 10–7) between rare genetic variants in ANKRD26 and ADP-induced platelet aggregation. One of the ANKRD26 SNVs - rs191015656, encoding a threonine to isoleucine substitution predicted to alter protein structure/function, was replicated in Europeans. Aggregation increases of ~20–50% were observed in heterozygotes in all cohorts. Novel genetic signals in ABCG1 and HCP5 were also associated with platelet aggregation to ADP in meta-analyses, although only results for HCP5 could be replicated. The SNV in HCP5 intersects epigenetic signatures in CD41+ megakaryocytes suggesting a new functional role in platelet biology for HCP5. This is the first study to use gene-based association methods from SNV array genotypes to identify rare variants related to platelet function. The molecular mechanisms and pathophysiological relevance for the identified genetic associations requires further study.

The burden of asthma varies between ancestries, but GWAS have so far focused on mainly European ancestry populations. Here, Daya et al. perform GWAS for asthma in 14,654 individuals of African ancestry and, besides confirming previously known loci, identify two potentially African ancestry-specific loci.

Background: Genetic factors that influence kidney traits have been understudied for low frequency and ancestry-specific variants. Methods: We combined whole genome sequencing (WGS) data from 23,732 participants from 10 NHLBI Trans-Omics for Precision Medicine (TOPMed) Program multi-ethnic studies to identify novel loci for estimated glomerular filtration rate (eGFR). Participants included European, African, East Asian, and Hispanic ancestries. We applied linear mixed models using a genetic relationship matrix estimated from the WGS data and adjusted for age, sex, study, and ethnicity. Findings: When testing single variants, we identified three novel loci driven by low frequency variants more commonly observed in non-European ancestry (PRKAA2, rs180996919, minor allele frequency [MAF] 0.04%, P = 6.1 × 10−11; METTL8, rs116951054, MAF 0.09%, P = 4.5 × 10−9; and MATK, rs539182790, MAF 0.05%, P = 3.4 × 10−9). We also replicated two known loci for common variants (rs2461702, MAF=0.49, P = 1.2 × 10−9, nearest gene GATM, and rs71147340, MAF=0.34, P = 3.3 × 10−9, CDK12). Testing aggregated variants within a gene identified the MAF gene. A statistical approach based on local ancestry helped to identify replication samples for ancestry-specific variants. Interpretation: This study highlights challenges in studying variants influencing kidney traits that are low frequency in populations and more common in non-European ancestry.

Human induced pluripotent stem cell (hiPSC) lines have previously been generated through the NHLBI sponsored NextGen program at nine individual study sites. Here, we examined the structural integrity of 506 hiPSC lines as determined by copy number variations (CNVs). We observed that 149 hiPSC lines acquired 258 CNVs relative to donor DNA. We identified six recurrent regions of CNVs on chromosomes 1, 2, 3, 16 and 20 that overlapped with cancer associated genes. Furthermore, the genes mapping to regions of acquired CNVs show an enrichment in cancer related biological processes (IL6 production) and signaling cascades (JNK cascade & NFκB cascade). The genomic region of instability on chr20 (chr20q11.2) includes transcriptomic signatures for cancer associated genes such as ID1, BCL2L1, TPX2, PDRG1 and HCK. Of these HCK shows statistically significant differential expression between carrier and non-carrier hiPSC lines. Overall, while a low level of genomic instability was observed in the NextGen generated hiPSC lines, the observation of structural instability in regions with known cancer associated genes substantiates the importance of systematic evaluation of genetic variations in hiPSCs before using them as disease/research models.

BackgroundGenome-wide association studies have identified multiple genomic loci associated with coronary artery disease, but most are common variants in non-coding regions that provide limited information on causal genes and etiology of the disease. To overcome the limited scope that common variants provide, we focused our investigation on low-frequency and rare sequence variations primarily residing in coding regions of the genome.Methods and resultsUsing samples of individuals of European ancestry from ten cohorts within the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, both cross-sectional and prospective analyses were conducted to examine associations between genetic variants and myocardial infarction (MI), coronary heart disease (CHD), and all-cause mortality following these events. For prevalent events, a total of 27,349 participants of European ancestry, including 1831 prevalent MI cases and 2518 prevalent CHD cases were used. For incident cases, a total of 55,736 participants of European ancestry were included (3,031 incident MI cases and 5,425 incident CHD cases). There were 1,860 all-cause deaths among the 3,751 MI and CHD cases from six cohorts that contributed to the analysis of all-cause mortality. Single variant and gene-based analyses were performed separately in each cohort and then meta-analyzed for each outcome. A low-frequency intronic variant (rs988583) in PLCL1 was significantly associated with prevalent MI (OR = 1.80, 95% confidence interval: 1.43, 2.27; P = 7.12 × 10-7). We conducted gene-based burden tests for genes with a cumulative minor allele count (cMAC) ≥ 5 and variants with minor allele frequency (MAF) < 5%. TMPRSS5 and LDLRAD1 were significantly associated with prevalent MI and CHD, respectively, and RC3H2 and ANGPTL4 were significantly associated with incident MI and CHD, respectively. No loci were significantly associated with all-cause mortality following a MI or CHD event.ConclusionThis study identified one known locus (ANGPTL4) and four new loci (PLCL1, RC3H2, TMPRSS5, and LDLRAD1) associated with cardiovascular disease risk that warrant further investigation.

Most genome-wide association and fine-mapping studies to date have been conducted in individuals of European descent, and genetic studies of populations of Hispanic/Latino and African ancestry are limited. In addition, these populations have more complex linkage disequilibrium structure. In order to better define the genetic architecture of these understudied populations, we leveraged >100,000 phased sequences available from deep-coverage whole genome sequencing through the multi-ethnic NHLBI Trans-Omics for Precision Medicine (TOPMed) program to impute genotypes into admixed African and Hispanic/Latino samples with genome-wide genotyping array data. We demonstrated that using TOPMed sequencing data as the imputation reference panel improves genotype imputation quality in these populations, which subsequently enhanced gene-mapping power for complex traits. For rare variants with minor allele frequency (MAF) < 0.5%, we observed a 2.3- to 6.1-fold increase in the number of well-imputed variants, with 11-34% improvement in average imputation quality, compared to the state-of-the-art 1000 Genomes Project Phase 3 and Haplotype Reference Consortium reference panels. Impressively, even for extremely rare variants with minor allele count 86%. Subsequent association analyses of TOPMed reference panel-imputed genotype data with hematological traits (hemoglobin (HGB), hematocrit (HCT), and white blood cell count (WBC)) in ~21,600 African-ancestry and ~21,700 Hispanic/Latino individuals identified associations with two rare variants in the HBB gene (rs33930165 with higher WBC [p = 8.8x10-15] in African populations, rs11549407 with lower HGB [p = 1.5x10-12] and HCT [p = 8.8x10-10] in Hispanics/Latinos). By comparison, neither variant would have been genome-wide significant if either 1000 Genomes Project Phase 3 or Haplotype Reference Consortium reference panels had been used for imputation. Our findings highlight the utility of the TOPMed imputation reference panel for identification of novel rare variant associations not previously detected in similarly sized genome-wide studies of under-represented African and Hispanic/Latino populations.

BackgroundFibrinogen is an essential hemostatic factor and cardiovascular disease risk factor. Early attempts at evaluating the causal effect of fibrinogen on coronary heart disease (CHD) and myocardial infraction (MI) using Mendelian randomization (MR) used single variant approaches, and did not take advantage of recent genome-wide association studies (GWAS) or multi-variant, pleiotropy robust MR methodologies.Methods and findingsWe evaluated evidence for a causal effect of fibrinogen on both CHD and MI using MR. We used both an allele score approach and pleiotropy robust MR models. The allele score was composed of 38 fibrinogen-associated variants from recent GWAS. Initial analyses using the allele score used a meta-analysis of 11 European-ancestry prospective cohorts, free of CHD and MI at baseline, to examine incidence CHD and MI. We also applied 2 sample MR methods with data from a prevalent CHD and MI GWAS. Results are given in terms of the hazard ratio (HR) or odds ratio (OR), depending on the study design, and associated 95% confidence interval (CI). In single variant analyses no causal effect of fibrinogen on CHD or MI was observed. In multi-variant analyses using incidence CHD cases and the allele score approach, the estimated causal effect (HR) of a 1 g/L higher fibrinogen concentration was 1.62 (CI = 1.12, 2.36) when using incident cases and the allele score approach. In 2 sample MR analyses that accounted for pleiotropy, the causal estimate (OR) was reduced to 1.18 (CI = 0.98, 1.42) and 1.09 (CI = 0.89, 1.33) in the 2 most precise (smallest CI) models, out of 4 models evaluated. In the 2 sample MR analyses for MI, there was only very weak evidence of a causal effect in only 1 out of 4 models.ConclusionsA small causal effect of fibrinogen on CHD is observed using multi-variant MR approaches which account for pleiotropy, but not single variant MR approaches. Taken together, results indicate that even with large sample sizes and multi-variant approaches MR analyses still cannot exclude the null when estimating the causal effect of fibrinogen on CHD, but that any potential causal effect is likely to be much smaller than observed in epidemiological studies.

The Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) aims to better understand population genetics of the African diaspora. Here, it uses deeply sequenced whole-genomes to describe the impact of admixture and potential disease burden of deleterious variants.

Personalized medicine requires accurate and ethnicity-optimized reference genome panels. Here, the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA) evaluates typical variant filters and existing genome databases against newly sequenced African-ancestry populations.

Omega-6 (n-6) and omega-3 (n-3) long (≥ 20 carbon) chain polyunsaturated fatty acids (LC-PUFAs) play a critical role in human health and disease. Biosynthesis of LC-PUFAs from dietary 18 carbon PUFAs in tissues such as the liver is highly associated with genetic variation within the fatty acid desaturase (FADS) gene cluster, containing FADS1 and FADS2 that encode the rate-limiting desaturation enzymes in the LC-PUFA biosynthesis pathway. However, the molecular mechanisms by which FADS genetic variants affect LC-PUFA biosynthesis, and in which tissues, are unclear. The current study examined associations between common single nucleotide polymorphisms (SNPs) within the FADS gene cluster and FADS1 and FADS2 gene expression in 44 different human tissues (sample sizes ranging 70-361) from the Genotype-Tissue Expression (GTEx) Project. FADS1 and FADS2 expression were detected in all 44 tissues. Significant cis-eQTLs (within 1 megabase of each gene, False Discovery Rate, FDR

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

Previously, we have described our feeder-free, xeno-free approach to generate megakaryocytes (MKs) in culture from human induced pluripotent stem cells (iPSCs). Here, we focus specifically on the integrity of these MKs using: (1) genotype discordance between parent cell DNA to iPSC cell DNA and onward to the differentiated MK DNA; (2) genomic structural integrity using copy number variation (CNV); and (3) transcriptomic signatures of the derived MK lines compared to the iPSC lines. We detected a very low rate of genotype discordance; estimates were 0.0001%-0.01%, well below the genotyping error rate for our assay (0.37%). No CNVs were generated in the iPSCs that were subsequently passed on to the MKs. Finally, we observed highly biologically relevant gene sets as being upregulated in MKs relative to the iPSCs: platelet activation, blood coagulation, megakaryocyte development, platelet formation, platelet degranulation, and platelet aggregation. These data strongly support the integrity of the derived MK lines.