A high-resolution transcriptomic and spatial atlas of cell types in the whole mouse brain.

The mammalian brain consists of millions to billions of cells that are organized into many cell types with specific spatial distribution patterns and structural and functional properties1-3. Here we report a comprehensive and high-resolution transcriptomic and spatial cell-type atlas for the whole adult mouse brain. The cell-type atlas was created by combining a single-cell RNA-sequencing (scRNA-seq) dataset of around 7 million cells profiled (approximately 4.0 million cells passing quality control), and a spatial transcriptomic dataset of approximately 4.3 million cells using multiplexed error-robust fluorescence in situ hybridization (MERFISH). The atlas is hierarchically organized into 4 nested levels of classification: 34 classes, 338 subclasses, 1,201 supertypes and 5,322 clusters. We present an online platform, Allen Brain Cell Atlas, to visualize the mouse whole-brain cell-type atlas along with the single-cell RNA-sequencing and MERFISH datasets. We systematically analysed the neuronal and non-neuronal cell types across the brain and identified a high degree of correspondence between transcriptomic identity and spatial specificity for each cell type. The results reveal unique features of cell-type organization in different brain regions-in particular, a dichotomy between the dorsal and ventral parts of the brain. The dorsal part contains relatively fewer yet highly divergent neuronal types, whereas the ventral part contains more numerous neuronal types that are more closely related to each other. Our study also uncovered extraordinary diversity and heterogeneity in neurotransmitter and neuropeptide expression and co-expression patterns in different cell types. Finally, we found that transcription factors are major determinants of cell-type classification and identified a combinatorial transcription factor code that defines cell types across all parts of the brain. The whole mouse brain transcriptomic and spatial cell-type atlas establishes a benchmark reference atlas and a foundational resource for integrative investigations of cellular and circuit function, development and evolution of the mammalian brain.

Personal omics profiling reveals dynamic molecular and medical phenotypes.

Personalized medicine is expected to benefit from combining genomic information with regular monitoring of physiological states by multiple high-throughput methods. Here, we present an integrative personal omics profile (iPOP), an analysis that combines genomic, transcriptomic, proteomic, metabolomic, and autoantibody profiles from a single individual over a 14 month period. Our iPOP analysis revealed various medical risks, including type 2 diabetes. It also uncovered extensive, dynamic changes in diverse molecular components and biological pathways across healthy and diseased conditions. Extremely high-coverage genomic and transcriptomic data, which provide the basis of our iPOP, revealed extensive heteroallelic changes during healthy and diseased states and an unexpected RNA editing mechanism. This study demonstrates that longitudinal iPOP can be used to interpret healthy and diseased states by connecting genomic information with additional dynamic omics activity.

The conformational cycle of prestin underlies outer-hair cell electromotility.

The voltage-dependent motor protein prestin (also known as SLC26A5) is responsible for the electromotive behaviour of outer-hair cells and underlies the cochlear amplifier1. Knockout or impairment of prestin causes severe hearing loss2-5. Despite the key role of prestin in hearing, the mechanism by which mammalian prestin senses voltage and transduces it into cellular-scale movements (electromotility) is poorly understood. Here we determined the structure of dolphin prestin in six distinct states using single-particle cryo-electron microscopy. Our structural and functional data suggest that prestin adopts a unique and complex set of states, tunable by the identity of bound anions (Cl- or SO42-). Salicylate, a drug that can cause reversible hearing loss, competes for the anion-binding site of prestin, and inhibits its function by immobilizing prestin in a new conformation. Our data suggest that the bound anion together with its coordinating charged residues and helical dipole act as a dynamic voltage sensor. An analysis of all of the anion-dependent conformations reveals how structural rearrangements in the voltage sensor are coupled to conformational transitions at the protein-membrane interface, suggesting a previously undescribed mechanism of area expansion. Visualization of the electromotility cycle of prestin distinguishes the protein from the closely related SLC26 anion transporters, highlighting the basis for evolutionary specialization of the mammalian cochlear amplifier at a high resolution.

Electromechanical coupling in the hyperpolarization-activated K+ channel KAT1.

