Synthesizing subject-matter expertise for variable selection in causal effect estimation: A case study.

BACKGROUND: Causal graphs are an important tool for covariate selection but there is limited applied research on how best to create them. Here, we used data from the Coronary Drug Project (CDP) trial to assess a range of approaches to directed acyclic graph (DAG) creation. We focused on the effect of adherence on mortality in the placebo arm, since the true causal effect is believed with a high degree of certainty. METHODS: We created DAGs for the effect of placebo adherence on mortality using different approaches for identifying variables and links to include or exclude. For each DAG, we identified minimal adjustment sets of covariates for estimating our causal effect of interest, and applied these to analyses of the CDP data. RESULTS: When we used only baseline covariate values to estimate the cumulative effect of placebo adherence on mortality, all adjustment sets performed similarly. The specific choice of covariates had minimal effect on these (biased) point estimates, but including non-confounding prognostic factors resulted in smaller variance estimates. When we additionally adjusted for time-varying covariates of adherence using inverse probability weighting, covariates identified from the DAG created by focusing on prognostic factors performed best. CONCLUSION: Theoretical advice on covariate selection suggests including prognostic factors that are not exposure predictors can reduce variance without increasing bias. In contrast, for exposure predictors that are not prognostic factors, inclusion may result in less bias control. Our results empirically confirm this advice. We recommend that hand-creating DAGs begin with identification of all potential outcome-prognostic factors.

Genomic and phenotypic correlates of mosaic loss of chromosome Y in blood.

Mosaic loss of Y (mLOY) is the most common somatic chromosomal alteration detected in human blood. The presence of mLOY is associated with altered blood cell counts and increased risk of Alzheimer's disease, solid tumors, and other age-related diseases. We sought to gain a better understanding of genetic drivers and associated phenotypes of mLOY through analyses of whole genome sequencing of a large set of genetically diverse males from the Trans-Omics for Precision Medicine (TOPMed) program. This approach enabled us to identify differences in mLOY frequencies across populations defined by genetic similarity, revealing a higher frequency of mLOY in the European American (EA) ancestry group compared to those of Hispanic American (HA), African American (AA), and East Asian (EAS) ancestry. Further, we identified two genes ( CFHR1 and LRP6 ) that harbor multiple rare, putatively deleterious variants associated with mLOY susceptibility, show that subsets of human hematopoietic stem cells are enriched for activity of mLOY susceptibility variants, and that certain alleles on chromosome Y are more likely to be lost than others.

Neighborhood Disadvantage and Risk of Heart Failure: The Reasons for Geographic and Racial Differences in Stroke (REGARDS) Study.

BACKGROUND: Heart failure (HF) affects >6 million US adults, with recent increases in HF hospitalizations. We aimed to investigate the association between neighborhood disadvantage and incident HF events and potential differences by diabetes status. METHODS: We included 23 645 participants from the REGARDS study (Reasons for Geographic and Racial Differences in Stroke), a prospective cohort of Black and White adults aged ≥45 years living in the continental United States (baseline 2005-2007). Neighborhood disadvantage was assessed using a Z score of 6 census tract variables (2000 US Census) and categorized as quartiles. Incident HF hospitalizations or HF-related deaths through 2017 were adjudicated. Multivariable-adjusted Cox regression was used to examine the association between neighborhood disadvantage and incident HF. Heterogeneity by diabetes was assessed using an interaction term. RESULTS: The mean age was 64.4 years, 39.5% were Black adults, 54.9% females, and 18.8% had diabetes. During a median follow-up of 10.7 years, there were 1125 incident HF events with an incidence rate of 3.3 (quartile 1), 4.7 (quartile 2), 5.2 (quartile 3), and 6.0 (quartile 4) per 1000 person-years. Compared to adults living in the most advantaged neighborhoods (quartile 1), those living in neighborhoods in quartiles 2, 3, and 4 (most disadvantaged) had 1.30 (95% CI, 1.06-1.60), 1.36 (95% CI, 1.11-1.66), and 1.45 (95% CI, 1.18-1.79) times greater hazard of incident HF even after accounting for known confounders. This association did not significantly differ by diabetes status (interaction P=0.59). For adults with diabetes, the adjusted incident HF hazards comparing those in quartile 4 versus quartile 1 was 1.34 (95% CI, 0.92-1.96), and it was 1.50 (95% CI, 1.16-1.94) for adults without diabetes. CONCLUSIONS: In this large contemporaneous prospective cohort, neighborhood disadvantage was associated with an increased risk of incident HF events. This increase in HF risk did not differ by diabetes status. Addressing social, economic, and structural factors at the neighborhood level may impact HF prevention.

