Uncovering methylation-dependent genetic effects on regulatory element function in diverse genomes.

A major goal in evolutionary biology and biomedicine is to understand the complex interactions between genetic variants, the epigenome, and gene expression. However, the causal relationships between these factors remain poorly understood. mSTARR-seq, a methylation-sensitive massively parallel reporter assay, is capable of identifying methylation-dependent regulatory activity at many thousands of genomic regions simultaneously, and allows for the testing of causal relationships between DNA methylation and gene expression on a region-by-region basis. Here, we developed a multiplexed mSTARR-seq protocol to assay naturally occurring human genetic variation from 25 individuals sampled from 10 localities in Europe and Africa. We identified 6,957 regulatory elements in either the unmethylated or methylated state, and this set was enriched for enhancer and promoter annotations, as expected. The expression of 58% of these regulatory elements was modulated by methylation, which was generally associated with decreased RNA expression. Within our set of regulatory elements, we used allele-specific expression analyses to identify 8,020 sites with genetic effects on gene regulation; further, we found that 42.3% of these genetic effects varied between methylated and unmethylated states. Sites exhibiting methylation-dependent genetic effects were enriched for GWAS and EWAS annotations, implicating them in human disease. Compared to datasets that assay DNA from a single European individual, our multiplexed assay uncovers dramatically more genetic effects and methylation-dependent genetic effects, highlighting the importance of including diverse individuals in assays which aim to understand gene regulatory processes.

Cost-effective solutions for high-throughput enzymatic DNA methylation sequencing.

Characterizing DNA methylation patterns is important for addressing key questions in evolutionary biology, geroscience, and medical genomics. While costs are decreasing, whole-genome DNA methylation profiling remains prohibitively expensive for most population-scale studies, creating a need for cost-effective, reduced representation approaches (i.e., assays that rely on microarrays, enzyme digests, or sequence capture to target a subset of the genome). Most common whole genome and reduced representation techniques rely on bisulfite conversion, which can damage DNA resulting in DNA loss and sequencing biases. Enzymatic methyl sequencing (EM-seq) was recently proposed to overcome these issues, but thorough benchmarking of EM-seq combined with cost-effective, reduced representation strategies has not yet been performed. To do so, we optimized Targeted Methylation Sequencing protocol (TMS)-which profiles ∼4 million CpG sites-for miniaturization, flexibility, and multispecies use at a cost of ∼$80. First, we tested modifications to increase throughput and reduce cost, including increasing multiplexing, decreasing DNA input, and using enzymatic rather than mechanical fragmentation to prepare DNA. Second, we compared our optimized TMS protocol to commonly used techniques, specifically the Infinium MethylationEPIC BeadChip (n=55 paired samples) and whole genome bisulfite sequencing (n=6 paired samples). In both cases, we found strong agreement between technologies (R² = 0.97 and 0.99, respectively). Third, we tested the optimized TMS protocol in three non-human primate species (rhesus macaques, geladas, and capuchins). We captured a high percentage (mean=77.1%) of targeted CpG sites and produced methylation level estimates that agreed with those generated from reduced representation bisulfite sequencing (R² = 0.98). Finally, we applied our protocol to profile age-associated DNA methylation variation in two subsistence-level populations-the Tsimane of lowland Bolivia and the Orang Asli of Peninsular Malaysia-and found age-methylation patterns that were strikingly similar to those reported in high income cohorts, despite known differences in age-health relationships between lifestyle contexts. Altogether, our optimized TMS protocol will enable cost-effective, population-scale studies of genome-wide DNA methylation levels across human and non-human primate species.

The built environment is more predictive of cardiometabolic health than other aspects of lifestyle in two rapidly transitioning Indigenous populations.