Voltage-gated potassium (Kv) channels coordinate electrical signalling and control cell volume by gating in response to membrane depolarization or hyperpolarization. However, although voltage-sensing domains transduce transmembrane electric field changes by a common mechanism involving the outward or inward translocation of gating charges1-3, the general determinants of channel gating polarity remain poorly understood4. Here we suggest a molecular mechanism for electromechanical coupling and gating polarity in non-domain-swapped Kv channels on the basis of the cryo-electron microscopy structure of KAT1, the hyperpolarization-activated Kv channel from Arabidopsis thaliana. KAT1 displays a depolarized voltage sensor, which interacts with a closed pore domain directly via two interfaces and indirectly via an intercalated phospholipid. Functional evaluation of KAT1 structure-guided mutants at the sensor-pore interfaces suggests a mechanism in which direct interaction between the sensor and the C-linker hairpin in the adjacent pore subunit is the primary determinant of gating polarity. We suggest that an inward motion of the S4 sensor helix of approximately 5-7 Å can underlie a direct-coupling mechanism, driving a conformational reorientation of the C-linker and ultimately opening the activation gate formed by the S6 intracellular bundle. This direct-coupling mechanism contrasts with allosteric mechanisms proposed for hyperpolarization-activated cyclic nucleotide-gated channels5, and may represent an unexpected link between depolarization- and hyperpolarization-activated channels.

IsoVis - a webserver for visualization and annotation of alternative RNA isoforms.

Genes commonly express multiple RNA products (RNA isoforms), which differ in exonic content and can have different functions. Making sense of the plethora of known and novel RNA isoforms being identified by transcriptomic approaches requires a user-friendly way to visualize gene isoforms and how they differ in exonic content, expression levels and potential functions. Here we introduce IsoVis, a freely available webserver that accepts user-supplied transcriptomic data and visualizes the expressed isoforms in a clear, intuitive manner. IsoVis contains numerous features, including the ability to visualize all RNA isoforms of a gene and their expression levels; the annotation of known isoforms from external databases; mapping of protein domains and features to exons, allowing changes to protein sequence and function between isoforms to be established; and extensive species compatibility. Datasets visualised on IsoVis remain private to the user, allowing analysis of sensitive data. IsoVis visualisations can be downloaded to create publication-ready figures. The IsoVis webserver enables researchers to perform isoform analyses without requiring programming skills, is free to use, and available at https://isomix.org/isovis/.

Estimating the environmental impacts of 57,000 food products.

Understanding and communicating the environmental impacts of food products is key to enabling transitions to environmentally sustainable food systems [El Bilali and Allahyari, Inf. Process. Agric. 5, 456-464 (2018)]. While previous analyses compared the impacts of food commodities such as fruits, wheat, and beef [Poore and Nemecek, Science 360, 987-992 (2018)], most food products contain numerous ingredients. However, because the amount of each ingredient in a product is often known only by the manufacturer, it has been difficult to assess their environmental impacts. Here, we develop an approach to overcome this limitation. It uses prior knowledge from ingredient lists to infer the composition of each ingredient, and then pairs this with environmental databases [Poore and Nemecek Science 360, 987-992 (2018); Gephart et al., Nature 597, 360-365 (2021)] to derive estimates of a food product's environmental impact across four indicators: greenhouse gas emissions, land use, water stress, and eutrophication potential. Using the approach on 57,000 products in the United Kingdom and Ireland shows food types have low (e.g., sugary beverages, fruits, breads), to intermediate (e.g., many desserts, pastries), to high environmental impacts (e.g., meat, fish, cheese). Incorporating NutriScore reveals more nutritious products are often more environmentally sustainable but there are exceptions to this trend, and foods consumers may view as substitutable can have markedly different impacts. Sensitivity analyses indicate the approach is robust to uncertainty in ingredient composition and in most cases sourcing. This approach provides a step toward enabling consumers, retailers, and policy makers to make informed decisions on the environmental impacts of food products.

Discovery of NRG1-VII: the myeloid-derived class of NRG1.

The growth factor Neuregulin-1 (NRG1) has pleiotropic roles in proliferation and differentiation of the stem cell niche in different tissues. It has been implicated in gut, brain and muscle development and repair. Six isoform classes of NRG1 and over 28 protein isoforms have been previously described. Here we report a new class of NRG1, designated NRG1-VII to denote that these NRG1 isoforms arise from a myeloid-specific transcriptional start site (TSS) previously uncharacterized. Long-read sequencing was used to identify eight high-confidence NRG1-VII transcripts. These transcripts presented major structural differences from one another, through the use of cassette exons and alternative stop codons. Expression of NRG1-VII was confirmed in primary human monocytes and tissue resident macrophages and induced pluripotent stem cell-derived macrophages (iPSC-derived macrophages). Isoform switching via cassette exon usage and alternate polyadenylation was apparent during monocyte maturation and macrophage differentiation. NRG1-VII is the major class expressed by the myeloid lineage, including tissue-resident macrophages. Analysis of public gene expression data indicates that monocytes and macrophages are a primary source of NRG1. The size and structure of class VII isoforms suggests that they may be more diffusible through tissues than other NRG1 classes. However, the specific roles of class VII variants in tissue homeostasis and repair have not yet been determined.

