CSF proteomic profiling with amyloid/tau positivity identifies distinctive sex-different alteration of multiple proteins involved in Alzheimer's disease.

In Alzheimer's disease (AD), the most common cause of dementia, females have higher prevalence and faster progression, but sex-specific molecular findings in AD are limited. Here, we comprehensively examined and validated 7,006 aptamers targeting 6,162 proteins in cerebral spinal fluid (CSF) from 2,077 amyloid/tau positive cases and controls to identify sex-specific proteomic signatures of AD. In discovery (N=1,766), we identified 330 male-specific and 121 female-specific proteomic alternations in CSF (FDR <0.05). These sex-specific proteins strongly predicted amyloid/tau positivity (AUC=0.98 in males; 0.99 in females), significantly higher than those with age, sex, and APOE-ε4 (AUC=0.85). The identified sex-specific proteins were well validated (r≥0.5) in the Stanford study (N=108) and Emory study (N=148). Biological follow-up of these proteins led to sex differences in cell-type specificity, pathways, interaction networks, and drug targets. Male-specific proteins, enriched in astrocytes and oligodendrocytes, were involved in postsynaptic and axon-genesis. The male network exhibited direct connections among 152 proteins and highlighted PTEN, NOTCH1, FYN, and MAPK8 as hubs. Drug target suggested melatonin (used for sleep-wake cycle regulation), nabumetone (used for pain), daunorubicin, and verteporfin for treating AD males. In contrast, female-specific proteins, enriched in neurons, were involved in phosphoserine residue binding including cytokine activities. The female network exhibits strong connections among 51 proteins and highlighted JUN and 14-3-3 proteins (YWHAG and YWHAZ) as hubs. Drug target suggested biperiden (for muscle control of Parkinson's disease), nimodipine (for cerebral vasospasm), quinostatin and ethaverine for treating AD females. Together, our findings provide mechanistic understanding of sex differences for AD risk and insights into clinically translatable interventions.

Merged Affinity Network Association Clustering: Joint multi-omic/clinical clustering to identify disease endotypes.

Although clinical and laboratory data have long been used to guide medical practice, this information is rarely integrated with multi-omic data to identify endotypes. We present Merged Affinity Network Association Clustering (MANAclust), a coding-free, automated pipeline enabling integration of categorical and numeric data spanning clinical and multi-omic profiles for unsupervised clustering to identify disease subsets. Using simulations and real-world data from The Cancer Genome Atlas, we demonstrate that MANAclust's feature selection algorithms are accurate and outperform competitors. We also apply MANAclust to a clinically and multi-omically phenotyped asthma cohort. MANAclust identifies clinically and molecularly distinct clusters, including heterogeneous groups of "healthy controls" and viral and allergy-driven subsets of asthmatic subjects. We also find that subjects with similar clinical presentations have disparate molecular profiles, highlighting the need for additional testing to uncover asthma endotypes. This work facilitates data-driven personalized medicine through integration of clinical parameters with multi-omics. MANAclust is freely available at https://bitbucket.org/scottyler892/manaclust/src/master/.

Whole-Exome Sequencing and hiPSC Cardiomyocyte Models Identify MYRIP, TRAPPC11, and SLC27A6 of Potential Importance to Left Ventricular Hypertrophy in an African Ancestry Population.

Background: Indices of left ventricular (LV) structure and geometry represent useful intermediate phenotypes related to LV hypertrophy (LVH), a predictor of cardiovascular (CV) disease (CVD) outcomes. Methods and Results: We conducted an exome-wide association study of LV mass (LVM) adjusted to height2.7, LV internal diastolic dimension (LVIDD), and relative wall thickness (RWT) among 1,364 participants of African ancestry (AAs) in the Hypertension Genetic Epidemiology Network (HyperGEN). Both single-variant and gene-based sequence kernel association tests were performed to examine whether common and rare coding variants contribute to variation in echocardiographic traits in AAs. We then used a data-driven procedure to prioritize and select genes for functional validation using a human induced pluripotent stem cell cardiomyocyte (hiPSC-CM) model. Three genes [myosin VIIA and Rab interacting protein (MYRIP), trafficking protein particle complex 11 (TRAPPC11), and solute carrier family 27 member 6 (SLC27A6)] were prioritized based on statistical significance, variant functional annotations, gene expression in the hiPSC-CM model, and prior biological evidence and were subsequently knocked down in the hiPSC-CM model. Expression profiling of hypertrophic gene markers in the knockdowns suggested a decrease in hypertrophic expression profiles. MYRIP knockdowns showed a significant decrease in atrial natriuretic factor (NPPA) and brain natriuretic peptide (NPPB) expression. Knockdowns of the heart long chain fatty acid (FA) transporter SLC27A6 resulted in downregulated caveolin 3 (CAV3) expression, which has been linked to hypertrophic phenotypes in animal models. Finally, TRAPPC11 knockdown was linked to deficient calcium handling. Conclusions: The three genes are biologically plausible candidates that provide new insight to hypertrophic pathways.