Background: Many subsistence-level and Indigenous societies around the world are rapidly experiencing urbanization, nutrition transition, and integration into market-economies, resulting in marked increases in cardiometabolic diseases. Determining the most potent and generalized drivers of changing health is essential for identifying vulnerable communities and creating effective policies to combat increased chronic disease risk across socio-environmental contexts. However, comparative tests of how different lifestyle features affect the health of populations undergoing lifestyle transitions remain rare, and require comparable, integrated anthropological and health data collected in diverse contexts. Methods: We developed nine scales to quantify different facets of lifestyle (e.g., urban infrastructure, market-integration, acculturation) in two Indigenous, transitioning subsistence populations currently undergoing rapid change in very different ecological and sociopolitical contexts: Turkana pastoralists of northwest Kenya (n = 3,692) and Orang Asli mixed subsistence groups of Peninsular Malaysia (n = 688). We tested the extent to which these lifestyle scales predicted 16 measures of cardiometabolic health and compared the generalizability of each scale across the two populations. We used factor analysis to decompose comprehensive lifestyle data into salient axes without supervision, sensitivity analyses to understand which components of the multidimensional scales were most important, and sex-stratified analyses to understand how facets of lifestyle variation differentially impacted cardiometabolic health among males and females. Findings: Cardiometabolic health was best predicted by measures that quantified urban infrastructure and market-derived material wealth compared to metrics encompassing diet, mobility, or acculturation, and these results were highly consistent across both populations and sexes. Factor analysis results were also highly consistent between the Turkana and Orang Asli and revealed that lifestyle variation decomposes into two distinct axes-the built environment and diet-which change at different paces and have different relationships with health. Interpretation: Our analysis of comparable data from Indigenous peoples in East Africa and Southeast Asia revealed a surprising amount of generalizability: in both contexts, measures of local infrastructure and built environment are consistently more predictive of cardiometabolic health than other facets of lifestyle that are seemingly more proximate to health, such as diet. We hypothesize that this is because the built environment impacts unmeasured proximate drivers like physical activity, increased stress, and broader access to market goods, and serves as a proxy for the duration of time that communities have been market-integrated.

Immune gene regulation is associated with age and environmental adversity in a nonhuman primate.

Phenotypic aging is ubiquitous across mammalian species, suggesting shared underlying mechanisms of aging. Aging is linked to molecular changes to DNA methylation and gene expression, and environmental factors, such as severe external challenges or adversities, can moderate these age-related changes. Yet, it remains unclear whether environmental adversities affect gene regulation via the same molecular pathways as chronological, or 'primary', aging. Investigating molecular aging in naturalistic animal populations can fill this gap by providing insight into shared molecular mechanisms of aging and the effects of a greater diversity of environmental adversities - particularly those that can be challenging to study in humans or laboratory organisms. Here, we characterised molecular aging - specifically, CpG methylation - in a sample of free-ranging rhesus macaques living off the coast of Puerto Rico (n samples = 571, n individuals = 499), which endured a major hurricane during our study. Age was associated with methylation at 78,661 sites (31% of all sites tested). Age-associated hypermethylation occurred more frequently in areas of active gene regulation, while hypomethylation was enriched in regions that show less activity in immune cells, suggesting these regions may become de-repressed in older individuals. Age-associated hypomethylation also co-occurred with increased chromatin accessibility while hypermethylation showed the opposite trend, hinting at a coordinated, multi-level loss of epigenetic stability during aging. We detected 32,048 CpG sites significantly associated with exposure to a hurricane, and these sites overlapped age-associated sites, most strongly in regulatory regions and most weakly in quiescent regions. Together, our results suggest that environmental adversity may contribute to aging-related molecular phenotypes in regions of active gene transcription, but that primary aging has specific signatures in non-regulatory regions.

Testing frameworks for early life effects: the developmental constraints and adaptive response hypotheses do not explain key fertility outcomes in wild female baboons.

In evolutionary ecology, two classes of explanations are frequently invoked to explain "early life effects" on adult outcomes. Developmental constraints (DC) explanations contend that costs of early adversity arise from limitations adversity places on optimal development. Adaptive response (AR) hypotheses propose that later life outcomes will be worse when early and adult environments are poorly "matched." Here, we use recently proposed mathematical definitions for these hypotheses and a quadratic-regression based approach to test the long-term consequences of variation in developmental environments on fertility in wild baboons. We evaluate whether low rainfall and/or dominance rank during development predict three female fertility measures in adulthood, and whether any observed relationships are consistent with DC and/or AR. Neither rainfall during development nor the difference between rainfall in development and adulthood predicted any fertility measures. Females who were low-ranking during development had an elevated risk of losing infants later in life, and greater change in rank between development and adulthood predicted greater risk of infant loss. However, both effects were statistically marginal and consistent with alternative explanations, including adult environmental quality effects. Consequently, our data do not provide compelling support for either of these common explanations for the evolution of early life effects.