Standardizing to specific target populations in distributed networks and multisite pharmacoepidemiologic studies.

Distributed network studies and multisite studies assess drug safety and effectiveness in diverse populations by pooling information. Targeting groups of clinical or policy interest (including specific sites or site combinations) and applying weights based on effect measure modifiers (EMMs) prior to pooling estimates within multisite studies may increase interpretability and improve precision. We simulated a 4-site study, standardized each site using inverse odds weights (IOWs) to resemble the 3 smallest sites or the smallest site, estimated IOW-weighted risk differences (RDs), and combined estimates with inverse variance weights (IVWs). We also created an artificial distributed network in the Clinical Practice Research Datalink (CPRD) Aurum consisting of 1 site for each geographic region. We compared metformin and sulfonylurea initiators with respect to mortality, targeting the smallest region. In the simulation, IOWs reduced differences between estimates and increased precision when targeting the 3 smallest sites or the smallest site. In the CPRD Aurum study, the IOW + IVW estimate was also more precise (smallest region: RD = 5.41% [95% CI, 1.03-9.79]; IOW + IVW estimate: RD = 3.25% [95% CI, 3.07-3.43]). When performing pharmacoepidemiologic research in distributed networks or multisite studies in the presence of EMMs, designation of target populations has the potential to improve estimate precision and interpretability. This article is part of a Special Collection on Pharmacoepidemiology.

Variable selection when estimating effects in external target populations.

External validity is an important part of epidemiologic research. To validly estimate effects in specific external target populations using a chosen effect measure (ie, "transport"), some methods require that one account for all effect measure modifiers (EMMs). However, little is known about how including other variables that are not EMMs (ie, non-EMMs) in adjustment sets affects estimates. Using simulations, we evaluated how inclusion of non-EMMs affected estimation of the transported risk difference (RD) by assessing the impacts of covariates that (1) differ (or not) between the trial and the target, (2) are associated with the outcome (or not), and (3) modify the RD (or not). We assessed variation and bias when covariates with each possible combination of these factors were used to transport RDs using outcome modeling or inverse odds weighting. Inclusion of variables that differed in distribution between the populations but were non-EMMs reduced precision, regardless of whether they were associated with the outcome. However, non-EMMs associated with selection did not amplify bias resulting from omission of necessary EMMs. Including all variables associated with the outcome may result in unnecessarily imprecise estimates when estimating treatment effects in external target populations.

Targeted single-cell RNA sequencing of transcription factors enhances the identification of cell types and trajectories.

Single-cell RNA sequencing (scRNA-seq) is a widely used method for identifying cell types and trajectories in biologically heterogeneous samples, but it is limited in its detection and quantification of lowly expressed genes. This results in missing important biological signals, such as the expression of key transcription factors (TFs) driving cellular differentiation. We show that targeted sequencing of ∼1000 TFs (scCapture-seq) in iPSC-derived neuronal cultures greatly improves the biological information garnered from scRNA-seq. Increased TF resolution enhanced cell type identification, developmental trajectories, and gene regulatory networks. This allowed us to resolve differences among neuronal populations, which were generated in two different laboratories using the same differentiation protocol. ScCapture-seq improved TF-gene regulatory network inference and thus identified divergent patterns of neurogenesis into either excitatory cortical neurons or inhibitory interneurons. Furthermore, scCapture-seq revealed a role for of retinoic acid signaling in the developmental divergence between these different neuronal populations. Our results show that TF targeting improves the characterization of human cellular models and allows identification of the essential differences between cellular populations, which would otherwise be missed in traditional scRNA-seq. scCapture-seq TF targeting represents a cost-effective enhancement of scRNA-seq, which could be broadly applied to improve scRNA-seq resolution.

Visualizing External Validity: Graphical Displays to Inform the Extension of Treatment Effects from Trials to Clinical Practice.