Network study of nasal transcriptome profiles reveals master regulator genes of asthma.

BACKGROUND: Nasal transcriptomics can provide an accessible window into asthma pathobiology. OBJECTIVE: Our goal was to move beyond gene signatures of asthma to identify master regulator genes that causally regulate genes associated with asthma phenotypes. METHODS: We recruited 156 children with severe persistent asthma and controls for nasal transcriptome profiling and applied network-based and probabilistic causal methods to identify severe asthma genes and their master regulators. We then took the same approach in an independent cohort of 190 adults with mild/moderate asthma and controls to identify mild/moderate asthma genes and their master regulators. Comparative analysis of the master regulator genes followed by validation testing in independent children with severe asthma (n = 21) and mild/moderate asthma (n = 154) was then performed. RESULTS: Nasal gene signatures for severe persistent asthma and for mild/moderate persistent asthma were identified; both were found to be enriched in coexpression network modules for ciliary function and inflammatory response. By applying probabilistic causal methods to these gene signatures and validation testing in independent cohorts, we identified (1) a master regulator gene common to asthma across severity and ages (FOXJ1); (2) master regulator genes of severe persistent asthma in children (LRRC23, TMEM231, CAPS, PTPRC, and FYB); and (3) master regulator genes of mild/moderate persistent asthma in children and adults (C1orf38 and FMNL1). The identified master regulators were statistically inferred to causally regulate the expression of downstream genes that modulate ciliary function and inflammatory response to influence asthma. CONCLUSION: The identified master regulator genes of asthma provide a novel path forward to further uncovering asthma mechanisms and therapy.

Endotoxin, food allergen sensitization, and food allergy: A complementary epidemiologic and experimental study.

BACKGROUND: Household endotoxin levels have been variably associated with risk for asthma and atopy. METHODS: We studied participants from the 2005-2006 National Health and Nutrition Examination Survey (NHANES, n = 6963), a large cohort representative of the US population (aged 1-84 years). We built logistic regression models to test for associations between house dust endotoxin and sensitization to specific foods (milk, egg, and peanut). To experimentally explore the detected epidemiologic associations, peripheral blood mononuclear cells (PBMCs) were collected from 21 children (aged 1-19 years) mono-food allergic (ie, sensitized and clinically reactive) to milk, egg, or peanut and nonallergic controls for stimulation with endotoxin and secreted cytokine measurement. For each food allergy, linear mixed-effects models were built to test the association between endotoxin stimulation and cytokine level. RESULTS: Among NHANES subjects, the geometric mean household endotoxin level was 15.5 EU/mg (GSE 0.5). Prevalence of food allergen sensitization (sIgE ≥ 0.35 kUA /L) varied by food: milk 5.7%, egg 4.0%, and peanut 7.9%. In models adjusted for potential confounders (age, race, country of birth, total people per household, US region, and history of wheezing in the past year), household endotoxin level was associated with sensitization to milk (OR 1.7, 95% CI 1.2-2.1) and egg (OR 1.4, 95% CI 1.01-1.9), but not peanut (OR 0.98, 95% CI 0.8-1.2). Interferon-γ levels of endotoxin-stimulated PBMCs from children allergic to milk or egg, but not peanut, were significantly lower compared to controls in linear mixed-effects models adjusted for repeated measures, experimental variables, age, and inter-individual variability (P-values .007, .018, and .058, respectively). CONCLUSION: Higher household endotoxin is associated with increased odds of milk and egg sensitization. Altered cytokine responsiveness to endotoxin is also observed in PBMCs from individuals with milk and egg allergy.

