Mixed Patterns of Intergenerational DNA Methylation Inheritance in Acropora.

For sessile organisms at high risk from climate change, phenotypic plasticity can be critical to rapid acclimation. Epigenetic markers like DNA methylation are hypothesized as mediators of plasticity; methylation is associated with the regulation of gene expression, can change in response to ecological cues, and is a proposed basis for the inheritance of acquired traits. Within reef-building corals, gene-body methylation (gbM) can change in response to ecological stressors. If coral DNA methylation is transmissible across generations, this could potentially facilitate rapid acclimation to environmental change. We investigated methylation heritability in Acropora, a stony reef-building coral. Two Acropora millepora and two Acropora selago adults were crossed, producing eight offspring crosses (four hybrid, two of each species). We used whole-genome bisulfite sequencing to identify methylated loci and allele-specific alignments to quantify per-locus inheritance. If methylation is heritable, differential methylation (DM) between the parents should equal DM between paired offspring alleles at a given locus. We found a mixture of heritable and nonheritable loci, with heritable portions ranging from 44% to 90% among crosses. gBM was more heritable than intergenic methylation, and most loci had a consistent degree of heritability between crosses (i.e. the deviation between parental and offspring DM were of similar magnitude and direction). Our results provide evidence that coral methylation can be inherited but that heritability is heterogenous throughout the genome. Future investigations into this heterogeneity and its phenotypic implications will be important to understanding the potential capability of intergenerational environmental acclimation in reef building corals.

Should I stay or should I go? Coral bleaching from the symbionts' perspective.

Coral bleaching, the stress-induced breakdown of coral-algal symbiosis, threatens reefs globally. Paradoxically, despite adverse fitness effects, corals bleach annually, even outside of abnormal temperatures. This generally occurs shortly after the once-per-year mass coral spawning. Here, we propose a hypothesis linking annual coral bleaching and the transmission of symbionts to the next generation of coral hosts. We developed a dynamic model with two symbiont growth strategies, and found that high sexual recruitment and low adult coral survivorship and growth favour bleaching susceptibility, while the reverse promotes bleaching resilience. Otherwise, unexplained trends in the Indo-Pacific align with our hypothesis, where reefs and coral taxa exhibiting higher recruitment are more bleaching susceptible. The results from our model caution against interpreting potential shifts towards more bleaching-resistant symbionts as evidence of climate adaptation-we predict such a shift could also occur in declining systems experiencing low recruitment rates, a common scenario on today's reefs.

Decrypting Corals: Does Regulatory Evolution Underlie Environmental Specialisation of Coral Cryptic Lineages?

A recent sequencing study has shown that two common Caribbean corals, Montastraea cavernosa and Siderastrea siderea, each consist of four genetically distinct lineages in the Florida Keys. These lineages are specialised to a certain depth and, to a lesser extent, to nearshore or offshore habitats. We hypothesised that the lineages' environmental specialisation is at least in part due to regulatory evolution, which would manifest as the emergence of groups of coregulated genes ('modules') demonstrating lineage-specific responses to different reef environments. Our hypothesis also predicted that genes belonging to such modules would show greater genetic divergence between lineages than other genes. Contrary to these expectations, responses of cryptic lineages to natural environmental variation were essentially the same at the genome-wide gene coexpression network level, with much fewer differences in gene expression between lineages compared to between habitats. Moreover, none of the identified coregulated gene expression modules exhibit elevated genetic divergence between lineages. Possible explanations of these unexpected results range from the leading role of algal symbionts and/or microbiome in adaptation to strong action of spatially varying selection equalising gene expression patterns despite different genetic background. We discuss how future studies could assist in discriminating between these possibilities.

Evidence of sweepstakes reproductive success in a broadcast-spawning coral and its implications for coral metapopulation persistence.

