USGS ScienceSearch

SEARCH · USGS Science

Results for “Genome”

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 55 records · Page 3Linked to original sources

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

Comparative genomics analyses support the reclassification of Bisgaard taxon 40 as Mergibacter gen. nov., with Mergibacter septicus sp. nov. as type species: Novel insights into the phylogeny and virulence factors of a Pasteurellaceae family member associated with mortality events in seabirds

The Pasteurellaceae family has been associated with fatal diseases in numerous avian species. Several new taxa within this family, including Bisgaard taxon 40, have been recently described in wild birds, but their genomic characteristics and pathogenicity are not well understood. We isolated Bisgaard taxon 40 from four species of seabirds, including one sampled during a mass, multi-species mortality event in Florida, United States. Here, we present a comprehensive phenotypic and genetic characterization of Bisgaard taxon 40 and comparative genomic analysis with reference strains from the Pasteurellaceae family, aiming at determining its phylogenetic position, antimicrobial susceptibility profile, and identifying putative virulence factors. In silico multilocus sequence-based and whole-genome-based phylogenetic analysis clustered all Bisgaard taxon 40 strains together on a distinct branch separated from the other members of the Pasteurellaceae family, indicating that Bisgaard taxon 40 could represent a new genus. These findings were further supported by protein similarity analyses using the concatenation of 31 conserved proteins and other taxonomic approaches such as the percentage of conserved protein test. Additionally, several putative virulence factors were identified, including those associated with adhesion (capsule, ompA , ompH ) and colonization ( exbD , fur , galU , galE , lpxA , lpxC , and kdsA ) of the host and a cytolethal distending toxin ( cdt ), which may have played a role in disease development leading to the mortality event. Considerably low minimum inhibitory concentrations (MICs) were found for all the drugs tested, in concordance with the absence of antimicrobial resistance genes in these genomes. The novel findings of this study highlight genomic and phenotypic characteristics of this bacterium, providing insights into genome evolution and pathogenicity. We propose a reclassification of these organisms within the Pasteurellaceae family, designated as Mergibacter gen. nov., with Mergibacter septicus sp. nov. as the type species. The type strain is Mergibacter septicus A25201 T (=DSM 112696).

Frontiers in Microbiology

Genomic evolution, recombination, and inter-strain diversity of chelonid alphaherpesvirus 5 from Florida and Hawaii green sea turtles with fibropapillomatosis

Chelonid alphaherpesvirus 5 (ChHV5) is a herpesvirus associated with fibropapillomatosis (FP) in sea turtles worldwide. Single-locus typing has previously shown differentiation between Atlantic and Pacific strains of this virus, with low variation within each geographic clade. However, a lack of multi-locus genomic sequence data hinders understanding of the rate and mechanisms of ChHV5 evolutionary divergence, as well as how these genomic changes may contribute to differences in disease manifestation. To assess genomic variation in ChHV5 among five Hawaii and three Florida green sea turtles, we used high-throughput short-read sequencing of long-range PCR products amplified from tumor tissue using primers designed from the single available ChHV5 reference genome from a Hawaii green sea turtle. This strategy recovered sequence data from both geographic regions for approximately 75% of the predicted ChHV5 coding sequences. The average nucleotide divergence between geographic populations was 1.5%; most of the substitutions were fixed differences between regions. Protein divergence was generally low (average 0.08%), and ranged between 0 and 5.3%. Several atypical genes originally identified and annotated in the reference genome were confirmed in ChHV5 genomes from both geographic locations. Unambiguous recombination events between geographic regions were identified, and clustering of private alleles suggests the prevalence of recombination in the evolutionary history of ChHV5. This study significantly increased the amount of sequence data available from ChHV5 strains, enabling informed selection of loci for future population genetic and natural history studies, and suggesting the (possibly latent) co-infection of individuals by well-differentiated geographic variants.

Florida, Hawaii

