Explore the vast range of titles available.

For many complex traits, gene regulation is likely to play a crucial mechanistic role. How the genetic architectures of complex traits vary between populations and subsequent effects on genetic prediction are not well understood, in part due to the historical paucity of GWAS in populations of non-European ancestry. We used data from the MESA (Multi-Ethnic Study of Atherosclerosis) cohort to characterize the genetic architecture of gene expression within and between diverse populations. Genotype and monocyte gene expression were available in individuals with African American (AFA, n = 233), Hispanic (HIS, n = 352), and European (CAU, n = 578) ancestry. We performed expression quantitative trait loci (eQTL) mapping in each population and show genetic correlation of gene expression depends on shared ancestry proportions. Using elastic net modeling with cross validation to optimize genotypic predictors of gene expression in each population, we show the genetic architecture of gene expression for most predictable genes is sparse. We found the best predicted gene in each population, TACSTD2 in AFA and CHURC1 in CAU and HIS, had similar prediction performance across populations with R2 > 0.8 in each population. However, we identified a subset of genes that are well-predicted in one population, but poorly predicted in another. We show these differences in predictive performance are due to allele frequency differences between populations. Using genotype weights trained in MESA to predict gene expression in independent populations showed that a training set with ancestry similar to the test set is better at predicting gene expression in test populations, demonstrating an urgent need for diverse population sampling in genomics. Our predictive models and performance statistics in diverse cohorts are made publicly available for use in transcriptome mapping methods at https://github.com/WheelerLab/DivPop.

Dilated cardiomyopathy (DCM) is an important cause of heart failure and the leading indication for heart transplantation. Many rare genetic variants have been associated with DCM, but common variant studies of the disease have yielded few associated loci. As structural changes in the heart are a defining feature of DCM, we report a genome-wide association study of cardiac magnetic resonance imaging (MRI)-derived left ventricular measurements in 36,041 UK Biobank participants, with replication in 2184 participants from the Multi-Ethnic Study of Atherosclerosis. We identify 45 previously unreported loci associated with cardiac structure and function, many near well-established genes for Mendelian cardiomyopathies. A polygenic score of MRI-derived left ventricular end systolic volume strongly associates with incident DCM in the general population. Even among carriers of TTN truncating mutations, this polygenic score influences the size and function of the human heart. These results further implicate common genetic polymorphisms in the pathogenesis of DCM. Structural changes to the left ventricle are characteristic of dilated cardiomyopathy (DCM), a disease for which many rare genetic variants are known. Here, Pirruccello et al. report GWAS of seven cardiac MRI measurements in the left ventricle and describe shared loci and polygenic association with DCM.

Mitochondrial DNA copy number (mtDNA-CN), a measure of the number of mitochondrial genomes per cell, is a minimally invasive proxy measure for mitochondrial function and has been associated with several aging-related diseases. Although quantitative real-time PCR (qPCR) is the current gold standard method for measuring mtDNA-CN, mtDNA-CN can also be measured from genotyping microarray probe intensities and DNA sequencing read counts. To conduct a comprehensive examination on the performance of these methods, we use known mtDNA-CN correlates (age, sex, white blood cell count, Duffy locus genotype, incident cardiovascular disease) to evaluate mtDNA-CN calculated from qPCR, two microarray platforms, as well as whole genome (WGS) and whole exome sequence (WES) data across 1,085 participants from the Atherosclerosis Risk in Communities (ARIC) study and 3,489 participants from the Multi-Ethnic Study of Atherosclerosis (MESA). We observe mtDNA-CN derived from WGS data is significantly more associated with known correlates compared to all other methods (p < 0.001). Additionally, mtDNA-CN measured from WGS is on average more significantly associated with traits by 5.6 orders of magnitude and has effect size estimates 5.8 times more extreme than the current gold standard of qPCR. We further investigated the role of DNA extraction method on mtDNA-CN estimate reproducibility and found mtDNA-CN estimated from cell lysate is significantly less variable than traditional phenol-chloroform-isoamyl alcohol (p = 5.44x10-4) and silica-based column selection (p = 2.82x10-7). In conclusion, we recommend the field moves towards more accurate methods for mtDNA-CN, as well as re-analyze trait associations as more WGS data becomes available from larger initiatives such as TOPMed.

Clonal hematopoiesis of indeterminate potential (CHIP) is associated with an increased risk of cardiovascular diseases (CVDs), putatively via inflammasome activation. We pursued an inflammatory gene modifier scan for CHIP-associated CVD risk among 424,651 UK Biobank participants. We identified CHIP using whole-exome sequencing data of blood DNA and modeled as a composite, considering all driver genes together, as well as separately for common drivers (DNMT3A, TET2, ASXL1, and JAK2). We developed predicted gene expression scores for 26 inflammasome-related genes and assessed how they modify CHIP-associated CVD risk. We identified IL1RAP as a potential key molecule for CHIP-associated CVD risk across genes and increased AIM2 gene expression leading to heightened JAK2- and ASXL1-associated CVD risk. We show that CRISPR-induced Asxl1-mutated murine macrophages had a particularly heightened inflammatory response to AIM2 agonism, associated with an increased DNA damage response, as well as increased IL-10 secretion, mirroring a CVD-protective effect of IL10 expression in ASXL1 CHIP. Our study supports the role of inflammasomes in CHIP-associated CVD and provides evidence to support gene-specific strategies to address CHIP-associated CVD risk.