Dual transcriptomic and epigenomic study of reaction severity in peanut-allergic children.

BACKGROUND: Unexpected allergic reactions to peanut are the most common cause of fatal food-related anaphylaxis. Mechanisms underlying the variable severity of peanut-allergic reactions remain unclear. OBJECTIVES: We sought to expand mechanistic understanding of reaction severity in peanut allergy. METHODS: We performed an integrated transcriptomic and epigenomic study of peanut-allergic children as they reacted in vivo during double-blind, placebo-controlled peanut challenges. We integrated whole-blood transcriptome and CD4+ T-cell epigenome profiles to identify molecular signatures of reaction severity (ie, how severely a peanut-allergic child reacts when exposed to peanut). A threshold-weighted reaction severity score was calculated for each subject based on symptoms experienced during peanut challenge and the eliciting dose. Through linear mixed effects modeling, network construction, and causal mediation analysis, we identified genes, CpGs, and their interactions that mediate reaction severity. Findings were replicated in an independent cohort. RESULTS: We identified 318 genes with changes in expression during the course of reaction associated with reaction severity, and 203 CpG sites with differential DNA methylation associated with reaction severity. After replicating these findings in an independent cohort, we constructed interaction networks with the identified peanut severity genes and CpGs. These analyses and leukocyte deconvolution highlighted neutrophil-mediated immunity. We identified NFKBIA and ARG1 as hubs in the networks and 3 groups of interacting key node CpGs and peanut severity genes encompassing immune response, chemotaxis, and regulation of macroautophagy. In addition, we found that gene expression of PHACTR1 and ZNF121 causally mediates the association between methylation at corresponding CpGs and reaction severity, suggesting that methylation may serve as an anchor upon which gene expression modulates reaction severity. CONCLUSIONS: Our findings enhance current mechanistic understanding of the genetic and epigenetic architecture of reaction severity in peanut allergy.

Genome-wide meta-analysis of SNP and antihypertensive medication interactions on left ventricular traits in African Americans.

BACKGROUND: Left ventricular (LV) hypertrophy affects up to 43% of African Americans (AAs). Antihypertensive treatment reduces LV mass (LVM). However, interindividual variation in LV traits in response to antihypertensive treatments exists. We hypothesized that genetic variants may modify the association of antihypertensive treatment class with LV traits measured by echocardiography. METHODS: We evaluated the main effects of the three most common antihypertensive treatments for AAs as well as the single nucleotide polymorphism (SNP)-by-drug interaction on LVM and relative wall thickness (RWT) in 2,068 participants across five community-based cohorts. Treatments included thiazide diuretics (TDs), angiotensin converting enzyme inhibitors (ACE-Is), and dihydropyridine calcium channel blockers (dCCBs) and were compared in a pairwise manner. We performed fixed effects inverse variance weighted meta-analyses of main effects of drugs and 2.5 million SNP-by-drug interaction estimates. RESULTS: We observed that dCCBs versus TDs were associated with higher LVM after adjusting for covariates (p = 0.001). We report three SNPs at a single locus on chromosome 20 that modified the association between RWT and treatment when comparing dCCBs to ACE-Is with consistent effects across cohorts (smallest p = 4.7 × 10-8 , minor allele frequency range 0.09-0.12). This locus has been linked to LV hypertrophy in a previous study. A marginally significant locus in BICD1 (rs326641) was validated in an external population. CONCLUSIONS: Our study identified one locus having genome-wide significant SNP-by-drug interaction effect on RWT among dCCB users in comparison to ACE-I users. Upon additional validation in future studies, our findings can enhance the precision of medical approaches in hypertension treatment.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals.

BACKGROUND: Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans. OBJECTIVE: The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes. METHODS: A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes. RESULTS: The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs. CONCLUSIONS: We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.

Epigenome-wide association study of metabolic syndrome in African-American adults.