BACKGROUND: In the presence of effect measure modification, estimates of treatment effects from randomized controlled trials may not be valid in clinical practice settings. The development and application of quantitative approaches for extending treatment effects from trials to clinical practice settings is an active area of research. METHODS: In this article, we provide researchers with a practical roadmap and four visualizations to assist in variable selection for models to extend treatment effects observed in trials to clinical practice settings and to assess model specification and performance. We apply this roadmap and visualizations to an example extending the effects of adjuvant chemotherapy (5-fluorouracil vs. plus oxaliplatin) for colon cancer from a trial population to a population of individuals treated in community oncology practices in the United States. RESULTS: The first visualization screens for potential effect measure modifiers to include in models extending trial treatment effects to clinical practice populations. The second visualization displays a measure of covariate overlap between the clinical practice populations and the trial population. The third and fourth visualizations highlight considerations for model specification and influential observations. The conceptual roadmap describes how the output from the visualizations helps interrogate the assumptions required to extend treatment effects from trials to target populations. CONCLUSIONS: The roadmap and visualizations can inform practical decisions required for quantitatively extending treatment effects from trials to clinical practice settings.

Interplay between gain and loss in arrays of nonlinear plasmonic nanoparticles: toward parametric downconversion and amplification.

With the help of a theoretical model and finite-difference time-domain (FDTD) simulations based on the hydrodynamic-Maxwell model, we examine the effect of difference-frequency generation (DFG) in an array of L-shaped metal nanoparticles (MNPs) characterized by intrinsic plasmonic nonlinearity. The outcomes of the calculations reveal the spectral interplay between gain and loss in the vicinity of the fundamental frequency of the localized surface plasmon resonances. Subsequently, we identify different array thicknesses and pumping regimes facilitating parametric amplification and spontaneous parametric downconversion. Our results suggest that the parametric amplification regime becomes feasible on a scale of hundreds of nanometers and spontaneous parametric downconversion on the scale of tens of nanometers, opening up new exciting opportunities for developing building blocks of photonic metasurfaces.

Options for keeping the food system within environmental limits.

The food system is a major driver of climate change, changes in land use, depletion of freshwater resources, and pollution of aquatic and terrestrial ecosystems through excessive nitrogen and phosphorus inputs. Here we show that between 2010 and 2050, as a result of expected changes in population and income levels, the environmental effects of the food system could increase by 50-90% in the absence of technological changes and dedicated mitigation measures, reaching levels that are beyond the planetary boundaries that define a safe operating space for humanity. We analyse several options for reducing the environmental effects of the food system, including dietary changes towards healthier, more plant-based diets, improvements in technologies and management, and reductions in food loss and waste. We find that no single measure is enough to keep these effects within all planetary boundaries simultaneously, and that a synergistic combination of measures will be needed to sufficiently mitigate the projected increase in environmental pressures.

Clarifying the causal contrast: An empirical example applying the prevalent new user study design.

PURPOSE: The prevalent new user design extends the active comparator new user design to include patients switching to a treatment of interest from a comparator. We examined the impact of adding "switchers" to incident new users on the estimated hazard ratio (HR) of hospitalized heart failure. METHODS: Using MarketScan claims data (2000-2014), we estimated HRs of hospitalized heart failure between patients initiating GLP-1 receptor agonists (GLP-1 RA) and sulfonylureas (SU). We considered three estimands: (1) the effect of incident new use; (2) the effect of switching; and (3) the effect of incident new use or switching, combining the two population. We used time-conditional propensity scores (TCPS) and time-stratified standardized morbidity ratio (SMR) weighting to adjust for confounding. RESULTS: We identified 76 179 GLP-1 RA new users, of which 12% were direct switchers (within 30 days) from SU. Among incident new users, GLP-1 RA was protective against heart failure (adjHRSMR = 0.74 [0.69, 0.80]). Among switchers, GLP-1 RA was not protective (adjHRSMR = 0.99 [0.83, 1.18]). Results in the combined population were largely driven by the incident new users, with GLP-1 RA having a protective effect (adjHRSMR = 0.77 [0.72, 0.83]). Results using TCPS were consistent with those estimated using SMR weighting. CONCLUSIONS: When analyses were conducted only among incident new users, GLP-1 RA had a protective effect. However, among switchers from SU to GLP-1 RA, the effect estimates substantially shifted toward the null. Combining patients with varying treatment histories can result in poor confounding control and camouflage important heterogeneity.

Harnessing complexity: Nonlinear optical phenomena in L-shapes, nanocrescents, and split-ring resonators.

