Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

BACKGROUND

The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP.

METHODS

We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses.

RESULTS

We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

CONCLUSIONS

We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

Refinement of a published gene-physical activity interaction impacting HDL-cholesterol: role of sex and lipoprotein subfractions

Large-scale gene-environment interaction (GxE) discovery efforts often involve compromises in the definition of outcomes and choice of covariates for the sake of data harmonization and statistical power. Consequently, refinement of exposures, covariates, outcomes, and population subsets may be helpful to establish often-elusive replication and evaluate potential clinical utility. Here, we used additional datasets, an expanded set of statistical models, and interrogation of lipoprotein metabolism via nuclear magnetic resonance (NMR)-based lipoprotein subfractions to refine a previously discovered GxE modifying the relationship between physical activity (PA) and HDL-cholesterol (HDL-C). This GxE was originally identified by Kilpeläinen et al., with the strongest cohort-specific signal coming from the Women's Genome Health Study (WGHS). We thus explored this GxE further in the WGHS (N = 23,294), with follow-up in the UK Biobank (UKB; N = 281,380), and the Multi-Ethnic Study of Atherosclerosis (MESA; N = 4,587). Self-reported PA (MET-hrs/wk), genotypes at rs295849 (nearest gene: LHX1), and NMR metabolomics data were available in all three cohorts. As originally reported, minor allele carriers of rs295849 in WGHS had a stronger positive association between PA and HDL-C (pint = 0.002). When testing a range of NMR metabolites (primarily lipoprotein and lipid subfractions) to refine the HDL-C outcome, we found a stronger interaction effect on medium-sized HDL particle concentrations (M-HDL-P; pint = 1.0×10-4) than HDL-C. Meta-regression revealed a systematically larger interaction effect in cohorts from the original meta-analysis with a greater fraction of women (p = 0.018). In the UKB, GxE effects were stronger both in women and using M-HDL-P as the outcome. In MESA, the primary interaction for HDL-C showed nominal significance (pint = 0.013), but without clear differences by sex and with a greater magnitude using large, rather than medium, HDL-P as an outcome. Towards reconciling these observations, further exploration leveraging NMR platform-specific HDL subfraction diameter annotations revealed modest agreement across all cohorts in the interaction affecting medium-to-large particles. Taken together, our work provides additional insights into a specific known gene-PA interaction while illustrating the importance of phenotype and model refinement towards understanding and replicating GxEs.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure

Objective

We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP).

Methods

We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses.

Results

We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1.

Conclusion

We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

A new Approach to Identify Gene-Environment Interactions and Reveal New Biological Insight in Complex traits

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the CHARGE Gene-Lifestyle Interactions Working Group has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identified and confirmed 5 loci (6 independent signals) interacting with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrated that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated contribution ranges from 1.76% to 14.05% of SNP heritability of serum lipids in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

BACKGROUND

Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery.

RESULTS

To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism.

CONCLUSIONS

Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

Rare coding variants in RCN3 are associated with blood pressure

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.

Rare coding variants in 35 genes associate with circulating lipid levels-A multi-ancestry analysis of 170,000 exomes

Large-scale gene sequencing studies for complex traits have the potential to identify causal genes with therapeutic implications. We performed gene-based association testing of blood lipid levels with rare (minor allele frequency < 1%) predicted damaging coding variation by using sequence data from >170,000 individuals from multiple ancestries: 97,493 European, 30,025 South Asian, 16,507 African, 16,440 Hispanic/Latino, 10,420 East Asian, and 1,182 Samoan. We identified 35 genes associated with circulating lipid levels; some of these genes have not been previously associated with lipid levels when using rare coding variation from population-based samples. We prioritize 32 genes in array-based genome-wide association study (GWAS) loci based on aggregations of rare coding variants; three (EVI5, SH2B3, and PLIN1) had no prior association of rare coding variants with lipid levels. Most of our associated genes showed evidence of association among multiple ancestries. Finally, we observed an enrichment of gene-based associations for low-density lipoprotein cholesterol drug target genes and for genes closest to GWAS index single-nucleotide polymorphisms (SNPs). Our results demonstrate that gene-based associations can be beneficial for drug target development and provide evidence that the gene closest to the array-based GWAS index SNP is often the functional gene for blood lipid levels.

Regular exercise and patterns of response across multiple cardiometabolic traits: the HERITAGE family study

OBJECTIVES

We investigated whether high responsiveness or low responsiveness to exercise training aggregates in the same individuals across seven cardiometabolic traits.

METHODS