Processes governing genetic diversity and adaptive potential in reef-building corals are of interest both for fundamental evolutionary biology and for reef conservation. Here, we investigated the possibility of "sweepstakes reproductive success" (SRS) in a broadcast spawning coral, Acropora hyacinthus, at Yap Island, Micronesia. SRS is an extreme yearly variation in the number of surviving offspring among parents. It is predicted to generate genetically differentiated, low-genetic-diversity recruit cohorts, containing close kin individuals. We have tested these predictions by comparing genetic composition of size classes (adults and juveniles) at several sites on the island of Yap. We did see the genome-wide dip in genetic diversity in juveniles compared to adults at two of the four sites; however, both adults and juveniles varied in genetic diversity across sites, and there was no detectable genetic structure among juveniles, which does not conform to the classical SRS scenario. Yet, we have identified a pair of juvenile siblings at the site where juveniles had the lowest genetic diversity compared to adults, an observation that is hard to explain without invoking SRS. While further support for SRS is needed to fully settle the issue, we show that incorporating SRS into the Indo-West Pacific coral metapopulation adaptation model had surprisingly little effect on mean rates of coral cover decline during warming. Still, SRS notably increases year-to-year variation in coral cover throughout the simulation.

Gene expression underlying variation in the survival skills of red drum larvae (Sciaenops ocellatus).

Mortality rates of marine fish larvae are incredibly high and can determine year-class strength. The major causes of larval mortality are predation and starvation, and the performance of larvae in survival skills that can mitigate this mortality (predator evasion, foraging) varies among individuals and cohorts, but the causes of the variation are not known. Transcriptomics can link gene expression variation to phenotypic variation at the whole-system level to investigate the molecular basis of behavioural variation. We used tag-based RNA-sequencing to examine the molecular basis of variation in predator evasion and routine swimming (trait related to foraging efficiency) in the larval red drum, Sciaenops ocellatus. We looked for functional gene networks in which interindividual variation would explain variation in larval behavioural performance. We identified co-expressed gene groups ("modules") associated with predator evasion traits and found enrichment of motor, neural and energy metabolism pathways. These functional associations and pattern of correlations between modules and traits suggest that energy availability and allocation were responsible for the magnitude of startle responses, while differential neural and motor activation were associated with differences in response latency.

Millennia-old coral holobiont DNA provides insight into future adaptive trajectories.

Ancient DNA (aDNA) has been applied to evolutionary questions across a wide variety of taxa. Here, for the first time, we utilized aDNA from millennia-old fossil coral fragments to gain new insights into a rapidly declining western Atlantic reef ecosystem. We sampled four Acropora palmata fragments (dated 4215 BCE to 1099 CE) obtained from two Florida Keys reef cores. From these samples, we established that it is possible both to sequence aDNA from reef cores and place the data in the context of modern-day genetic variation. We recovered varying amounts of nuclear DNA exhibiting the characteristic signatures of aDNA from the A. palmata fragments. To describe the holobiont sensu lato, which plays a crucial role in reef health, we utilized metagenome-assembled genomes as a reference to identify a large additional proportion of ancient microbial DNA from the samples. The samples shared many common microbes with modern-day coral holobionts from the same region, suggesting remarkable holobiont stability over time. Despite efforts, we were unable to recover ancient Symbiodiniaceae reads from the samples. Comparing the ancient A. palmata data to whole-genome sequencing data from living acroporids, we found that while slightly distinct, ancient samples were most closely related to individuals of their own species. Together, these results provide a proof-of-principle showing that it is possible to carry out direct analysis of coral holobiont change over time, which lays a foundation for studying the impacts of environmental stress and evolutionary constraints.

Comparative neurotranscriptomics reveal widespread species differences associated with bonding.