We conduct systematic studies of the optical characteristics of plasmonic nanoparticles that exhibit C2v symmetry. In particular, we analyze three distinct geometric configurations: an L-type shape, a crescent, and a split-ring resonator shaped like the Greek letter π. Optical properties are examined using the finite-difference time-domain method. It is demonstrated that all three shapes exhibit two prominent plasmon modes associated with the two axes of symmetry. This is in addition to a wide range of resonances observed at high frequencies corresponding to quadrupole modes and peaks due to sharp corners. Next, to facilitate nonlinear analysis, we employ a semiclassical hydrodynamic model, where the electron pressure term is explicitly accounted for. This model goes beyond the standard Drude description and enables capturing nonlocal and nonlinear effects. Employing this model enables us to rigorously examine the second-order angular resolved nonlinear optical response of these nanoparticles in each of the three configurations. Two pumping regimes are considered, namely, continuous wave (CW) and pulsed excitations. For CW pumping, we explore the properties of the second harmonic generation (SHG). Polarization and angle-resolved SHG spectra are obtained, revealing strong dependence on the nanoparticle geometry and incident wave polarization. The C2v symmetry is shown to play a key role in determining the polarization states and selection rules of the SHG signal. For pulsed excitations, we discuss the phenomenon of broadband terahertz (THz) generation induced by the difference-frequency generation . It is shown that the THz emission spectra exhibit unique features attributed to the plasmonic resonances and symmetry of the nanoparticles. The polarization of the generated THz waves is also examined, revealing interesting patterns tied to the nanoparticle geometry. To gain deeper insight, we propose an analytical theory that agrees very well with the numerical experiments. The theory shows that the physical origin of the THz radiation is the mixing of various frequency components of the fundamental pulse by the second-order nonlinear susceptibility. An expression for the far-field THz intensity is derived in terms of the incident pulse parameters and the nonlinear response tensor of the nanoparticle. The results presented in this work offer new insights into the linear and nonlinear optical properties of nanoparticles with C2v symmetry. The demonstrated strong SHG response and efficient broadband THz generation hold great promise for applications in nonlinear spectroscopy, nanophotonics, and optoelectronics. The proposed theoretical framework also provides a valuable tool for understanding and predicting the nonlinear behavior of other related nanostructures.

Accurate expression quantification from nanopore direct RNA sequencing with NanoCount.

Accurately quantifying gene and isoform expression changes is essential to understanding cell functions, differentiation and disease. Sequencing full-length native RNAs using long-read direct RNA sequencing (DRS) has the potential to overcome many limitations of short and long-read sequencing methods that require RNA fragmentation, cDNA synthesis or PCR. However, there are a lack of tools specifically designed for DRS and its ability to identify differential expression in complex organisms is poorly characterised. We developed NanoCount for fast, accurate transcript isoform quantification in DRS and demonstrate it outperforms similar methods. Using synthetic controls and human SH-SY5Y cell differentiation into neuron-like cells, we show that DRS accurately quantifies RNA expression and identifies differential expression of genes and isoforms. Differential expression of 231 genes, 333 isoforms, plus 27 isoform switches were detected between undifferentiated and differentiated SH-SY5Y cells and samples clustered by differentiation state at the gene and isoform level. Genes upregulated in neuron-like cells were associated with neurogenesis. NanoCount quantification of thousands of novel isoforms discovered with DRS likewise enabled identification of their differential expression. Our results demonstrate enhanced DRS isoform quantification with NanoCount and establish the ability of DRS to identify biologically relevant differential expression of genes and isoforms.

Air quality-related health damages of food.

Agriculture is a major contributor to air pollution, the largest environmental risk factor for mortality in the United States and worldwide. It is largely unknown, however, how individual foods or entire diets affect human health via poor air quality. We show how food production negatively impacts human health by increasing atmospheric fine particulate matter (PM2.5), and we identify ways to reduce these negative impacts of agriculture. We quantify the air quality-related health damages attributable to 95 agricultural commodities and 67 final food products, which encompass >99% of agricultural production in the United States. Agricultural production in the United States results in 17,900 annual air quality-related deaths, 15,900 of which are from food production. Of those, 80% are attributable to animal-based foods, both directly from animal production and indirectly from growing animal feed. On-farm interventions can reduce PM2.5-related mortality by 50%, including improved livestock waste management and fertilizer application practices that reduce emissions of ammonia, a secondary PM2.5 precursor, and improved crop and animal production practices that reduce primary PM2.5 emissions from tillage, field burning, livestock dust, and machinery. Dietary shifts toward more plant-based foods that maintain protein intake and other nutritional needs could reduce agricultural air quality-related mortality by 68 to 83%. In sum, improved livestock and fertilization practices, and dietary shifts could greatly decrease the health impacts of agriculture caused by its contribution to reduced air quality.