A total of 564 adults (29.2% black, 53.7% female) from the HERITAGE family study completed a 20-week endurance training programme (at 55%-75% of participants' maximal oxygen uptake (VO₂max)) with VO₂max, per cent body fat, visceral adipose tissue, fasting levels of insulin, high-density lipoprotein cholesterol, small low-density lipoprotein particles and inflammatory marker GlycA measured before and after training. For each exercise response trait, we created ethnicity-specific, sex-specific and generation-specific quintiles. High responses were defined as those within the 20th percentile representing the favourable end of the response trait distribution, low responses were defined as the 20th percentile from the least favourable end, and the remaining were labelled as average responses.

RESULTS

Only one individual had universally high or low responses for all seven cardiometabolic traits. Almost half (49%) of the cohort had at least one high response and one low response across the seven traits. About 24% had at least one high response but no low responses, 24% had one or more low responses but no high responses, and 2.5% had average responses across all traits.

CONCLUSIONS

Interindividual variation in exercise responses was evident in all the traits we investigated, and responsiveness did not aggregate consistently in the same individuals. While adherence to an exercise prescription is known to produce health benefits, targeted risk factors may not improve.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Population sequencing data reveal a compendium of mutational processes in the human germ line

Biological mechanisms underlying human germline mutations remain largely unknown. We statistically decompose variation in the rate and spectra of mutations along the genome using volume-regularized nonnegative matrix factorization. The analysis of a sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. We provide a biological interpretation for seven of these processes. We associate one process with bulky DNA lesions that are resolved asymmetrically with respect to transcription and replication. Two processes track direction of replication fork and replication timing, respectively. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions and a mutagenic effect of long interspersed nuclear elements. We localize a mutagenic process specific to oocytes from population sequencing data. This process appears transcriptionally asymmetric.

Lifestyle Risk Score: handling missingness of individual lifestyle components in meta-analysis of gene-by-lifestyle interactions

Recent studies consider lifestyle risk score (LRS), an aggregation of multiple lifestyle exposures, in identifying association of gene-lifestyle interaction with disease traits. However, not all cohorts have data on all lifestyle factors, leading to increased heterogeneity in the environmental exposure in collaborative meta-analyses. We compared and evaluated four approaches (Naïve, Safe, Complete and Moderator Approaches) to handle the missingness in LRS-stratified meta-analyses under various scenarios. Compared to "benchmark" results with all lifestyle factors available for all cohorts, the Complete Approach, which included only cohorts with all lifestyle components, was underpowered due to lower sample size, and the Naïve Approach, which utilized all available data and ignored the missingness, was slightly inflated. The Safe Approach, which used all data in LRS-exposed group and only included cohorts with all lifestyle factors available in the LRS-unexposed group, and the Moderator Approach, which handled missingness via moderator meta-regression, were both slightly conservative and yielded almost identical p values. We also evaluated the performance of the Safe Approach under different scenarios. We observed that the larger the proportion of cohorts without missingness included, the more accurate the results compared to "benchmark" results. In conclusion, we generally recommend the Safe Approach, a straightforward and non-inflated approach, to handle heterogeneity among cohorts in the LRS based genome-wide interaction meta-analyses.

Chromosome Xq23 is associated with lower atherogenic lipid concentrations and favorable cardiometabolic indices

Autosomal genetic analyses of blood lipids have yielded key insights for coronary heart disease (CHD). However, X chromosome genetic variation is understudied for blood lipids in large sample sizes. We now analyze genetic and blood lipid data in a high-coverage whole X chromosome sequencing study of 65,322 multi-ancestry participants and perform replication among 456,893 European participants. Common alleles on chromosome Xq23 are strongly associated with reduced total cholesterol, LDL cholesterol, and triglycerides (min P = 8.5 × 10-72), with similar effects for males and females. Chromosome Xq23 lipid-lowering alleles are associated with reduced odds for CHD among 42,545 cases and 591,247 controls (P = 1.7 × 10-4), and reduced odds for diabetes mellitus type 2 among 54,095 cases and 573,885 controls (P = 1.4 × 10-5). Although we observe an association with increased BMI, waist-to-hip ratio adjusted for BMI is reduced, bioimpedance analyses indicate increased gluteofemoral fat, and abdominal MRI analyses indicate reduced visceral adiposity. Co-localization analyses strongly correlate increased CHRDL1 gene expression, particularly in adipose tissue, with reduced concentrations of blood lipids.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Deriving stratified effects from joint models investigating gene-environment interactions

Background

Models including an interaction term and performing a joint test of SNP and/or interaction effect are often used to discover Gene-Environment (GxE) interactions. When the environmental exposure is a binary variable, analyses from exposure-stratified models which consist of estimating genetic effect in unexposed and exposed individuals separately can be of interest. In large-scale consortia focusing on GxE interactions in which only the joint test has been performed, it may be challenging to get summary statistics from both exposure-stratified and marginal (i.e not accounting for interaction) models.

Results

In this work, we developed a simple framework to estimate summary statistics in each stratum of a binary exposure and in the marginal model using summary statistics from the “joint” model. We performed simulation studies to assess our estimators’ accuracy and examined potential sources of bias, such as correlation between genotype and exposure and differing phenotypic variances within exposure strata. Results from these simulations highlight the high theoretical accuracy of our estimators and yield insights into the impact of potential sources of bias. We then applied our methods to real data and demonstrate our estimators’ retained accuracy after filtering SNPs by sample size to mitigate potential bias.

Conclusions

These analyses demonstrated the accuracy of our method in estimating both stratified and marginal summary statistics from a joint model of gene-environment interaction. In addition to facilitating the interpretation of GxE screenings, this work could be used to guide further functional analyses. We provide a user-friendly Python script to apply this strategy to real datasets. The Python script and documentation are available at https://gitlab.pasteur.fr/statistical-genetics/j2s.

Smoking-by-genotype interaction in type 2 diabetes risk and fasting glucose

Smoking is a potentially causal behavioral risk factor for type 2 diabetes (T2D), but not all smokers develop T2D. It is unknown whether genetic factors partially explain this variation. We performed genome-environment-wide interaction studies to identify loci exhibiting potential interaction with baseline smoking status (ever vs. never) on incident T2D and fasting glucose (FG). Analyses were performed in participants of European (EA) and African ancestry (AA) separately. Discovery analyses were conducted using genotype data from the 50,000-single-nucleotide polymorphism (SNP) ITMAT-Broad-CARe (IBC) array in 5 cohorts from from the Candidate Gene Association Resource Consortium (n = 23,189). Replication was performed in up to 16 studies from the Cohorts for Heart Aging Research in Genomic Epidemiology Consortium (n = 74,584). In meta-analysis of discovery and replication estimates, 5 SNPs met at least one criterion for potential interaction with smoking on incident T2D at p<1x10-7 (adjusted for multiple hypothesis-testing with the IBC array). Two SNPs had significant joint effects in the overall model and significant main effects only in one smoking stratum: rs140637 (FBN1) in AA individuals had a significant main effect only among smokers, and rs1444261 (closest gene C2orf63) in EA individuals had a significant main effect only among nonsmokers. Three additional SNPs were identified as having potential interaction by exhibiting a significant main effects only in smokers: rs1801232 (CUBN) in AA individuals, rs12243326 (TCF7L2) in EA individuals, and rs4132670 (TCF7L2) in EA individuals. No SNP met significance for potential interaction with smoking on baseline FG. The identification of these loci provides evidence for genetic interactions with smoking exposure that may explain some of the heterogeneity in the association between smoking and T2D.

Identifying blood pressure loci whose effects are modulated by multiple lifestyle exposures

Although multiple lifestyle exposures simultaneously impact blood pressure (BP) and cardiovascular health, most analysis so far has considered each single lifestyle exposure (e.g., smoking) at a time. Here, we exploit gene-multiple lifestyle exposure interactions to find novel BP loci. For each of 6,254 Framingham Heart Study participants, we computed lifestyle risk score (LRS) value by aggregating the risk of four lifestyle exposures (smoking, alcohol, education, and physical activity) on BP. Using the LRS, we performed genome-wide gene-environment interaction analysis in systolic and diastolic BP using the joint 2 degree of freedom (DF) and 1 DF interaction tests. We identified one genome-wide significant (p < 5 × 10^-8 ) and 11 suggestive (p < 1 × 10^-6 ) loci. Gene-environment analysis using single lifestyle exposures identified only one of the 12 loci. Nine of the 12 BP loci detected were novel. Loci detected by the LRS were located within or nearby genes with biologically plausible roles in the pathophysiology of hypertension, including KALRN, VIPR2, SNX1, and DAPK2. Our results suggest that simultaneous consideration of multiple lifestyle exposures in gene-environment interaction analysis can identify additional loci missed by single lifestyle approaches.

Genome-wide meta-analysis of variant-by-diuretic interactions as modulators of lipid traits in persons of European and African ancestry

Hypertension (HTN) is a significant risk factor for cardiovascular morbidity and mortality. Metabolic abnormalities, including adverse cholesterol and triglycerides (TG) profiles, are frequent comorbid findings with HTN and contribute to cardiovascular disease. Diuretics, which are used to treat HTN and heart failure, have been associated with worsening of fasting lipid concentrations. Genome-wide meta-analyses with 39,710 European-ancestry (EA) individuals and 9,925 African-ancestry (AA) individuals were performed to identify genetic variants that modify the effect of loop or thiazide diuretic use on blood lipid concentrations. Both longitudinal and cross-sectional data were used to compute cohort-specific interaction results, which were then combined through meta-analysis in each ancestry. These ancestry-specific results were further combined through trans-ancestry meta-analysis. Analysis of EA data identified two genome-wide significant (p < 5×10⁻⁸) loci with single nucleotide variant (SNV)-loop diuretic interaction on TG concentrations (including COL11A1). Analysis of AA data identified one genome-wide significant locus adjacent to BMP2 with SNV-loop diuretic interaction on TG concentrations. Trans-ancestry analysis strengthened evidence of association for SNV-loop diuretic interaction at two loci (KIAA1217 and BAALC). There were few significant SNV-thiazide diuretic interaction associations on TG concentrations and for either diuretic on cholesterol concentrations. Several promising loci were identified that may implicate biologic pathways that contribute to adverse metabolic side effects from diuretic therapy.

VarExp: estimating variance explained by genome-wide GxE summary statistics

Summary

Many genome-wide association studies and genome-wide screening for gene-environment (GxE) interactions have been performed to elucidate the underlying mechanisms of human traits and diseases. When the analyzed outcome is quantitative, the overall contribution of identified genetic variants to the outcome is often expressed as the percentage of phenotypic variance explained. This is commonly done using individual-level genotype data but it is challenging when results are derived through meta-analyses. Here, we present R package, 'VarExp', that allows for the estimation of the percentage of phenotypic variance explained using summary statistics only. It allows for a range of models to be evaluated, including marginal genetic effects, GxE interaction effects and both effects jointly. Its implementation integrates all recent methodological developments and does not need external data to be uploaded by users.

Availability and implementation

The R package is available at https://gitlab.pasteur.fr/statistical-genetics/VarExp.git.

Supplementary information

Supplementary data are available at Bioinformatics online.

The Promise of Selecting Individuals from the Extremes of Exposure in the Analysis of Gene-Physical Activity Interactions

BACKGROUND

Dichotomization using the lower quartile as cutoff is commonly used for harmonizing heterogeneous physical activity (PA) measures across studies. However, this may create misclassification and hinder discovery of new loci.

OBJECTIVES

This study aimed to evaluate the performance of selecting individuals from the extremes of the exposure (SIEE) as an alternative approach to reduce such misclassification.

METHOD

For systolic and diastolic blood pressure in the Framingham Heart Study, we performed a genome-wide association study with gene-PA interaction analysis using three PA variables derived by SIEE and two other dichotomization approaches. We compared number of loci detected and overlap with loci found using a quantitative PA variable. In addition, we performed simulation studies to assess bias, false discovery rates (FDR), and power under synergistic/antagonistic genetic effects in exposure groups and in the presence/absence of measurement error.

RESULTS

In the empirical analysis, SIEE's performance was neither the best nor the worst. In most simulation scenarios, SIEE was consistently outperformed in terms of FDR and power. Particularly, in a scenario characterized by antagonistic effects and measurement error, SIEE had the least bias and highest power.

CONCLUSION

SIEE's promise appears limited to detecting loci with antagonistic effects. Further studies are needed to evaluate SIEE's full advantage.

Leveraging linkage evidence to identify low-frequency and rare variants on 16p13 associated with blood pressure using TOPMed whole genome sequencing data

In this study, we investigated low-frequency and rare variants associated with blood pressure (BP) by focusing on a linkage region on chromosome 16p13. We used whole genome sequencing (WGS) data obtained through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program on 395 Cleveland Family Study (CFS) European Americans (CFS-EA). By analyzing functional coding variants and non-coding rare variants with CADD score > 10 residing within the chromosomal region in families with linkage evidence, we observed 25 genes with nominal statistical evidence (burden or SKAT p < 0.05). One of the genes is RBFOX1, an evolutionarily conserved RNA-binding protein that regulates tissue-specific alternative splicing that we previously reported to be associated with BP using exome array data in CFS. After follow-up analysis of the 25 genes in ten independent TOPMed studies with individuals of European, African, and East Asian ancestry, and Hispanics (N = 29,988), we identified variants in SLX4 (p = 2.19 × 10^-4) to be significantly associated with BP traits when accounting for multiple testing. We also replicated the associations previously reported for RBFOX1 (p = 0.007). Follow-up analysis with GTEx eQTL data shows SLX4 variants are associated with gene expression in coronary artery, multiple brain tissues, and right atrial appendage of the heart. Our study demonstrates that linkage analysis of family data can provide an efficient approach for detecting rare variants associated with complex traits in WGS data.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Effects of regular endurance exercise on GlycA: Combined analysis of 14 exercise interventions

BACKGROUND AND AIMS

GlycA is a relatively new biomarker for inflammation as well as cardiometabolic disease risk. However, the effect of exercise on GlycA is largely unknown. Therefore, the purpose of this study was to examine the effects of regular exercise on the inflammatory marker GlycA across seven studies and 14 exercise interventions.

METHODS

Nuclear magnetic resonance spectroscopy, specifically signal amplitudes originating from the N-acetyl methyl group protons of the N-acetylglucosamine residues on the glycan branches of glycoproteins, was used to quantify GlycA concentrations. GlycA was measured before and after completion of an exercise intervention in 1568 individuals across seven studies and 14 exercise interventions. Random effects inverse variance weighting models were used to pool effects across interventions.

RESULTS

Combined analysis of unadjusted data showed that regular exercise significantly (p = 2 × 10-6) reduced plasma GlycA (-8.26 ± 1.8 μmol/L). This reduction remained significant (-9.12 ± 1.9 μmol/L, p = 1.22 × 10-6) following adjustment for age, sex, race, baseline BMI, and baseline GlycA. Changes in GlycA were correlated with changes in traditional inflammatory markers, C-reactive protein, interleukin-6, and fibrinogen, however, these correlations were relatively weak (range r: 0.21-0.38, p < 0.0001).

CONCLUSIONS

Regular exercise significantly reduced plasma GlycA across 14 different exercise interventions despite differences in exercise programs and study populations. The current study provides a greater understanding of the use of exercise as a potential therapy for the reduction of systemic inflammation. Further research is needed to understand the mechanisms behind the exercise-related reductions in GlycA.

The genetic underpinnings of variation in ages at menarche and natural menopause among women from the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) Study: A trans-ethnic meta-analysis