High-quality, chromosome-level reference genomes of the viviparous Caribbean skinks Spondylurus nitidus and S. culebrae

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae . Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

Genome Biology and Evolution

The roles of antimicrobial resistance, phage diversity, isolation source, and selection in shaping the genomic architecture of Bacillus anthracis

Bacillus anthracis, the causative agent of anthrax disease, is a worldwide threat to livestock, wildlife and public health. While analyses of genetic data from across the globe have increased our understanding of this bacterium’s population genomic structure, the influence of selective pressures on this successful pathogen is not well understood. In this study, we investigate the effects of antimicrobial resistance, phage diversity, geography and isolation source in shaping population genomic structure. We also identify a suite of candidate genes potentially under selection, driving patterns of diversity across 356 globally extant B. anthracis genomes. We report ten antimicrobial resistance genes and 11 different prophage sequences, resulting in the first large-scale documentation of these genetic anomalies for this pathogen. Results of random forest classification suggest genomic structure may be driven by a combination of antimicrobial resistance, geography and isolation source, specific to the population cluster examined. We found strong evidence that a recombination event linked to a gene involved in protein synthesis may be responsible for phenotypic differences between comparatively disparate populations. We also offer a list of genes for further examination of B. anthracis evolution, based on high-impact single nucleotide polymorphisms (SNPs) and clustered mutations. The information presented here sheds new light on the factors driving genomic structure in this notorious pathogen and may act as a road map for future studies aimed at understanding functional differences in terms of B. anthracis biogeography, virulence and evolution.

Microbial Genomics

Comparative genomic analyses and a novel linkage map for cisco (Coregonus artedi) provides insight into chromosomal evolution and rediploidization across salmonids

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco ( Coregonus artedi ), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

G3 Genes, Genomes, Genetics

Improving conservation policy with genomics: A guide to integrating adaptive potential into U.S. Endangered Species Act decisions for conservation practitioners and geneticists

Rapid environmental change makes adaptive potential—the capacity of populations to evolve genetically based changes in response to selection—more important than ever for long-term persistence of at-risk species. At the same time, advances in genomics provide unprecedented power to test for and quantify adaptive potential, enabling consideration of adaptive potential in estimates of extinction risk and laws protecting endangered species. The U.S. Endangered Species Act (ESA) is one of the most powerful environmental laws in the world, but so far, the full potential of genomics in ESA listing and recovery decisions has not been realized by the federal agencies responsible for implementing the ESA or by conservation geneticists. The goal of our paper is to chart a path forward for integrating genomics into ESA decision making to facilitate full consideration of adaptive potential in evaluating long-term risk of extinction. For policy makers, managers, and other conservation practitioners, we outline why adaptive potential is important for population persistence and what genomic tools are available for quantifying it. For conservation geneticists, we discuss how federal agencies can integrate information on the effect of adaptive potential on extinction risk—and the related uncertainty—into decisions, and suggest next steps for advancing understanding of the effect of adaptive potential on extinction risk. The mechanisms and consequences of adaptation are incredibly complex, and we may never have a complete understanding of adaptive potential for any organism. Nevertheless, we argue that the best available evidence regarding adaptive potential should be incorporated by federal agencies into modeling and decision making processes now, while at the same time conserving genome-wide variation and striving for a deeper understanding of adaptive potential and its effects on population persistence to improve decision-making into the future.

Conservation Genetics

Interspecies transmission and limited persistence of low pathogenic avian influenza genomes among Alaska dabbling ducks