Background: The high prevalence of obesity among US adults has resulted in significant increases in associated metabolic disorders such as diabetes, dyslipidemia, and high blood pressure. Together, these disorders constitute metabolic syndrome, a clinically defined condition highly prevalent among African-Americans. Identifying epigenetic alterations associated with metabolic syndrome may provide additional information regarding etiology beyond current evidence from genome-wide association studies. Methods: Data on metabolic syndrome and DNA methylation was assessed on 614 African-Americans from the Hypertension Genetic Epidemiology Network (HyperGEN) study. Metabolic syndrome was defined using the joint harmonized criteria, and DNA methylation was assessed using the Illumina HumanMethylation450K Bead Chip assay on DNA extracted from buffy coat. Linear mixed effects regression models were used to examine the association between CpG methylation at > 450,000 CpG sites and metabolic syndrome adjusted for study covariates. Replication using DNA from a separate sample of 69 African-Americans, as well as meta-analysis combining both cohorts, was conducted. Results: Two differentially methylated CpG sites in the IGF2BP1 gene on chromosome 17 (cg06638433; p value = 3.10 × 10- 7) and the ABCG1 gene on chromosome 21 (cg06500161; p value = 2.60 × 10- 8) were identified. Results for the ABCG1 gene remained statistically significant in the replication dataset and meta-analysis. Conclusion: Metabolic syndrome was consistently associated with increased methylation in the ABCG1 gene in the discovery and replication datasets, a gene that encodes a protein in the ATP-binding cassette transporter family and is involved in intra- and extra-cellular signaling and lipid transport.

Association of Methylation Signals With Incident Coronary Heart Disease in an Epigenome-Wide Assessment of Circulating Tumor Necrosis Factor α.

Importance: Tumor necrosis factor α (TNF-α) is a proinflammatory cytokine with manifold consequences for mammalian pathophysiology, including cardiovascular disease. A deeper understanding of TNF-α biology may enhance treatment precision. Objective: To conduct an epigenome-wide analysis of blood-derived DNA methylation and TNF-α levels and to assess the clinical relevance of findings. Design, Setting, and Participants: This meta-analysis assessed epigenome-wide associations in circulating TNF-α concentrations from 5 cohort studies and 1 interventional trial, with replication in 3 additional cohort studies. Follow-up analyses investigated associations of identified methylation loci with gene expression and incident coronary heart disease; this meta-analysis included 11 461 participants who experienced 1895 coronary events. Exposures: Circulating TNF-α concentration. Main Outcomes and Measures: DNA methylation at approximately 450 000 loci, neighboring DNA sequence variation, gene expression, and incident coronary heart disease. Results: The discovery cohort included 4794 participants, and the replication study included 816 participants (overall mean [SD] age, 60.7 [8.5] years). In the discovery stage, circulating TNF-α levels were associated with methylation of 7 cytosine-phosphate-guanine (CpG) sites, 3 of which were located in or near DTX3L-PARP9 at cg00959259 (ß [SE] = -0.01 [0.003]; P = 7.36 × 10-8), cg08122652 (ß [SE] = -0.008 [0.002]; P = 2.24 × 10-7), and cg22930808(ß [SE] = -0.01 [0.002]; P = 6.92 × 10-8); NLRC5 at cg16411857 (ß [SE] = -0.01 [0.002]; P = 2.14 × 10-13) and cg07839457 (ß [SE] = -0.02 [0.003]; P = 6.31 × 10-10); or ABO, at cg13683939 (ß [SE] = 0.04 [0.008]; P = 1.42 × 10-7) and cg24267699 (ß [SE] = -0.009 [0.002]; P = 1.67 × 10-7), after accounting for multiple testing. Of these, negative associations between TNF-α concentration and methylation of 2 loci in NLRC5 and 1 in DTX3L-14 PARP9 were replicated. Replicated TNF-α-linked CpG sites were associated with 9% to 19% decreased risk of incident coronary heart disease per 10% higher methylation per CpG site (cg16411857: hazard ratio [HR], 0.86; 95% CI, 0.78-1.95; P = .003; cg07839457: HR, 0.89; 95% CI, 0.80-0.94; P = 3.1 × 10-5; cg00959259: HR, 0.91; 95% CI, 0.84-0.97; P = .002; cg08122652: HR, 0.81; 95% CI, 0.74-0.89; P = 2.0 × 10-5). Conclusions and Relevance: We identified and replicated novel epigenetic correlates of circulating TNF-α concentration in blood samples and linked these loci to coronary heart disease risk, opening opportunities for validation and therapeutic applications.

Whole exome analyses to examine the impact of rare variants on left ventricular traits in African American participants from the HyperGEN and GENOA studies.