Current knowledge of the genetic architecture of key reproductive events across the female life course is largely based on association studies of European descent women. The relevance of known loci for age at menarche (AAM) and age at natural menopause (ANM) in diverse populations remains unclear. We investigated 32 AAM and 14 ANM previously-identified loci and sought to identify novel loci in a trans-ethnic array-wide study of 196,483 SNPs on the MetaboChip (Illumina, Inc.). A total of 45,364 women of diverse ancestries (African, Hispanic/Latina, Asian American and American Indian/Alaskan Native) in the Population Architecture using Genomics and Epidemiology (PAGE) Study were included in cross-sectional analyses of AAM and ANM. Within each study we conducted a linear regression of SNP associations with self-reported or medical record-derived AAM or ANM (in years), adjusting for birth year, population stratification, and center/region, as appropriate, and meta-analyzed results across studies using multiple meta-analytic techniques. For both AAM and ANM, we observed more directionally consistent associations with the previously reported risk alleles than expected by chance (p-values_binomial≤0.01). Eight densely genotyped reproductive loci generalized significantly to at least one non-European population. We identified one trans-ethnic array-wide SNP association with AAM and two significant associations with ANM, which have not been described previously. Additionally, we observed evidence of independent secondary signals at three of six AAM trans-ethnic loci. Our findings support the transferability of reproductive trait loci discovered in European women to women of other race/ethnicities and indicate the presence of additional trans-ethnic associations both at both novel and established loci. These findings suggest the benefit of including diverse populations in future studies of the genetic architecture of female growth and development.

