Whole-Genome Sequencing Analysis of Human Metabolome in Multi-Ethnic Populations

Circulating metabolite levels may reflect the state of the human organism in health and disease, however, the genetic architecture of metabolites is not fully understood. We have performed a whole-genome sequencing association analysis of both common and rare variants in up to 11,840 multi-ethnic participants from five studies with up to 1666 circulating metabolites. We have discovered 1985 novel variant-metabolite associations, and validated 761 locus-metabolite associations reported previously. Seventy-nine novel variant-metabolite associations have been replicated, including three genetic loci located on the X chromosome that have demonstrated its involvement in metabolic regulation. Gene-based analysis have provided further support for seven metabolite-replicated loci pairs and their biologically plausible genes. Among those novel replicated variant-metabolite pairs, follow-up analyses have revealed that 26 metabolites have colocalized with 21 tissues, seven metabolite-disease outcome associations have been putatively causal, and 7 metabolites might be regulated by plasma protein levels. Our results have depicted the genetic contribution to circulating metabolite levels, providing additional insights into understanding human disease.

A multi-ethnic polygenic risk score is associated with hypertension prevalence and progression throughout adulthood

In a multi-stage analysis of 52,436 individuals aged 17-90 across diverse cohorts and biobanks, we train, test, and evaluate a polygenic risk score (PRS) for hypertension risk and progression. The PRS is trained using genome-wide association studies (GWAS) for systolic, diastolic blood pressure, and hypertension, respectively. For each trait, PRS is selected by optimizing the coefficient of variation (CV) across estimated effect sizes from multiple potential PRS using the same GWAS, after which the 3 trait-specific PRSs are combined via an unweighted sum called "PRSsum", forming the HTN-PRS. The HTN-PRS is associated with both prevalent and incident hypertension at 4-6 years of follow up. This association is further confirmed in age-stratified analysis. In an independent biobank of 40,201 individuals, the HTN-PRS is confirmed to be predictive of increased risk for coronary artery disease, ischemic stroke, type 2 diabetes, and chronic kidney disease.

Host and gut microbial tryptophan metabolism and type 2 diabetes: an integrative analysis of host genetics, diet, gut microbiome and circulating metabolites in cohort studies

OBJECTIVE

Tryptophan can be catabolised to various metabolites through host kynurenine and microbial indole pathways. We aimed to examine relationships of host and microbial tryptophan metabolites with incident type 2 diabetes (T2D), host genetics, diet and gut microbiota.

METHOD

We analysed associations between circulating levels of 11 tryptophan metabolites and incident T2D in 9180 participants of diverse racial/ethnic backgrounds from five cohorts. We examined host genome-wide variants, dietary intake and gut microbiome associated with these metabolites.

RESULTS

Tryptophan, four kynurenine-pathway metabolites (kynurenine, kynurenate, xanthurenate and quinolinate) and indolelactate were positively associated with T2D risk, while indolepropionate was inversely associated with T2D risk. We identified multiple host genetic variants, dietary factors, gut bacteria and their potential interplay associated with these T2D-relaetd metabolites. Intakes of fibre-rich foods, but not protein/tryptophan-rich foods, were the dietary factors most strongly associated with tryptophan metabolites. The fibre-indolepropionate association was partially explained by indolepropionate-associated gut bacteria, mostly fibre-using Firmicutes. We identified a novel association between a host functional LCT variant (determining lactase persistence) and serum indolepropionate, which might be related to a host gene-diet interaction on gut Bifidobacterium, a probiotic bacterium significantly associated with indolepropionate independent of other fibre-related bacteria. Higher milk intake was associated with higher levels of gut Bifidobacterium and serum indolepropionate only among genetically lactase non-persistent individuals.

CONCLUSION

Higher milk intake among lactase non-persistent individuals, and higher fibre intake were associated with a favourable profile of circulating tryptophan metabolites for T2D, potentially through the host-microbial cross-talk shifting tryptophan metabolism toward gut microbial indolepropionate production.

Exome sequence association study of levels and longitudinal change of cardiovascular risk factor phenotypes in European Americans and African Americans from the Atherosclerosis Risk in Communities Study

