A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

Multi-ancestry transcriptome-wide association analyses yield insights into tobacco use biology and drug repurposing

Most transcriptome-wide association studies (TWASs) so far focus on European ancestry and lack diversity. To overcome this limitation, we aggregated genome-wide association study (GWAS) summary statistics, whole-genome sequences and expression quantitative trait locus (eQTL) data from diverse ancestries. We developed a new approach, TESLA (multi-ancestry integrative study using an optimal linear combination of association statistics), to integrate an eQTL dataset with a multi-ancestry GWAS. By exploiting shared phenotypic effects between ancestries and accommodating potential effect heterogeneities, TESLA improves power over other TWAS methods. When applied to tobacco use phenotypes, TESLA identified 273 new genes, up to 55% more compared with alternative TWAS methods. These hits and subsequent fine mapping using TESLA point to target genes with biological relevance. In silico drug-repurposing analyses highlight several drugs with known efficacy, including dextromethorphan and galantamine, and new drugs such as muscle relaxants that may be repurposed for treating nicotine addiction.

Improving Specialty Care Access via Telemedicine

Introduction: Telehealth is a potential solution to persistent disparities in health and health care access by eliminating structural barriers to care. However, its adoption in urban underserved settings has been limited and remains poorly characterized. Methods: This is a prospective cohort study of patients receiving telemedicine (TM) consultation for specialty care of diabetes, hypertension, and/or kidney disease with a Federally Qualified Health Center (FQHC) as the originating site and an academic medical center (AMC) multispecialty group practice as the distant site in an urban setting. Primary data were collected onsite at a local FQHC and an urban AMC between March 2017 and March 2020, before the COVID-19 pandemic. Clinical outcomes of study participants were compared with matched controls (CON) from a sister FQHC site who were referred for traditional in-person specialty visits at the AMC. No-show rates for study participants were calculated and compared to their no-show rates for standard (STD) in-person specialty visits at the AMC during the study period. A patient satisfaction questionnaire was administered at the end of each TM visit. Results: Visit attendance data were analyzed for 104 patients (834 visits). The no-show rate was 15%. The adjusted odds ratio for no-show for TM versus STD visits was 1.03 [0.66-1.63], p = 0.87. There were no significant differences between TM and CON groups in the change from pre- to intervention periods for mean arterial pressure (p = 0.26), serum creatinine (p = 0.90), or estimated glomerular filtration rate (p = 0.56). The reduction in hemoglobin A1c was significant at a trend level (p = 0.053). Patients indicated high overall satisfaction with TM. Discussion: The study demonstrated improved glycemic control and equivalent outcomes in TM management of hypertension and kidney disease with excellent patient satisfaction. This supports ongoing efforts to increase the availability of TM to improve access to care for urban underserved populations.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Serial 7-Day Electrocardiogram Patch Screening for AF in High-Risk Older Women by the CHARGE-AF Score

BACKGROUND

Asymptomatic atrial fibrillation (AF) is associated with an increased risk of stroke. The yield of serial electrocardiographic (ECG) screening for AF is unknown.

OBJECTIVES

The aim of this study was to determine the frequency of AF detected by serial, 7-day ECG patch screenings in older women identified as having an elevated risk of AF according to the CHARGE (Cohorts for Heart and Aging Research in Genomic Epidemiology)-AF clinical prediction score.

METHODS