Family History of Coronary Heart Disease and Hemostatic Variables in Middle-Aged Adults

Individuals with a family history of coronary heart disease (CHD) may be predisposed to atherothrombosis. To investigate this hypothesis, a family CHD risk score was computed for approximately 13,000 men and women aged 45 to 64; hemostatic variables (fibrinogen, factor VIIc, factor VIIIc, von Willebrand factor, antithrombin III, protein C) were also measured in plasma. After adjustment for age and ethnicity, there was a statistically significant, positive association between the family risk score and four of the six hemostatic variables (fibrinogen, factor VIIc, factor VIIIc, von Willebrand factor) in women and all six hemostatic variables in men. In general, these associations were weak and substantially attenuated after adjustment for constitutional, lifestyle, and biochemical covariates. These results indicate that mean levels of selected hemostatic variables, like traditional CHD risk factors, are higher in individuals with a family history of heart disease.

Multiancestry Study of Gene-Lifestyle Interactions for Cardiovascular Traits in 610 475 Individuals From 124 Cohorts: Design and Rationale

BACKGROUND

Several consortia have pursued genome-wide association studies for identifying novel genetic loci for blood pressure, lipids, hypertension, etc. They demonstrated the power of collaborative research through meta-analysis of study-specific results.

METHODS AND RESULTS