The reassortment and geographic distribution of low pathogenic avian influenza (LPAI) virus genes are well documented, but little is known about the persistence of intact LPAI genomes among species and locations. To examine persistence of entire LPAI genome constellations in Alaska, we calculated the genetic identities among 161 full-genome LPAI viruses isolated across 4 years from five species of duck: northern pintail ( Anas acuta ), mallard ( Anas platyrhynchos ), American green-winged teal ( Anas crecca ), northern shoveler ( Anas clypeata ) and American wigeon ( Anas americana ). Based on pairwise genetic distance, highly similar LPAI genomes (>99% identity) were observed within and between species and across a range of geographic distances (up to and >1000 km), but most often between isolates collected 0–10 km apart. Highly similar viruses were detected between years, suggesting inter-annual persistence, but these were rare in our data set with the majority occurring within 0–9 days of sampling. These results identify LPAI transmission pathways in the context of species, space and time, an initial perspective into the extent of regional virus distribution and persistence, and insight into why no completely Eurasian genomes have ever been detected in Alaska. Such information will be useful in forecasting the movement of foreign-origin avian influenza strains should they be introduced to North America.

Alaska

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5–10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4–5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

Proceedings of the National Academy of Sciences

Reference genome of an iconic lizard in western North America, Blainville’s horned lizard Phrynosoma blainvillii

Genome assemblies are increasingly being used to identify adaptive genetic variation that can help prioritize the population management of protected species. This approach may be particularly relevant to species like Blainville’s horned lizard, Phrynosoma blainvillii , due to its specialized diet on noxious harvester ants, numerous adaptative traits for avoiding predation (e.g. cranial horns, dorsoventrally compressed body, cryptic coloration, and blood squirting from the orbital sinuses), and status as Species of Special Concern in California. Rangewide decline since the early 20th century, the basis of its conservation status, has been driven mainly by habitat conversion, over-collecting, and invasion of a non-native ant that displaces its native ant prey base. Here, we report on a scaffold-level genome assembly for P. blainvillii as part of the California Conservation Genomics Project (CCGP), produced using Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology. The de novo assembly has 78 scaffolds, a total length of ~2.21 Gb, a scaffold N50 length of ~352 Mb, and BUSCO score of 97.4%. This is the second species of Phrynosoma for which a reference genome has been assembled and represents a considerable improvement in terms of contiguity and completeness. Combined with the landscape genomics data being compiled by the CCGP, this assembly will help strategize efforts to maintain and/or restore local genetic diversity, where interventions like genetic rescue, translocation, and strategic land preservation may be the only means by which P. blainvillii and other low-vagility species can survive in the fragmented habitats of California.

California

Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis

Spiny lizards (genus Sceloporus ) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis , is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis , as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis , the species status of the California endemic island fence lizard ( S. becki ), and the spectacular radiation of Sceloporus lizards.

California

Development of high-throughput genomic resources to inform white-tailed deer population and disease management

White-tailed deer ( Odocoileus virginianus ) are the most abundant and widespread cervid in North America. Genetic data are used as a tool to monitor populations and make management decisions for this game species. However, the development and use of genomic tools that can generate a set of markers suitable for longitudinal genomic data collection, whether for management purposes or to study the demographic and evolutionary processes of widely distributed species, have been challenging. This is mainly due to the cost required to fully implement and interpret the data produced. Here, we generated whole genome resequencing data for 44 free-ranging deer from three regions in their central and eastern North American range and identified over 89 million single nucleotide polymorphisms (SNPs). We used a subset of these SNPs to develop two nested SNP tools, a high-density array (702,183 SNPs) and a medium-density array (72,723 SNPs) to support deer and chronic wasting disease (CWD) management and research. SNPs were selected to ensure an even distribution across scaffolds of the reference genome and include SNPs associated with CWD susceptibility. Using genotyping results for 469 deer from 15 states in the US and Mexico generated by the high-density array and 1335 deer from 18 states generated by the medium-density array, we assessed genotyping success across different populations and explored some insights into population structure. These genomic tools offer a standard set of markers that will enable researchers and managers to address important questions related to white-tailed deer and CWD management. Our SNP arrays also offer the opportunity to examine aspects of white-tailed deer ecology and evolutionary history that were previously difficult to address.

Molecular Ecology Resources

