Associations of plasma neurofilament light chain with cognition and neuroimaging measures in community-dwelling early old age men

BACKGROUND

Plasma neurofilament light chain (NfL) is a promising biomarker of neurodegeneration with potential clinical utility in monitoring the progression of neurodegenerative diseases. However, the cross-sectional associations of plasma NfL with measures of cognition and brain have been inconsistent in community-dwelling populations.

METHODS

We examined these associations in a large community-dwelling sample of early old age men (N = 969, mean age = 67.57 years, range = 61-73 years), who are either cognitively unimpaired (CU) or with mild cognitive impairment (MCI). Specifically, we investigated five cognitive domains (executive function, episodic memory, verbal fluency, processing speed, visual-spatial ability), as well as neuroimaging measures of gray and white matter.

RESULTS

After adjusting for age, health status, and young adult general cognitive ability, plasma NfL level was only significantly associated with processing speed and white matter hyperintensity (WMH) volume, but not with other cognitive or neuroimaging measures. The association with processing speed was driven by individuals with MCI, as it was not detected in CU individuals.

CONCLUSIONS

These results suggest that in early old age men without dementia, plasma NfL does not appear to be sensitive to cross-sectional individual differences in most domains of cognition or neuroimaging measures of gray and white matter. The revealed plasma NfL associations were limited to WMH for all participants and processing speed only within the MCI cohort. Importantly, considering cognitive status in community-based samples will better inform the interpretation of the relationships of plasma NfL with cognition and brain and may help resolve mixed findings in the literature.

Using Instrumental Variables to Measure Causation over Time in Cross-Lagged Panel Models

Cross-lagged panel models (CLPMs) are commonly used to estimate causal influences between two variables with repeated assessments. The lagged effects in a CLPM depend on the time interval between assessments, eventually becoming undetectable at longer intervals. To address this limitation, we incorporate instrumental variables (IVs) into the CLPM with two study waves and two variables. Doing so enables estimation of both the lagged (i.e., "distal") effects and the bidirectional cross-sectional (i.e., "proximal") effects at each wave. The distal effects reflect Granger-causal influences across time, which decay with increasing time intervals. The proximal effects capture causal influences that accrue over time and can help infer causality when the distal effects become undetectable at longer intervals. Significant proximal effects, with a negligible distal effect, would imply that the time interval is too long to estimate a lagged effect at that time interval using the standard CLPM. Through simulations and an empirical application, we demonstrate the impact of time intervals on causal inference in the CLPM and present modeling strategies to detect causal influences regardless of the time interval in a study. Furthermore, to motivate empirical applications of the proposed model, we highlight the utility and limitations of using genetic variables as IVs in large-scale panel studies.

Changes in affect longitudinally mediate associations between emotion regulation strategy use and disordered eating

BACKGROUND

Trait-level emotion regulation (ER) difficulties are associated with eating disorders (EDs) transdiagnostically. However, little research has examined whether within-person fluctuations in ER longitudinally predict ED behaviors in daily life or the mechanisms of ER effects. Investigating daily ER could help us better understand why people experience ED behaviors at a given time. We examined whether day-to-day changes in adaptive (e.g., cognitive reappraisal) and maladaptive (e.g., rumination) ER longitudinally predicted core ED behaviors (binge eating, purging, dieting) and whether changes in affect mediated effects.

METHOD

Female participants (N = 688) ages 15-30 from the Michigan State University Twin Registry reported their adaptive and maladaptive ER use, negative affect (NA), positive affect (PA), binge eating, purging, and dieting on 49 consecutive days. Using structural equation modeling, we examined whether within-person fluctuations in ER predicted same- and next-day ED behaviors and whether changes in affect mediated longitudinal ER effects.

RESULTS

Greater maladaptive ER predicted increased likelihood of same-day binge eating and next-day binge eating and purging. The association between maladaptive ER and next-day binge eating and purging was mediated by increased next-day NA. In contrast, dieting was more closely related to changes in PA. Adaptive ER did not predict reduced likelihood of any ED behavior.

CONCLUSIONS

Maladaptive ER may longitudinally increase risk for binge eating and purging by amplifying NA. Interventions focused on decreasing maladaptive ER and subsequent NA might help disrupt binge eating-purging cycles. Conversely, results add to evidence that PA fluctuations may play a unique role in maintaining restrictive behaviors.

PUBLIC SIGNIFICANCE

Little is known about how daily changes in emotion regulation may impact disordered eating. We found that maladaptive emotion regulation (e.g., rumination) was associated with a higher likelihood of binge eating and purging on the next day because it predicted increased next-day negative affect. In contrast, dieting was more closely tied to fluctuations in positive affect. Targeting daily emotion regulation and affective processes may help disrupt cycles of disordered eating.

The heritability of blood-based biomarkers related to risk of Alzheimer's disease in a population-based sample of early old-age men

INTRODUCTION

Despite their increased application, the heritability of Alzheimer's disease (AD)-related blood-based biomarkers remains unexplored.

METHODS

Plasma amyloid beta 40 (Aβ40), Aβ42, the Aβ42/40 ratio, total tau (t-tau), and neurofilament light (NfL) data came from 1035 men 60 to 73 years of age (μ = 67.0, SD = 2.6). Twin models were used to calculate heritability and the genetic and environmental correlations between them.

RESULTS

Additive genetics explained 44% to 52% of Aβ42, Aβ40, t-tau, and NfL. The Aβ42/40 ratio was not heritable. Aβ40 and Aβ42 were genetically near identical (rg  = 0.94). Both Aβ40 and Aβ42 were genetically correlated with NfL (rg  = 0.35 to 0.38), but genetically unrelated to t-tau.

DISCUSSION

Except for Aβ42/40, plasma biomarkers are heritable. Aβ40 and Aβ42 share mostly the same genetic influences, whereas genetic influences on plasma t-tau and NfL are largely unique in early old-age men. The absence of genetic associations between the Aβs and t-tau is not consistent with the amyloid cascade hypothesis.

Power, measurement error, and pleiotropy robustness in twin-design extensions to Mendelian Randomization

Mendelian Randomization (MR) has become an important tool for causal inference in the health sciences. It takes advantage of the random segregation of alleles to control for background confounding factors. In brief, the method works by using genetic variants as instrumental variables, but it depends on the assumption of exclusion restriction, i.e., that the variants affect the outcome exclusively via the exposure variable. Equivalently, the assumption states that there is no horizontal pleiotropy from the variant to the outcome. This assumption is unlikely to hold in nature, so several extensions to MR have been developed to increase its robustness against horizontal pleiotropy, though not eliminating the problem entirely (Sanderson et al. 2022). The Direction of Causation (DoC) model, which affords information from the cross-twin cross-trait correlations to estimate causal paths, was extended with polygenic scores to explicitly model horizontal pleiotropy and a causal path (MR-DoC, Minică et al 2018). MR-DoC was further extended to accommodate bidirectional causation (MR-DoC2 ; Castro-de-Araujo et al. 2023). In the present paper, we compared the power of the DoC model, MR-DoC, and MR-DoC2. We investigated the effect of phenotypic measurement error and the effect of misspecification of unshared (individual-specific) environmental factors on the parameter estimates.

Higher cortical thickness/volume in Alzheimer's-related regions: protective factor or risk factor?

Some evidence suggests a biphasic pattern of changes in cortical thickness wherein higher, rather than lower, thickness is associated with very early Alzheimer's disease (AD) pathology. We examined whether integrating information from AD brain signatures based on mean diffusivity (MD) can aid in the interpretation of cortical thickness/volume as a risk factor for future AD-related changes. Participants were 572 men in the Vietnam Era Twin Study of Aging who were cognitively unimpaired at baseline (mean age = 56 years; range = 51-60). Individuals with both high thickness/volume signatures and high MD signatures at baseline had lower cortical thickness/volume in AD signature regions and lower episodic memory performance 12 years later compared to those with high thickness/volume and low MD signatures at baseline. Groups did not differ in level of young adult cognitive reserve. Our findings are in line with a biphasic model in which increased cortical thickness may precede future decline and establish the value of examining cortical MD alongside cortical thickness to identify subgroups with differential risk for poorer brain and cognitive outcomes.

Genetic and Environmental Influences on Structural and Diffusion-Based Alzheimer's Disease Neuroimaging Signatures Across Midlife and Early Old Age

BACKGROUND

Composite scores of magnetic resonance imaging-derived metrics in brain regions associated with Alzheimer's disease (AD), commonly termed AD signatures, have been developed to distinguish early AD-related atrophy from normal age-associated changes. Diffusion-based gray matter signatures may be more sensitive to early AD-related changes compared with thickness/volume-based signatures, demonstrating their potential clinical utility. The timing of early (i.e., midlife) changes in AD signatures from different modalities and whether diffusion- and thickness/volume-based signatures each capture unique AD-related phenotypic or genetic information remains unknown.