DNA methylation signatures of early-life adversity are exposure-dependent in wild baboons.

The early-life environment can profoundly shape the trajectory of an animal's life, even years or decades later. One mechanism proposed to contribute to these early-life effects is DNA methylation. However, the frequency and functional importance of DNA methylation in shaping early-life effects on adult outcomes is poorly understood, especially in natural populations. Here, we integrate prospectively collected data on fitness-associated variation in the early environment with DNA methylation estimates at 477,270 CpG sites in 256 wild baboons. We find highly heterogeneous relationships between the early-life environment and DNA methylation in adulthood: aspects of the environment linked to resource limitation (e.g., low-quality habitat, early-life drought) are associated with many more CpG sites than other types of environmental stressors (e.g., low maternal social status). Sites associated with early resource limitation are enriched in gene bodies and putative enhancers, suggesting they are functionally relevant. Indeed, by deploying a baboon-specific, massively parallel reporter assay, we show that a subset of windows containing these sites are capable of regulatory activity, and that, for 88% of early drought-associated sites in these regulatory windows, enhancer activity is DNA methylation-dependent. Together, our results support the idea that DNA methylation patterns contain a persistent signature of the early-life environment. However, they also indicate that not all environmental exposures leave an equivalent mark and suggest that socioenvironmental variation at the time of sampling is more likely to be functionally important. Thus, multiple mechanisms must converge to explain early-life effects on fitness-related traits.

Diet composition impacts eQTL discovery across multiple tissues in baboons.

Understanding how genetic variation impacts gene expression is a major goal of genomics; however, only a fraction of disease-associated loci have been demonstrated to impact gene expression when cells are in an unperturbed "steady state." In this issue of Cell Genomics, Lin et al.1 investigate how exposure to a particular cellular context (i.e., a high-cholesterol, high-fat diet) can enhance our ability to identify new regulatory variants through longitudinal sampling of three tissue types in the baboon.

DNA methylation-environment interactions in the human genome.

Previously, we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here, we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Chromatin activity identifies differential gene regulation across human ancestries.

BACKGROUND: Current evidence suggests that cis-regulatory elements controlling gene expression may be the predominant target of natural selection in humans and other species. Detecting selection acting on these elements is critical to understanding evolution but remains challenging because we do not know which mutations will affect gene regulation. RESULTS: To address this, we devise an approach to search for lineage-specific selection on three critical steps in transcriptional regulation: chromatin activity, transcription factor binding, and chromosomal looping. Applying this approach to lymphoblastoid cells from 831 individuals of either European or African descent, we find strong signals of differential chromatin activity linked to gene expression differences between ancestries in numerous contexts, but no evidence of functional differences in chromosomal looping. Moreover, we show that enhancers rather than promoters display the strongest signs of selection associated with sites of differential transcription factor binding. CONCLUSIONS: Overall, our study indicates that some cis-regulatory adaptation may be more easily detected at the level of chromatin than DNA sequence. This work provides a vast resource of genomic interaction data from diverse human populations and establishes a novel selection test that will benefit future study of regulatory evolution in humans and other species.

Applying an evolutionary mismatch framework to understand disease susceptibility.

Noncommunicable diseases (NCDs) are on the rise worldwide. Obesity, cardiovascular disease, and type 2 diabetes are among a long list of "lifestyle" diseases that were rare throughout human history but are now common. The evolutionary mismatch hypothesis posits that humans evolved in environments that radically differ from those we currently experience; consequently, traits that were once advantageous may now be "mismatched" and disease causing. At the genetic level, this hypothesis predicts that loci with a history of selection will exhibit "genotype by environment" (GxE) interactions, with different health effects in "ancestral" versus "modern" environments. To identify such loci, we advocate for combining genomic tools in partnership with subsistence-level groups experiencing rapid lifestyle change. In these populations, comparisons of individuals falling on opposite extremes of the "matched" to "mismatched" spectrum are uniquely possible. More broadly, the work we propose will inform our understanding of environmental and genetic risk factors for NCDs across diverse ancestries and cultures.