BACKGROUND: Pair bonding with a reproductive partner is rare among mammals but is an important feature of human social behavior. Decades of research on monogamous prairie voles (Microtus ochrogaster), along with comparative studies using the related non-bonding meadow vole (M. pennsylvanicus), have revealed many of the neural and molecular mechanisms necessary for pair-bond formation in that species. However, these studies have largely focused on just a few neuromodulatory systems. To test the hypothesis that neural gene expression differences underlie differential capacities to bond, we performed RNA-sequencing on tissue from three brain regions important for bonding and other social behaviors across bond-forming prairie voles and non-bonding meadow voles. We examined gene expression in the amygdala, hypothalamus, and combined ventral pallidum/nucleus accumbens in virgins and at three time points after mating to understand species differences in gene expression at baseline, in response to mating, and during bond formation. RESULTS: We first identified species and brain region as the factors most strongly associated with gene expression in our samples. Next, we found gene categories related to cell structure, translation, and metabolism that differed in expression across species in virgins, as well as categories associated with cell structure, synaptic and neuroendocrine signaling, and transcription and translation that varied among the focal regions in our study. Additionally, we identified genes that were differentially expressed across species after mating in each of our regions of interest. These include genes involved in regulating transcription, neuron structure, and synaptic plasticity. Finally, we identified modules of co-regulated genes that were strongly correlated with brain region in both species, and modules that were correlated with post-mating time points in prairie voles but not meadow voles. CONCLUSIONS: These results reinforce the importance of pre-mating differences that confer the ability to form pair bonds in prairie voles but not promiscuous species such as meadow voles. Gene ontology analysis supports the hypothesis that pair-bond formation involves transcriptional regulation, and changes in neuronal structure. Together, our results expand knowledge of the genes involved in the pair bonding process and open new avenues of research in the molecular mechanisms of bond formation.

Fantastic Beasts and How To Sequence Them: Ecological Genomics for Obscure Model Organisms.

The application of genomic approaches to 'obscure model organisms' (OMOs), meaning species with no prior genomic resources, enables increasingly sophisticated studies of the genomic basis of evolution, acclimatization, and adaptation in real ecological contexts. I consider here ecological questions that can be addressed using OMOs, and indicate optimal sequencing and data-handling solutions for each case. With this I hope to promote the diversity of OMO-based projects that would capitalize on the peculiarities of the natural history of OMOs and could feasibly be completed within the scope of a single PhD thesis.

How mitonuclear discordance and geographic variation have confounded species boundaries in a widely studied snake.

As DNA sequencing technologies and methods for delimiting species with genomic data become more accessible and numerous, researchers have more tools than ever to investigate questions in systematics and phylogeography. However, easy access to sophisticated computational tools is not without its drawbacks. Choosing the right approach for one's question can be challenging when presented with multitudinous options, some of which fail to distinguish between species and intraspecific population structure. Here, we employ a methodology that emphasizes intensive geographic sampling, particularly at contact zones between populations, with a focus on differentiating intraspecific genetic clusters from species in the Pantherophis guttatus complex, a group of North American ratsnakes. Using a mitochondrial marker as well as ddRADseq data, we find evidence of mitonuclear discordance which has contributed to historical confusion about the relationships within this group. Additionally, we identify geographically and genetically structured populations within the species Pantherophis emoryi that are congruent with previously described morphological variation. Importantly, we find that these structured populations within P. emoryi are highly admixed throughout the range of the species and show no evidence of any reproductive isolation. Our data support a revision of the taxonomy of this group, and we recognize two species within the complex and three subspecies within P. emoryi. This study illustrates the importance of thorough sampling of contact zones and consideration of gene flow when delimiting species in widespread complexes containing parapatric lineages.

Role of gene body methylation in acclimatization and adaptation in a basal metazoan.