Epigenetic age in African American adolescents with type 2 diabetes: A cross-sectional case-control study protocol.

Background and Aims: Type 2 diabetes (T2D) is a disease caused by a relative insulin deficiency compared to the significant insulin requirement needed by the body to achieve glycemic control. T2D in adolescence appears to be increasing in prevalence over the past several decades, necessitating studies to understand for the onset of the disease to occur early in the lifespan. Given the high burden of disease, specifically in young African American adolescents, our study chose to focus initially on feasibility of recruitment of this population. Methods: Data was collected at a single study center at Children's of Alabama. The protocol was completed as part of routine care or at a study visit. The study team was able to leverage the Electronic Medical Record to prescreen eligible patients to discuss the study. A variety of times of day were utilized to improve likely of success with reaching potential participants. Inclusion criteria for patients with T2D was focused on the adolescent population (ages 12-18 years), with no history of an obesity syndrome. DNA methylation age will be calculated using the EPIC 850K array. Statistical analysis will be done using linear regression analysis, adjusting for covariates. Conclusions: This study's aim was to screen and enroll young African American adolescents for a study investigating epigenetic aging and T2D. Our study found that more direct contact (face-to-face- or phone call) improve success of recruitment. Leveraging the electronic medical record also helped improve success with pre-screening participants. Challenges included recruiting participants who might come from long distances to a tertiary care center. Consolidating appointments helped improve the success of reaching these participants. Other challenges included frequent address changes and changed phone numbers. Close attention to the barriers as well as the successes will aid in understanding effective strategies for this important population.

Neighborhood Environment and Epigenetic Age: A Scoping Review.

BACKGROUND: Interest in how the neighborhood environment impacts age-related health conditions has been increasing for decades. Epigenetic changes are environmentally derived modifications to the genome that alter the way genes function-thus altering health status. Epigenetic age, a biomarker for biological age, has been shown to be a useful predictor of several age-related health conditions. Consequently, its relation to the neighborhood environment has been the focus of a growing body of literature. OBJECTIVE: We aimed to describe the scope of the evidence on the relationship between neighborhood environmental characteristics and epigenetic age. METHODS: Using scoping review following methods established by Arksey and O'Malley, we first defined our research questions and searched the literature in PubMed, PsycINFO, and EMBASE. Next, we selected the literature to be included, and finally, we analyzed and summarized the information. RESULTS: Nine articles met the inclusion criteria. Most studies examined deprivation as the neighborhood characteristic of interest. While all studies were observational in design, the articles included diverse participants, including men and women, adults and children, and multiple ethnicities. Results demonstrated a relationship between the neighborhood environment and epigenetic age, whether the characteristic of interest is socioeconomic or physical. CONCLUSIONS: Overall, studies concluded there was a relationship between neighborhood characteristics and epigenetic age, whether the characteristic of interest was socioeconomic or physical. However, findings varied based on how the neighborhood characteristic and/or epigenetic age was measured. Furthermore, a paucity of investigations on physical characteristics was noticeable and warrants increased attention.

WHOLE GENOME SEQUENCING ANALYSIS OF BODY MASS INDEX IDENTIFIES NOVEL AFRICAN ANCESTRY-SPECIFIC RISK ALLELE.

Obesity is a major public health crisis associated with high mortality rates. Previous genome-wide association studies (GWAS) investigating body mass index (BMI) have largely relied on imputed data from European individuals. This study leveraged whole-genome sequencing (WGS) data from 88,873 participants from the Trans-Omics for Precision Medicine (TOPMed) Program, of which 51% were of non-European population groups. We discovered 18 BMI-associated signals (P < 5 × 10-9). Notably, we identified and replicated a novel low frequency single nucleotide polymorphism (SNP) in MTMR3 that was common in individuals of African descent. Using a diverse study population, we further identified two novel secondary signals in known BMI loci and pinpointed two likely causal variants in the POC5 and DMD loci. Our work demonstrates the benefits of combining WGS and diverse cohorts in expanding current catalog of variants and genes confer risk for obesity, bringing us one step closer to personalized medicine.