Cardiovascular disease (CVD) is responsible for 31% of all deaths worldwide. Among CVD risk factors are age, race, increased systolic blood pressure (BP), and dyslipidemia. Both BP and blood lipids levels change with age, with a dose-dependent relationship between the cumulative exposure to hyperlipidemia and the risk of CVD. We performed an exome sequence association study using longitudinal data with up to 7805 European Americans (EAs) and 3171 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study. We assessed associations of common (minor allele frequency > 5%) nonsynonymous and splice-site variants and gene-based sets of rare variants with levels and with longitudinal change of seven CVD risk factor phenotypes (BP traits: systolic BP, diastolic BP, pulse pressure; lipids traits: triglycerides, total cholesterol, high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C]). Furthermore, we investigated the relationship of the identified variants and genes with select CVD endpoints. We identified two novel genes: DCLK3 associated with the change of HDL-C levels in AAs and RAB7L1 associated with the change of LDL-C levels in EAs. RAB7L1 is further associated with an increased risk of heart failure in ARIC EAs. Investigation of the contribution of genetic factors to the longitudinal change of CVD risk factor phenotypes promotes our understanding of the etiology of CVD outcomes, stressing the importance of incorporating the longitudinal structure of the cohort data in future analyses.

Genome-wide association study of serum metabolites in the African American Study of Kidney Disease and Hypertension

The genome-wide association study (GWAS) is a powerful means to study genetic determinants of disease traits and generate insights into disease pathophysiology. To date, few GWAS of circulating metabolite levels have been performed in African Americans with chronic kidney disease. Hypothesizing that novel genetic-metabolite associations may be identified in a unique population of African Americans with a lower glomerular filtration rate (GFR), we conducted a GWAS of 652 serum metabolites in 619 participants (mean measured glomerular filtration rate 45 mL/min/1.73m2) in the African American Study of Kidney Disease and Hypertension, a clinical trial of blood pressure lowering and antihypertensive medication in African Americans with chronic kidney disease. We identified 42 significant variant metabolite associations. Twenty associations had been previously identified in published GWAS, and eleven novel associations were replicated in a separate cohort of 818 African Americans with genetic and metabolomic data from the Atherosclerosis Risk in Communities Study. The replicated novel variant-metabolite associations comprised eight metabolites and eleven distinct genomic loci. Nine of the replicated associations represented clear enzyme-metabolite interactions, with high expression in the kidneys as well as the liver. Three loci (ACY1, ACY3, and NAT8) were associated with a common pool of metabolites, acetylated amino acids, but with different individual affinities. Thus, extensive metabolite profiling in an African American population with chronic kidney disease aided identification of novel genome-wide metabolite associations, providing clues about substrate specificity and the key roles of enzymes in modulating systemic levels of metabolites.

Discovery of rare variants associated with blood pressure regulation through meta-analysis of 1.3 million individuals

Genetic studies of blood pressure (BP) to date have mainly analyzed common variants (minor allele frequency > 0.05). In a meta-analysis of up to ~1.3 million participants, we discovered 106 new BP-associated genomic regions and 87 rare (minor allele frequency ≤ 0.01) variant BP associations (P < 5 × 10-8), of which 32 were in new BP-associated loci and 55 were independent BP-associated single-nucleotide variants within known BP-associated regions. Average effects of rare variants (44% coding) were ~8 times larger than common variant effects and indicate potential candidate causal genes at new and known loci (for example, GATA5 and PLCB3). BP-associated variants (including rare and common) were enriched in regions of active chromatin in fetal tissues, potentially linking fetal development with BP regulation in later life. Multivariable Mendelian randomization suggested possible inverse effects of elevated systolic and diastolic BP on large artery stroke. Our study demonstrates the utility of rare-variant analyses for identifying candidate genes and the results highlight potential therapeutic targets.

A Genome-wide Association Study Discovers 46 Loci of the Human Metabolome in the Hispanic Community Health Study/Study of Latinos

Variation in levels of the human metabolome reflect changes in homeostasis, providing a window into health and disease. The genetic impact on circulating metabolites in Hispanics, a population with high cardiometabolic disease burden, is largely unknown. We conducted genome-wide association analyses on 640 circulating metabolites in 3,926 Hispanic Community Health Study/Study of Latinos participants. The estimated heritability for 640 metabolites ranged between 0%-54% with a median at 2.5%. We discovered 46 variant-metabolite pairs (p value < 1.2 × 10-10, minor allele frequency ≥ 1%, proportion of variance explained [PEV] mean = 3.4%, PEVrange = 1%-22%) with generalized effects in two population-based studies and confirmed 301 known locus-metabolite associations. Half of the identified variants with generalized effect were located in genes, including five nonsynonymous variants. We identified co-localization with the expression quantitative trait loci at 105 discovered and 151 known loci-metabolites sets. rs5855544, upstream of SLC51A, was associated with higher levels of three steroid sulfates and co-localized with expression levels of SLC51A in several tissues. Mendelian randomization (MR) analysis identified several metabolites associated with coronary heart disease (CHD) and type 2 diabetes. For example, two variants located in or near CYP4F2 (rs2108622 and rs79400241, respectively), involved in vitamin E metabolism, were associated with the levels of octadecanedioate and vitamin E metabolites (gamma-CEHC and gamma-CEHC glucuronide); MR analysis showed that genetically high levels of these metabolites were associated with lower odds of CHD. Our findings document the genetic architecture of circulating metabolites in an underrepresented Hispanic/Latino community, shedding light on disease etiology.

