Geology ReportsSearch

SEARCH · Geology Reports

Results for “Molecular Ecology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8Linked to original sources

Loss of dendritic connectivity in southern California's urban riverscape facilitates decline of an endemic freshwater fish

Life history adaptations and spatial configuration of metapopulation networks allow certain species to persist in extreme fluctuating environments, yet long-term stability within these systems relies on the maintenance of linkage habitat. Degradation of such linkages in urban riverscapes can disrupt this dynamic in aquatic species, leading to increased extinction debt in local populations experiencing environment-related demographic flux. We used microsatellites and mtDNA to examine the effects of collapsed network structure in the endemic Santa Ana sucker Catostomus santaanae of southern California, a threatened species affected by natural flood-drought cycles, ‘boom-and-bust’ demography, hybridization, and presumed artificial transplantation. Our results show a predominance of drift-mediated processes in shaping population structure, and that reverse mechanisms for counterbalancing the genetic effects of these phenomena have dissipated with the collapse of dendritic connectivity. We use approximate Bayesian models to support two cases of artificial transplantation, and provide evidence that one of the invaded systems better represents the historic processes that maintained genetic variation within watersheds than any remaining drainages where C. santaanae is considered native. We further show that a stable dry gap in the northern range is preventing genetic dilution of pure C. santaanae persisting upstream of a hybrid assemblage involving a non-native sucker, and that local accumulation of genetic variation in the same drainage is influenced by position within the network. This work has important implications for declining species that have historically relied on dendritic metapopulation networks to maintain source-sink dynamics in phasic environments, but no longer possess this capacity in urban-converted landscapes.

California

Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern United States

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in an invasive riparian shrub, saltcedar ( Tamarix L.) by analysing 1,997 genomewide SNP markers for 259 individuals from 25 populations collected throughout the southwestern United States. Our results revealed a broad-scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large-scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers and founder effects during range expansion.

Molecular Ecology

Adaptive population divergence and directional gene flow across steep elevational gradients in a climate‐sensitive mammal

The American pika is a thermally sensitive, alpine lagomorph species. Recent climate-associated population extirpations and genetic signatures of reduced population sizes range-wide indicate the viability of this species is sensitive to climate change. To test for potential adaptive responses to climate stress, we sampled pikas along two elevational gradients (each ~470 to 1640 m) and employed three outlier detection methods, BAYESCAN, LFMM, and BAYPASS, to scan for genotype-environment associations in samples genotyped at 30,763 SNP loci. We resolved 173 loci with robust evidence of natural selection detected by either two independent analyses or replicated in both transects. A BLASTN search of these outlier loci revealed several genes associated with metabolic function and oxygen transport, indicating natural selection from thermal stress and hypoxia. We also found evidence of directional gene flow primarily downslope from large high-elevation populations and reduced gene flow at outlier loci, a pattern suggesting potential impediments to the upward elevational movement of adaptive alleles in response to contemporary climate change. Finally, we documented evidence of reduced genetic diversity associated the south-facing transect and an increase in corticosterone stress levels associated with inbreeding. This study suggests the American pika is already undergoing climate-associated natural selection at multiple genomic regions. Further analysis is needed to determine if the rate of climate adaptation in the American pika and other thermally sensitive species will be able to keep pace with rapidly changing climate conditions.

Washington

Regional variation in drivers of connectivity for two frog species (Rana pretiosa and R. luteiventris) from the U.S. Pacific Northwest

Comparative landscape genetics has uncovered high levels of variability in which landscape factors affect connectivity among species and regions. However, the relative importance of species traits versus environmental variation for predicting landscape patterns of connectivity is unresolved. We provide evidence from a landscape genetics study of two sister taxa of frogs, the Oregon spotted frog ( Rana pretiosa ) and the Columbia spotted frog ( Rana luteiventris ) in Oregon and Idaho, USA. Rana pretiosa is relatively more dependent on moisture for dispersal than R. luteiventris , so if species traits influence connectivity, we predicted that connectivity among R. pretiosa populations would be more positively associated with moisture than R. luteiventris . However, if environmental differences are important drivers of gene flow, we predicted that connectivity would be more positively related to moisture in arid regions. We tested these predictions using eight microsatellite loci and gravity models in two R. pretiosa regions and four R. luteiventris regions ( n = 1,168 frogs). In R. pretiosa , but not R. luteiventris , connectivity was positively related to mean annual precipitation, supporting our first prediction. In contrast, connectivity was not more positively related to moisture in more arid regions. Various temperature metrics were important predictors for both species and in all regions, but the directionality of their effects varied. Therefore, the pattern of variation in drivers of connectivity was consistent with predictions based on species traits rather than on environmental variation.