The Gene-Lifestyle Interactions Working Group was formed to facilitate the first large, concerted, multiancestry study to systematically evaluate gene-lifestyle interactions. In stage 1, genome-wide interaction analysis is performed in 53 cohorts with a total of 149 684 individuals from multiple ancestries. In stage 2 involving an additional 71 cohorts with 460 791 individuals from multiple ancestries, focused analysis is performed for a subset of the most promising variants from stage 1. In all, the study involves up to 610 475 individuals. Current focus is on cardiovascular traits including blood pressure and lipids, and lifestyle factors including smoking, alcohol, education (as a surrogate for socioeconomic status), physical activity, psychosocial variables, and sleep. The total sample sizes vary among projects because of missing data. Large-scale gene-lifestyle or more generally gene-environment interaction (G×E) meta-analysis studies can be cumbersome and challenging. This article describes the design and some of the approaches pursued in the interaction projects.

CONCLUSIONS

The Gene-Lifestyle Interactions Working Group provides an excellent framework for understanding the lifestyle context of genetic effects and to identify novel trait loci through analysis of interactions. An important and novel feature of our study is that the gene-lifestyle interaction (G×E) results may improve our knowledge about the underlying mechanisms for novel and already known trait loci.

Adropin: An endocrine link between the biological clock and cholesterol homeostasis

Objective

Identify determinants of plasma adropin concentrations, a secreted peptide translated from the Energy Homeostasis Associated (ENHO) gene linked to metabolic control and vascular function.

Methods

Associations between plasma adropin concentrations, demographics (sex, age, BMI) and circulating biomarkers of lipid and glucose metabolism were assessed in plasma obtained after an overnight fast in humans. The regulation of adropin expression was then assessed in silico, in cultured human cells, and in animal models.

Results

In humans, plasma adropin concentrations are inversely related to atherogenic LDL-cholesterol (LDL-C) levels in men (n = 349), but not in women (n = 401). Analysis of hepatic Enho expression in male mice suggests control by the biological clock. Expression is rhythmic, peaking during maximal food consumption in the dark correlating with transcriptional activation by RORα/γ. The nadir in the light phase coincides with the rest phase and repression by Rev-erb. Plasma adropin concentrations in nonhuman primates (rhesus monkeys) also exhibit peaks coinciding with feeding times (07:00 h, 15:00 h). The ROR inverse agonists SR1001 and the 7-oxygenated sterols 7-β-hydroxysterol and 7-ketocholesterol, or the Rev-erb agonist SR9009, suppress ENHO expression in cultured human HepG2 cells. Consumption of high-cholesterol diets suppress expression of the adropin transcript in mouse liver. However, adropin over expression does not prevent hypercholesterolemia resulting from a high cholesterol diet and/or LDL receptor mutations.

Conclusions

In humans, associations between plasma adropin concentrations and LDL-C suggest a link with hepatic lipid metabolism. Mouse studies suggest that the relationship between adropin and cholesterol metabolism is unidirectional, and predominantly involves suppression of adropin expression by cholesterol and 7-oxygenated sterols. Sensing of fatty acids, cholesterol and oxysterols by the RORα/γ ligand-binding domain suggests a plausible functional link between adropin expression and cellular lipid metabolism. Furthermore, the nuclear receptors RORα/γ and Rev-erb may couple adropin synthesis with circadian rhythms in carbohydrate and lipid metabolism.

•

In male humans, plasma adropin concentrations are inversely related to low-density circulating cholesterol (LDL-C) levels.

•

Adropin expression is regulated by core elements of the biological clock (RORA/G, Rev-Erb).

•

Sterol-sensing by the ROR ligand-binding domain provides a plausible link between adropin expression and lipid metabolism.

•

In mouse liver, adropin expression is rhythmic and suppressed by exogenous (dietary) cholesterol.

A Review of the Genetics of Hypertension with a Focus on Gene-Environment Interactions

PURPOSE OF REVIEW

Here, we discuss the interpretation and modeling of gene-environment interactions in hypertension-related phenotypes, with a focus on the necessary assumptions and possible challenges.

RECENT FINDINGS

Recently, small cohort studies have discovered several novel genetic variants associated with hypertension-related phenotypes through modeling gene-environment interactions. Several consortia-based meta-analytic efforts have uncovered many novel genetic variants in hypertension without modeling interaction terms, giving promise to future meta-analytic efforts that incorporate gene-environment interactions. Heritability studies and genome-wide association studies have established that hypertension, a prevalent cardiovascular disease, has a genetic component that may be modulated by the environment (such as lifestyle factors). This review includes a discussion of known genetic associations for hypertension/blood pressure, including those resulting from the incorporation of gene-environmental interaction modeling.

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.

Single-trait and multi-trait genome-wide association analyses identify novel loci for blood pressure in African-ancestry populations