Comparison of anadromous and landlocked Atlantic salmon genomes reveals signatures of parallel and relaxed selection across the northern hemisphere

Most Atlantic salmon ( Salmo salar L.) populations follow an anadromous life cycle, spending early life in freshwater, migrating to the sea for feeding and returning to rivers to spawn. At the end of the last ice age ~10,000 years ago, several populations of Atlantic salmon became landlocked. Comparing their genomes to their anadromous counterparts can help identify genetic variation related to either freshwater residency or anadromy. The objective of this study was to identify consistently divergent loci between anadromous and landlocked Atlantic salmon strains throughout their geographical distribution, with the long‐term aim of identifying traits relevant for salmon aquaculture, including fresh and seawater growth, omega‐3 metabolism, smoltification and disease resistance. We used a Pool‐seq approach (n=10‐40 individuals per population) to sequence the genomes of twelve anadromous and six landlocked Atlantic salmon populations covering a large part of the northern hemisphere and conducted a genome‐wide association study to identify genomic regions having been under different selection pressure in landlocked and anadromous strains. A total of 28 genomic regions were identified, and included cadm1 on Chr 13, and ppargc1a on Chr 18. Seven of the regions additionally displayed consistently reduced heterozygosity in fish obtained from landlocked populations, including the genes gpr132 , cdca4 and sertad2 on Chr 15. We also found 16 regions, including igf1 on Chr 17, which consistently display reduced heterozygosity in the anadromous populations compared to the freshwater populations, indicating relaxed selection on traits associated with anadromy in landlocked salmon. In conclusion, we have identified 37 regions which may harbor genetic variation relevant for improving fish welfare and quality in the salmon farming industry and for understanding life history traits in fish.

Evolutionary Applications

Using landscape genomics to delineate future adaptive potential for climate change in the Yosemite toad (Anaxyrus canorus)

An essential goal in conservation biology is delineating population units that maximize the probability of species persisting into the future and adapting to future environmental change. However, future-facing conservation concerns are often addressed using retrospective patterns that could be irrelevant. We recommend a novel landscape genomics framework for delineating future “Geminate Evolutionary Units” (GEUs) in a focal species: (1) identify loci under environmental selection, (2) model and map adaptive conservation units that may spawn future lineages, (3) forecast relative selection pressures on each future lineage, and (4) estimate their fitness and likelihood of persistence using geo-genomic simulations. Using this process, we delineated conservation units for the Yosemite toad ( Anaxyrus canorus ), a U.S. federally threatened species that is highly vulnerable to climate change. We used a genome-wide dataset, redundancy analysis, and Bayesian association methods to identify 24 candidate loci responding to climatic selection ( R 2 ranging from 0.09 to 0.52), after controlling for demographic structure. Candidate loci included genes such as MAP3K5, involved in cellular response to environmental change. We then forecasted future genomic response to climate change using the multivariate machine learning algorithm Gradient Forests. Based on all available evidence, we found three GEUs in Yosemite National Park, reflecting contrasting adaptive optima: YF-North (high winter snowpack with moderate summer rainfall), YF-East (low to moderate snowpack with high summer rainfall), and YF-Low-Elevation (low snowpack and rainfall). Simulations under the RCP 8.5 climate change scenario suggest that the species will decline by 29% over 90 years, but the highly diverse YF-East lineage will be least impacted for two reasons: (1) geographically it will be sheltered from the largest climatic selection pressures, and (2) its standing genetic diversity will promote a faster adaptive response. Our approach provides a comprehensive strategy for protecting imperiled non-model species with genomic data alone and has wide applicability to other declining species.

California

Whole‐genome resequencing reveals persistence of forest‐associated mammals in Late Pleistocene refugia along North America’s North Pacific Coast