Gene body methylation (GBM) has been hypothesized to modulate responses to environmental change, including transgenerational plasticity, but the evidence thus far has been lacking. Here we show that coral fragments reciprocally transplanted between two distant reefs respond predominantly by increase or decrease in genome-wide GBM disparity: The range of methylation levels between lowly and highly methylated genes becomes either wider or narrower. Remarkably, at a broad functional level this simple adjustment correlated very well with gene expression change, reflecting a shifting balance between expressions of environmentally responsive and housekeeping genes. In our experiment, corals in a lower-quality habitat up-regulated genes involved in environmental responses, while corals in a higher-quality habitat invested more in housekeeping genes. Transplanted fragments showing closer GBM match to local corals attained higher fitness characteristics, which supports GBM's role in acclimatization. Fixed differences in GBM between populations did not align with plastic GBM changes and were mostly observed in genes with elevated FST, which suggests that they arose predominantly through genetic divergence. However, we cannot completely rule out transgenerational inheritance of acquired GBM states.

Potential and limits for rapid genetic adaptation to warming in a Great Barrier Reef coral.

Can genetic adaptation in reef-building corals keep pace with the current rate of sea surface warming? Here we combine population genomics, biophysical modeling, and evolutionary simulations to predict future adaptation of the common coral Acropora millepora on the Great Barrier Reef (GBR). Genomics-derived migration rates were high (0.1-1% of immigrants per generation across half the latitudinal range of the GBR) and closely matched the biophysical model of larval dispersal. Both genetic and biophysical models indicated the prevalence of southward migration along the GBR that would facilitate the spread of heat-tolerant alleles to higher latitudes as the climate warms. We developed an individual-based metapopulation model of polygenic adaptation and parameterized it with population sizes and migration rates derived from the genomic analysis. We find that high migration rates do not disrupt local thermal adaptation, and that the resulting standing genetic variation should be sufficient to fuel rapid region-wide adaptation of A. millepora populations to gradual warming over the next 20-50 coral generations (100-250 years). Further adaptation based on novel mutations might also be possible, but this depends on the currently unknown genetic parameters underlying coral thermal tolerance and the rate of warming realized. Despite this capacity for adaptation, our model predicts that coral populations would become increasingly sensitive to random thermal fluctuations such as ENSO cycles or heat waves, which corresponds well with the recent increase in frequency of catastrophic coral bleaching events.

CRISPR/Cas9-mediated genome editing in a reef-building coral.

Reef-building corals are critically important species that are threatened by anthropogenic stresses including climate change. In attempts to understand corals' responses to stress and other aspects of their biology, numerous genomic and transcriptomic studies have been performed, generating a variety of hypotheses about the roles of particular genes and molecular pathways. However, it has not generally been possible to test these hypotheses rigorously because of the lack of genetic tools for corals. Here, we demonstrate efficient genome editing using the CRISPR/Cas9 system in the coral Acropora millepora We targeted the genes encoding fibroblast growth factor 1a (FGF1a), green fluorescent protein (GFP), and red fluorescent protein (RFP). After microinjecting CRISPR/Cas9 ribonucleoprotein complexes into fertilized eggs, we detected induced mutations in the targeted genes using changes in restriction-fragment length, Sanger sequencing, and high-throughput Illumina sequencing. We observed mutations in ∼50% of individuals screened, and the proportions of wild-type and various mutant gene copies in these individuals indicated that mutation induction continued for at least several cell cycles after injection. Although multiple paralogous genes encoding green fluorescent proteins are present in A. millepora, appropriate design of the guide RNA allowed us to induce mutations simultaneously in more than one paralog. Because A. millepora larvae can be induced to settle and begin colony formation in the laboratory, CRISPR/Cas9-based gene editing should allow rigorous tests of gene function in both larval and adult corals.

Cladocopium community divergence in two Acropora coral hosts across multiple spatial scales.