METHODS

Our validated thickness/volume signature, our novel mean diffusivity (MD) signature, and a magnetic resonance imaging-derived measure of brain age were used in biometrical analyses to examine genetic and environmental influences on the measures as well as phenotypic and genetic relationships between measures over 12 years. Participants were 736 men from 3 waves of the Vietnam Era Twin Study of Aging (VETSA) (baseline/wave 1: mean age [years] = 56.1, SD = 2.6, range = 51.1-60.2). Subsequent waves occurred at approximately 5.7-year intervals.

RESULTS

MD and thickness/volume signatures were highly heritable (56%-72%). Baseline MD signatures predicted thickness/volume signatures over a decade later, but baseline thickness/volume signatures showed a significantly weaker relationship with future MD signatures. AD signatures and brain age were correlated, but each measure captured unique phenotypic and genetic variance.

CONCLUSIONS

Cortical MD and thickness/volume AD signatures are heritable, and each signature captures unique variance that is also not explained by brain age. Moreover, results are in line with changes in MD emerging before changes in cortical thickness, underscoring the utility of MD as a very early predictor of AD risk.

Effects of smoking on genome-wide DNA methylation profiles: A study of discordant and concordant monozygotic twin pairs

Background

Smoking-associated DNA methylation levels identified through epigenome-wide association studies (EWASs) are generally ascribed to smoking-reactive mechanisms, but the contribution of a shared genetic predisposition to smoking and DNA methylation levels is typically not accounted for.

Methods

We exploited a strong within-family design, that is, the discordant monozygotic twin design, to study reactiveness of DNA methylation in blood cells to smoking and reversibility of methylation patterns upon quitting smoking. Illumina HumanMethylation450 BeadChip data were available for 769 monozygotic twin pairs (mean age = 36 years, range = 18-78, 70% female), including pairs discordant or concordant for current or former smoking.

Results

In pairs discordant for current smoking, 13 differentially methylated CpGs were found between current smoking twins and their genetically identical co-twin who never smoked. Top sites include multiple CpGs in CACNA1D and GNG12, which encode subunits of a calcium voltage-gated channel and G protein, respectively. These proteins interact with the nicotinic acetylcholine receptor, suggesting that methylation levels at these CpGs might be reactive to nicotine exposure. All 13 CpGs have been previously associated with smoking in unrelated individuals and data from monozygotic pairs discordant for former smoking indicated that methylation patterns are to a large extent reversible upon smoking cessation. We further showed that differences in smoking level exposure for monozygotic twins who are both current smokers but differ in the number of cigarettes they smoke are reflected in their DNA methylation profiles.

Conclusions

In conclusion, by analysing data from monozygotic twins, we robustly demonstrate that DNA methylation level in human blood cells is reactive to cigarette smoking.

Funding

We acknowledge funding from the National Institute on Drug Abuse grant DA049867, the Netherlands Organization for Scientific Research (NWO): Biobanking and Biomolecular Research Infrastructure (BBMRI-NL, NWO 184.033.111) and the BBRMI-NL-financed BIOS Consortium (NWO 184.021.007), NWO Large Scale infrastructures X-Omics (184.034.019), Genotype/phenotype database for behaviour genetic and genetic epidemiological studies (ZonMw Middelgroot 911-09-032); Netherlands Twin Registry Repository: researching the interplay between genome and environment (NWO-Groot 480-15-001/674); the Avera Institute, Sioux Falls (USA), and the National Institutes of Health (NIH R01 HD042157-01A1, MH081802, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995); epigenetic data were generated at the Human Genomics Facility (HuGe-F) at ErasmusMC Rotterdam. Cotinine assaying was sponsored by the Neuroscience Campus Amsterdam. DIB acknowledges the Royal Netherlands Academy of Science Professor Award (PAH/6635).

Environment-by-PGS Interaction in the Classical Twin Design: An Application to Childhood Anxiety and Negative Affect

One type of genotype-environment interaction occurs when genetic effects on a phenotype are moderated by an environment; or when environmental effects on a phenotype are moderated by genes. Here we outline these types of genotype-environment interaction models, and propose a test of genotype-environment interaction based on the classical twin design, which includes observed genetic variables (polygenic scores: PGSs) that account for part of the genetic variance of the phenotype. We introduce environment-by-PGS interaction and the results of a simulation study to address statistical power and parameter recovery. Next, we apply the model to empirical data on anxiety and negative affect in children. The power to detect environment-by-PGS interaction depends on the heritability of the phenotype, and the strength of the PGS. The simulation results indicate that under realistic conditions of sample size, heritability and strength of the interaction, the environment-by-PGS model is a viable approach to detect genotype-environment interaction. In 7-year-old children, we defined two PGS based on the largest genetic association studies for 2 traits that are genetically correlated to childhood anxiety and negative affect, namely major depression (MDD) and intelligence (IQ). We find that common environmental influences on negative affect are amplified for children with a lower IQ-PGS.

Researching COVID to enhance recovery (RECOVER) pediatric study protocol: Rationale, objectives and design

Importance

The prevalence, pathophysiology, and long-term outcomes of COVID-19 (post-acute sequelae of SARS-CoV-2 [PASC] or "Long COVID") in children and young adults remain unknown. Studies must address the urgent need to define PASC, its mechanisms, and potential treatment targets in children and young adults.

Observations

We describe the protocol for the Pediatric Observational Cohort Study of the NIH's RE searching COV ID to E nhance R ecovery (RECOVER) Initiative. RECOVER-Pediatrics is an observational meta-cohort study of caregiver-child pairs (birth through 17 years) and young adults (18 through 25 years), recruited from more than 100 sites across the US. This report focuses on two of five cohorts that comprise RECOVER-Pediatrics: 1) a de novo RECOVER prospective cohort of children and young adults with and without previous or current infection; and 2) an extant cohort derived from the Adolescent Brain Cognitive Development (ABCD) study ( n =10,000). The de novo cohort incorporates three tiers of data collection: 1) remote baseline assessments (Tier 1, n=6000); 2) longitudinal follow-up for up to 4 years (Tier 2, n=6000); and 3) a subset of participants, primarily the most severely affected by PASC, who will undergo deep phenotyping to explore PASC pathophysiology (Tier 3, n=600). Youth enrolled in the ABCD study participate in Tier 1. The pediatric protocol was developed as a collaborative partnership of investigators, patients, researchers, clinicians, community partners, and federal partners, intentionally promoting inclusivity and diversity. The protocol is adaptive to facilitate responses to emerging science.

Conclusions and Relevance

RECOVER-Pediatrics seeks to characterize the clinical course, underlying mechanisms, and long-term effects of PASC from birth through 25 years old. RECOVER-Pediatrics is designed to elucidate the epidemiology, four-year clinical course, and sociodemographic correlates of pediatric PASC. The data and biosamples will allow examination of mechanistic hypotheses and biomarkers, thus providing insights into potential therapeutic interventions.

Clinical Trialsgov Identifier

Clinical Trial Registration: http://www.clinicaltrials.gov . Unique identifier: NCT05172011.

A Traitlike Dimension of Subjective Memory Concern Over 30 Years Among Adult Male Twins

Importance

Subjective memory concern has long been considered a state-related indicator of impending cognitive decline or dementia. The possibility that subjective memory concern may itself be a heritable trait is largely ignored, yet such an association would substantially confound its use in clinical or research settings.

Objective

To assess the heritability and traitlike dimensions of subjective memory concern and its clinical correlates.

Design, Setting, and Participants

This longitudinal twin cohort study was conducted from 1967 to 2019 among male adults with a mean (SD) age of 37.75 (2.52) years to follow-up at mean ages of 56.15 (2.72), 61.50 (2.43), and 67.35 (2.57) years (hereafter, 38, 56, 62, and 67 years, respectively) in the Vietnam Era Twin Study of Aging. The study included a national community-dwelling sample with health, education, and lifestyle characteristics comparable to a general sample of US men in this age cohort. Participants were monozygotic and dizygotic twins randomly recruited from the Vietnam Era Twin Registry. Data were analyzed from May 2021 to December 2022.

Main Outcomes and Measures

Measures included subjective memory concern at 4 time points; objective memory, depressive symptoms, and anxiety at the last 3 time points; negative emotionality (trait neuroticism) at age 56 years; polygenic risk scores (PRSs) for neuroticism, depression, and Alzheimer disease; APOE genotype; and parental history of dementia. Primary outcomes were heritability and correlations between subjective memory concern and other measures.