Epigenetic age acceleration mediates the relationship between neighborhood deprivation and pain severity in adults with or at risk for knee osteoarthritis pain.

An estimated 250 million people worldwide suffer from knee osteoarthritis (KOA), with older adults having greater risk. Like other age-related diseases, residents of high-deprivation neighborhoods experience worse KOA pain outcomes compared to their more affluent neighbors. The purpose of this study was to examine the relationship between neighborhood deprivation and pain severity in KOA and the influence of epigenetic age acceleration (EpAA) on that relationship. The sample of 128 participants was mostly female (60.9%), approximately half non-Hispanic Black (49.2%), and had a mean age of 58 years. Spearman bivariate correlations revealed that pain severity positively correlated with EpAA (ρ = 0.47, p ≤ 0.001) and neighborhood deprivation (ρ = 0.25, p = 0.004). We found a positive significant relationship between neighborhood deprivation and EpAA (ρ = 0.47, p ≤ 0.001). Results indicate a mediating relationship between neighborhood deprivation (predictor), EpAA (mediator), and pain severity (outcome variable). There was a significant indirect effect of neighborhood deprivation on pain severity through EpAA, as the mediator accounted for a moderate portion of the total effect, PM = 0.44. Epigenetic age acceleration may act as a mechanism through which neighborhood deprivation leads to worse KOA pain outcomes and may play a role in the well-documented relationship between the neighborhood of residence and age-related diseases.

Still Exhausted: The Role of Residual Caregiving Fatigue on Women in Medicine and Science Across the Pipeline.

Understanding the impact of caregiving responsibilities on women in medicine is crucial for ensuring a healthy and intact workforce, as caregiving responsibilities have the potential to affect the careers of women in health care along the entire pipeline, from students and trainees to physicians, physician-scientists, and biomedical researchers.

The CTSA Diversity, Equity, Inclusion, and Accessibility (DEIA) Task Force's recommendations for the CTSA program consortium.

The Clinical and Translational Science Award (CTSA) Program recognizes that advancing diversity, equity, inclusion, and accessibility (DEIA) requires moving beyond statements of commitment to transformative actions. In 2021, the CTSA Program created a Task Force (TF) to initiate work in support of structural and transformational initiatives that advance DEIA for the consortium and its individual hubs. We describe the process of forming the expertise-driven (DEIA) TF and our activities to date. We 1) developed and adopted the DEIA Learning Systems Framework to guide our approach; 2) defined a set of recommendations across four focus areas (Institutional; Programmatic; Community-Centered; and Social, Cultural, Environmental); and 3) designed and disseminated a survey to capture the CTSA Program's baseline demographic, community, infrastructural, and leadership diversity. The CTSA Consortium also elevated the TF to a standing Committee to extend our understanding, development, and implementation of DEIA approaches to translational and clinical science. These initial steps provide a foundation for collectively fostering environment that support DEIA across the research continuum.

Investigating Gene-Diet Interactions Impacting the Association Between Macronutrient Intake and Glycemic Traits.