Molecular Ecology

Parallel signatures of selection at genomic islands of divergence and the major histocompatibility complex in ecotypes of sockeye salmon across Alaska

Understanding the genetic mechanisms that facilitate adaptive radiation is an important component of evolutionary biology. Here, we genotyped 82 neutral SNPs, seven SNPs in islands of divergence identified in a previous study (island SNPs), and a region of the major histocompatibility complex (MHC) in 32 populations of sockeye salmon to investigate whether conserved genes and genomic regions are involved in adaptive radiation. Populations representing three ecotypes were sampled from seven drainages with differing habitats and colonization histories spanning a range of 2,000 km. We found strong signatures of parallel selection across drainages at the island SNPs and MHC, suggesting that the same loci undergo divergent selection during adaptive radiation. However, patterns of differentiation at most island SNPs and the MHC were not associated with ecotypes, suggesting that these loci are responding differently to a mosaic of selective pressures. Our study provides some of the first evidence that conserved genomic islands may be involved in adaptive divergence of salmon populations. Additionally, our data provide further support for the hypothesis that sockeye salmon inhabiting rivers unconnected to lakes harbour similar genetic diversity across large distances, are likely the ancestral form of the species, and have repeatedly recolonized lake systems as they have become available after glacial recession. Finally, our results highlight the value and importance of validating outlier loci by screening additional populations and regions, a practice that will hopefully become more common in the future.

Alaska

Satellite tracking of gulls and genomic characterization of fecal bacteria reveals environmentally mediated acquisition and dispersal of antimicrobial resistant Escherichia coli on the Kenai Peninsula, Alaska

Gulls (Larus spp.) have frequently been reported to carry Escherichia coli exhibiting antimicrobial resistance (AMR E. coli); however, the pathways governing the acquisition and dispersal of such bacteria are not well-described. We equipped 17 landfill-foraging gulls with satellite transmitters and collected gull fecal samples longitudinally from four locations on the Kenai Peninsula, Alaska to assess: 1) gull attendance and transitions between sites, 2) spatiotemporal prevalence of fecally-shed AMR E. coli, and 3) genomic relatedness of AMR E. coli isolates among sites. We also sampled Pacific salmon (Oncorhynchus spp.) harvested as part of personal-use dipnet fisheries at two sites to assess potential contamination with AMR E. coli. Among our study sites, marked gulls most commonly occupied the lower Kenai River (61% of site locations) followed by the Soldotna landfill (11%), lower Kasilof River (5%), and upper Kenai River (<1%). Gulls primarily moved between the Soldotna landfill and the lower Kenai River (94% of transitions among sites), which were also the two locations with the highest prevalence of AMR E. coli. There was relatively high spatial and temporal variability in AMR E. coli prevalence in gull feces and there was no evidence of contamination on salmon harvested in personal-use fisheries. We identified E. coli sequence types and AMR genes of clinical importance, with some isolates possessing genes associated with resistance to as many as eight antibiotic classes. Our findings suggest that gulls acquire AMR E. coli at habitats with anthropogenic inputs and subsequent movements may represent pathways through which AMR is dispersed.

Alaska

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still don’t know how extensive variation in outcomes of hybridization is across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under assumptions of random mating and no selection against hybrids. Since this implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modeled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Wyoming

Urbanization reduces genetic connectivity in bobcats (Lynx rufus) at both intra- and interpopulation spatial scales

Urbanization is a major factor driving habitat fragmentation and connectivity loss in wildlife. However, the impacts of urbanization on connectivity can vary among species and even populations due to differences in local landscape characteristics, and our ability to detect these relationships may depend on the spatial scale at which they are measured. Bobcats ( Lynx rufus ) are relatively sensitive to urbanization and the status of bobcat populations is an important indicator of connectivity in urban coastal southern California. We genotyped 271 bobcats at 13,520 SNP loci to conduct a replicated landscape resistance analysis in five genetically distinct populations. We tested urban and natural factors potentially influencing individual connectivity in each population separately, as well as study–wide. Overall, landscape genomic effects were most frequently detected at the study–wide spatial scale, with urban land cover (measured as impervious surface) having negative effects and topographic roughness having positive effects on gene flow. The negative effect of urban land cover on connectivity was also evident when populations were analyzed separately despite varying substantially in spatial area and the proportion of urban development, confirming a pervasive impact of urbanization largely independent of spatial scale. The effect of urban development was strongest in one population where stream habitat had been lost to development, suggesting that riparian corridors may help mitigate reduced connectivity in urbanizing areas. Our results demonstrate the importance of replicating landscape genetic analyses across populations and considering how landscape genetic effects may vary with spatial scale and local landscape structure.