Increased paternal corticosterone exposure influences offspring behaviour and expression of urinary pheromones.

BACKGROUND: Studies have shown that paternal stress prior to conception can influence the innate behaviours of their offspring. The evolutionary impacts of such intergenerational effects are therefore of considerable interest. Our group previously showed in a model of daily stress that glucocorticoid treatment of adult male mouse breeders prior to conception leads to increased anxiety-related behaviours in male offspring. Here, we aimed to understand the transgenerational effects of paternal stress exposure on the social behaviour of progeny and its potential influence on reproductive success. RESULTS: We assessed social parameters including social reward, male attractiveness and social dominance, in the offspring (F1) and grand-offspring (F2). We report that paternal corticosterone treatment was associated with increased display of subordination towards other male mice. Those mice were unexpectedly more attractive to female mice while expressing reduced levels of the key rodent pheromone Darcin, contrary to its conventional role in driving female attraction. We investigated the epigenetic regulation of major urinary protein (Mup) expression by performing the first Oxford Nanopore direct methylation of sperm DNA in a mouse model of stress, but found no differences in Mup genes that could be attributed to corticosterone-treatment. Furthermore, no overt differences of the prefrontal cortex transcriptome were found in F1 offspring, implying that peripheral mechanisms are likely contributing to the phenotypic differences. Interestingly, no phenotypic differences were observed in the F2 grand-offspring. CONCLUSIONS: Overall, our findings highlight the potential of moderate paternal stress to affect intergenerational (mal)adaptive responses, informing future studies of adaptiveness in rodents, humans and other species.

Use of Next-Generation Sequencing in a State-Wide Strategy of HIV-1 Surveillance: Impact of the SARS-COV-2 Pandemic on HIV-1 Diagnosis and Transmission.

BACKGROUND: The ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic posed an unpreceded threat to the management of other pandemics such as human immunodeficiency virus-1 (HIV-1) in the United States. The full impact of the SARS-CoV-2 pandemic on the HIV-1 pandemic needs to be evaluated. METHODS: All individuals with newly reported HIV-1 diagnoses from NC State Laboratory of Public Health were enrolled in this prospective observational study, 2018-2021. We used a sequencing-based recency assay to identify recent HIV-1 infections and to determine the days postinfection (DPI) for each person at the time of diagnosis. RESULTS: Sequencing used diagnostic serum samples from 814 individuals with new HIV-1 diagnoses spanning this 4-year period. Characteristics of individuals diagnosed in 2020 differed from those in other years. People of color diagnosed in 2021 were on average 6 months delayed in their diagnosis compared to those diagnosed in 2020. There was a trend that genetic networks were more known for individuals diagnosed in 2021. We observed no major integrase resistance mutations over the course of the study. CONCLUSIONS: SARS-CoV-2 pandemic may contribute to the spread of HIV-1. Public health resources need to focus on restoring HIV-1 testing and interrupting active, ongoing, transmission.

How Choice of Effect Measure Influences Minimally Sufficient Adjustment Sets for External Validity.

Epidemiologic researchers generalizing or transporting effect estimates from a study to a target population must account for effect-measure modifiers (EMMs) on the scale of interest. However, little attention is paid to how the EMMs required may vary depending on the mathematical nuances of each effect measure. We defined 2 types of EMMs: a marginal EMM, where the effect on the scale of interest differs across levels of a variable, and a conditional EMM, where the effect differs conditional on other variables associated with the outcome. These types define 3 classes of variables: class 1 (conditional EMM), class 2 (marginal but not conditional EMM), and class 3 (neither marginal nor conditional EMM). Class 1 variables are necessary to achieve a valid estimate of the RD in a target population, while an RR requires class 1 and class 2 and an OR requires classes 1, 2, and 3 (i.e., all variables associated with the outcome). This does not mean that fewer variables are required for an externally valid RD (because variables may not modify effects on all scales), but it does suggest that researchers should consider the scale of the effect measure when identifying an EMM necessary for an externally valid treatment effect estimate.