Few studies have demonstrated reproducible gene-diet interactions (GDIs) impacting metabolic disease risk factors, likely due in part to measurement error in dietary intake estimation and insufficient capture of rare genetic variation. We aimed to identify GDIs across the genetic frequency spectrum impacting the macronutrient-glycemia relationship in genetically and culturally diverse cohorts. We analyzed 33,187 participants free of diabetes from 10 National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program cohorts with whole-genome sequencing, self-reported diet, and glycemic trait data. We fit cohort-specific, multivariable-adjusted linear mixed models for the effect of diet, modeled as an isocaloric substitution of carbohydrate for fat, and its interactions with common and rare variants genome-wide. In main effect meta-analyses, participants consuming more carbohydrate had modestly lower glycemic trait values (e.g., for glycated hemoglobin [HbA1c], -0.013% HbA1c/250 kcal substitution). In GDI meta-analyses, a common African ancestry-enriched variant (rs79762542) reached study-wide significance and replicated in the UK Biobank cohort, indicating a negative carbohydrate-HbA1c association among major allele homozygotes only. Simulations revealed that >150,000 samples may be necessary to identify similar macronutrient GDIs under realistic assumptions about effect size and measurement error. These results generate hypotheses for further exploration of modifiable metabolic disease risk in additional cohorts with African ancestry. ARTICLE HIGHLIGHTS: We aimed to identify genetic modifiers of the dietary macronutrient-glycemia relationship using whole-genome sequence data from 10 Trans-Omics for Precision Medicine program cohorts. Substitution models indicated a modest reduction in glycemia associated with an increase in dietary carbohydrate at the expense of fat. Genome-wide interaction analysis identified one African ancestry-enriched variant near the FRAS1 gene that may interact with macronutrient intake to influence hemoglobin A1c. Simulation-based power calculations accounting for measurement error suggested that substantially larger sample sizes may be necessary to discover further gene-macronutrient interactions.

Obesity and Overweight: Probing Causes, Consequences, and Novel Therapeutic Approaches Through the American Heart Association's Strategically Focused Research Network.

As the worldwide prevalence of overweight and obesity continues to rise, so too does the urgency to fully understand mediating mechanisms, to discover new targets for safe and effective therapeutic intervention, and to identify biomarkers to track obesity and the success of weight loss interventions. In 2016, the American Heart Association sought applications for a Strategically Focused Research Network (SFRN) on Obesity. In 2017, 4 centers were named, including Johns Hopkins University School of Medicine, New York University Grossman School of Medicine, University of Alabama at Birmingham, and Vanderbilt University Medical Center. These 4 centers were convened to study mechanisms and therapeutic targets in obesity, to train a talented cadre of American Heart Association SFRN-designated fellows, and to initiate and sustain effective and enduring collaborations within the individual centers and throughout the SFRN networks. This review summarizes the central themes, major findings, successful training of highly motivated and productive fellows, and the innovative collaborations and studies forged through this SFRN on Obesity. Leveraging expertise in in vitro and cellular model assays, animal models, and humans, the work of these 4 centers has made a significant impact in the field of obesity, opening doors to important discoveries, and the identification of a future generation of obesity-focused investigators and next-step clinical trials. The creation of the SFRN on Obesity for these 4 centers is but the beginning of innovative science and, importantly, the birth of new collaborations and research partnerships to propel the field forward.

Whole-exome sequence analysis of anthropometric traits illustrates challenges in identifying effects of rare genetic variants.

Anthropometric traits, measuring body size and shape, are highly heritable and significant clinical risk factors for cardiometabolic disorders. These traits have been extensively studied in genome-wide association studies (GWASs), with hundreds of genome-wide significant loci identified. We performed a whole-exome sequence analysis of the genetics of height, body mass index (BMI) and waist/hip ratio (WHR). We meta-analyzed single-variant and gene-based associations of whole-exome sequence variation with height, BMI, and WHR in up to 22,004 individuals, and we assessed replication of our findings in up to 16,418 individuals from 10 independent cohorts from Trans-Omics for Precision Medicine (TOPMed). We identified four trait associations with single-nucleotide variants (SNVs; two for height and two for BMI) and replicated the LECT2 gene association with height. Our expression quantitative trait locus (eQTL) analysis within previously reported GWAS loci implicated CEP63 and RFT1 as potential functional genes for known height loci. We further assessed enrichment of SNVs, which were monogenic or syndromic variants within loci associated with our three traits. This led to the significant enrichment results for height, whereas we observed no Bonferroni-corrected significance for all SNVs. With a sample size of ∼20,000 whole-exome sequences in our discovery dataset, our findings demonstrate the importance of genomic sequencing in genetic association studies, yet they also illustrate the challenges in identifying effects of rare genetic variants.

Epigenetic Age Acceleration in Mothers and Offspring 4-10 Years after a Pregnancy Complicated by Gestational Diabetes and Obesity.