Results

The sample included 1555 male adults examined at age 38 years, 520 at age 56 years (due to late introduction of subjective memory concern questions), 1199 at age 62 years, and 1192 at age 67 years. Phenotypically, subjective memory concerns were relatively stable over time. At age 56 years, subjective memory concern had larger correlations with depressive symptoms (r, 0.32; 95% CI, 0.21 to 0.42), anxiety (r, 0.36; 95% CI, 0.18 to 0.51), and neuroticism (r, 0.34; 95% CI, 0.26 to 0.41) than with objective memory (r, -0.24; 95% CI, -0.33 to -0.13). Phenotypic results were similar at ages 62 and 67 years. A best-fitting autoregressive twin model indicated that genetic influences on subjective memory concern accumulated and persisted over time (h2 = 0.26-0.34 from age 38-67 years). At age 56 years, genetic influences for subjective memory concern were moderately correlated with genetic influences for anxiety (r, 0.36; 95% CI, 0.18 to 0.51), negative emotionality (r, 0.51; 95% CI, 0.44-0.57), and depressive symptoms (r, 0.20; 95% CI, 0.10 to 0.29) as well as objective memory (r, -0.22; 95% CI, -0.30 to -0.14). Similar genetic correlations were seen at ages 62 and 67 years. The neuroticism PRS was associated with subjective memory concern at age 38 years (r, 0.10; 95% CI, 0.03. to 0.18) and age 67 years (r, 0.09; 95% CI, 0.01 to 0.16). Subjective memory concern was not associated with any Alzheimer disease risk measures.

Conclusions and Relevance

This cohort study found stable genetic influences underlying subjective memory concern dating back to age 38 years. Subjective memory concern had larger correlations with affect-related measures than with memory-related measures. Improving the utility of subjective memory concern as an indicator of impending cognitive decline and dementia may depend on isolating its statelike component from its traitlike component.

Associations between resting state functional brain connectivity and childhood anhedonia: A reproduction and replication study

BACKGROUND

Previously, a study using a sample of the Adolescent Brain Cognitive Development (ABCD)® study from the earlier 1.0 release found differences in several resting state functional MRI (rsfMRI) brain connectivity measures associated with children reporting anhedonia. Here, we aim to reproduce, replicate, and extend the previous findings using data from the later ABCD study 4.0 release, which includes a significantly larger sample.

METHODS

To reproduce and replicate the previous authors' findings, we analyzed data from the ABCD 1.0 release (n = 2437), from an independent subsample from the newer ABCD 4.0 release (excluding individuals from the 1.0 release) (n = 6456), and from the full ABCD 4.0 release sample (n = 8866). Additionally, we assessed whether using a multiple linear regression approach could improve replicability by controlling for the effects of comorbid psychiatric conditions and sociodemographic covariates.

RESULTS

While the previously reported associations were reproducible, effect sizes for most rsfMRI measures were drastically reduced in replication analyses (including for both t-tests and multiple linear regressions) using the ABCD 4.0 (excluding 1.0) sample. However, 2 new rsfMRI measures (the Auditory vs. Right Putamen and the Retrosplenial-Temporal vs. Right-Thalamus-Proper measures) exhibited replicable associations with anhedonia and stable, albeit small, effect sizes across the ABCD samples, even after accounting for sociodemographic covariates and comorbid psychiatric conditions using a multiple linear regression approach.

CONCLUSION

The most statistically significant associations between anhedonia and rsfMRI connectivity measures found in the ABCD 1.0 sample tended to be non-replicable and inflated. Contrastingly, replicable associations exhibited smaller effects with less statistical significance in the ABCD 1.0 sample. Multiple linear regressions helped assess the specificity of these findings and control the effects of confounding covariates.

Genetic and Cultural Transmission of Alcohol Use Disorder in Swedish Twin Pedigrees

OBJECTIVE

Using Swedish nationwide registry data, we investigated the contribution of genetic and environmental risk factors to the etiology of alcohol use disorder (AUD) by extended twin pedigree modeling.

METHOD

AUD was defined using public inpatient, outpatient, prescription, and criminal records. Three-generational pedigrees were selected for index individuals born between 1980 and 1990, obtained from the national twin and genealogical registers, whose parents were twins. Relatives of the twins included in the pedigrees were their parents, siblings, spouses, and children. Genetic structural equation modeling was applied to the population-based data on AUD, using OpenMx, with age used as a covariate.

RESULTS

Analyses including up to 162,469 individuals in 18,971 pedigrees estimated AUD prevalence at 5%-12% in men and 2%-5% in women. Results indicated substantial heritability (about 50%-60%), of which a portion upwards of 5% was attributable to the consequences of assortative mating. Contributions of shared environmental factors to AUD, which represent a mix of within- and cross-generational effects, appeared to be moderate (about 10%-20%). Unique environment accounted for the remaining variance (about 20%-30%). Sex differences in the magnitude of the variance components suggested higher heritability in men and correspondingly higher shared environmental contributions in women.

CONCLUSIONS

Using objective registry data, we found that AUD is highly heritable. Furthermore, shared environmental factors contributed significantly to the liability of AUD in both men and women.

Alcohol use and cognitive aging in middle-aged men: The Vietnam Era Twin Study of Aging

OBJECTIVE

To determine associations of alcohol use with cognitive aging among middle-aged men.

METHOD

1,608 male twins (mean 57 years at baseline) participated in up to three visits over 12 years, from 2003-2007 to 2016-2019. Participants were classified into six groups based on current and past self-reported alcohol use: lifetime abstainers, former drinkers, very light (1-4 drinks in past 14 days), light (5-14 drinks), moderate (15-28 drinks), and at-risk drinkers (>28 drinks in past 14 days). Linear mixed-effects regressions modeled cognitive trajectories by alcohol group, with time-based models evaluating rate of decline as a function of baseline alcohol use, and age-based models evaluating age-related differences in performance by current alcohol use. Analyses used standardized cognitive domain factor scores and adjusted for sociodemographic and health-related factors.

RESULTS

Performance decreased over time in all domains. Relative to very light drinkers, former drinkers showed worse verbal fluency performance, by -0.21 SD (95% CI -0.35, -0.07), and at-risk drinkers showed faster working memory decline, by 0.14 SD (95% CI 0.02, -0.20) per decade. There was no evidence of protective associations of light/moderate drinking on rate of decline. In age-based models, light drinkers displayed better memory performance at advanced ages than very light drinkers (+0.14 SD; 95% CI 0.02, 0.20 per 10-years older age); likely attributable to residual confounding or reverse association.

CONCLUSIONS

Alcohol consumption showed minimal associations with cognitive aging among middle-aged men. Stronger associations of alcohol with cognitive aging may become apparent at older ages, when cognitive abilities decline more rapidly.

Genetic and Environmental Variation in Continuous Phenotypes in the ABCD Study®

Twin studies yield valuable insights into the sources of variation, covariation and causation in human traits. The ABCD Study® (abcdstudy.org) was designed to take advantage of four universities known for their twin research, neuroimaging, population-based sampling, and expertise in genetic epidemiology so that representative twin studies could be performed. In this paper we use the twin data to: (i) provide initial estimates of heritability for the wide range of phenotypes assessed in the ABCD Study using a consistent direct variance estimation approach, assuring that both data and methodology are sound; and (ii) provide an online resource for researchers that can serve as a reference point for future behavior genetic studies of this publicly available dataset. Data were analyzed from 772 pairs of twins aged 9-10 years at study inception, with zygosity determined using genotypic data, recruited and assessed at four twin hub sites. The online tool provides twin correlations and both standardized and unstandardized estimates of additive genetic, and environmental variation for 14,500 continuously distributed phenotypic features, including: structural and functional neuroimaging, neurocognition, personality, psychopathology, substance use propensity, physical, and environmental trait variables. The estimates were obtained using an unconstrained variance approach, so they can be incorporated directly into meta-analyses without upwardly biasing aggregate estimates. The results indicated broad consistency with prior literature where available and provided novel estimates for phenotypes without prior twin studies or those assessed at different ages. Effects of site, self-identified race/ethnicity, age and sex were statistically controlled. Results from genetic modeling of all 53,172 continuous variables, including 38,672 functional MRI variables, will be accessible via the user-friendly open-access web interface we have established, and will be updated as new data are released from the ABCD Study. This paper provides an overview of the initial results from the twin study embedded within the ABCD Study, an introduction to the primary research domains in the ABCD study and twin methodology, and an evaluation of the initial findings with a focus on data quality and suitability for future behavior genetic studies using the ABCD dataset. The broad introductory material is provided in recognition of the multidisciplinary appeal of the ABCD Study. While this paper focuses on univariate analyses, we emphasize the opportunities for multivariate, developmental and causal analyses, as well as those evaluating heterogeneity by key moderators such as sex, demographic factors and genetic background.