Aim Numerous glacial refugia have been hypothesized along North America's North Pacific Coast that may have increased divergence of refugial taxa, leading to elevated endemism and subsequently clustered hybrid zones following deglaciation. The locations and community composition of these ice‐free areas remains controversial, but whole‐genome sequences now enable detailed analysis of the demographic and evolutionary histories of refugial taxa. Here, we use genomic data to test spatial and temporal processes of diversification among martens with respect to the Coastal Refugium Hypothesis, to understand the role of climate cycling in shaping diversity across complex landscapes. Location North America and North Pacific Coast archipelagos. Taxon North American martens ( Martes ). Methods Short‐read whole‐genome resequencing data were generated for 11 martens: four M. americana , four M. caurina , two hybrids, and one outgroup ( Martes zibellina ). Sampling was representative of known genetic clades within New World martens, including sampling within insular and continental hybrid zones and along the North Pacific Coast (five island populations). ADMIXTURE , F‐statistics, and D ‐statistics (ABBA‐BABA) were used to identify introgression and infer directionality. Heterozygosity densities, estimated via PSMC, were used to characterize historical demography at and below the species level to infer refugial and colonization processes. Results Forest‐associated Pacific martens ( M. caurina ) are divided into distinct insular and continental clades consistent with the Coastal Refugium Hypothesis. There was no evidence of introgression on islands that received historical translocations of American pine martens ( M. americana ), but introgression was detected in two active zones of secondary contact: one insular and one continental. Only early‐generational hybrids were identified across multiple hybrid zones, a pattern consistent with potential genetic swamping of M. caurina by M. americana . Main conclusions Despite an incomplete fossil record, genomic evidence supports the persistence of forest‐associated martens, likely the insular Pacific marten lineage, along the western edges of the Alexander Archipelago during the Last Glacial Maximum. This discovery informs our understanding of refugial paleoenvironments, critical to interpreting refugial timing, duration, and community composition. Genomic reevaluations of other taxa along North America's North Pacific Coast may yield new and deeper perspectives on the history of refugial forest communities and the role of dynamic climate shifts in shaping high‐latitude diversity across complex insular landscapes.

Alaska, British Columbia, California, Oregon, Wash

Combining ecological and genomic diversity surveys to inform conservation and restoration of an endangered wetland plant, soft salty bird’s-beak (Chloropyron molle ssp. molle)

Emergent tidal wetlands are declining globally as a result of sea level rise and land use change. This habitat loss can keenly affect rare plant species within wetlands, and may require restoration to meet species recovery goals related to retaining populations throughout species' ranges. Soft salty bird’s-beak ( Chloropyron molle ssp. molle ) is a federally- and state-endangered hemi-parasitic plant that occurs at the upper marsh transition zone in the San Francisco Bay–Delta, California, USA. We combined field surveys to document habitat associations and trends in abundance with genomic surveys to understand patterns of genetic structure in this rare endemic. We found that C . molle ssp. molle persisted at nine previously occupied marsh sites, although four sites (Hill Slough, MOTCO East, Fagan Marsh, and Joice Island) were smaller in population size than when surveyed in the 1990s. Additionally, twelve sites contained plots with suitable but unoccupied habitat that could be further assessed for restoration. Genomic analysis of over 40,000 single-nucleotide polymorphisms (SNPs) and 253 individuals grouped C . molle ssp. molle into six to seven regional genetic clusters with isolation by distance, and confirmed that C . molle ssp. molle is genetically distinct from adjacent populations of its closest relative ( C . molle ssp. hispidum ). The western-most C . molle ssp. molle sites of Point Pinole and Fagan Marsh were the most genetically and geographically isolated and had the lowest genome-wide diversity. Heterozygosity in sets of genes associated with tidal elevation, salinity, and annual and summer precipitation varied independently across populations. Overall, these genomic patterns indicate that selecting donor sites with similar environmental conditions and utilizing composite seeding approaches from multiple sites could allow for local adaptation to a range of possible environmental conditions. This comprehensive survey of habitat and genomic patterns can allow for the development of restoration actions and build climate-adaptation planning to help prevent the loss of a rare plant.