A known association exists between exposure to gestational diabetes mellitus (GDM) and epigenetic age acceleration (EAA) in GDM-exposed offspring compared to those without GDM exposure. This association has not been assessed previously in mothers with pregnancies complicated by GDM. A total of 137 mother-child dyads with an index pregnancy 4−10 years before study enrollment were included. Clinical data and whole blood samples were collected and quantified to obtain DNA methylation (DNAm) estimates using the Illumina MethylEPIC 850K array in mothers and offspring. DNAm age and age acceleration were evaluated using the Horvath and Hannum clocks. Multivariable linear regression models were performed to determine the association between EAA and leptin, high-density lipoprotein cholesterol (HDL-C), fasting glucose, fasting insulin, and HOMA-IR. Mothers with a GDM and non-GDM pregnancy had strong correlations between chronological age and DNAm age (r > 0.70). Offspring of GDM mothers had moderate to strong correlations, whereas offspring of non-GDM mothers had moderate correlations between chronological age and DNAm age. Association analyses revealed a significant association between EAA and fasting insulin in offspring (FDR < 0.05), while HDL-C was the only metabolic marker significantly associated with EAA in mothers (FDR < 0.05). Mothers in the GDM group had a higher predicted epigenetic age and age acceleration than mothers in the non-GDM group. The association between EAA with elevated fasting insulin in offspring and elevated HDL-C in mothers suggests possible biomarkers that can better elucidate the effects of exposure to a GDM pregnancy and future cardiometabolic outcomes.

SMOC2 gene interacts with APOL1 in the development of end-stage kidney disease: A genome-wide association study.

Background: Some but not all African-Americans (AA) who carry APOL1 nephropathy risk variants (APOL1) develop kidney failure (end-stage kidney disease, ESKD). To identify genetic modifiers, we assessed gene-gene interactions in a large prospective cohort of the REasons for Geographic and Racial Differences in Stroke (REGARDS) study. Methods: Genotypes from 8,074 AA participants were obtained from Illumina Infinium Multi-Ethnic AMR/AFR Extended BeadChip. We compared 388 incident ESKD cases with 7,686 non-ESKD controls, using a two-locus interaction approach. Logistic regression was used to examine the effect of APOL1 risk status (using recessive and additive models), single nucleotide polymorphism (SNP), and APOL1*SNP interaction on incident ESKD, adjusting for age, sex, and ancestry. APOL1 *SNP interactions that met the threshold of 1.0 × 10-5 were replicated in the Genetics of Hypertension Associated Treatment (GenHAT) study (626 ESKD cases and 6,165 controls). In a sensitivity analysis, models were additionally adjusted for diabetes status. We conducted additional replication in the BioVU study. Results: Two APOL1 risk alleles prevalence (recessive model) was similar in the REGARDS and GenHAT studies. Only one APOL1-SNP interaction, for rs7067944 on chromosome 10, ~10 KB from the PCAT5 gene met the genome-wide statistical threshold (P interaction = 3.4 × 10-8), but this interaction was not replicated in the GenHAT study. Among other relevant top findings (with P interaction < 1.0 × 10-5), a variant (rs2181251) near SMOC2 on chromosome six interacted with APOL1 risk status (additive) on ESKD outcomes (REGARDS study, P interaction =5.3 × 10-6) but the association was not replicated (GenHAT study, P interaction = 0.07, BioVU study, P interaction = 0.53). The association with the locus near SMOC2 persisted further in stratified analyses. Among those who inherited ≥1 alternate allele of rs2181251, APOL1 was associated with an increased risk of incident ESKD (OR [95%CI] = 2.27[1.53, 3.37]) but APOL1 was not associated with ESKD in the absence of the alternate allele (OR [95%CI] = 1.34[0.96, 1.85]) in the REGARDS study. The associations were consistent after adjusting for diabetes. Conclusion: In a large genome-wide association study of AAs, a locus SMOC2 exhibited a significant interaction with the APOL1 locus. SMOC2 contributes to the progression of fibrosis after kidney injury and the interaction with APOL1 variants may contribute to an explanation for why only some APOLI high-risk individuals develop ESKD.

Mother-child cardiometabolic health 4-10 years after pregnancy complicated by obesity with and without gestational diabetes.