MR-DoC2: Bidirectional Causal Modeling with Instrumental Variables and Data from Relatives

Establishing causality is an essential step towards developing interventions for psychiatric disorders, substance use and many other conditions. While randomized controlled trials (RCTs) are considered the gold standard for causal inference, they are unethical in many scenarios. Mendelian randomization (MR) can be used in such cases, but importantly both RCTs and MR assume unidirectional causality. In this paper, we developed a new model, MRDoC2, that can be used to identify bidirectional causation in the presence of confounding due to both familial and non-familial sources. Our model extends the MRDoC model (Minică et al. in Behav Genet 48:337-349,  https://doi.org/10.1007/s10519-018-9904-4 , 2018), by simultaneously including risk scores for each trait. Furthermore, the power to detect causal effects in MRDoC2 does not require the phenotypes to have different additive genetic or shared environmental sources of variance, as is the case in the direction of causation twin model (Heath et al. in Behav Genet 23:29-50,  https://doi.org/10.1007/BF01067552 , 1993).

A cotwin control study of associations between financial hardship and binge eating phenotypes during COVID-19

BACKGROUND

COVID-19 was associated with significant financial hardship and increased binge eating (BE). However, it is largely unknown whether financial stressors contributed to BE during the pandemic. We used a longitudinal, cotwin control design that controls for genetic/environmental confounds by comparing twins in the same family to examine whether financial hardship during COVID-19 was associated with BE.

METHODS

Female twins (N = 158; M_age  = 22.13) from the Michigan State University Twin Registry rated financial stressors (e.g., inability to afford necessities) daily for 49 consecutive days during COVID-19. We first examined whether financial hardship was associated with BE phenotypes across the full sample. We then examined whether cotwins who differed on financial hardship also differed in BE.

RESULTS

Participants who experienced greater mean financial hardship across the study had significantly greater dimensional BE symptoms, and participants who experienced greater financial hardship on a given day reported significantly more emotional eating that day. These results were replicated in cotwin control analyses. Twins who experienced more financial hardship than their cotwin across the study reported greater dimensional BE symptoms than their cotwin, and participants who experienced more financial hardship than their cotwin on a given day reported greater emotional eating that day. Results were identical when restricting analyses to monozygotic twins, suggesting associations were not due to genetic confounds.

CONCLUSIONS

Results suggest that BE-related symptoms may be elevated in women who experienced financial hardship during COVID-19 independent of potential genetic/environmental confounds. However, additional research in larger samples is needed.

PUBLIC SIGNIFICANCE

Little is known regarding how financial difficulties during the COVID-19 pandemic may have contributed to increased binge eating (BE). We found preliminary evidence that financial hardship during COVID-19 may be associated with greater rates of BE-related symptoms even when comparing twins from the same family. While additional research is needed, results suggest that people who experienced financial hardship during COVID-19 may be at increased risk for BE.

Association of neurofilament light chain with renal function: mechanisms and clinical implications

BACKGROUND

Blood-based neurofilament light chain (NfL) is a promising biomarker of neurodegeneration across multiple neurodegenerative diseases. However, blood-based NfL is highly associated with renal function in older adults, which leads to the concern that blood-based NfL levels may be influenced by renal function, rather than neurodegeneration alone. Despite growing interest in using blood-based NfL as a biomarker of neurodegeneration in research and clinical practices, whether renal function should always be accounted for in these settings remains unclear. Moreover, the mechanisms underlying this association between blood-based measures of NfL and renal function remain elusive. In this study, we first evaluated the effect of renal function on the associations of plasma NfL with other measures of neurodegeneration. We then examined the extent of genetic and environmental contributions to the association between plasma NfL and renal function.

METHODS

In a sample of 393 adults (mean age=75.22 years, range=54-90), we examined the associations of plasma NfL with cerebrospinal fluid (CSF) NfL and brain volumetric measures before and after adjusting for levels of serum creatinine (an index of renal function). In an independent sample of 969 men (mean age=67.57 years, range=61-73) that include monozygotic and dizygotic twin pairs, we replicated the same analyses and leveraged biometrical twin modeling to examine the genetic and environmental influences on the plasma NfL and creatinine association.

RESULTS

Plasma NfL's associations with cerebrospinal fluid NfL and brain volumetric measures did not meaningfully change after adjusting for creatinine levels. Both plasma NfL and creatinine were significantly heritable (h2=0.54 and 0.60, respectively). Their phenotypic correlation (r=0.38) was moderately explained by shared genetic influences (genetic correlation=0.46) and unique environmental influences (unique environmental correlation=0.27).

CONCLUSIONS

Adjusting for renal function is unnecessary when assessing associations between plasma NfL and other measures of neurodegeneration but is necessary if plasma NfL is compared to a cutoff for classifying neurodegeneration-positive versus neurodegeneration-negative individuals. Blood-based measures of NfL and renal function are heritable and share common genetic influences.

A daily diary study of emotion regulation as a moderator of negative affect-binge eating associations

BACKGROUND

While negative affect (NA) typically increases risk for binge eating, the ultimate impact of NA may depend on a person's ability to regulate their emotions. In this daily, longitudinal study, we examined whether emotion regulation (ER) modified the strength of NA-dysregulated eating associations.

METHODS

Women (N = 311) from the Michigan State University Twin Registry first reported dimensional binge eating symptoms and broad ER difficulties (e.g., limited emotional awareness, difficulty controlling emotional impulses). Participants then rated use of adaptive (cognitive reappraisal, social sharing, situation modification, and acceptance) and maladaptive (rumination, expressive suppression, and self-criticism) ER strategies, emotional eating (EE), objective binge eating (OBE), and NA once daily for 49 consecutive days.

RESULTS

There were several main effects of ER on binge-eating pathology in both between-person (i.e., comparing women who differed on average) and within-person (i.e., examining fluctuations in variables day-to-day) analyses. Between-person, greater broad ER difficulties, greater maladaptive strategy use, and lower adaptive strategy use were all associated with greater binge-eating pathology. Within-person, greater maladaptive strategy use was associated with greater odds of OBE on that day and on the following day. However, neither broad ER difficulties nor use of specific strategies moderated associations between NA and dysregulated eating in between- or within-person analyses.

CONCLUSIONS

While ER is independently associated with risk for dysregulated eating, it may not fully mitigate the impact of NA. Additional strategies (e.g., decreasing environmental stressors and increasing social support) may be needed to minimize NA and its impact on dysregulated eating.

PUBLIC SIGNIFICANCE

Negative affect (NA; e.g., sadness, guilt) increases dysregulated eating risk. Because NA is sometimes unavoidable, we examined whether emotion regulation (ER; i.e., how a person responds to their emotions) might impact whether NA leads to dysregulated eating. Although more effective ER was associated with less dysregulated eating overall, ER did not impact the association between NA and dysregulated eating. Other approaches may therefore be needed to mitigate NA-dysregulated eating associations.

Associations between depression and cardiometabolic health: A 27-year longitudinal study

BACKGROUND

Clarifying the relationship between depression symptoms and cardiometabolic and related health could clarify risk factors and treatment targets. The objective of this study was to assess whether depression symptoms in midlife are associated with the subsequent onset of cardiometabolic health problems.

METHODS

The study sample comprised 787 male twin veterans with polygenic risk score data who participated in the Harvard Twin Study of Substance Abuse ('baseline') and the longitudinal Vietnam Era Twin Study of Aging ('follow-up'). Depression symptoms were assessed at baseline [mean age 41.42 years (s.d. = 2.34)] using the Diagnostic Interview Schedule, Version III, Revised. The onset of eight cardiometabolic conditions (atrial fibrillation, diabetes, erectile dysfunction, hypercholesterolemia, hypertension, myocardial infarction, sleep apnea, and stroke) was assessed via self-reported doctor diagnosis at follow-up [mean age 67.59 years (s.d. = 2.41)].

RESULTS

Total depression symptoms were longitudinally associated with incident diabetes (OR 1.29, 95% CI 1.07-1.57), erectile dysfunction (OR 1.32, 95% CI 1.10-1.59), hypercholesterolemia (OR 1.26, 95% CI 1.04-1.53), and sleep apnea (OR 1.40, 95% CI 1.13-1.74) over 27 years after controlling for age, alcohol consumption, smoking, body mass index, C-reactive protein, and polygenic risk for specific health conditions. In sensitivity analyses that excluded somatic depression symptoms, only the association with sleep apnea remained significant (OR 1.32, 95% CI 1.09-1.60).

CONCLUSIONS

A history of depression symptoms by early midlife is associated with an elevated risk for subsequent development of several self-reported health conditions. When isolated, non-somatic depression symptoms are associated with incident self-reported sleep apnea. Depression symptom history may be a predictor or marker of cardiometabolic risk over decades.

The Impact of Genes and Environment on Brain Ageing in Males Aged 51 to 72 Years

Magnetic resonance imaging data are being used in statistical models to predicted brain ageing (PBA) and as biomarkers for neurodegenerative diseases such as Alzheimer's Disease. Despite their increasing application, the genetic and environmental etiology of global PBA indices is unknown. Likewise, the degree to which genetic influences in PBA are longitudinally stable and how PBA changes over time are also unknown. We analyzed data from 734 men from the Vietnam Era Twin Study of Aging with repeated MRI assessments between the ages 51-72 years. Biometrical genetic analyses "twin models" revealed significant and highly correlated estimates of additive genetic heritability ranging from 59 to 75%. Multivariate longitudinal modeling revealed that covariation between PBA at different timepoints could be explained by a single latent factor with 73% heritability. Our results suggest that genetic influences on PBA are detectable in midlife or earlier, are longitudinally very stable, and are largely explained by common genetic influences.

Long-term associations of cigarette smoking in early mid-life with predicted brain aging from mid- to late life

BACKGROUND AND AIMS

Smoking is associated with increased risk for brain aging/atrophy and dementia. Few studies have examined early associations with brain aging. This study aimed to measure whether adult men with a history of heavier smoking in early mid-life would have older than predicted brain age 16-28 years later.

DESIGN

Prospective cohort observational study, utilizing smoking pack years data from average age 40 (early mid-life) predicting predicted brain age difference scores (PBAD) at average ages 56, 62 (later mid-life) and 68 years (early old age). Early mid-life alcohol use was also evaluated.

SETTING

Population-based United States sample.

PARTICIPANTS/CASES

Participants were male twins of predominantly European ancestry who served in the United States military between 1965 and 1975. Structural magnetic resonance imaging (MRI) began at average age 56. Subsequent study waves included most baseline participants; attrition replacement subjects were added at later waves.

MEASUREMENTS

Self-reported smoking information was used to calculate pack years smoked at ages 40, 56, 62, and 68. MRIs were processed with the Brain-Age Regression Analysis and Computation Utility software (BARACUS) program to create PBAD scores (chronological age-predicted brain age) acquired at average ages 56 (n = 493; 2002-08), 62 (n = 408; 2009-14) and 68 (n = 499; 2016-19).

FINDINGS

In structural equation modeling, age 40 pack years predicted more advanced age 56 PBAD [β = -0.144, P = 0.012, 95% confidence interval (CI) = -0.257, -0.032]. Age 40 pack years did not additionally predict PBAD at later ages. Age 40 alcohol consumption, but not a smoking × alcohol interaction, predicted more advanced PBAD at age 56 (β = -0.166, P = 0.001, 95% CI = -0.261, -0.070) with additional influences at age 62 (β = -0.115, P = 0.005, 95% CI = -0.195, -0.036). Age 40 alcohol did not predict age 68 PBAD. Within-twin-pair analyses suggested some genetic mechanism partially underlying effects of alcohol, but not smoking, on PBAD.

CONCLUSIONS

Heavier smoking and alcohol consumption by age 40 appears to predict advanced brain aging by age 56 in men.

Paradoxical cognitive trajectories in men from earlier to later adulthood

Because longitudinal studies of aging typically lack cognitive data from earlier ages, it is unclear how general cognitive ability (GCA) changes throughout the life course. In 1173 Vietnam Era Twin Study of Aging (VETSA) participants, we assessed young adult GCA at average age 20 and current GCA at 3 VETSA assessments beginning at average age 56. The same GCA index was used throughout. Higher young adult GCA and better GCA maintenance were associated with stronger specific cognitive abilities from age 51 to 73. Given equivalent GCA at age 56, individuals who had higher age 20 GCA outperformed those whose GCA remained stable in terms of memory, executive function, and working memory abilities from age 51 to 73. Thus, paradoxically, despite poorer maintenance of GCA, high young adult GCA still conferred benefits. Advanced predicted brain age and the combination of elevated vascular burden and APOE-ε4 status were associated with poorer maintenance of GCA. These findings highlight the importance of distinguishing between peak and current GCA for greater understanding of cognitive aging.

Associations Between Traumatic Stress, Brain Volumes and Post-traumatic Stress Disorder Symptoms in Children: Data from the ABCD Study

Reduced volumes in brain regions of interest (ROIs), primarily from adult samples, are associated with posttraumatic stress disorder (PTSD). We extended this work to children using data from the Adolescent Brain Cognitive Development (ABCD) Study® (N = 11,848; M_age = 9.92). Structural equation modeling and an elastic-net (EN) machine-learning approach were used to identify potential effects of traumatic events (TEs) on PTSD symptoms (PTSDsx) directly, and indirectly via the volumes 300 subcortical and cortical ROIs. We then estimated the genetic and environmental variation in the phenotypes. TEs were directly associated with PTSDsx (r = 0.92) in children, but their indirect effects (r < 0.0004)-via the volumes of EN-identified subcortical and cortical ROIs-were negligible at this age. Additive genetic factors explained a modest proportion of the variance in TEs (23.4%) and PTSDsx (21.3%), and accounted for most of the variance of EN-identified volumes of four of the five subcortical (52.4-61.8%) three of the nine cortical ROIs (46.4-53.3%) and cerebral white matter in the left hemisphere (57.4%). Environmental factors explained most of the variance in TEs (C = 61.6%, E = 15.1%), PTSDsx (residual-C = 18.4%, residual-E = 21.8%), right lateral ventricle (C = 15.2%, E = 43.1%) and six of the nine EN-identified cortical ROIs (C = 4.0-13.6%, E = 56.7-74.8%). There is negligible evidence that the volumes of brain ROIs are associated with the indirect effects of TEs on PTSDsx at this age. Overall, environmental factors accounted for more of the variation in TEs and PTSDsx. Whereas additive genetic factors accounted for most of the variability in the volumes of a minority of cortical and in most of subcortical ROIs.

Passive Sensing of Preteens' Smartphone Use: An Adolescent Brain Cognitive Development (ABCD) Cohort Substudy

BACKGROUND

Concerns abound regarding childhood smartphone use, but studies to date have largely relied on self-reported screen use. Self-reporting of screen use is known to be misreported by pediatric samples and their parents, limiting the accurate determination of the impact of screen use on social, emotional, and cognitive development. Thus, a more passive, objective measurement of smartphone screen use among children is needed.

OBJECTIVE

This study aims to passively sense smartphone screen use by time and types of apps used in a pilot sample of children and to assess the feasibility of passive sensing in a larger longitudinal sample.

METHODS

The Adolescent Brain Cognitive Development (ABCD) study used passive, objective phone app methods for assessing smartphone screen use over 4 weeks in 2019-2020 in a subsample of 67 participants (aged 11-12 years; 31/67, 46% female; 23/67, 34% White). Children and their parents both reported average smartphone screen use before and after the study period, and they completed a questionnaire regarding the acceptability of the study protocol. Descriptive statistics for smartphone screen use, app use, and protocol feasibility and acceptability were reviewed. Analyses of variance were run to assess differences in categorical app use by demographics. Self-report and parent report were correlated with passive sensing data.

RESULTS

Self-report of smartphone screen use was partly consistent with objective measurement (r=0.49), although objective data indicated that children used their phones more than they reported. Passive sensing revealed the most common types of apps used were for streaming (mean 1 hour 57 minutes per day, SD 1 hour 32 minutes), communication (mean 48 minutes per day, SD 1 hour 17 minutes), gaming (mean 41 minutes per day, SD 41 minutes), and social media (mean 36 minutes per day, SD 1 hour 7 minutes). Passive sensing of smartphone screen use was generally acceptable to children (43/62, 69%) and parents (53/62, 85%).

CONCLUSIONS

The results of passive, objective sensing suggest that children use their phones more than they self-report. Therefore, use of more robust methods for objective data collection is necessary and feasible in pediatric samples. These data may then more accurately reflect the impact of smartphone screen use on behavioral and emotional functioning. Accordingly, the ABCD study is implementing a passive sensing protocol in the full ABCD cohort. Taken together, passive assessment with a phone app provided objective, low-burden, novel, informative data about preteen smartphone screen use.

Lifestyle and the aging brain: interactive effects of modifiable lifestyle behaviors and cognitive ability in men from midlife to old age

We examined the influence of lifestyle on brain aging after nearly 30 years, and tested the hypothesis that young adult general cognitive ability (GCA) would moderate these effects. In the community-dwelling Vietnam Era Twin Study of Aging (VETSA), 431 largely non-Hispanic white men completed a test of GCA at mean age 20. We created a modifiable lifestyle behavior composite from data collected at mean age 40. During VETSA, MRI-based measures at mean age 68 included predicted brain age difference (PBAD), Alzheimer's disease (AD) brain signature, and abnormal white matter scores. There were significant main effects of young adult GCA and lifestyle on PBAD and the AD signature (ps ≤ 0.012), and a GCA-by-lifestyle interaction on both (ps ≤ 0.006). Regardless of GCA level, having more favorable lifestyle behaviors predicted less advanced brain age and less AD-like brain aging. Unfavorable lifestyles predicted advanced brain aging in those with lower age 20 GCA, but did not affect brain aging in those with higher age 20 GCA. Targeting early lifestyle modification may promote dementia risk reduction, especially among lower reserve individuals.

Vaginal microbiome Lactobacillus crispatus is heritable among European American women

The diversity and dominant bacterial taxa in the vagina are reported to be influenced by multiple intrinsic and extrinsic factors, including but not limited to pregnancy, contraceptive use, pathogenic states, socioeconomic status, and ancestry. However, the extent to which host genetic factors influence variation in the vaginal microbiota is unclear. We used a biometrical genetic approach to determine whether host genetic factors contribute to inter-individual differences in taxa from a sample of 332 twins who self-identified as being of African (44 pairs) or European ancestry (122 pairs). Lactobacillus crispatus, a major determinant of vaginal health, was identified as heritable among European American women (narrow-sense heritability = 34.7%, P-value = 0.018). Heritability of L. crispatus is consistent with the reduced prevalence of adverse reproductive disorders, including bacterial vaginosis and preterm birth, among women of European ancestry.

Wright et al. apply biometric genetic approach to identify the extent to which host genetic factors influence species-level variation in the vaginal microbiota. Their study suggests that Lactobacillus crispatus, a major determinant of vaginal health, is heritable among European American women, consistent with the reduced prevalence of adverse reproductive disorders in them.

Baseline brain function in the preadolescents of the ABCD Study

The Adolescent Brain Cognitive Development (ABCD) Study® is a 10-year longitudinal study of children recruited at ages 9 and 10. A battery of neuroimaging tasks are administered biennially to track neurodevelopment and identify individual differences in brain function. This study reports activation patterns from functional MRI (fMRI) tasks completed at baseline, which were designed to measure cognitive impulse control with a stop signal task (SST; N = 5,547), reward anticipation and receipt with a monetary incentive delay (MID) task (N = 6,657) and working memory and emotion reactivity with an emotional N-back (EN-back) task (N = 6,009). Further, we report the spatial reproducibility of activation patterns by assessing between-group vertex/voxelwise correlations of blood oxygen level-dependent (BOLD) activation. Analyses reveal robust brain activations that are consistent with the published literature, vary across fMRI tasks/contrasts and slightly correlate with individual behavioral performance on the tasks. These results establish the preadolescent brain function baseline, guide interpretation of cross-sectional analyses and will enable the investigation of longitudinal changes during adolescent development.

12-year prediction of mild cognitive impairment aided by Alzheimer's brain signatures at mean age 56

Neuroimaging signatures based on composite scores of cortical thickness and hippocampal volume predict progression from mild cognitive impairment to Alzheimer's disease. However, little is known about the ability of these signatures among cognitively normal adults to predict progression to mild cognitive impairment. Towards that end, a signature sensitive to microstructural changes that may predate macrostructural atrophy should be useful. We hypothesized that: (i) a validated MRI-derived Alzheimer's disease signature based on cortical thickness and hippocampal volume in cognitively normal middle-aged adults would predict progression to mild cognitive impairment; and (ii) a novel grey matter mean diffusivity signature would be a better predictor than the thickness/volume signature. This cohort study was part of the Vietnam Era Twin Study of Aging. Concurrent analyses compared cognitively normal and mild cognitive impairment groups at each of three study waves (ns = 246-367). Predictive analyses included 169 cognitively normal men at baseline (age = 56.1, range = 51-60). Our previously published thickness/volume signature derived from independent data, a novel mean diffusivity signature using the same regions and weights as the thickness/volume signature, age, and an Alzheimer's disease polygenic risk score were used to predict incident mild cognitive impairment an average of 12 years after baseline (follow-up age = 67.2, range = 61-71). Additional analyses adjusted for predicted brain age difference scores (chronological age minus predicted brain age) to determine if signatures were Alzheimer-related and not simply ageing-related. In concurrent analyses, individuals with mild cognitive impairment had higher (worse) mean diffusivity signature scores than cognitively normal participants, but thickness/volume signature scores did not differ between groups. In predictive analyses, age and polygenic risk score yielded an area under the curve of 0.74 (sensitivity = 80.00%; specificity = 65.10%). Prediction was significantly improved with addition of the mean diffusivity signature (area under the curve = 0.83; sensitivity = 85.00%; specificity = 77.85%; P = 0.007), but not with addition of the thickness/volume signature. A model including both signatures did not improve prediction over a model with only the mean diffusivity signature. Results held up after adjusting for predicted brain age difference scores. The novel mean diffusivity signature was limited by being yoked to the thickness/volume signature weightings. An independently derived mean diffusivity signature may thus provide even stronger prediction. The young age of the sample at baseline is particularly notable. Given that the brain signatures were examined when participants were only in their 50 s, our results suggest a promising step towards improving very early identification of Alzheimer's disease risk and the potential value of mean diffusivity and/or multimodal brain signatures.

Genetically Informed Regression Analysis: Application to Aggression Prediction by Inattention and Hyperactivity in Children and Adults

We present a procedure to simultaneously fit a genetic covariance structure model and a regression model to multivariate data from mono- and dizygotic twin pairs to test for the prediction of a dependent trait by multiple correlated predictors. We applied the model to aggressive behavior as an outcome trait and investigated the prediction of aggression from inattention (InA) and hyperactivity (HA) in two age groups. Predictions were examined in twins with an average age of 10 years (11,345 pairs), and in adult twins with an average age of 30 years (7433 pairs). All phenotypes were assessed by the same, but age-appropriate, instruments in children and adults. Because of the different genetic architecture of aggression, InA and HA, a model was fitted to these data that specified additive and non-additive genetic factors (A and D) plus common and unique environmental (C and E) influences. Given appropriate identifying constraints, this ADCE model is identified in trivariate data. We obtained different results for the prediction of aggression in children, where HA was the more important predictor, and in adults, where InA was the more important predictor. In children, about 36% of the total aggression variance was explained by the genetic and environmental components of HA and InA. Most of this was explained by the genetic components of HA and InA, i.e., 29.7%, with 22.6% due to the genetic component of HA. In adults, about 21% of the aggression variance was explained. Most was this was again explained by the genetic components of InA and HA (16.2%), with 8.6% due to the genetic component of InA.

Combining Structural-Equation Modeling with Genomic-Relatedness-Matrix Restricted Maximum Likelihood in OpenMx

There is a long history of fitting biometrical structural-equation models (SEMs) in the pregenomic behavioral-genetics literature of twin, family, and adoption studies. Recently, a method has emerged for estimating biometrical variance-covariance components based not upon the expected degree of genetic resemblance among relatives, but upon the observed degree of genetic resemblance among unrelated individuals for whom genome-wide genotypes are available-genomic-relatedness-matrix restricted maximum-likelihood (GREML). However, most existing GREML software is concerned with quickly and efficiently estimating heritability coefficients, genetic correlations, and so on, rather than with allowing the user to fit SEMs to multitrait samples of genotyped participants. We therefore introduce a feature in the OpenMx package, "mxGREML", designed to fit the biometrical SEMs from the pregenomic era in present-day genomic study designs. We explain the additional functionality this new feature has brought to OpenMx, and how the new functionality works. We provide an illustrative example of its use. We discuss the feature's current limitations, and our plans for its further development.

Best Practices for Binary and Ordinal Data Analyses

The measurement of many human traits, states, and disorders begins with a set of items on a questionnaire. The response format for these questions is often simply binary (e.g., yes/no) or ordered (e.g., high, medium or low). During data analysis, these items are frequently summed or used to estimate factor scores. In clinical applications, such assessments are often non-normally distributed in the general population because many respondents are unaffected, and therefore asymptomatic. As a result, in many cases these measures violate the statistical assumptions required for subsequent analyses. To reduce the influence of the non-normality and quasi-continuous assessment, variables are frequently recoded into binary (affected-unaffected) or ordinal (mild-moderate-severe) diagnoses. Ordinal data therefore present challenges at multiple levels of analysis. Categorizing continuous variables into ordered categories typically results in a loss of statistical power, which represents an incentive to the data analyst to assume that the data are normally distributed, even when they are not. Despite prior zeitgeists suggesting that, e.g., variables with more than 10 ordered categories may be regarded as continuous and analyzed as if they were, we show via simulation studies that this is not generally the case. In particular, using Pearson product-moment correlations instead of maximum likelihood estimates of polychoric correlations biases the estimated correlations towards zero. This bias is especially severe when a plurality of the observations fall into a single observed category, such as a score of zero. By contrast, estimating the ordinal correlation by maximum likelihood yields no estimation bias, although standard errors are (appropriately) larger. We also illustrate how odds ratios depend critically on the proportion or prevalence of affected individuals in the population, and therefore are sub-optimal for studies where comparisons of association metrics are needed. Finally, we extend these analyses to the classical twin model and demonstrate that treating binary data as continuous will underestimate genetic and common environmental variance components, and overestimate unique environment (residual) variance. These biases increase as prevalence declines. While modeling ordinal data appropriately may be more computationally intensive and time consuming, failing to do so will likely yield biased correlations and biased parameter estimates from modeling them.

Clarifying the Genetic Influences on Nicotine Dependence and Quantity of Use in Cigarette Smokers

Nicotine dependence and smoking quantity are both robustly associated with the CHRNA5-A3-B4 gene cluster in the 15q25 region, and SNP rs16969968 in particular. The purpose of this paper is to use structural equation modeling techniques (SEM) to disentangle the complex pattern of relationships between rs16969968, nicotine quantity (as measured by the number of cigarettes an individual smokes per day; CPD) and nicotine dependence (as measured by the Fagerström Test for Nicotine Dependence; FTND). CPD is an indicator, but also a potential cause, of FTND, complicating the interpretation of associations between these constructs and requires a more detailed investigation than standard GWAS or general linear regression models can provide. FTND items and genotypes were collected in four samples, with a combined sample size of 5,373 respondents. A mega-analysis was conducted using a multiple group SEM approach to test competing hypotheses regarding the relationships between the SNP rs16969968, FTND and CPD. In the best fitting model, the FTND items loaded onto two correlated factors. The first, labeled "maintenance," assesses the motivation to maintain constant levels of nicotine through out the day. The second was labeled "urgency" as its items concern the urgency to restore nicotine levels after abstinence. We focus our attention on the "maintenance" factor, of which CPD was an indicator. The best fitting model included a negative feedback loop between the Maintenance factor and CPD. Accordingly, the motivation to maintain higher levels of nicotine increased the quantity of nicotine consumed, which subsequently decreases the maintenance motivation. The fact that the Maintenance-CPD feedback model fits the data best implies that there are at least two biological pathways that lead from rs16969968 to smoking behaviors. The model is consistent with a supply and demand system, which allows individuals to achieve a homeostatic equilibrium for their nicotine concentration.

Periventricular and deep abnormal white matter differ in associations with cognitive performance at midlife

Objective: Abnormal white matter (AWM) on magnetic resonance imaging is associated with cognitive performance in older adults. We explored cognitive associations with AWM during late-midlife. Method: Participants were community-dwelling men (n = 242; M = 61.90 years; range = 56-66). Linear-mixed effects regression models examined associations of total, periventricular, and deep AWM with cognitive performance, controlling for multiple comparisons. Models considering specific cognitive domains controlled for current general cognitive ability (GCA). We hypothesized that total AWM would be associated with worse processing speed, executive function, and current GCA; deep AWM would correlate with GCA and periventricular AWM would relate to specific cognitive abilities. We also assessed the potential influence of cognitive reserve by examining a moderation effect of early life (mean age of 20) cognition. Results: Greater total and deep AWM were associated with poorer current GCA. Periventricular AWM was associated with worse executive function, working memory, and episodic memory. When periventricular and deep AWM were modeled simultaneously, both retained their respective significant associations with cognitive performance. Cognitive reserve did not moderate associations. Conclusions: Our findings suggest that AWM contributes to poorer cognitive function in late-midlife. Examining only total AWM may obscure the potential differential impact of regional AWM. Separating total AWM into subtypes while controlling for current GCA revealed a dissociation in relationships with cognitive performance; deep AWM was associated with nonspecific cognitive ability whereas periventricular AWM was associated with specific frontal-related abilities and memory. Management of vascular or other risk factors that may increase the risk of AWM should begin during or before early late-midlife. (PsycInfo Database Record (c) 2021 APA, all rights reserved).

GW-SEM 2.0: Efficient, Flexible, and Accessible Multivariate GWAS

Most genome-wide association study (GWAS) analyses test the association between single-nucleotide polymorphisms (SNPs) and a single trait or outcome. While valuable second-step analyses of these associations (e.g., calculating genetic correlations between traits) are common, single-step multivariate analyses of GWAS data are rarely performed. This is unfortunate because multivariate analyses can reveal information which is irrevocably obscured in multi-step analysis. One simple example is the distinction between variance common to a set of measures, and variance specific to each. Neither GWAS of sum- or factor-scores, nor GWAS of the individual measures will deliver a clean picture of loci associated with each measure's specific variance. While multivariate GWAS opens up a broad new landscape of feasible and informative analyses, its adoption has been slow, likely due to the heavy computational demands and difficulties specifying models it requires. Here we describe GW-SEM 2.0, which is designed to simplify model specification and overcome the inherent computational challenges associated with multivariate GWAS. In addition, GW-SEM 2.0 allows users to accurately model ordinal items, which are common in behavioral and psychological research, within a GWAS context. This new release enhances computational efficiency, allows users to select the fit function that is appropriate for their analyses, expands compatibility with standard genomic data formats, and outputs results for seamless reading into other standard post-GWAS processing software. To demonstrate GW-SEM's utility, we conducted (1) a series of GWAS using three substance use frequency items from data in the UK Biobank, (2) a timing study for several predefined GWAS functions, and (3) a Type I Error rate study. Our multivariate GWAS analyses emphasize the utility of GW-SEM for identifying novel patterns of associations that vary considerably between genomic loci for specific substances, highlighting the importance of differentiating between substance-specific use behaviors and polysubstance use. The timing studies demonstrate that the analyses take a reasonable amount of time and show the cost of including additional items. The Type I Error rate study demonstrates that hypothesis tests for genetic associations with latent variable models follow the hypothesized uniform distribution. Taken together, we suggest that GW-SEM may provide substantially deeper insights into the underlying genomic architecture for multivariate behavioral and psychological systems than is currently possible with standard GWAS methods. The current release of GW-SEM 2.0 is available on CRAN (stable release) and GitHub (beta release), and tutorials are available on our github wiki ( https://jpritikin.github.io/gwsem/ ).

MRI-assessed locus coeruleus integrity is heritable and associated with multiple cognitive domains, mild cognitive impairment, and daytime dysfunction

Introduction

The locus coeruleus (LC) undergoes extensive neurodegeneration in early Alzheimer's disease (AD). The LC is implicated in regulating the sleep–wake cycle, modulating cognitive function, and AD progression.

Methods

Participants were 481 men (ages 62 to 71.7) from the Vietnam Era Twin Study of Aging. LC structural integrity was indexed by neuromelanin‐sensitive magnetic resonance imaging (MRI) contrast‐to‐noise ratio (LC_CNR). We examined LC_CNR, cognition, amnestic mild cognitive impairment (aMCI), and daytime dysfunction.

Results

Heritability of LC_CNR was .48. Participants with aMCI showed greater daytime dysfunction. Lower LC_CNR was associated with poorer episodic memory, general verbal fluency, semantic fluency, and processing speed, as well as increased odds of aMCI and greater daytime dysfunction.

Discussion

Reduced LC integrity is associated with widespread differences across cognitive domains, daytime sleep‐related dysfunction, and risk for aMCI. These findings in late‐middle‐aged adults highlight the potential of MRI‐based measures of LC integrity in early identification of AD risk.

The Augmented Classical Twin Design: Incorporating Genome-Wide Identity by Descent Sharing Into Twin Studies in Order to Model Violations of the Equal Environments Assumption

The Classical Twin Method (CTM) compares the similarity of monozygotic (MZ) twins with that of dizygotic (DZ) twins to make inferences about the relative importance of genes and environment in the etiology of individual differences. The design has been applied to thousands of traits across the biomedical, behavioral and social sciences and is arguably the most widely used natural experiment known to science. The fundamental assumption of the CTM is that trait relevant environmental covariation within MZ pairs is the same as that found within DZ pairs, so that zygosity differences in within-pair variance must be due to genetic factors uncontaminated by the environment. This equal environments assumption (EEA) has been, and still is hotly contested, and has been mentioned as a possible contributing factor to the missing heritability conundrum. In this manuscript, we introduce a new model for testing the EEA, which we call the Augmented Classical Twin Design which uses identity by descent (IBD) sharing between DZ twin pairs to estimate separate environmental variance components for MZ and DZ twin pairs, and provides a test of whether these are equal. We show through simulation that given large samples of DZ twin pairs, the model provides unbiased estimates of variance components and valid tests of the EEA under strong assumptions (e.g. no epistatic variance, IBD sharing in DZ twins estimated accurately etc.) which may not hold in reality. Sample sizes in excess of 50,000 DZ twin pairs with genome-wide genetic data are likely to be required in order to detect substantial violations of the EEA with moderate power. Consequently, we recommend that the Augmented Classical Twin Design only be applied to datasets with very large numbers of DZ twin pairs (> 50,000 DZ twin pairs), and given the strong assumptions relating to the absence of epistatic variance, appropriate caution be exercised regarding interpretation of the results.

Modeling Etiology of Smoking During Pregnancy in Swedish Twins, Full-, and Half-Siblings, Reared Together and Apart

INTRODUCTION

Using Swedish nationwide registry data, we investigated the contribution of genetic and environmental risk factors to the etiology of smoking status across stages of pregnancy with increasing degrees of social and psychological pressure to reduce or quit smoking, by twin and sibling modeling.

AIMS AND METHODS

Smoking status was available before, and during early and late pregnancy from the Medical Birth Register. Twin, full-, and half-sibling pairs, both reared together and apart, born between 1960 and 1990 were obtained from national twin and genealogical registers. Genetic structural equation modeling in OpenMx was applied to the population-based data to estimate shared genetic and/or environmental covariance across stages of pregnancy, accounting for maternal birth cohort and age at pregnancy.

RESULTS

Analyses of paired data on 258 749 individuals suggested that risk factors for smoking status changed across stages of pregnancy. Results predicted substantial heritability (60-70%) and moderate contributions of shared environmental factors (10-15%) for smoking status. Whilst the same shared environmental factors were amplified from before pregnancy to late pregnancy, new primarily unique environmental factors explained ~10% of the variance during early pregnancy which was carried forward to late pregnancy.

CONCLUSIONS

Using registry data on women across pregnancy, we replicated that smoking status is highly heritable. Furthermore, we found support for increased impact of shared environmental factors during pregnancy of factors already present prior to pregnancy, and an independent set of mostly new unique environmental factors that may be triggered by increased social pressure to reduce or quit smoking during pregnancy.

IMPLICATIONS

As new factors partially explain smoking status during pregnancy and the effects of familial factors increase across pregnancy, efforts to prevent or reduce smoking during pregnancy should receive continued attention, with a focus on both the individual and the family unit.

Using Multimodel Inference/Model Averaging to Model Causes of Covariation Between Variables in Twins

OBJECTIVE

To explore and apply multimodel inference to test the relative contributions of latent genetic, environmental and direct causal factors to the covariation between two variables with data from the classical twin design by estimating model-averaged parameters.

METHODS

Behavior genetics is concerned with understanding the causes of variation in phenotypes and the causes of covariation between two or more phenotypes. Two variables may correlate as a result of genetic, shared environmental or unique environmental factors contributing to variation in both variables. Two variables may also correlate because one or both directly cause variation in the other. Furthermore, covariation may result from any combination of these sources, leading to 25 different identified structural equation models. OpenMx was used to fit all these models to account for covariation between two variables collected in twins. Multimodel inference and model averaging were used to summarize the key sources of covariation, and estimate the magnitude of these causes of covariance. Extensions of these models to test heterogeneity by sex are discussed.

RESULTS

We illustrate the application of multimodel inference by fitting a comprehensive set of bivariate models to twin data from the Virginia Twin Study of Psychiatric and Substance Use Disorders. Analyses of body mass index and tobacco consumption data show sufficient power to reject distinct models, and to estimate the contribution of each of the five potential sources of covariation, irrespective of selecting the best fitting model. Discrimination between models on sample size, type of variable (continuous versus binary or ordinal measures) and the effect size of sources of variance and covariance.

CONCLUSIONS

We introduce multimodel inference and model averaging approaches to the behavior genetics community, in the context of testing models for the causes of covariation between traits in term of genetic, environmental and causal explanations.

Expanding the genetic architecture of nicotine dependence and its shared genetics with multiple traits

Cigarette smoking is the leading cause of preventable morbidity and mortality. Genetic variation contributes to initiation, regular smoking, nicotine dependence, and cessation. We present a Fagerström Test for Nicotine Dependence (FTND)-based genome-wide association study in 58,000 European or African ancestry smokers. We observe five genome-wide significant loci, including previously unreported loci MAGI2/GNAI1 (rs2714700) and TENM2 (rs1862416), and extend loci reported for other smoking traits to nicotine dependence. Using the heaviness of smoking index from UK Biobank (N = 33,791), rs2714700 is consistently associated; rs1862416 is not associated, likely reflecting nicotine dependence features not captured by the heaviness of smoking index. Both variants influence nearby gene expression (rs2714700/MAGI2-AS3 in hippocampus; rs1862416/TENM2 in lung), and expression of genes spanning nicotine dependence-associated variants is enriched in cerebellum. Nicotine dependence (SNP-based heritability = 8.6%) is genetically correlated with 18 other smoking traits (rg = 0.40-1.09) and co-morbidities. Our results highlight nicotine dependence-specific loci, emphasizing the FTND as a composite phenotype that expands genetic knowledge of smoking.

Increased rates of eating disorders and their symptoms in women with major depressive disorder and anxiety disorders

BACKGROUND

Individuals with eating disorders (EDs) have increased rates of major depressive disorder (MDD) and anxiety disorders. Yet, few studies have investigated rates of EDs and their symptoms in individuals presenting with MDD/anxiety disorders. Identifying potential disordered eating in people with MDD/anxiety disorders is important because even subclinical disordered eating is associated with reduced quality of life, and undiagnosed eating pathology may hinder treatment progress for both MDD/anxiety disorders and comorbid EDs.

METHOD

We compared rates of EDs (anorexia nervosa, bulimia nervosa, binge-eating disorder, and other specified feeding and eating disorders) and their symptoms in 130 women with, and 405 women without, lifetime MDD or an anxiety disorder (generalized anxiety disorder, obsessive-compulsive disorder, social phobia, specific phobia, panic disorder, agoraphobia, or post-traumatic stress disorder) recruited from the population-based Michigan State University Twin Registry. Lifetime ED and MDD/anxiety diagnoses, and lifetime clinically significant disordered eating behaviors (e.g., binge eating, excessive exercise) were assessed using the Structured Clinical Interview for DSM-IV (SCID).

RESULTS

Among participants with lifetime MDD or any anxiety disorder, 13% met criteria for a lifetime ED and 39% reported engaging in at least one lifetime clinically significant disordered eating behavior (e.g., binge eating) on the SCID. In contrast, only 3% of participants without a history of MDD/an anxiety disorder met criteria for a lifetime ED, and only 11% reported lifetime clinically significant disordered eating behavior.

DISCUSSION

Our findings suggest that women with MDD and anxiety disorders have elevated rates of EDs, and it is therefore imperative to screen for disordered eating in these populations.

Mendelian randomization study of maternal influences on birthweight and future cardiometabolic risk in the HUNT cohort

There is a robust observational relationship between lower birthweight and higher risk of cardiometabolic disease in later life. The Developmental Origins of Health and Disease (DOHaD) hypothesis posits that adverse environmental factors in utero increase future risk of cardiometabolic disease. Here, we explore if a genetic risk score (GRS) of maternal SNPs associated with offspring birthweight is also associated with offspring cardiometabolic risk factors, after controlling for offspring GRS, in up to 26,057 mother-offspring pairs (and 19,792 father-offspring pairs) from the Nord-Trøndelag Health (HUNT) Study. We find little evidence for a maternal (or paternal) genetic effect of birthweight associated variants on offspring cardiometabolic risk factors after adjusting for offspring GRS. In contrast, offspring GRS is strongly related to many cardiometabolic risk factors, even after conditioning on maternal GRS. Our results suggest that the maternal intrauterine environment, as proxied by maternal SNPs that influence offspring birthweight, is unlikely to be a major determinant of adverse cardiometabolic outcomes in population based samples of individuals.