Early life adversity has sex-dependent effects on survival across the lifespan in rhesus macaques.

Exposure to adversity during early life is linked to lasting detrimental effects on evolutionary fitness across many taxa. However, due to the challenges of collecting longitudinal data, especially in species where one sex disperses, direct evidence from long-lived species remains relatively scarce. Here we test the effects of early life adversity on male and female longevity in a free-ranging population of rhesus macaques (Macaca mulatta) at Cayo Santiago, Puerto Rico. We leveraged six decades of data to quantify the relative importance of ten forms of early life adversity for 6,599 macaques (3,230 male, 3,369 female), with a smaller sample size (N=299) for one form of adversity (maternal social isolation) which required high-resolution behavioral data. We found that individuals who experienced more early life adversity died earlier than those who experienced less adversity. Mortality risk was highest during early life, defined as birth to four years old, suggesting acute survival effects of adversity, but heightened mortality risk was also present in macaques who survived to adulthood. Females and males were affected differently by some forms of adversity, and these differences might be driven by varying energetic demands, female philopatry, and male dispersal. By leveraging data on thousands of macaques collected over decades, our results show that the fitness consequences of early life adversity are not uniform across individuals but vary as a function of the type of adversity, timing, and social context, and thus contribute to our limited but growing understanding of the evolution of early life sensitivities in long-lived species.

Apolipoprotein-ε4 is associated with higher fecundity in a natural fertility population.

In many populations, the apolipoprotein-ε4 (APOE-ε4) allele increases the risk for several chronic diseases of aging, including dementia and cardiovascular disease; despite these harmful effects at later ages, the APOE-ε4 allele remains prevalent. We assess the impact of APOE-ε4 on fertility and its proximate determinants (age at first reproduction, interbirth interval) among the Tsimane, a natural fertility population of forager-horticulturalists. Among 795 women aged 13 to 90 (20% APOE-ε4 carriers), those with at least one APOE-ε4 allele had 0.3 to 0.5 more children than (ε3/ε3) homozygotes, while those with two APOE-ε4 alleles gave birth to 1.4 to 2.1 more children. APOE-ε4 carriers achieve higher fertility by beginning reproduction 0.8 years earlier and having a 0.23-year shorter interbirth interval. Our findings add to a growing body of literature suggesting a need for studies of populations living in ancestrally relevant environments to assess how alleles that are deleterious in sedentary urban environments may have been maintained by selection throughout human evolutionary history.

Does vision therapy for visual information processing improve academic performance? A randomised clinical trial.

PURPOSE: To determine whether a typical vision therapy (VT) programme designed to improve visual information processing (VIP) skills is effective in improving these skills and/or academic performance. METHODS: We used a double-blind, randomised clinical trial to compare VIP VT to placebo training. Participating schools referred a sample of 579 early primary school children identified as being within the lower third of their class for literacy. From the referred sample, we identified 247 children eligible to participate (passed visions and auditory processing screening, and VIP performance <34th percentile), 94 of whom participated. Matching IQ, school grade and sex was achieved by sorting hierarchically on these values and then alternately allocating to VT or placebo groups. Both programmes ran for 10 weeks and consisted of 33 h working at home and 4 h working in office. The VT programme was indicative of that employed in Australian paediatric optometry practices, with the placebo programme containing similar activities, except targeting skills within a child's competencies and with specific VIP development activities removed. The main outcome measures were score change on three standardised educational tests (reading comprehension, spelling and mathematics) and six VIP tests, both immediately post-intervention (PI) and 6 months later. RESULTS: Sixty-nine children completed the programmes. The VT programme produced no significant improvement in the three educational tests or in five of the six VIP tests compared to the control. The VT programme improved visual sequential memory (VSM) by a moderate amount compared to the control (Cohen's d = 0.57 and 0.52, immediately PI and at 6 months, respectively: p < 0.03 and p < 0.02). CONCLUSIONS: The VIP and academic performance benefits from a VT programme were largely identical to those from a control programme, both immediately and 6-month PI. Placebo effects and general effects such as improvements in executive function and/or regression-to-the-mean could be mistaken for specific programme effectiveness.