Hypertension is a leading cause of global disease, mortality, and disability. While individuals of African descent suffer a disproportionate burden of hypertension and its complications, they have been underrepresented in genetic studies. To identify novel susceptibility loci for blood pressure and hypertension in people of African ancestry, we performed both single and multiple-trait genome-wide association analyses. We analyzed 21 genome-wide association studies comprised of 31,968 individuals of African ancestry, and validated our results with additional 54,395 individuals from multi-ethnic studies. These analyses identified nine loci with eleven independent variants which reached genome-wide significance (P < 1.25×10-8) for either systolic and diastolic blood pressure, hypertension, or for combined traits. Single-trait analyses identified two loci (TARID/TCF21 and LLPH/TMBIM4) and multiple-trait analyses identified one novel locus (FRMD3) for blood pressure. At these three loci, as well as at GRP20/CDH17, associated variants had alleles common only in African-ancestry populations. Functional annotation showed enrichment for genes expressed in immune and kidney cells, as well as in heart and vascular cells/tissues. Experiments driven by these findings and using angiotensin-II induced hypertension in mice showed altered kidney mRNA expression of six genes, suggesting their potential role in hypertension. Our study provides new evidence for genes related to hypertension susceptibility, and the need to study African-ancestry populations in order to identify biologic factors contributing to hypertension.

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Discovery and fine-mapping of adiposity loci using high density imputation of genome-wide association studies in individuals of African ancestry: African Ancestry Anthropometry Genetics Consortium

Genome-wide association studies (GWAS) have identified >300 loci associated with measures of adiposity including body mass index (BMI) and waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), but few have been identified through screening of the African ancestry genomes. We performed large scale meta-analyses and replications in up to 52,895 individuals for BMI and up to 23,095 individuals for WHR_adjBMI from the African Ancestry Anthropometry Genetics Consortium (AAAGC) using 1000 Genomes phase 1 imputed GWAS to improve coverage of both common and low frequency variants in the low linkage disequilibrium African ancestry genomes. In the sex-combined analyses, we identified one novel locus (TCF7L2/HABP2) for WHR_adjBMI and eight previously established loci at P < 5×10⁻⁸: seven for BMI, and one for WHR_adjBMI in African ancestry individuals. An additional novel locus (SPRYD7/DLEU2) was identified for WHR_adjBMI when combined with European GWAS. In the sex-stratified analyses, we identified three novel loci for BMI (INTS10/LPL and MLC1 in men, IRX4/IRX2 in women) and four for WHR_adjBMI (SSX2IP, CASC8, PDE3B and ZDHHC1/HSD11B2 in women) in individuals of African ancestry or both African and European ancestry. For four of the novel variants, the minor allele frequency was low (<5%). In the trans-ethnic fine mapping of 47 BMI loci and 27 WHR_adjBMI loci that were locus-wide significant (P < 0.05 adjusted for effective number of variants per locus) from the African ancestry sex-combined and sex-stratified analyses, 26 BMI loci and 17 WHR_adjBMI loci contained ≤ 20 variants in the credible sets that jointly account for 99% posterior probability of driving the associations. The lead variants in 13 of these loci had a high probability of being causal. As compared to our previous HapMap imputed GWAS for BMI and WHR_adjBMI including up to 71,412 and 27,350 African ancestry individuals, respectively, our results suggest that 1000 Genomes imputation showed modest improvement in identifying GWAS loci including low frequency variants. Trans-ethnic meta-analyses further improved fine mapping of putative causal variants in loci shared between the African and European ancestry populations.

Genome-wide association studies (GWAS) have identified >300 genetic regions that influence body size and shape as measured by body mass index (BMI) and waist-to-hip ratio (WHR), respectively, but few have been identified in populations of African ancestry. We conducted large scale high coverage GWAS and replication of these traits in 52,895 and 23,095 individuals of African ancestry, respectively, followed by additional replication in European populations. We identified 10 genome-wide significant loci in all individuals, and an additional seven loci by analyzing men and women separately. We combined African and European ancestry GWAS and were able to narrow down 43 out of 74 African ancestry associated genetic regions to contain small number of putative causal variants. Our results highlight the improvement of applying high density genome coverage and combining multiple ancestries in the identification and refinement of location of genetic regions associated with adiposity traits.

Phenomapping for the Identification of Hypertensive Patients with the Myocardial Substrate for Heart Failure with Preserved Ejection Fraction