Objective: Obesity in pregnancy and gestational diabetes (GDM) increase cardiometabolic disease risk but are difficult to disentangle. This study aimed to test the hypothesis that 4-10 years after a pregnancy complicated by overweight/obesity and GDM (OB-GDM), women and children would have greater adiposity and poorer cardiometabolic health than those with overweight/obesity (OB) or normal weight (NW) and no GDM during the index pregnancy. Methods: In this cross-sectional study, mother-child dyads were stratified into three groups based on maternal health status during pregnancy (OB-GDM = 67; OB = 76; NW = 76). Weight, height, waist and hip circumferences, and blood pressure were measured, along with fasting glucose, insulin, HbA1c, lipids, adipokines, and cytokines. Results: Women in the OB and OB-GDM groups had greater current adiposity and poorer cardiometabolic health outcomes than those in the NW group (p < 0.05). After adjusting for current adiposity, women in the OB-GDM group had higher HbA1c, glucose, HOMA-IR and triglycerides than NW and OB groups (p < 0.05). Among children, adiposity was greater in the OB-GDM versus NW group (p < 0.05), but other indices of cardiometabolic health did not differ. Conclusions: Poor cardiometabolic health in women with prior GDM is independent of current adiposity. Although greater adiposity among children exposed to GDM is evident at 4-10 years, differences in cardiometabolic health may not emerge until later.

Whole genome sequence analysis of blood lipid levels in >66,000 individuals.

Blood lipids are heritable modifiable causal factors for coronary artery disease. Despite well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing (WGS), partly due to limited sample sizes, ancestral diversity, and interpretation of clinical significance. Among 66,329 ancestrally diverse (56% non-European) participants, we associate 428M variants from deep-coverage WGS with lipid levels; ~400M variants were not assessed in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with blood lipids through analysis of common and rare coding variants. We discover several associated rare non-coding variants, largely at Mendelian lipid genes. Notably, we observe rare LDLR intronic variants associated with markedly increased LDL-C, similar to rare LDLR exonic variants. In conclusion, we conducted a systematic whole genome scan for blood lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.

Recommendations on the use and reporting of race, ethnicity, and ancestry in genetic research: Experiences from the NHLBI TOPMed program.

How race, ethnicity, and ancestry are used in genomic research has wide-ranging implications for how research is translated into clinical care and incorporated into public understanding. Correlation between race and genetic ancestry contributes to unresolved complexity for the scientific community, as illustrated by heterogeneous definitions and applications of these variables. Here, we offer commentary and recommendations on the use of race, ethnicity, and ancestry across the arc of genetic research, including data harmonization, analysis, and reporting. While informed by our experiences as researchers affiliated with the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, these recommendations are applicable to basic and translational genomic research in diverse populations with genome-wide data. Moving forward, considerable collaborative effort will be required to ensure that race, ethnicity, and ancestry are described and used appropriately to generate scientific knowledge that yields broad and equitable benefit.

Genetic Contributors of Efficacy and Adverse Metabolic Effects of Chlorthalidone in African Americans from the Genetics of Hypertension Associated Treatments (GenHAT) Study.

Hypertension is a leading risk factor for cardiovascular disease mortality. African Americans (AAs) have the highest prevalence of hypertension in the United States, and to alleviate the burden of hypertension in this population, better control of blood pressure (BP) is needed. Previous studies have shown considerable interpersonal differences in BP response to antihypertensive treatment, suggesting a genetic component. Utilizing data from 4297 AA participants randomized to chlorthalidone from the Genetics of Hypertension Associated Treatments (GenHAT) study, we aimed to identify variants associated with the efficacy of chlorthalidone. An additional aim was to find variants that contributed to changes in fasting glucose (FG) in these individuals. We performed genome-wide association analyses on the change of systolic and diastolic BP (SBP and DBP) over six months and FG levels over 24 months of treatment. We sought replication in the International Consortia of Pharmacogenomics Studies. We identified eight variants statistically associated with BP response and nine variants associated with FG response. One suggestive LINC02211-CDH9 intergenic variant was marginally replicated with the same direction of effect. Given the impact of hypertension in AAs, this study implies that understanding the genetic background for BP control and glucose changes during chlorthalidone treatment may help prevent adverse cardiovascular events in this population.