Many broadly-dispersing corals acquire their algal symbionts (Symbiodiniaceae) "horizontally" from their environment upon recruitment. Horizontal transmission could promote coral fitness across diverse environments provided that corals can associate with divergent algae across their range and that these symbionts exhibit reduced dispersal potential. Here we quantified community divergence of Cladocopium algal symbionts in two coral host species (Acropora hyacinthus, Acropora digitifera) across two spatial scales (reefs on the same island, and between islands) across the Micronesian archipelago using microsatellites. We find that both hosts associated with a variety of multilocus genotypes (MLG) within two genetically distinct Cladocopium lineages (C40, C21), confirming that Acropora coral hosts associate with a range of Cladocopium symbionts across this region. Both C40 and C21 included multiple asexual lineages bearing identical MLGs, many of which spanned host species, reef sites within islands, and even different islands. Both C40 and C21 exhibited moderate host specialization and divergence across islands. In addition, within every island, algal symbiont communities were significantly clustered by both host species and reef site, highlighting that coral-associated Cladocopium communities are structured across small spatial scales and within hosts on the same reef. This is in stark contrast to their coral hosts, which never exhibited significant genetic divergence between reefs on the same island. These results support the view that horizontal transmission could improve local fitness for broadly dispersing Acropora coral species.

Estimating the potential for coral adaptation to global warming across the Indo-West Pacific.

The potential of reef-building corals to adapt to increasing sea-surface temperatures is often debated but has rarely been comprehensively modeled on a region-wide scale. We used individual-based simulations to model adaptation to warming in a coral metapopulation comprising 680 reefs and representing the whole of the Central Indo-West Pacific. Encouragingly, some reefs-most notably Vietnam, Japan, Taiwan, New Caledonia and the southern half of the Great Barrier Reef-exhibited high capacity for adaptation and, in our model, maintained coral cover even under a rapid "business-as-usual" warming scenario throughout the modeled period (200 years). Higher resilience of these reefs was observed under all tested parameter settings except the models prohibiting selection and/or migration during warming. At the same time, the majority of reefs in the region tended to collapse within the first 100 years of warming. The adaptive potential (odds of maintaining high coral cover) of a given reef could be predicted based on two metrics: the reef's present-day temperature, and the proportion of recruits immigrating from warmer locations. The latter metric explains the most variation in adaptive potential, and significantly correlates with actual coral cover changes observed throughout the region between the 1970s and the early 2000s. These findings will help prioritize coral conservation efforts and plan assisted gene flow interventions to boost the adaptive potential of specific coral populations.

Positive genetic associations among fitness traits support evolvability of a reef-building coral under multiple stressors.

Climate change threatens organisms in a variety of interactive ways that requires simultaneous adaptation of multiple traits. Predicting evolutionary responses requires an understanding of the potential for interactions among stressors and the genetic variance and covariance among fitness-related traits that may reinforce or constrain an adaptive response. Here we investigate the capacity of Acropora millepora, a reef-building coral, to adapt to multiple environmental stressors: rising sea surface temperature, ocean acidification, and increased prevalence of infectious diseases. We measured growth rates (weight gain), coral color (a proxy for Symbiodiniaceae density), and survival, in addition to nine physiological indicators of coral and algal health in 40 coral genets exposed to each of these three stressors singly and combined. Individual stressors resulted in predicted responses (e.g., corals developed lesions after bacterial challenge and bleached under thermal stress). However, corals did not suffer substantially more when all three stressors were combined. Nor were trade-offs observed between tolerances to different stressors; instead, individuals performing well under one stressor also tended to perform well under every other stressor. An analysis of genetic correlations between traits revealed positive covariances, suggesting that selection to multiple stressors will reinforce rather than constrain the simultaneous evolution of traits related to holobiont health (e.g., weight gain and algal density). These findings support the potential for rapid coral adaptation under climate change and emphasize the importance of accounting for corals' adaptive capacity when predicting the future of coral reefs.

Role of host genetics and heat-tolerant algal symbionts in sustaining populations of the endangered coral Orbicella faveolata in the Florida Keys with ocean warming.