California

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow ( Corvus caurinus ) and American Crow ( Corvus brachyrhynchos ) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%‐divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late‐generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human‐modified habitats.

Washington

Robust geographical determinants of infection prevalence and a contrasting latitudinal diversity gradient for haemosporidian parasites in Western Palearctic birds

Identifying robust environmental predictors of infection probability is central to forecasting and mitigating the ongoing impacts of climate change on vector‐borne disease threats. We applied phylogenetic hierarchical models to a data set of 2,171 Western Palearctic individual birds from 47 species to determine how climate and landscape variation influence infection probability for three genera of haemosporidian blood parasites ( Haemoproteus , Leucocytozoon , and Plasmodium ). Our comparative models found compelling evidence that birds in areas with higher vegetation density (captured by the normalized difference vegetation index [NDVI]) had higher likelihoods of carrying parasite infection. Magnitudes of this relationship were remarkably similar across parasite genera considering that these parasites use different arthropod vectors and are widely presumed to be epidemiologically distinct. However, we also uncovered key differences among genera that highlighted complexities in their climate responses. In particular, prevalences of Haemoproteus and Plasmodium showed strong but contrasting relationships with winter temperatures, supporting mounting evidence that winter warming is a key environmental filter impacting the dynamics of host‐parasite interactions. Parasite phylogenetic community diversities demonstrated a clear but contrasting latitudinal gradient, with Haemoproteus diversity increasing towards the equator and Leucocytozoon diversity increasing towards the poles. Haemoproteus diversity also increased in regions with higher vegetation density, supporting our evidence that summer vegetation density is important for structuring the distributions of these parasites. Ongoing variation in winter temperatures and vegetation characteristics will probably have far‐reaching consequences for the transmission and spread of vector‐borne diseases.

Molecular Ecology

Does harvest affect genetic diversity in grey wolves?

Harvest can affect vital rates such as reproduction and survival, but also genetic measures of individual and population health. Grey wolves ( Canis lupus ) live and breed in groups, and effective population size is a small fraction of total abundance. As a result, genetic diversity of wolves may be particularly sensitive to harvest. We evaluated how harvest affected genetic diversity and relatedness in wolves. We hypothesized that harvest would (a) reduce relatedness of individuals within groups in a subpopulation but increase relatedness of individuals between groups due to increased local immigration, (b) increase individual heterozygosity and average allelic richness across groups in subpopulations and (c) add new alleles to a subpopulation and decrease the number of private alleles in subpopulations due to an increase in breeding opportunities for unrelated individuals. We found harvest had no effect on observed heterozygosity of individuals or allelic richness at loci within subpopulations but was associated with a small, biologically insignificant effect on within-group relatedness values in grey wolves. Harvest was, however, positively associated with increased relatedness of individuals between groups and a net gain (+16) of alleles into groups in subpopulations monitored since harvest began, although the number of private alleles in subpopulations overall declined. Harvest likely created opportunities for wolves to immigrate into nearby groups and breed, thereby making groups in subpopulations more related over time. Harvest appears to affect genetic diversity in wolves at the group and population levels, but its effects are less apparent at the individual level given the population sizes we studied.

Idaho

Harnessing landscape genomics to identify future climate resilient genotypes in a desert annual