Postmenopausal women with a 5-year predicted risk of new-onset AF ≥5% according to CHARGE-AF were recruited from the ongoing WHISH (Women's Health Initiative Strong and Healthy) randomized trial of a physical activity intervention. Participants with AF at baseline by self-report or medical records review were excluded. Screening with 7-day ECG patch monitors was performed at baseline, 6 months, and 12 months from study enrollment.

RESULTS

On baseline monitoring, 2.5% of the cohort had AF detected, increasing to 3.7% by 6 months and 4.9% cumulatively by 12 months. Yield of patch screening was higher among participants with a higher (≥10%) CHARGE-AF score: 4.2% had AF detected at baseline, 5.9% at 6 months, and 7.2% at 12 months. Most participants with patch-identified AF never had a clinical diagnosis of AF (36 of 46 [78%]).

CONCLUSIONS

Older women with an elevated CHARGE-AF score had a high prevalence of AF on 7-day ECG patch screening. Serial screening over 12 months substantially increased the detection of AF. These data can be useful in helping identify high-risk participants for enrollment in future studies of the management of asymptomatic AF.(Women's Health Initiative Silent Atrial Fibrillation Recording Study [WHISH STAR]; NCT05366803.).

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Insights From a Large-Scale Whole-Genome Sequencing Study of Systolic Blood Pressure, Diastolic Blood Pressure, and Hypertension

BACKGROUND

The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure.

METHODS

We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants.

RESULTS

Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings (P<5×10^-8). Among them, a rare intergenic variant at novel locus, LOC100506274, was associated with lower systolic blood pressure in stage-1 (beta [SE]=-32.6 [6.0]; P=4.99×10^-8) but not stage-2 analysis (P=0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=-0.36 [0.07]; P=4.18×10^-7) and attained genome-wide significance in stage-2 (beta [SE]=-0.29 [0.03]; P=7.28×10^-23). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings (P<1×10^-6 and P<1×10^-4, respectively).

DISCUSSION

We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

Enabling ultra-low-voltage switching in BaTiO3

Single crystals of BaTiO₃ exhibit small switching fields and energies, but thin-film performance is considerably worse, thus precluding their use in next-generation devices. Here, we demonstrate high-quality BaTiO₃ thin films with nearly bulk-like properties. Thickness scaling provides access to the coercive voltages (<100 mV) and fields (<10 kV cm^-1) required for future applications and results in a switching energy of <2 J cm^-3 (corresponding to <2 aJ per bit in a 10 × 10 × 10 nm³ device). While reduction in film thickness reduces coercive voltage, it does so at the expense of remanent polarization. Depolarization fields impact polar state stability in thicker films but fortunately suppress the coercive field, thus driving a deviation from Janovec-Kay-Dunn scaling and enabling a constant coercive field for films <150 nm in thickness. Switching studies reveal fast speeds (switching times of ~2 ns for 25-nm-thick films with 5-µm-diameter capacitors) and a pathway to subnanosecond switching. Finally, integration of BaTiO₃ thin films onto silicon substrates is shown. We also discuss what remains to be demonstrated to enable the use of these materials for next-generation devices.

The association of walking pace and incident heart failure and subtypes among postmenopausal women

BACKGROUND

To investigate the association between walking pace and the risk of heart failure (HF) and HF sub-types.

METHODS

We examined associations of self-reported walking pace with risk of incident HF and HF subtypes of preserved (HFpEF) and reduced (HFrEF) ejection fractions, among 25,183 postmenopausal women, ages 50-79 years. At enrollment into the Women's Health Initiative cohort in 1993-1998, this subset of women was free of HF, cancer, or the inability to walk one block, with self-reported information on walking pace and walking duration. Multivariable Cox regression was used to examine associations of walking pace (casual <2 mph [referent], average 2-3 mph, and fast >3 mph) with incident HF. We also examined the joint association of walking pace and duration with incident HF.

RESULTS

There were 1455 incident adjudicated acute decompensated HF hospitalization cases during a median of 16.9 years of follow-up. There was a strong inverse association between walking pace and overall risk of HF (HR = 0.73, 95% CI [0.65, 0.83] for average vs. casual walking; HR = 0.66, 95%CI [0.56, 0.78] for fast vs. casual walking). There were similar associations of walking pace with HFpEF (HR = 0.73, 95%CI [0.62, 0.86] average vs. casual; HR = 0.63, 95%CI [0.50, 0.80] for fast vs. casual) and with HFrEF (HR = 0.72, 95%CI [0.57, 0.91] for average vs. casual; HR = 0.74, 95%CI [0.54, 0.99] for fast vs. casual). The risk of HF associated with fast walking with less than 1 h/week walking duration was comparable with the risk of HF among casual and average walkers with more than 2 h/week walking duration.

CONCLUSION

Walking pace was inversely associated with risks of overall HF, HFpEF, and HFrEF in postmenopausal women. Whether interventions to increase the walking pace in older adults will reduce HF risk and whether fast pace will compensate for the short duration of walking warrants further study.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Sugar-Sweetened Beverage Consumption May Modify Associations Between Genetic Variants in the CHREBP (Carbohydrate Responsive Element Binding Protein) Locus and HDL-C (High-Density Lipoprotein Cholesterol) and Triglyceride Concentrations

Supplemental Digital Content is available in the text.

Background:

ChREBP (carbohydrate responsive element binding protein) is a transcription factor that responds to sugar consumption. Sugar-sweetened beverage (SSB) consumption and genetic variants in the CHREBP locus have separately been linked to HDL-C (high-density lipoprotein cholesterol) and triglyceride concentrations. We hypothesized that SSB consumption would modify the association between genetic variants in the CHREBP locus and dyslipidemia.

Methods:

Data from 11 cohorts from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (N=63 599) and the UK Biobank (N=59 220) were used to quantify associations of SSB consumption, genetic variants, and their interaction on HDL-C and triglyceride concentrations using linear regression models. A total of 1606 single nucleotide polymorphisms within or near CHREBP were considered. SSB consumption was estimated from validated questionnaires, and participants were grouped by their estimated intake.

Results:

In a meta-analysis, rs71556729 was significantly associated with higher HDL-C concentrations only among the highest SSB consumers (β, 2.12 [95% CI, 1.16–3.07] mg/dL per allele; P<0.0001), but not significantly among the lowest SSB consumers (P=0.81; P_Diff <0.0001). Similar results were observed for 2 additional variants (rs35709627 and rs71556736). For triglyceride, rs55673514 was positively associated with triglyceride concentrations only among the highest SSB consumers (β, 0.06 [95% CI, 0.02–0.09] ln-mg/dL per allele, P=0.001) but not the lowest SSB consumers (P=0.84; P_Diff=0.0005).

Conclusions:

Our results identified genetic variants in the CHREBP locus that may protect against SSB-associated reductions in HDL-C and other variants that may exacerbate SSB-associated increases in triglyceride concentrations.

Registration:

URL: https://www.clinicaltrials.gov; Unique identifier: NCT00005133, NCT00005121, NCT00005487, and NCT00000479.

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.

When the At-Risk Do Not Develop Heart Failure: Understanding Positive Deviance Among Postmenopausal African American and Hispanic Women

BACKGROUND

African American and Hispanic postmenopausal women have the highest risk for heart failure compared with other races, but heart failure prevalence is lower than expected in some national cohorts. It is unknown whether psychosocial factors are associated with lower risk of incident heart failure hospitalization among high-risk postmenopausal minority women.

METHODS AND RESULTS

Using the Women's Health Initiative Study, African American and US Hispanic women were classified as high-risk for incident heart failure hospitalization with 1 or more traditional heart failure risk factors and the highest tertile heart failure genetic risk scores. Positive psychosocial factors (optimism, social support, religion) and negative psychosocial factors (living alone, social strain, depressive symptoms) were measured using validated survey instruments at baseline. Adjusted subdistribution hazard ratios of developing heart failure hospitalization were determined with death as a competing risk. Positive deviance indicated not developing incident heart failure hospitalization with 1 or more risk factors and the highest tertile for genetic risk. Among 7986 African American women (mean follow-up of 16 years), 27.0% demonstrated positive deviance. Among high-risk African American women, optimism was associated with modestly reduced risk of heart failure hospitalization (subdistribution hazard ratio 0.94, 95% confidence interval 0.91-0.99), and social strain was associated with modestly increased risk of heart failure hospitalization (subdistribution hazard ratio 1.07, 95% confidence interval 1.02-1.12) in the initial models; however, no psychosocial factors were associated with heart failure hospitalization in fully adjusted analyses. Among 3341 Hispanic women, 25.1% demonstrated positive deviance. Among high-risk Hispanic women, living alone was associated with increased risk of heart failure hospitalization (subdistribution hazard ratio 1.97, 95% confidence interval 1.06-3.63) in unadjusted analyses; however, no psychosocial factors were associated with heart failure hospitalization in fully adjusted analyses.

CONCLUSIONS

Among postmenopausal African American and Hispanic women, a significant proportion remained free from heart failure hospitalization despite having the highest genetic risk profile and 1 or more traditional risk factors. No observed psychosocial factors were associated with incident heart failure hospitalization in high-risk African Americans and Hispanics. Additional investigation is needed to understand protective factors among high-risk African American and Hispanic women.

Local negative permittivity and topological phase transition in polar skyrmions

Topological solitons such as magnetic skyrmions have drawn attention as stable quasi-particle-like objects. The recent discovery of polar vortices and skyrmions in ferroelectric oxide superlattices has opened up new vistas to explore topology, emergent phenomena and approaches for manipulating such features with electric fields. Using macroscopic dielectric measurements, coupled with direct scanning convergent beam electron diffraction imaging on the atomic scale, theoretical phase-field simulations and second-principles calculations, we demonstrate that polar skyrmions in (PbTiO₃)_n/(SrTiO₃)_n superlattices are distinguished by a sheath of negative permittivity at the periphery of each skyrmion. This enhances the effective dielectric permittivity compared with the individual SrTiO₃ and PbTiO₃ layers. Moreover, the response of these topologically protected structures to electric field and temperature shows a reversible phase transition from the skyrmion state to a trivial uniform ferroelectric state, accompanied by large tunability of the dielectric permittivity. Pulsed switching measurements show a time-dependent evolution and recovery of the skyrmion state (and macroscopic dielectric response). The interrelationship between topological and dielectric properties presents an opportunity to simultaneously manipulate both by a single, and easily controlled, stimulus, the applied electric field.

Multi-Ancestry Genome-wide Association Study Accounting for Gene-Psychosocial Factor Interactions Identifies Novel Loci for Blood Pressure Traits

Psychological and social factors are known to influence blood pressure (BP) and risk of hypertension and associated cardiovascular diseases. To identify novel BP loci, we carried out genome-wide association meta-analyses of systolic, diastolic, pulse, and mean arterial BP taking into account the interaction effects of genetic variants with three psychosocial factors: depressive symptoms, anxiety symptoms, and social support. Analyses were performed using a two-stage design in a sample of up to 128,894 adults from 5 ancestry groups. In the combined meta-analyses of Stages 1 and 2, we identified 59 loci (p value <5e-8), including nine novel BP loci. The novel associations were observed mostly with pulse pressure, with fewer observed with mean arterial pressure. Five novel loci were identified in African ancestry, and all but one showed patterns of interaction with at least one psychosocial factor. Functional annotation of the novel loci supports a major role for genes implicated in the immune response (PLCL2), synaptic function and neurotransmission (LIN7A, PFIA2), as well as genes previously implicated in neuropsychiatric or stress-related disorders (FSTL5, CHODL). These findings underscore the importance of considering psychological and social factors in gene discovery for BP, especially in non-European populations.

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

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

Diet quality indices and risk of metabolic syndrome among postmenopausal women of Mexican ethnic descent in the Women's Health Initiative Observational Study

BACKGROUND

The prevalence of metabolic syndrome is higher among minority populations, including individuals of Mexican ethnic descent. Whether alignment to healthy dietary patterns is associated with lower risk of metabolic syndrome in this population is largely unknown.

OBJECTIVE

To prospectively evaluate the associations between a priori diet quality scores and risk of metabolic syndrome and its components among postmenopausal women of Mexican ethnic descent.

METHODS

A total of 334 women of Mexican ethnic descent who participated in the Women's Health Initiative (WHI) observational study without metabolic syndrome or diabetes at baseline (1993-1998) were included. Baseline diets were scored with the Alternate Mediterranean Diet (aMED), the Dietary Approaches to Stop Hypertension (DASH), the Healthy Eating Index (HEI-2010), the Mediterranean Diet Score (MDS), and the traditional Mexican Diet (MexD) score. Multivariable linear and logistic regression models were used to test the associations between baseline diet quality and risk of metabolic syndrome and its individual components at follow-up (2012-2013).

RESULTS

Approximately 16% of women met the criteria for metabolic syndrome at follow-up. None of the diet quality indices were associated with risk of metabolic syndrome. However, higher vs lower DASH scores were associated with lower waist circumference (85.2 vs 88.0 cm) and glucose concentrations (90.0 vs 95.1 mg/dL), and higher HDL cholesterol (62.6 vs 59.0 mg/dL), while higher vs lower HEI-2010 scores were associated with lower waist circumference (83.9 vs 88.1 cm), triglycerides (103 vs 117 mg/dL) and glucose concentrations (89.5 vs 94.4 mg/dL), and higher HDL cholesterol levels (63.9 vs 58.5 mg/dL).

CONCLUSIONS

Diet quality was not associated with risk of metabolic syndrome in this population. However, the results suggest that alignment to DASH and HEI-2010 recommendations may be beneficial for reducing some individual components of metabolic syndrome among postmenopausal women of Mexican descent.

Walking Volume and Speed Are Inversely Associated With Incidence of Treated Hypertension in Postmenopausal Women

Few studies have evaluated hypertension incidence in relation to walking, which is a common physical activity among adults. We examined the association between walking and hypertension incidence in 83 435 postmenopausal women who at baseline were aged 50 to 79 years, without known hypertension, heart failure, coronary heart disease, or stroke, and reported the ability to walk at least one block without assistance. Walking volume (metabolic equivalent hours per week) and speed (miles per hour) were assessed by questionnaire. Incident physician-diagnosed hypertension treated with medication was ascertained through annual questionnaires. During a mean 11-year follow-up, 38 230 hypertension cases were identified. After adjustment for covariates including nonwalking activities, a significant inverse association with hypertension was observed across categories of baseline walking volume (0 [referent], >0-3.5, 3.6-7.5, and >7.5 metabolic equivalent hours per week), hazard ratio: 1.00 (referent), 0.98, 0.95, 0.89; trend P<0.001. Faster walking speeds (<2, 2-3, 3-4, and >4 miles per hour) also were associated with lower hypertension risk, hazard ratio: 1.00 (referent), 1.07, 0.95, 0.86, 0.79; trend P<0.001. Further adjustment for walking duration (h/wk) had little impact on the association for walking speed (hazard ratio: 1.00 [referent], 1.08, 0.96, 0.86, 0.77; trend P<0.001). Significant inverse associations for walking volume and speed persisted after additional control for baseline blood pressure. Results for time-varying walking were comparable to those for baseline exposures. This study showed that walking at guideline-recommended volumes (>7.5 metabolic equivalent hours per week) and at faster speeds (≥2 miles per hour) is associated with lower hypertension risk in postmenopausal women. Walking should be encouraged as part of hypertension prevention in older adults.

Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce 'annotation principal components', multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

Role of Rare and Low-Frequency Variants in Gene-Alcohol Interactions on Plasma Lipid Levels

BACKGROUND

Alcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in gene by alcohol consumption interactions associated with fasting plasma lipid levels.

METHODS

In the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, fasting plasma triglycerides and high- and low-density lipoprotein cholesterol were measured in 34 153 individuals with European ancestry from 5 discovery studies and 32 277 individuals from 6 replication studies. Rare and low-frequency functional protein-coding variants (minor allele frequency, ≤5%) measured by an exome array were aggregated by genes and evaluated by a gene-environment interaction test and a joint test of genetic main and gene-environment interaction effects. Two dichotomous self-reported alcohol consumption variables, current drinker, defined as any recurrent drinking behavior, and regular drinker, defined as the subset of current drinkers who consume at least 2 drinks per week, were considered.

RESULTS

We discovered and replicated 21 gene-lipid associations at 13 known lipid loci through the joint test. Eight loci (PCSK9, LPA, LPL, LIPG, ANGPTL4, APOB, APOC3, and CD300LG) remained significant after conditioning on the common index single-nucleotide polymorphism identified by previous genome-wide association studies, suggesting an independent role for rare and low-frequency variants at these loci. One significant gene-alcohol interaction on triglycerides in a novel locus was significantly discovered (P=6.65×10^-6 for the interaction test) and replicated at nominal significance level (P=0.013) in SMC5.

CONCLUSIONS

In conclusion, this study applied new gene-based statistical approaches and suggested that rare and low-frequency genetic variants interacted with alcohol consumption on lipid levels.

Relationship Between Dietary Magnesium Intake and Incident Heart Failure Among Older Women: The WHI

Background

Women represent a large proportion of the growing heart failure (HF) epidemic, yet data are lacking regarding optimal dietary and lifestyle prevention strategies for them. Specifically, the association between magnesium intake and HF in a multiracial cohort of women is uncertain.

Methods and Results

We included 97 725 postmenopausal women from the WHI (Women's Health Initiative) observational studies and placebo arms of the hormone trial. Magnesium intake was measured at baseline by a 122‐item validated food‐frequency questionnaire and stratified into quartiles based on diet only, total intake (diet with supplements), and residual intake (calibration by total energy). Incident hospitalized HF (2153 events, median follow‐up 8.1 years) was adjudicated by medical record abstraction. In Cox proportional hazards models, we evaluated the association between magnesium intake and HF adjusting for potential confounders. Analyses were repeated on a subcohort (n=18 745; median‐follow‐up, 13.2 years) for whom HF cases were subclassified into preserved ejection fraction (526 events), reduced ejection fraction (291 events) or unknown (168 events). Most women were white (85%) with a mean age of 63. Compared with the highest quartile of magnesium intake, women in the lowest quartile had an increased risk of incident HF, with adjusted hazard ratios of 1.32 (95% CI, 1.02–1.71) for diet only (P trend=0.03), 1.26 (95% CI, 1.03–1.56) for total intake, and 1.31 (95% CI, 1.02–1.67) for residual intake. Results did not significantly vary by race. Subcohort analyses showed low residual magnesium intake was associated with HF with reduced ejection fraction (hazard ratio, 1.81, lowest versus highest quartile; 95% CI, 1.08–3.05) but not HF with preserved ejection fraction.

Conclusions

Low magnesium intake in a multiracial cohort of postmenopausal women was associated with a higher risk of incident HF, especially HF with reduced ejection fraction.