Identifying which factors lead to coral bleaching resistance is a priority given the global decline of coral reefs with ocean warming. During the second year of back-to-back bleaching events in the Florida Keys in 2014 and 2015, we characterized key environmental and biological factors associated with bleaching resilience in the threatened reef-building coral Orbicella faveolata. Ten reefs (five inshore, five offshore, 179 corals total) were sampled during bleaching (September 2015) and recovery (May 2016). Corals were genotyped with 2bRAD and profiled for algal symbiont abundance and type. O. faveolata at the inshore sites, despite higher temperatures, demonstrated significantly higher bleaching resistance and better recovery compared to offshore. The thermotolerant Durusdinium trenchii (formerly Symbiondinium trenchii) was the dominant endosymbiont type region-wide during initial (78.0% of corals sampled) and final (77.2%) sampling; >90% of the nonbleached corals were dominated by D. trenchii. 2bRAD host genotyping found no genetic structure among reefs, but inshore sites showed a high level of clonality. While none of the measured environmental parameters were correlated with bleaching, 71% of variation in bleaching resistance and 73% of variation in the proportion of D. trenchii was attributable to differences between genets, highlighting the leading role of genetics in shaping natural bleaching patterns. Notably, D. trenchii was rarely dominant in O. faveolata from the Florida Keys in previous studies, even during bleaching. The region-wide high abundance of D. trenchii was likely driven by repeated bleaching associated with the two warmest years on record for the Florida Keys (2014 and 2015). On inshore reefs in the Upper Florida Keys, O. faveolata was most abundant, had the highest bleaching resistance, and contained the most corals dominated by D. trenchii, illustrating a causal link between heat tolerance and ecosystem resilience with global change.

Considerations for maximizing the adaptive potential of restored coral populations in the western Atlantic.

Active coral restoration typically involves two interventions: crossing gametes to facilitate sexual larval propagation; and fragmenting, growing, and outplanting adult colonies to enhance asexual propagation. From an evolutionary perspective, the goal of these efforts is to establish self-sustaining, sexually reproducing coral populations that have sufficient genetic and phenotypic variation to adapt to changing environments. Here, we provide concrete guidelines to help restoration practitioners meet this goal for most Caribbean species of interest. To enable the persistence of coral populations exposed to severe selection pressure from many stressors, a mixed provenance strategy is suggested: genetically unique colonies (genets) should be sourced both locally as well as from more distant, environmentally distinct sites. Sourcing three to four genets per reef along environmental gradients should be sufficient to capture a majority of intraspecies genetic diversity. It is best for practitioners to propagate genets with one or more phenotypic traits that are predicted to be valuable in the future, such as low partial mortality, high wound healing rate, high skeletal growth rate, bleaching resilience, infectious disease resilience, and high sexual reproductive output. Some effort should also be reserved for underperforming genets because colonies that grow poorly in nurseries sometimes thrive once returned to the reef and may harbor genetic variants with as yet unrecognized value. Outplants should be clustered in groups of four to six genets to enable successful fertilization upon maturation. Current evidence indicates that translocating genets among distant reefs is unlikely to be problematic from a population genetic perspective but will likely provide substantial adaptive benefits. Similarly, inbreeding depression is not a concern given that current practices only raise first-generation offspring. Thus, proceeding with the proposed management strategies even in the absence of a detailed population genetic analysis of the focal species at sites targeted for restoration is the best course of action. These basic guidelines should help maximize the adaptive potential of reef-building corals facing a rapidly changing environment.

Molecular characterization of larval development from fertilization to metamorphosis in a reef-building coral.