Local adaptation features critically in shaping species responses to changing environments, complicating efforts to revegetate degraded areas. Rapid climate change poses an additional challenge that could reduce fitness of even locally sourced seeds in restoration. Predictive restoration strategies that apply seeds with favourable adaptations to future climate may promote long‐term resilience. Landscape genomics is increasingly used to assess spatial patterns in local adaption and may represent a cost‐efficient approach for identifying future‐adapted genotypes. To demonstrate such an approach, we genotyped 760 plants from 64 Mojave Desert populations of the desert annual Plantago ovata . Genome scans on 5,960 SNPs identified 184 potentially adaptive loci related to climate and satellite vegetation metrics. Causal modelling indicated that variation in potentially adaptive loci was not confounded by isolation by distance or isolation by habitat resistance. A generalized dissimilarity model (GDM) attributed spatial turnover in potentially adaptive loci to temperature, precipitation and NDVI amplitude, a measure of vegetation green‐up potential. By integrating a species distribution model (SDM), we find evidence that summer maximum temperature may both constrain the range of P. ovata and drive adaptive divergence in populations exposed to higher temperatures. Within the species’ current range, warm‐adapted genotypes are predicted to experience a fivefold expansion in climate niche by midcentury and could harbour key adaptations to cope with future climate. We recommend eight seed transfer zones and project each zone into its relative position in future climate. Prioritizing seed collection efforts on genotypes with expanding future habitat represents a promising strategy for restoration practitioners to address rapidly changing climates.

California, Nevada

Nuclear eDNA estimates population allele frequencies and abundance in experimental mesocosms

Advances in environmental DNA (eDNA) methodologies have led to improvements in the ability to detect species and communities in aquatic environments, yet the majority of studies emphasize biological diversity at the species level by targeting variable sites within the mitochondrial genome. Here, we demonstrate that eDNA approaches also have the capacity to detect intraspecific diversity in the nuclear genome, allowing for assessments of population-level genetic diversity and estimates of the number of genetic contributors in a sample. Using a panel of microsatellite loci, we evaluated intraspecific genetic diversity in the round goby (Neogobius melanostomus) using eDNA samples from experimental mesocosms. First, we tested the similarity between eDNA and individual tissue-based estimates of allele frequencies. Subsequently, we used a likelihood-based DNA mixture framework to estimate the number of unique genetic contributors in mesocosm eDNA samples and in simulated mixtures of alleles. Allele frequencies from eDNA accurately reflected allele frequencies from genotyped round goby tissue samples, indicating nuclear markers can be reliably amplified from water samples under controlled conditions. DNA mixture analyses were able to estimate the number of genetic contributors from eDNA samples and simulated mixtures of DNA from up to 58 individuals, with the degree of positive or negative bias dependent on the filtering scheme of low-frequency alleles. This study is the first to document the application of eDNA and multiple amplicon-based methods to obtain intraspecific nuclear genetic information and estimate the absolute abundance of a species in mesocosms. With proper validation, this approach has the potential to advance non-invasive survey methods to characterize populations and broadens the application of eDNA methodologies to inform population-level management objectives.

New York

The role of neutral and adaptive genomic variation in population diversification and speciation in two ground squirrel species of conservation concern

Understanding the neutral (demographic) and adaptive processes leading to the differentiation of species and populations is a critical component of evolutionary and conservation biology. In this context, recently diverged taxa represent a unique opportunity to study the process of genetic differentiation. Northern and southern Idaho ground squirrels ( Urocitellus brunneus —NIDGS, and U . endemicus —SIDGS, respectively) are a recently diverged pair of sister species that have undergone dramatic declines in the last 50 years and are currently found in metapopulations across restricted spatial areas with distinct environmental pressures. Here we genotyped single-nucleotide polymorphisms (SNPs) from buccal swabs with restriction site-associated DNA sequencing (RADseq). With these data we evaluated neutral genetic structure at both the inter- and intraspecific level, and identified putatively adaptive SNPs using population structure outlier detection and genotype–environment association (GEA) analyses. At the interspecific level, we detected a clear separation between NIDGS and SIDGS, and evidence for adaptive differentiation putatively linked to torpor patterns. At the intraspecific level, we found evidence of both neutral and adaptive differentiation. For NIDGS, elevation appears to be the main driver of adaptive differentiation, while neutral variation patterns match and expand information on the low connectivity between some populations identified in previous studies using microsatellite markers. For SIDGS, neutral substructure generally reflected natural geographical barriers, while adaptive variation reflected differences in land cover and temperature, as well as elevation. These results clearly highlight the roles of neutral and adaptive processes for understanding the complexity of the processes leading to species and population differentiation, which can have important conservation implications in susceptible and threatened species.

Idaho

Gene flow influences the genomic architecture of local adaptation in six riverine fish species