Best practices for genotype imputation from low-coverage sequencing data in natural populations.

Monitoring genetic diversity in wild populations is a central goal of ecological and evolutionary genetics and is critical for conservation biology. However, genetic studies of nonmodel organisms generally lack access to species-specific genotyping methods (e.g. array-based genotyping) and must instead use sequencing-based approaches. Although costs are decreasing, high-coverage whole-genome sequencing (WGS), which produces the highest confidence genotypes, remains expensive. More economical reduced representation sequencing approaches fail to capture much of the genome, which can hinder downstream inference. Low-coverage WGS combined with imputation using a high-confidence reference panel is a cost-effective alternative, but the accuracy of genotyping using low-coverage WGS and imputation in nonmodel populations is still largely uncharacterized. Here, we empirically tested the accuracy of low-coverage sequencing (0.1-10×) and imputation in two natural populations, one with a large (n = 741) reference panel, rhesus macaques (Macaca mulatta), and one with a smaller (n = 68) reference panel, gelada monkeys (Theropithecus gelada). Using samples sequenced to coverage as low as 0.5×, we could impute genotypes at >95% of the sites in the reference panel with high accuracy (median r2 ≥ 0.92). We show that low-coverage imputed genotypes can reliably calculate genetic relatedness and population structure. Based on these data, we also provide best practices and recommendations for researchers who wish to deploy this approach in other populations, with all code available on GitHub (https://github.com/mwatowich/LoCSI-for-non-model-species). Our results endorse accurate and effective genotype imputation from low-coverage sequencing, enabling the cost-effective generation of population-scale genetic datasets necessary for tackling many pressing challenges of wildlife conservation.

Genetic variants in melanogenesis proteins TYRP1 and TYR are associated with the golden rhesus macaque phenotype.

Nonhuman primates (NHPs) are vital translational research models due to their high genetic, physiological, and anatomical homology with humans. The "golden" rhesus macaque (Macaca mulatta) phenotype is a naturally occurring, inherited trait with a visually distinct pigmentation pattern resulting in light blonde colored fur. Retinal imaging also reveals consistent hypopigmentation and occasional foveal hypoplasia. Here, we describe the use of genome-wide association in 2 distinct NHP populations to identify candidate variants in genes linked to the golden phenotype. Two missense variants were identified in the Tyrosinase-related protein 1 gene (Asp343Gly and Leu415Pro) that segregate with the phenotype. An additional and distinct association was also found with a Tyrosinase variant (His256Gln), indicating the light-colored fur phenotype can result from multiple genetic mechanisms. The implicated genes are related through their contribution to the melanogenesis pathway. Variants in these 2 genes are known to cause pigmentation phenotypes in other species and to be associated with oculocutaneous albinism in humans. The novel associations presented in this study will permit further investigations into the role these proteins and variants play in the melanogenesis pathway and model the effects of genetic hypopigmentation and altered melanogenesis in a naturally occurring nonhuman primate model.

Characterizing the landscape of gene expression variance in humans.

Gene expression variance has been linked to organismal function and fitness but remains a commonly neglected aspect of molecular research. As a result, we lack a comprehensive understanding of the patterns of transcriptional variance across genes, and how this variance is linked to context-specific gene regulation and gene function. Here, we use 57 large publicly available RNA-seq data sets to investigate the landscape of gene expression variance. These studies cover a wide range of tissues and allowed us to assess if there are consistently more or less variable genes across tissues and data sets and what mechanisms drive these patterns. We show that gene expression variance is broadly similar across tissues and studies, indicating that the pattern of transcriptional variance is consistent. We use this similarity to create both global and within-tissue rankings of variation, which we use to show that function, sequence variation, and gene regulatory signatures contribute to gene expression variance. Low-variance genes are associated with fundamental cell processes and have lower levels of genetic polymorphisms, have higher gene-gene connectivity, and tend to be associated with chromatin states associated with transcription. In contrast, high-variance genes are enriched for genes involved in immune response, environmentally responsive genes, immediate early genes, and are associated with higher levels of polymorphisms. These results show that the pattern of transcriptional variance is not noise. Instead, it is a consistent gene trait that seems to be functionally constrained in human populations. Furthermore, this commonly neglected aspect of molecular phenotypic variation harbors important information to understand complex traits and disease.

DNA methylation-environment interactions in the human genome.

Previously we showed that a massively parallel reporter assay, mSTARR-seq, could be used to simultaneously test for both enhancer-like activity and DNA methylation-dependent enhancer activity for millions of loci in a single experiment (Lea et al., 2018). Here we apply mSTARR-seq to query nearly the entire human genome, including almost all CpG sites profiled either on the commonly used Illumina Infinium MethylationEPIC array or via reduced representation bisulfite sequencing. We show that fragments containing these sites are enriched for regulatory capacity, and that methylation-dependent regulatory activity is in turn sensitive to the cellular environment. In particular, regulatory responses to interferon alpha (IFNA) stimulation are strongly attenuated by methyl marks, indicating widespread DNA methylation-environment interactions. In agreement, methylation-dependent responses to IFNA identified via mSTARR-seq predict methylation-dependent transcriptional responses to challenge with influenza virus in human macrophages. Our observations support the idea that pre-existing DNA methylation patterns can influence the response to subsequent environmental exposures-one of the tenets of biological embedding. However, we also find that, on average, sites previously associated with early life adversity are not more likely to functionally influence gene regulation than expected by chance.

Natural Animal Populations as Model Systems for Understanding Early Life Adversity Effects on Aging.

Adverse experiences in early life are associated with aging-related disease risk and mortality across many species. In humans, confounding factors, as well as the difficulty of directly measuring experiences and outcomes from birth till death, make it challenging to identify how early life adversity impacts aging and health. These challenges can be mitigated, in part, through the study of non-human animals, which are exposed to parallel forms of adversity and can age similarly to humans. Furthermore, studying the links between early life adversity and aging in natural populations of non-human animals provides an excellent opportunity to better understand the social and ecological pressures that shaped the evolution of early life sensitivities. Here, we highlight ongoing and future research directions that we believe will most effectively contribute to our understanding of the evolution of early life sensitivities and their repercussions.

DNA methylation signatures of early life adversity are exposure-dependent in wild baboons.

The early life environment can profoundly shape the trajectory of an animal's life, even years or decades later. One mechanism proposed to contribute to these early life effects is DNA methylation. However, the frequency and functional importance of DNA methylation in shaping early life effects on adult outcomes is poorly understood, especially in natural populations. Here, we integrate prospectively collected data on fitness-associated variation in the early environment with DNA methylation estimates at 477,270 CpG sites in 256 wild baboons. We find highly heterogeneous relationships between the early life environment and DNA methylation in adulthood: aspects of the environment linked to resource limitation (e.g., low-quality habitat, early life drought) are associated with many more CpG sites than other types of environmental stressors (e.g., low maternal social status). Sites associated with early resource limitation are enriched in gene bodies and putative enhancers, suggesting they are functionally relevant. Indeed, by deploying a baboon-specific, massively parallel reporter assay, we show that a subset of windows containing these sites are capable of regulatory activity, and that, for 88% of early drought-associated sites in these regulatory windows, enhancer activity is DNA methylation-dependent. Together, our results support the idea that DNA methylation patterns contain a persistent signature of the early life environment. However, they also indicate that not all environmental exposures leave an equivalent mark and suggest that socioenvironmental variation at the time of sampling is more likely to be functionally important. Thus, multiple mechanisms must converge to explain early life effects on fitness-related traits.

Leveraging massively parallel reporter assays for evolutionary questions.

A long-standing goal of evolutionary biology is to decode how gene regulation contributes to organismal diversity. Doing so is challenging because it is hard to predict function from non-coding sequence and to perform molecular research with non-model taxa. Massively parallel reporter assays (MPRAs) enable the testing of thousands to millions of sequences for regulatory activity simultaneously. Here, we discuss the execution, advantages, and limitations of MPRAs, with a focus on evolutionary questions. We propose solutions for extending MPRAs to rare taxa and those with limited genomic resources, and we underscore MPRA's broad potential for driving genome-scale, functional studies across organisms.