We sought to evaluate whether unbiased machine learning of dense phenotypic data ("phenomapping") could identify distinct hypertension subgroups that are associated with the myocardial substrate (i.e., abnormal cardiac mechanics) for heart failure with preserved ejection fraction (HFpEF). In the HyperGEN study, a population- and family-based study of hypertension, we studied 1273 hypertensive patients utilizing clinical, laboratory, and conventional echocardiographic phenotyping of the study participants. We used machine learning analysis of 47 continuous phenotypic variables to identify mutually exclusive groups constituting a novel classification of hypertension. The phenomapping analysis classified study participants into 2 distinct groups that differed markedly in clinical characteristics, cardiac structure/function, and indices of cardiac mechanics (e.g., phenogroup #2 had a decreased absolute longitudinal strain [12.8 ± 4.1 vs. 14.6 ± 3.5%] even after adjustment for traditional comorbidities [p < 0.001]). The 2 hypertension phenogroups may represent distinct subtypes that may benefit from targeted therapies for the prevention of HFpEF.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Variant Discovery and Fine Mapping of Genetic Loci Associated with Blood Pressure Traits in Hispanics and African Americans

Despite the substantial burden of hypertension in US minority populations, few genetic studies of blood pressure have been conducted in Hispanics and African Americans, and it is unclear whether many of the established loci identified in European-descent populations contribute to blood pressure variation in non-European descent populations. Using the Metabochip array, we sought to characterize the genetic architecture of previously identified blood pressure loci, and identify novel cardiometabolic variants related to systolic and diastolic blood pressure in a multi-ethnic US population including Hispanics (n = 19,706) and African Americans (n = 18,744). Several known blood pressure loci replicated in African Americans and Hispanics. Fourteen variants in three loci (KCNK3, FGF5, ATXN2-SH2B3) were significantly associated with blood pressure in Hispanics. The most significant diastolic blood pressure variant identified in our analysis, rs2586886/KCNK3 (P = 5.2 x 10⁻⁹), also replicated in independent Hispanic and European-descent samples. African American and trans-ethnic meta-analysis data identified novel variants in the FGF5, ULK4 and HOXA-EVX1 loci, which have not been previously associated with blood pressure traits. Our identification and independent replication of variants in KCNK3, a gene implicated in primary hyperaldosteronism, as well as a variant in HOTTIP (HOXA-EVX1) suggest that further work to clarify the roles of these genes may be warranted. Overall, our findings suggest that loci identified in European descent populations also contribute to blood pressure variation in diverse populations including Hispanics and African Americans—populations that are understudied for hypertension genetic risk factors.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

Familial resemblance of adiposity-related parameters: results from a health check-up population in Taiwan

The mechanisms of obesity is still unclear, however, genetic and environmental factors are two major possible causes of obesity. The purpose of this study was to assess the degree of familial resemblance of different obesity-related parameters in a health check-up population in Taiwan. We measured body mass index (BMI), waist-to-hip ratio (WHR) and percentage of body fat (BFAT) anthropometrics in 1724 members of 431 families participating in the MJ Health Screening program. Each family contributed four members, i.e. father, mother, son, and daughter. All the participants were free from coronary heart disease, hypertension, diabetes, dyslipidemia, and generally in good health. The degree of familial aggregation was measured by heritability that was calculated based on age-adjusted familial (parent-offspring, sibling, spouse) correlations. The maximal heritability estimates were 39, 30 and 35% for BMI, WHR and BFAT, respectively. For WHR, the correlation between spouses was not significant and the heritability appears to be predominantly due to genetic causes. Furthermore, for BMI and BFAT, the spouse correlations were 0.08 and 0.11, respectively. The heritabilities for BMI and BFAT were mostly explained by genetic factors and familial environmental factors such as dietary habits or physical activity. The familial resemblance of various obesity-related parameters was moderate in a health check-up population in Taiwan. While the heritability for WHR appears to be mainly due to genetic factors, the familial resemblance for BMI and percentage of BFAT may involve both genetic and familial environmental factors.

New loci for body fat percentage reveal link between adiposity and cardiometabolic disease risk

To increase our understanding of the genetic basis of adiposity and its links to cardiometabolic disease risk, we conducted a genome-wide association meta-analysis of body fat percentage (BF%) in up to 100,716 individuals. Twelve loci reached genome-wide significance (P<5 × 10(-8)), of which eight were previously associated with increased overall adiposity (BMI, BF%) and four (in or near COBLL1/GRB14, IGF2BP1, PLA2G6, CRTC1) were novel associations with BF%. Seven loci showed a larger effect on BF% than on BMI, suggestive of a primary association with adiposity, while five loci showed larger effects on BMI than on BF%, suggesting association with both fat and lean mass. In particular, the loci more strongly associated with BF% showed distinct cross-phenotype association signatures with a range of cardiometabolic traits revealing new insights in the link between adiposity and disease risk.

Genome-wide meta-analysis uncovers novel loci influencing circulating leptin levels

Leptin is an adipocyte-secreted hormone, the circulating levels of which correlate closely with overall adiposity. Although rare mutations in the leptin (LEP) gene are well known to cause leptin deficiency and severe obesity, no common loci regulating circulating leptin levels have been uncovered. Therefore, we performed a genome-wide association study (GWAS) of circulating leptin levels from 32,161 individuals and followed up loci reaching P<10(-6) in 19,979 additional individuals. We identify five loci robustly associated (P<5 × 10(-8)) with leptin levels in/near LEP, SLC32A1, GCKR, CCNL1 and FTO. Although the association of the FTO obesity locus with leptin levels is abolished by adjustment for BMI, associations of the four other loci are independent of adiposity. The GCKR locus was found associated with multiple metabolic traits in previous GWAS and the CCNL1 locus with birth weight. Knockdown experiments in mouse adipose tissue explants show convincing evidence for adipogenin, a regulator of adipocyte differentiation, as the novel causal gene in the SLC32A1 locus influencing leptin levels. Our findings provide novel insights into the regulation of leptin production by adipose tissue and open new avenues for examining the influence of variation in leptin levels on adiposity and metabolic health.

A Genome-wide study of blood pressure in African Americans accounting for gene-smoking interaction

Cigarette smoking has been shown to be a health hazard. In addition to being considered a negative lifestyle behavior, studies have shown that cigarette smoking has been linked to genetic underpinnings of hypertension. Because African Americans have the highest incidence and prevalence of hypertension, we examined the joint effect of genetics and cigarette smoking on health among this understudied population. The sample included African Americans from the genome wide association studies of HyperGEN (N = 1083, discovery sample) and GENOA (N = 1427, replication sample), both part of the FBPP. Results suggested that 2 SNPs located on chromosomes 14 (NEDD8; rs11158609; raw p = 9.80 × 10⁻⁹, genomic control-adjusted p = 2.09 × 10⁻⁷) and 17 (TTYH2; rs8078051; raw p = 6.28 × 10⁻⁸, genomic control-adjusted p = 9.65 × 10⁻⁷) were associated with SBP including the genetic interaction with cigarette smoking. These two SNPs were not associated with SBP in a main genetic effect only model. This study advances knowledge in the area of main and joint effects of genetics and cigarette smoking on hypertension among African Americans and offers a model to the reader for assessing these risks. More research is required to determine how these genes play a role in expression of hypertension.

The effects of exercise on the lipoprotein subclass profile: A meta-analysis of 10 interventions

OBJECTIVE

The goal was to examine lipoprotein subclass responses to regular exercise as measured in 10 exercise interventions derived from six cohorts.

METHODS

Nuclear magnetic resonance spectroscopy was used to quantify average particle size, total and subclass concentrations of very low-density lipoprotein, low-density lipoprotein, and high-density lipoprotein particles (VLDL-P, LDL-P, and HDL-P, respectively) before and after an exercise intervention in 1555 adults from six studies, encompassing 10 distinct exercise programs: APOE (N = 106), DREW (N = 385), GERS (N = 79), HERITAGE (N = 715), STRRIDE I (N = 168) and II (N = 102). Random-effects meta-analyses were performed to evaluate the overall estimate of mean change across the unadjusted and adjusted mean change values from each exercise group.

RESULTS

Meta-analysis of unadjusted data showed that regular exercise induced significant decreases in the concentration of large VLDL-P, small LDL-P, and medium HDL-P and mean VLDL-P size, with significant increases in the concentration of large LDL-P and large HDL-P and mean LDL-P size. These changes remained significant in meta-analysis with adjustment for age, sex, race, baseline body mass index, and baseline trait value.

CONCLUSIONS

Despite differences in exercise programs and study populations, regular exercise produced putatively beneficial changes in the lipoprotein subclass profile across 10 exercise interventions. Further research is needed to examine how exercise-induced changes in lipoprotein subclasses may be associated with (concomitant changes in) cardiovascular disease risk.

Genomic and transcriptomic predictors of triglyceride response to regular exercise

AIM

We performed genome-wide and transcriptome-wide profiling to identify genes and single nucleotide polymorphisms (SNPs) associated with the response of triglycerides (TG) to exercise training.

METHODS

Plasma TG levels were measured before and after a 20-week endurance training programme in 478 white participants from the HERITAGE Family Study. Illumina HumanCNV370-Quad v3.0 BeadChips were genotyped using the Illumina BeadStation 500GX platform. Affymetrix HG-U133+2 arrays were used to quantitate gene expression levels from baseline muscle biopsies of a subset of participants (N=52). Genome-wide association study (GWAS) analysis was performed using MERLIN, while transcriptomic predictor models were developed using the R-package GALGO.

RESULTS

The GWAS results showed that eight SNPs were associated with TG training-response (ΔTG) at p<9.9×10(-6), while another 31 SNPs showed p values <1×10(-4). In multivariate regression models, the top 10 SNPs explained 32.0% of the variance in ΔTG, while conditional heritability analysis showed that four SNPs statistically accounted for all of the heritability of ΔTG. A molecular signature based on the baseline expression of 11 genes predicted 27% of ΔTG in HERITAGE, which was validated in an independent study. A composite SNP score based on the top four SNPs, each from the genomic and transcriptomic analyses, was the strongest predictor of ΔTG (R(2)=0.14, p=3.0×10(-68)).

CONCLUSIONS

Our results indicate that skeletal muscle transcript abundance at 11 genes and SNPs at a number of loci contribute to TG response to exercise training. Combining data from genomics and transcriptomics analyses identified a SNP-based gene signature that should be further tested in independent samples.

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Aggregate blood pressure responses to serial dietary sodium and potassium intervention: defining responses using independent component analysis

BACKGROUND

Hypertension is a complex trait that often co-occurs with other conditions such as obesity and is affected by genetic and environmental factors. Aggregate indices such as principal components among these variables and their responses to environmental interventions may represent novel information that is potentially useful for genetic studies.

RESULTS

In this study of families participating in the Genetic Epidemiology Network of Salt Sensitivity (GenSalt) Study, blood pressure (BP) responses to dietary sodium interventions are explored. Independent component analysis (ICA) was applied to 20 variables indexing obesity and BP measured at baseline and during low sodium, high sodium and high sodium plus potassium dietary intervention periods. A "heat map" protocol that classifies subjects based on risk for hypertension is used to interpret the extracted components. ICA and heat map suggest four components best describe the data: (1) systolic hypertension, (2) general hypertension, (3) response to sodium intervention and (4) obesity. The largest heritabilities are for the systolic (64%) and general hypertension (56%) components. There is a pattern of higher heritability for the component response to intervention (40-42%) as compared to those for the traditional intervention responses computed as delta scores (24%-40%).

CONCLUSIONS

In summary, the present study provides intermediate phenotypes that are heritable. Using these derived components may prove useful in gene discovery applications.

Influence of Smoking Status and Intensity on Discovery of Blood Pressure Loci Through Gene-Smoking Interactions

BACKGROUND

Genetic variation accounts for approximately 30% of blood pressure (BP) variability but most of that variability has not been attributed to specific variants. Interactions between genes and BP-associated factors may explain some "missing heritability." Cigarette smoking increases BP after short-term exposure and decreases BP with longer exposure. Gene-smoking interactions have discovered novel BP loci, but the contribution of smoking status and intensity to gene discovery is unknown.

METHODS

We analyzed gene-smoking intensity interactions for association with systolic BP (SBP) in three subgroups from the Framingham Heart Study: current smokers only (N = 1,057), current and former smokers ("ever smokers," N = 3,374), and all subjects (N = 6,710). We used three smoking intensity variables defined at cutoffs of 10, 15, and 20 cigarettes per day (CPD). We evaluated the 1 degree-of-freedom (df) interaction and 2df joint test using generalized estimating equations.

RESULTS

Analysis of current smokers using a CPD cutoff of 10 produced two loci associated with SBP. The rs9399633 minor allele was associated with increased SBP (5 mmHg) in heavy smokers (CPD > 10) but decreased SBP (7 mmHg) in light smokers (CPD ≤ 10). The rs11717948 minor allele was associated with decreased SBP (8 mmHg) in light smokers but decreased SBP (2 mmHg) in heavy smokers. Across all nine analyses, 19 additional loci reached P < 1 × 10(-6).

DISCUSSION

Analysis of current smokers may have the highest power to detect gene-smoking interactions, despite the reduced sample size. Associations of loci near SASH1 and KLHL6/KLHL24 with SBP may be modulated by tobacco smoking.