Understanding how gene flow influences adaptive divergence is important for predicting adaptive responses. Theoretical studies suggest that when gene flow is high, clustering of adaptive genes in fewer genomic regions would protect adaptive alleles from recombination and thus be selected for, but few studies have tested it with empirical data. Here, we used restriction site-associated sequencing to generate genomic data for six fish species with contrasting life histories from six reaches of the Upper Mississippi River System, USA. We used four differentiation-based outlier tests and three genotype–environment association analyses to define neutral single nucleotide polymorphisms (SNPs) and outlier SNPs that were putatively under selection. We then examined the distribution of outlier SNPs along the genome and investigated whether these SNPs were found in genomic islands of differentiation and inversions. We found that gene flow varied among species, and outlier SNPs were clustered more tightly in species with higher gene flow. The two species with the highest overall F ST (0.0303–0.0720) and therefore lowest gene flow showed little evidence of clusters of outlier SNPs, with outlier SNPs in these species spreading uniformly across the genome. In contrast, nearly all outlier SNPs in the species with the lowest F ST (0.0003) were found in a single large putative inversion. Two other species with intermediate gene flow ( F ST ~ 0.0025–0.0050) also showed clustered genomic architectures, with most islands of differentiation clustered on a few chromosomes. Our results provide important empirical evidence to support the hypothesis that increasingly clustered architecture of local adaptation is associated with high gene flow.

Molecular Ecology

Scale-dependent influence of the sagebrush community on genetic connectivity of the sagebrush obligate Gunnison sage-grouse

Habitat fragmentation and degradation impacts an organism's ability to navigate the landscape, ultimately resulting in decreased gene flow and increased extinction risk. Understanding how landscape composition impacts gene flow (i.e., connectivity) and interacts with scale is essential to conservation decision-making. We used a landscape genetics approach implementing a recently developed statistical model based on the generalized Wishart probability distribution to identify the primary landscape features affecting gene flow and estimate the degree to which each component influences connectivity for Gunnison sage-grouse ( Centrocercus minimus ). We were interested in two spatial scales: among distinct populations rangewide and among leks (i.e., breeding grounds) within the largest population, Gunnison Basin. Populations and leks are nested within a landscape fragmented by rough terrain and anthropogenic features, although requisite sagebrush habitat is more contiguous within populations. Our best fit models for each scale confirm the importance of sagebrush habitat in connectivity, although the important sagebrush characteristics differ. For Gunnison Basin, taller shrubs and higher quality nesting habitat were the primary drivers of connectivity, while more sagebrush cover and less conifer cover facilitated connectivity rangewide. Our findings support previous assumptions that Gunnison sage-grouse range contraction is largely the result of habitat loss and degradation. Importantly, we report direct estimates of resistance for landscape components that can be used to create resistance surfaces for prioritization of specific locations for conservation or management (i.e., habitat preservation, restoration, or development) or as we demonstrated, can be combined with simulation techniques to predict impacts to connectivity from potential management actions.

Colorado, Utah

Hybrid enrichment of adaptive variation revealed by genotype-environment associations in montane sedges

The role of hybridization in diversification is complex and may result in many possible outcomes. Not only can hybridization produce new lineages, but those lineages may contain unique combinations of adaptive genetic variation derived from parental taxa that allow hybrid-origin lineages to occupy unique environmental space relative to one (or both) parents. We document such a case of hybridization between two sedge species, Carex nova and Carex nelsonii (Cyperaceae), that occupy partially overlapping environmental space in the southern Rocky Mountains, USA. In the region hypothesized to be the origin of the hybrid lineage, one parental taxon ( C. nelsonii ) is at the edge of its environmental tolerance. Hybrid-origin individuals display mixed ancestry between the parental taxa – of nearly 7,000 unlinked loci sampled, almost 30% showed evidence of excess ancestry from one parental lineage – approximately half displayed a genomic background skewed towards one parent, and half skewed towards the other. To test whether excess ancestry loci may have conferred an adaptive advantage to the hybrid-origin lineage, we conducted genotype-environment association analyses on different combinations of loci – with and without excess ancestry – and with multiple contrasts between the hybrids and parental taxa. Loci with skewed ancestry showed significant environmental associations distinguishing the hybrid lineage from one parent ( C. nelsonii ), whereas loci with relatively equal representation of parental ancestries showed no such environmental associations. Moreover, the overwhelming majority of candidate adaptive loci with respect to environmental gradients also had excess ancestry from a parental lineage, implying these loci have facilitated the persistence of the hybrid lineage in an environment unsuitable to at least one parent .

Arizona, Colorado, Idaho, Kansas, Montana, Nebrask

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 leverage 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 - 1099 CE) obtained from two Florida Keys reef cores. From these samples, we established that it is possible both to sequence ancient DNA 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.

Molecular Ecology