California

Discovery and genomic characterization of a novel hepadnavirus from asymptomatic anadromous alewife (Alosa pseudoharengus)

The alewife ( Alosa pseudoharengus) is an anadromous herring that inhabits waters of northeastern North America. This prey species is a critical forage for piscivorous birds, mammals, and fishes in estuarine and oceanic ecosystems. During a discovery project tailored to identify potentially emerging pathogens of this species, we obtained the full genome of a novel hepadnavirus (ApHBV) from clinically normal alewives collected from the Maurice River, Great Egg Harbor River, and Delaware River in New Jersey, USA during 2015–2018. This previously undescribed hepadnavirus contained a circular DNA genome of 3146 nucleotides. Phylogenetic analysis of the polymerase protein placed this virus in the clade of metahepadnaviruses (family: Hepadnaviridae ; genus: Metahepadnavirus ). There was no evidence of pathology in the internal organs of infected fish and virions were not observed in liver tissues by electron microscopy. We developed a Taqman-based quantitative (qPCR) assay and screened 182 individuals collected between 2015 and 2018 and detected additional qPCR positives (n = 6). An additional complete genome was obtained in 2018 and it has 99.4% genome nucleotide identity to the first virus. Single-nucleotide polymorphisms were observed between the two genomes, including 7/9 and 12/8 synonymous vs nonsynonymous mutations across the polymerase and surface proteins, respectively. While there was no evidence that this virus was associated with disease in this species, alewives are migratory interjurisdictional fishes of management concern. Identification of microbial agents using de novo sequencing and other advanced technologies is a critical aspect of understanding disease ecology for informed population management.

New Jersey

Genome resequencing clarifies phylogeny and reveals patterns of selection in the toxicogenomics model Pimephales promelas

Background The fathead minnow ( Pimephales promelas ) is a model species for toxicological research. A high-quality genome reference sequence is available, and genomic methods are increasingly used in toxicological studies of the species. However, phylogenetic relationships within the genus remain incompletely known and little population-genomic data are available for fathead minnow despite the potential effects of genetic background on toxicological responses. On the other hand, a wealth of extant samples is stored in museum collections that in principle allow fine-scale analysis of contemporary and historical genetic variation. Methods Here we use short-read shotgun resequencing to investigate sequence variation among and within Pimephales species. At the genus level, our objectives were to resolve phylogenetic relationships and identify genes with signatures of positive diversifying selection. At the species level, our objective was to evaluate the utility of archived-sample resequencing for detecting selective sweeps within fathead minnow, applied to a population introduced to the San Juan River of the southwestern United States sometime prior to 1950. Results We recovered well-supported but discordant phylogenetic topologies for nuclear and mitochondrial sequences that we hypothesize arose from mitochondrial transfer among species. The nuclear tree supported bluntnose minnow ( P. notatus ) as sister to fathead minnow, with the slim minnow ( P. tenellus ) and bullhead minnow ( P. vigilax ) more closely related to each other. Using multiple methods, we identified 11 genes that have diversified under positive selection within the genus. Within the San Juan River population, we identified selective-sweep regions overlapping several sets of related genes, including both genes that encode the giant sarcomere protein titin and the two genes encoding the MTORC1 complex, a key metabolic regulator. We also observed elevated polymorphism and reduced differentation among populations (F ST ) in genomic regions containing certain immune-gene clusters, similar to what has been reported in other taxa. Collectively, our data clarify evolutionary relationships and selective pressures within the genus and establish museum archives as a fruitful resource for characterizing genomic variation. We anticipate that large-scale resequencing will enable the detection of genetic variants associated with environmental toxicants such as heavy metals, high salinity, estrogens, and agrichemicals, which could be exploited as efficient biomarkers of exposure in natural populations.

PeerJ