BACKGROUND: Molecular mechanisms underlying coral larval competence, the ability of larvae to respond to settlement cues, determine their dispersal potential and are potential targets of natural selection. Here, we profiled competence, fluorescence and genome-wide gene expression in embryos and larvae of the reef-building coral Acropora millepora daily throughout 12 days post-fertilization. RESULTS: Gene expression associated with competence was positively correlated with transcriptomic response to the natural settlement cue, confirming that mature coral larvae are "primed" for settlement. Rise of competence through development was accompanied by up-regulation of sensory and signal transduction genes such as ion channels, genes involved in neuropeptide signaling, and G-protein coupled receptor (GPCRs). A drug screen targeting components of GPCR signaling pathways confirmed a role in larval settlement behavior and metamorphosis. CONCLUSIONS: These results gives insight into the molecular complexity underlying these transitions and reveals receptors and pathways that, if altered by changing environments, could affect dispersal capabilities of reef-building corals. In addition, this dataset provides a toolkit for asking broad questions about sensory capacity in multicellular animals and the evolution of development.

Contrasting effects of Symbiodinium identity on coral host transcriptional profiles across latitudes.

Reef-building corals can increase their resistance to heat-induced bleaching through adaptation and acclimatization and/or by associating with a more thermo-tolerant strain of algal symbiont (Symbiodinium sp.). Here, we show that these two adaptive pathways interact. We collected Acropora millepora corals from two contrasting thermal environments on the Great Barrier Reef: cooler, mid-latitude Orpheus Island, where all corals hosted a heat-sensitive clade C Symbiodinium, and warmer, low-latitude Wilkie Island, where corals hosted either a clade C or a more thermo-tolerant clade D. Corals were kept in a benign common garden to reveal differences in baseline gene expression, reflecting prior adaptation/long-term acclimatization. Model-based analysis identified gene expression differences between Wilkie and Orpheus corals that were negatively correlated with previously described transcriptome-wide signatures of heat stress, signifying generally elevated thermotolerance of Wilkie corals. Yet, model-free analyses of gene expression revealed that Wilkie corals hosting clade C were distinct from Wilkie corals hosting clade D, whereas Orpheus corals were more variable. Wilkie corals hosting clade C symbionts exhibited unique functional signatures, including downregulation of histone proteins and ion channels and upregulation of chaperones and RNA processing genes, putatively representing constitutive "frontloading" of stress response genes. Furthermore, clade C Symbiodinium exhibited constitutive expression differences between Wilkie and Orpheus, indicative of contrasting life history strategies. Our results demonstrate that hosting alternative Symbiodinium types is associated with different pathways of local adaptation for the coral host. These interactions could play a significant role in setting the direction of genetic adaptation to global warming in the two symbiotic partners.

Complex selection on a regulator of social cognition: Evidence of balancing selection, regulatory interactions and population differentiation in the prairie vole Avpr1a locus.

Adaptive variation in social behaviour depends upon standing genetic variation, but we know little about how evolutionary forces shape genetic diversity relevant to brain and behaviour. In prairie voles (Microtus ochrogaster), variants at the Avpr1a locus predict expression of the vasopressin 1a receptor in the retrosplenial cortex (RSC), a brain region that mediates spatial and contextual memory; cortical V1aR abundance in turn predicts diversity in space use and sexual fidelity in the field. To examine the potential contributions of adaptive and neutral forces to variation at the Avpr1a locus, we explore sequence diversity at the Avpr1a locus and throughout the genome in two populations of wild prairie voles. First, we refine results demonstrating balancing selection at the locus by comparing the frequency spectrum of variants at the locus to a random sample of the genome. Next, we find that the four single nucleotide polymorphisms that predict high V1aR expression in the RSC are in stronger linkage disequilibrium than expected by chance despite high recombination among intervening variants, suggesting that epistatic selection maintains their association. Analysis of population structure and a haplotype network for two populations revealed that this excessive LD was unlikely to be due to admixture alone. Furthermore, the two populations differed considerably in the region shown to be a regulator of V1aR expression despite the extremely low levels of genomewide genetic differentiation. Together, our data suggest that complex selection on Avpr1a locus favours specific combinations of regulatory polymorphisms, maintains the resulting alleles at population-specific frequencies, and may contribute to unique patterns of spatial cognition and sexual fidelity among populations.