Unraveling the functional role of the orphan solute carrier, SLC22A24 in the transport of steroid conjugates through metabolomic and genome-wide association studies

Variation in steroid hormone levels has wide implications for health and disease. The genes encoding the proteins involved in steroid disposition represent key determinants of interindividual variation in steroid levels and ultimately, their effects. Beginning with metabolomic data from genome-wide association studies (GWAS), we observed that genetic variants in the orphan transporter, SLC22A24 were significantly associated with levels of androsterone glucuronide and etiocholanolone glucuronide (sentinel SNPs p-value <1x10-30). In cells over-expressing human or various mammalian orthologs of SLC22A24, we showed that steroid conjugates and bile acids were substrates of the transporter. Phylogenetic, genomic, and transcriptomic analyses suggested that SLC22A24 has a specialized role in the kidney and appears to function in the reabsorption of organic anions, and in particular, anionic steroids. Phenome-wide analysis showed that functional variants of SLC22A24 are associated with human disease such as cardiovascular diseases and acne, which have been linked to dysregulated steroid metabolism. Collectively, these functional genomic studies reveal a previously uncharacterized protein involved in steroid homeostasis, opening up new possibilities for SLC22A24 as a pharmacological target for regulating steroid levels.

Remnant-Like Particle Cholesterol, Low-Density Lipoprotein Triglycerides, and Incident Cardiovascular Disease

BACKGROUND

Hypertriglyceridemia is associated with increased remnant-like particle cholesterol (RLP-C) and triglycerides in low-density lipoprotein (LDL-TG). Recent studies have focused on atherogenicity of RLP-C, with few data on LDL-TG.

OBJECTIVES

The aim of this study was to examine associations of RLP-C and LDL-TG with incident cardiovascular disease (CVD) events and genetic variants in the ARIC (Atherosclerosis Risk In Communities) study.

METHODS

Fasting plasma RLP-C and LDL-TG levels were measured in 9,334 men and women without prevalent CVD. Participants were followed for incident CVD events (coronary heart disease and ischemic stroke) for up to 16 years. Associations between LDL-TG and RLP-C levels and genetic variants were assessed by whole-exome sequencing using single-variant analysis for common variants and gene-based burden tests for rare variants; both an unbiased and a candidate gene approach were explored.

RESULTS

RLP-C and LDL-TG levels were correlated with triglyceride levels (r = 0.85 and r = 0.64, p < 0.0001). In minimally adjusted analyses, RLP-C and LDL-TG were associated with CVD risk, but in models adjusted for traditional risk factors including lipids, only LDL-TG was associated with incident CHD (hazard ratio: 1.28; 95% confidence interval: 1.10 to 1.50) and stroke (hazard ratio: 1.47; 95% confidence interval: 1.13 to 1.92). A common APOE variant, rs7412, had the strongest association with LDL-TG and RLP-C (p < 5 × 10-8).

CONCLUSIONS

RLP-C and LDL-TG levels were predictive of CVD and associated with APOE variants. LDL-TG may represent a marker of dysfunctional remnant lipoprotein metabolism associated with increased CVD risk. Further research is needed to determine whether LDL-TG plays a causal role in CVD and may be a target for therapy.

Whole-genome sequencing study of serum peptide levels: the Atherosclerosis Risk in Communities study

Oligopeptides are important markers of protein metabolism, as they are cleaved from larger polypeptides and proteins. Genetic association studies may help elucidate their origin and function. In 1,552 European Americans and 1,872 African Americans of the Atherosclerosis Risk in Communities study, we performed whole-genome and whole-exome sequencing and measured serum levels of 25 peptides. Common variants (minor allele frequency > 5%) were analysed individually. We grouped low-frequency variants (minor allele frequency ≤ 5%) by a genome-wide sliding window using region-based aggregate tests. Furthermore, low-frequency regulatory variants were grouped by gene, as were functional coding variants. All analyses were performed separately in each ancestry group and then meta-analysed. We identified 22 common variant associations with peptide levels (P-value < 4.2 × 10-10), including 16 novel gene-peptide pairs. Notably, variants in kinin-kallikrein genes KNG1, F12, KLKB1, and ACE were associated with several different peptides. Variants in KLKB1 and ACE were associated with a fragment of complement component 3f. Both common variants and low-frequency coding variants in CPN1 were associated with a fibrinogen cleavage peptide. Four sliding windows were significantly associated with peptide levels (P-value < 4.2 × 10-10). Our results highlight the importance of the kinin-kallikrein system in the regulation of serum peptide levels, strengthen the evidence for a broad link between the kinin-kallikrein and complement systems, and suggest a role of CPN1 in the conversion of fibrinogen to fibrin.