Left ventricular (LV) hypertrophy, highest in prevalence among African Americans, is an established risk factor heart failure. Several genome wide association studies have identified common variants associated with LV-related quantitative-traits in African Americans. To date, however, the effect of rare variants on these traits has not been extensively studied, especially in minority groups. We therefore investigated the association between rare variants and LV traits among 1,934 African Americans using exome chip data from the Hypertension Genetic Epidemiology Network (HyperGEN) study, with replication in 1,090 African American from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. We used single-variant analyses and gene-based tests to investigate the association between 86,927 variants and six structural and functional LV traits including LV mass, LV internal dimension-diastole, relative wall thickness, left atrial dimension (LAD), fractional shortening (FS), and the ratio of LV early-to-late transmitral velocity (E/A ratio). Only rare variants (MAF <1% and <5%) were considered in gene-based analyses. In gene-based analyses, we found a statistically significant association between potassium voltage-gated channel subfamily H member 4 (KCNH4) and E/A ratio (P=8.7*10-8 using a burden test). Endonuclease G (ENDOG) was associated with LAD using the Madsen Browning weighted burden (MB) test (P=1.4*10-7). Neither gene result was replicated in GENOA, but the direction of effect of single variants in common was comparable. G protein-coupled receptor 55 (GPR55) was marginally associated with LAD in HyperGEN (P=3.2*10-5 using the MB test) and E/A ratio in GENOA, but with opposing directions of association for variants in common (P=0.03 for the MB test). No single variant was statistically significantly associated with any trait after correcting for multiple testing. The findings in this study highlight the potential cumulative contributions of rare variants to LV traits which, if validated, could improve our understanding of heart failure in African Americans.

Association of DNA Methylation at CPT1A Locus with Metabolic Syndrome in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) Study.

In this study, we conducted an epigenome-wide association study of metabolic syndrome (MetS) among 846 participants of European descent in the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN). DNA was isolated from CD4+ T cells and methylation at ~470,000 cytosine-phosphate-guanine dinucleotide (CpG) pairs was assayed using the Illumina Infinium HumanMethylation450 BeadChip. We modeled the percentage methylation at individual CpGs as a function of MetS using linear mixed models. A Bonferroni-corrected P-value of 1.1 x 10(-7) was considered significant. Methylation at two CpG sites in CPT1A on chromosome 11 was significantly associated with MetS (P for cg00574958 = 2.6x10(-14) and P for cg17058475 = 1.2x10(-9)). Significant associations were replicated in both European and African ancestry participants of the Bogalusa Heart Study. Our findings suggest that methylation in CPT1A is a promising epigenetic marker for MetS risk which could become useful as a treatment target in the future.

The effects of angiotensinogen gene polymorphisms on cardiovascular disease outcomes during antihypertensive treatment in the GenHAT study.

Previous studies have reported that risk of cardiovascular morbidity and mortality substantially increases in hypertensive patients, especially among those with inadequate blood pressure control. Two common antihypertensive drug classes including thiazide diuretics and angiotensin-converting enzyme (ACE) inhibitors affect different enzymes in the renin-angiotensin-aldosterone system (RAAS). In the RAAS, angiotensinogen is converted into angiotensin II which increases blood pressure through vasoconstriction. Using a case-only design with 3448 high-risk hypertensive individuals from the Genetics of Hypertension Associated Treatment (GenHAT) study, we examined whether seven single nucleotide polymorphisms (SNPs) in the angiotensinogen gene (AGT) interact with three classes of antihypertensive drugs including chlorthalidone (a thiazide diuretic), lisinopril (an ACE inhibitor), and amlodipine (a calcium channel blocker) to modify the risk of incident coronary heart disease (CHD) and heart failure (HF) among Caucasian and African American participants, separately. We found no gene by treatment interactions to be statistically significant after correction for multiple testing. However, some suggestive results were found. African American participants with the minor allele of rs11122576 had over two-fold higher risk of CHD when using chlorthalidone compared to using amlodipine, or lisinopril compared to amlodipine (p = 0.006 and p = 0.01, respectively). Other marginal associations are also reported among both race groups. The findings reported here suggest that rs11122576 could contribute to future personalization of antihypertensive treatment among African Americans though more studies are needed.