The intestinal microflora, typically equated with bacteria, influences diseases such as obesity and inflammatory bowel disease. Here, we show that the mammalian gut contains a rich fungal community that interacts with the immune system through the innate immune receptor Dectin-1. Mice lacking Dectin-1 exhibited increased susceptibility to chemically induced colitis, which was the result of altered responses to indigenous fungi. In humans, we identified a polymorphism in the gene for Dectin-1 (CLEC7A) that is strongly linked to a severe form of ulcerative colitis. Together, our findings reveal a eukaryotic fungal community in the gut (the \"mycobiome\") that coexists with bacteria and substantially expands the repertoire of organisms interacting with the intestinal immune system to influence health and disease.

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia. Mitochondrial DNA is known to exhibit heterogeneity of variants, even within a single cell. Here, the authors assessed this heteroplasmy of mitochondrial DNA within the UK Biobank cohort and showed that the presence of heteroplasmy and a functional score generated from heteroplasmic SNVs were associated with all-cause mortality and certain cancers.

Integrative approaches that simultaneously model multi-omics data have gained increasing popularity because they provide holistic system biology views of multiple or all components in a biological system of interest. Canonical correlation analysis (CCA) is a correlation-based integrative method designed to extract latent features shared between multiple assays by finding the linear combinations of features–referred to as canonical variables (CVs)–within each assay that achieve maximal across-assay correlation. Although widely acknowledged as a powerful approach for multi-omics data, CCA has not been systematically applied to multi-omics data in large cohort studies, which has only recently become available. Here, we adapted sparse multiple CCA (SMCCA), a widely-used derivative of CCA, to proteomics and methylomics data from the Multi-Ethnic Study of Atherosclerosis (MESA) and Jackson Heart Study (JHS). To tackle challenges encountered when applying SMCCA to MESA and JHS, our adaptations include the incorporation of the Gram-Schmidt (GS) algorithm with SMCCA to improve orthogonality among CVs, and the development of Sparse Supervised Multiple CCA (SSMCCA) to allow supervised integration analysis for more than two assays. Effective application of SMCCA to the two real datasets reveals important findings. Applying our SMCCA-GS to MESA and JHS, we identified strong associations between blood cell counts and protein abundance, suggesting that adjustment of blood cell composition should be considered in protein-based association studies. Importantly, CVs obtained from two independent cohorts also demonstrate transferability across the cohorts. For example, proteomic CVs learned from JHS, when transferred to MESA, explain similar amounts of blood cell count phenotypic variance in MESA, explaining 39.0% ~ 50.0% variation in JHS and 38.9% ~ 49.1% in MESA. Similar transferability was observed for other omics-CV-trait pairs. This suggests that biologically meaningful and cohort-agnostic variation is captured by CVs. We anticipate that applying our SMCCA-GS and SSMCCA on various cohorts would help identify cohort-agnostic biologically meaningful relationships between multi-omics data and phenotypic traits.

Angiotensinogen plays an essential role in maintaining circulatory homeostasis. AGT rs5050(T > G) has been identified as a regulator of the transcription of AGT mRNA, with differential expression between sexes. We sought to determine if rs5050(T > G), an estrogen response element, modifies the relationship between estrogen and angiotensinogen levels. rs5050(T > G) was genotyped, and plasma angiotensinogen levels were measured in 4,831 MESA participants, including postmenopausal women, on hormone therapy (n = 709) or not (n = 1,551), and 2,581 men. Linear regression models were employed to determine the associations of angiotensinogen with rs5050(T > G) allele dosage; and to evaluate whether rs5050(T > G) modifies the association between estradiol and angiotensinogen, with a main effect term and interaction term between rs5050(T > G)*estradiol. Estimated marginal means (EMMs) were used to further evaluate the effect of estradiol on angiotensinogen across different rs5050 alleles (T > G). rs5050TT had the highest median levels of angiotensinogen, followed by TG and GG. Adjusted main effect model showed positive associations between estradiol and angiotensinogen, with each rs5050T allele associated with 0.329 SD higher log-angiotensinogen levels (CI 95% 0.293, 0.365). The interaction rs5050(T > G)*estradiol was not significant, with EMMs exhibiting overlapping slope confidence intervals across genotypes. The proportion of the variance in angiotensinogen explained by modeling increases from 47.9% to 51.6% when including rs5050(T > G) or interation rs5050(T > G)*estradiol in the model. rs5050(T > G) is associated with circulating angiotensinogen levels, but rs5050(T > G) alleles do not influence the relationship between estradiol and angiotensinogen. This suggests that estrogen's effect on angiotensinogen regulation occurs independently of rs5050(T > G), despite its location within an estrogen-responsive element.

Biological age may be estimated by proteomic aging clocks (PACs). Previous published PACs were constructed either in smaller studies or mainly in white individuals, and they used proteomic measures from only one-time point. In this study, we created de novo PACs and compared their performance to published PACs at 2 different time points in the Atherosclerosis Risk in Communities (ARIC) study of white and black participants (around 75% white and 25% black). A total of 4,712 plasma proteins were measured using SomaScan in blood samples collected in 1990 to 1992 from 11,761 midlife participants (aged 46 to 70 years) and in 2011 to 2013 from 5,183 late-life participants (aged 66 to 90 years). The de novo ARIC PACs were constructed by training them against chronological age using elastic net regression in two-thirds of healthy participants in midlife and late life and validated in the remaining one-third of healthy participants at the corresponding time point. We also computed 3 published PACs. We estimated age acceleration for each PAC as residuals after regressing each PAC on chronological age. We also calculated the change in age acceleration from midlife to late life. We examined the associations of age acceleration and change in age acceleration with mortality through 2019 from all-cause, cardiovascular disease (CVD), cancer, and lower respiratory disease (LRD) using Cox proportional hazards regression in participants (irrespective of health) after excluding the training set. The model was adjusted for chronological age, smoking, body mass index (BMI), and other confounders. We externally validated the midlife PAC using the Multi-Ethnic Study of Atherosclerosis (MESA) Exam 1 data. The ARIC PACs had a slightly stronger correlation with chronological age than published PACs in healthy participants at each time point. Associations with mortality were similar for the ARIC PACs and published PACs. For late-life and midlife age acceleration for the ARIC PACs, respectively, hazard ratios (HRs) per 1 standard deviation were 1.65 and 1.38 (both p < 0.001) for all-cause mortality, 1.37 and 1.20 (both p < 0.001) for CVD mortality, 1.21 (p = 0.028) and 1.04 (p = 0.280) for cancer mortality, and 1.68 and 1.36 (both p < 0.001) for LRD mortality. For the change in age acceleration, HRs for all-cause, CVD, and LRD mortality were comparable to the HRs for late-life age acceleration. The association between the change in age acceleration and cancer mortality was not significant. The external validation of the midlife PAC in MESA showed significant associations with mortality, as observed for midlife participants in ARIC. The main limitation is that our PACs were constructed in midlife and late-life participants. It is unknown whether these PACs could be applied to young individuals. In this longitudinal study, we found that the ARIC PACs and published PACs were similarly associated with an increased risk of mortality. These findings suggested that PACs show promise as biomarkers of biological age. PACs may be serve as tools to predict mortality and evaluate the effect of anti-aging lifestyle and therapeutic interventions.

Background Statins are widely prescribed to lower plasma low-density lipoprotein cholesterol levels. Though statins reduce cardiovascular disease risk overall, statin efficacy varies, and some people experience adverse side effects while on statin treatment. Statins also have pleiotropic effects not directly related to their cholesterol-lowering properties, but the mechanisms are not well understood. To identify potential genetic modulators of clinical statin response, we looked for genetic variants associated with statin-induced changes in gene expression (differential eQTLs or deQTLs) in lymphoblastoid cell lines (LCLs) derived from participants of the Cholesterol and Pharmacogenetics (CAP) 40 mg/day 6-week simvastatin clinical trial. We exposed CAP LCLs to 2 μM simvastatin or control buffer for 24 h and performed polyA-selected, strand-specific RNA-seq. Statin-induced changes in gene expression from 259 European ancestry or 153 African American ancestry LCLs were adjusted for potential confounders prior to association with genotyped and imputed genetic variants within 1 Mb of each gene’s transcription start site. Results From the deQTL meta-analysis of the two ancestral populations, we identified significant cis-deQTLs for 15 genes ( TBC1D4 , MDGA1 , CHI3L2 , OAS1 , GATM , ASNSD1 , GLUL , TDRD12 , PPIP5K2 , OAS3 , SERPINB1 , ANKDD1A , DTD1 , CYFIP2 , and GSDME ), eight of which were significant in at least one of the ancestry subsets alone. We also conducted eQTL analyses of the endogenous (control-treated), statin-treated, and average of endogenous and statin-treated LCL gene expression levels. We identified eQTLs for approximately 6000 genes in each of the three (endogenous, statin-treated, and average) eQTL meta-analyses, with smaller numbers identified in the ancestral subsets alone. Conclusions Several of the genes in which we identified deQTLs have functions in human health and disease, such as defense from viruses, glucose regulation, and response to chemotherapy drugs. This suggests that DNA variation may play a role in statin effects on various health outcomes. These findings could prove useful to future studies aiming to assess benefit versus risk of statin treatment using individual genetic profiles.