Whole genome sequence analysis of serum amino acid levels

Background

Blood levels of amino acids are important biomarkers of disease and are influenced by synthesis, protein degradation, and gene–environment interactions. Whole genome sequence analysis of amino acid levels may establish a paradigm for analyzing quantitative risk factors.

Results

In a discovery cohort of 1872 African Americans and a replication cohort of 1552 European Americans we sequenced exons and whole genomes and measured serum levels of 70 amino acids. Rare and low-frequency variants (minor allele frequency ≤5%) were analyzed by three types of aggregating motifs defined by gene exons, regulatory regions, or genome-wide sliding windows. Common variants (minor allele frequency >5%) were analyzed individually. Over all four analysis strategies, 14 gene–amino acid associations were identified and replicated. The 14 loci accounted for an average of 1.8% of the variance in amino acid levels, which ranged from 0.4 to 9.7%. Among the identified locus–amino acid pairs, four are novel and six have been reported to underlie known Mendelian conditions. These results suggest that there may be substantial genetic effects on amino acid levels in the general population that may underlie inborn errors of metabolism. We also identify a predicted promoter variant in AGA (the gene that encodes aspartylglucosaminidase) that is significantly associated with asparagine levels, with an effect that is independent of any observed coding variants.

Conclusions

These data provide insights into genetic influences on circulating amino acid levels by integrating -omic technologies in a multi-ethnic population. The results also help establish a paradigm for whole genome sequence analysis of quantitative traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-016-1106-x) contains supplementary material, which is available to authorized users.

Abstract 3076: p53 homologue p63 trans-activates epidermal growth factor (EGF) receptor and 14-3-3α and thus promotes tumorigenesis and chemoresistance in pancreatic cancer

Pancreatic cancer is the 4th leading cause of cancer-related mortality in the United States. Despite these statistics, progress in therapy of pancreatic cancer has been slow. p53 tumor suppressor pathway is altered in 50-75% of pancreatic cancer tumors. However, the role of p53 homologs p63 and p73 has not been well studied. We previously reported that p63 is expressed in pancreatic cancer cell lines and contributes to their motility and invasion. Here we tested expression of p63 in resected tumors and investigated its role as a transcription factor in pancreatic cancer. Expression of p63 was determined by immunohistochemistry in 20 resected tumor samples (15 - adenocarcinoma, 5 - squamous) using pan-p63 (clone 4A4) and ΔNp63-specific antibodies; and by immunoblotting, real-time quantitative and semi-quantitative PCR in 5 pancreatic cancer cell lines. p63 levels were manipulated employing transient overexpression and lentiviral approaches. We studied cell doubling times, motility in a wound assay and invasiveness in Matrigel chambers in the presence or absence of epidermal growth factor (EGF). Functional reporter and chromatin immunoprecipitation assays were performed to study the effects of p63 on selected promoters. We found that the p53 homolog p63 was expressed in resected pancreatic cancer samples. Pancreatic tumors with squamous differentiation showed strong diffuse expression of p63, whereas adenocarcinomas without squamous component had a focal expression of p63 (ΔNp63). In pancreatic cancer cell lines, p63 protein levels were high in 3/5 tested cell lines and ΔNp63α was the predominantly expressed variant. By contrast, TAp63 was not detected while p73 was barely detectable in 2/5 cell lines. Ectopic expression of ΔNp63α in PANC-1 cells (which had low endogenous p63) and shRNA-mediated knockdown in T3M4 cells (with high endogenous p63) indicated that ΔNp63α promoted anchorage-dependent and independent growth, motility and invasion, and enhanced resistance to cisplatin-induced apoptosis. Motogenic effects of ΔNp63α were augmented in presence of EGF. ΔNp63α augmented the EGF receptor expression in pancreatic cancer cells via transcriptional regulation and potentiated EGF-mediated activation of ERK, Akt and JNK signaling. 14-3-3α transcription was also positively regulated by ΔNp63α and we have previously shown that 14-3-3α contributes to chemoresistance in pancreatic cancer. Thus, p53 homolog ΔNp63α is expressed in pancreatic cancer, where it enhances tumor progression and chemoresistance via trans-activation of EGF receptor and 14-3-3α.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 3076. doi:1538-7445.AM2012-3076