Geology ReportsSearch

SEARCH · Geology Reports

Results for “Genomics”

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 19 recordsLinked to original sources

Rapid microsatellite identification from Illumina paired-end genomic sequencing in two birds and a snake

Identification of microsatellites, or simple sequence repeats (SSRs), can be a time-consuming and costly investment requiring enrichment, cloning, and sequencing of candidate loci. Recently, however, high throughput sequencing (with or without prior enrichment for specific SSR loci) has been utilized to identify SSR loci. The direct "Seq-to-SSR" approach has an advantage over enrichment-based strategies in that it does not require a priori selection of particular motifs, or prior knowledge of genomic SSR content. It has been more expensive per SSR locus recovered, however, particularly for genomes with few SSR loci, such as bird genomes. The longer but relatively more expensive 454 reads have been preferred over less expensive Illumina reads. Here, we use Illumina paired-end sequence data to identify potentially amplifiable SSR loci (PALs) from a snake (the Burmese python, Python molurus bivittatus ), and directly compare these results to those from 454 data. We also compare the python results to results from Illumina sequencing of two bird genomes (Gunnison Sage-grouse, Centrocercus minimus , and Clark's Nutcracker, Nucifraga columbiana ), which have considerably fewer SSRs than the python. We show that direct Illumina Seq-to-SSR can identify and characterize thousands of potentially amplifiable SSR loci for as little as $10 per sample – a fraction of the cost of 454 sequencing. Given that Illumina Seq-to-SSR is effective, inexpensive, and reliable even for species such as birds that have few SSR loci, it seems that there are now few situations for which prior hybridization is justifiable.

PLoS ONE

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (~3.5–5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus , and the Clark's Nutcracker, Nucifraga columbiana ). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.

PLoS ONE

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads ( Anaxyrus [Bufo] boreas ) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis ( Bd ). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity ( π = 0.00034–0.00040) and effective population sizes ( N e = 8.9–38.4) were low, likely due to post-Pleistocene founder effects and Bd -related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd .

Colorado

Genomic and epidemiological investigations reveal chromosomal integration of the acipenserid herpesvirus 3 genome in Lake Sturgeon Acipenser fulvescens

DNA sequence from a new alloherpesvirus named acipenserid herpesvirus 3 (AciHV-3) was found in sturgeon species that are vulnerable to decline globally. A study was undertaken to develop a better understanding of the virus genome and to develop diagnostic tools to support an epidemiological investigation. A 184,426 bp genome was assembled from PacBio HiFi sequences generated with DNA from a Lake Sturgeon Acipenser fulvescens gonad cell line. The AciHV-3 genome was contiguous with host chromosomal DNA and was structured with telomere-like terminal direct repeat regions, five internal direct repeat regions and a U region that included intact open reading frames encoding alloherpesvirus core proteins. Diagnostic testing conducted with a newly developed and analytically validated qPCR assay established the ubiquitous presence and high titer of AciHV-3 DNA in somatic and germline tissues from wild Lake Sturgeon in the Hudson Bay drainage basin. Phylogenetic reconstructions confirm that the monophyletic AciHV-3 lineage shares a common ancestor with AciHV-1 and that AciHV-3 taxa cluster according to their sturgeon host. The same genotype of AciHV-3 is found in disjunctive Lake Sturgeon populations within and among drainage basins. The results support the hypotheses that AciHV-3 has established latency through germline chromosomal integration, is vertically transmitted via a Mendelian pattern of inheritance, is evolving in a manner consistent with a replication competent virus and has co-evolved with its host reaching genetic fixation in Lake Sturgeon populations in central Canada.

Hudson Bay drainage basin

Comparative genomics and genomic epidemiology of mycobacterium avium subsp. paratuberculosis strains

Two phenotypically distinct strains of Mycobacterium avium subsp. paratuberculosis (MAP) were recognized in the 1930s but it was not until the introduction of restriction endonuclease analysis (REA) in the mid-1980s that these two strains, MAP-C and MAP-S, could be distinguished genetically. Since then, a plethora of molecular typing techniques has been applied to MAP isolates (reviewed by Li et al. 2016; Fawzy et al., 2018) and a complex nomenclature for MAP strains has evolved. Currently, the most widely used genotyping method is Mycobacterial Interspersed Repetitive Units – Variable-Number Tandem Repeats (MIRU-VNTR). However, it has limited discriminatory power within the major lineages and does not always accurately reflect genetic relatedness since the repeat sequences are subject to homoplasy (Ahlstrom et al., 2015; Bryant et al., 2016). Whole genome sequencing (WGS) supplies the ultimate resolution and has revolutionized MAP research. It has enabled determination of single nucleotide polymorphism (SNP) level diversity, clarified phylogenetic relationships between divergent lineages and closely related strains and spawned the development of novel genotyping methods based on informative canonical SNPs (Leão et al. 2016; Ahlstrom et al. 2016b). This chapter presents an overview of comparative genomics and epidemiology of MAP strains and also highlights the role that WGS has played in increasing our understanding of MAP strain diversity.

Book chapter

Genomes & islands & evolution: Oh my!

A central question in evolutionary biology is how lineages quickly diversify to occupy different ecological niches, along with determining genomic factors that facilitate evolutionary change. Isolated, oceanic archipelagos are famous for adaptive radiations characterized by endemic, species-rich clades with substantial ecological variation, yet genome resources key to determining eco-evo processes are generally lacking. Here I present a comparison of the number of genome reference assemblies available (as of May 31, 2023) for three major eukaryotic lineages, briefly describe genome sequencing and benchmarking strategies, and highlight as a case study a genome assembly project for Bidens hawaiensis (Koʻokoʻolau, Asteraceae or Compositae; Coreopsidae), a member of a hexaploid Hawaiian plant adaptive radiation. The total number of plant genome references (1,394) was found to substantially lag the total number of genome references for animal (6,003) and fungi (4,400). Improvements to the quality of de novo assembled genomes are fueled by second- and third-generation long-read sequencing advancements, among other sequencing approaches. In conjunction, strategies to improve genome contiguity include optical maps, Hi-C chromatin capture, or trio binning. Continual improvements to genome sequencing and assembly algorithms have brought within reach telomere-to-telomere genome assemblies, albeit this level of sequencing has to date only been achieved in a few cases. With improvements in sequencing techniques and per-base pair costs that continue to trend downward, the number of high-quality genomes is anticipated to continue to increase, leading to the filling in of taxonomic gaps and sampling of groups of taxa from under sampled geographic areas. Increasing the number of plant genome resources available for the study of island endemism could help to shed light on genome-phenome relationships and genome characteristics that have produced the stunning biological diversity that we now observe across the globe.

Capitulum

Using advanced population genomics to better understand the relationship between offshore and spawning habitat use for Atlantic Sturgeon

Atlantic Sturgeon ( Acipenser oxyrinchus oxyrinchus ) are a large-bodied anadromous fish that historically supported important fisheries along the east coast of the United States. Following years of overharvest and habitat degradation, populations experienced severe declines. In 2012, the National Marine Fisheries Service listed Atlantic Sturgeon under the Endangered Species Act (ESA; 61 FR 4722). Their listing named five Distinct Population Segments (DPSs), predicated on genetic groups composed of geographically proximate populations. Federal management of Atlantic Sturgeon presents challenges, as sturgeon from each of the five DPSs mix extensively in coastal and marine habitats yet take and recovery progress must be evaluated separately for each unit. Genetic assignment testing based on mitochondrial and microsatellite markers allows individuals to be assigned back to their natal river and DPS. However, this approach is not perfect and some individuals may be incorrectly assigned. Recent advances in genomics offer the potential of a higher resolution approach to genetic assignment testing, and thus may reduce uncertainty associated with assignment testing. In addition, genomics allows a greater number of markers to be examined from across a broader portion of the sturgeon genome, thus may provide an enhanced perspective of population structure for the species, and potentially allow other previously intractable questions to be addressed (Bernatchez et al. 2017, Supple and Shapiro 2018). We used next-generation sequencing to develop a draft genome for Atlantic Sturgeon and identify single nucleotide polymorphisms (SNPs) that could be used to resolve the natal river and DPS of individual Atlantic Sturgeon. We identified 1,210 candidate SNPs within the nuclear genome as well as 49 SNPs within the mitochondrial genome. After filtering and review, we selected 161 nuclear SNPs and 39 mitochondrial SNPs for further testing and evaluation. We used genotyping-in-thousands by sequencing (GT-seq) to simultaneously sequence nuclear SNP loci, mitochondrial SNP loci, and the existing panel of twelve microsatellite loci. This effort required a pilot sequencing run on a single sturgeon sample to test marker amplification and refine primer strengths, followed by a series of sequencing runs to generate baseline data for 288 individuals representing nine populations of Atlantic Sturgeon in four DPSs. Using baseline data from the nine populations, we ran a series of genomic analyses to characterize diversity within and among populations, providing a benchmark for this species using the new SNP markers. Allelic richness was similar for all populations, although there was a general trend of more northern population containing greater levels of allelic richness. Interestingly, we observed linkage disequilibrium among many pairs of loci within many populations. This might be the result of physical linkage but could also suggest these populations are recovering from genetic bottlenecks and/or are effectively small, leading to specific haplotypes to be favored by chance. Pairwise differentiation among populations varied among the populations ( F ST range: 0.010-0.098) and was significantly correlated ( r = 0.771; P < 0.001) to pairwise F ST observed using microsatellite markers). Population clustering and ordination techniques using the new genomic data both support an overall population structure that is similar to the current DPS management units (which were developed primarily based on microsatellite genetic data). Overall, this suggests that existing microsatellite markers and the panel of SNP markers developed in this study provide similar information about the populations structure and ecology of Atlantic Sturgeon. Given the observed differences in allele frequencies among populations, our genomic baseline supports previous assertations that Atlantic Sturgeon show natal homing, despite mixing extensively in marine waters during non-breeding periods. Lower levels of differentiation between populations in the South Atlantic DPS suggest that populations in this region may have greater levels of gene flow relative to their more northerly conspecifics, which has also previously been suggested based on microsatellite data. The observed differentiation among populations provides the necessary foundation for determining the natal river and DPS of Atlantic Sturgeon using assignment testing. We tested the utility of our new genomic baseline for resolving the population and DPS of Atlantic Sturgeon. Our nuclear SNP markers showed utility for identifying the origin of unknown Atlantic Sturgeon samples, as 86.5% were assigned to the correct DPS and 66.3% were assigned to the correct natal river. However, since this study was funded the Conservation Genetics and Genomics Laboratory at Leetown Science Center has made significant improvements to their microsatellite genetic baseline, which now performs more effectively than our new genomic approach (the genetic baseline includes 12 populations and 5 DPSs, and correctly assigns 95.8% of individuals to DPS and 84.9% of individuals to their natal population using 12 microsatellite loci). We conducted an ad hoc exploration of how additional microsatellite or nuclear SNP loci may further improve the accuracy of assignment testing. We found that additional microsatellite markers are likely to result in greater improvements in assignment efficiency than additional nuclear SNPs. However, a much larger number of SNP loci (which if identified could be sequenced using other methods that are now available; e.g., the RAD-capture approach published by Ali et al. 2016) could produce assignment efficiencies that are greater than what is currently feasible using microsatellites. In the absence of further research and development of additional SNP markers for Atlantic Sturgeon (possibly using an approach other than GT-seq), the existing microsatellite loci are the most effective means available to determine the natal river and DPS of Atlantic Sturgeon encountered in offshore waters. Because our new genomic markers were less effective than the existing panel of 12 microsatellite markers, we chose to use the existing microsatellite markers to assign Atlantic Sturgeon captured in another BOEM-funded study (cooperative agreement M16AC00003; Monitoring endangered Atlantic Sturgeon and commercial finfish habitat use offshore New York) following consultation with our project officer. Using this approach, we genotyped and assigned 186 Atlantic Sturgeon captured in coastal waters off the Rockaway Peninsula, New York. The vast majority of these sturgeon were assigned to the New York Bight DPS (94.62%), and most appear to belong to the Hudson River population (87.10%) with smaller contributions from the Delaware River population (7.53%). Smaller contributions (2.15%) were observed from six other populations, including those from the James, York, Kennebec, Ogeechee, and Edisto rivers. Although most of the fish we assigned were assigned to the nearest spawning rivers (Hudson and Delaware), the contributions from distant rivers is consistent with the propensity of this species to move long distances and form mixed stock aggregations along the continental shelf. This finding indicates that spawning populations (and their corresponding DPS) from distant locations may potentially be impacted by offshore activities. In fact, activities in this region of the New York Bight could negatively impact Atlantic Sturgeon population from at least four different DPSs. Genetic or genomic assignment testing remains an essential tool to characterize potential impacts to Atlantic Sturgeon populations and should be applied more broadly to better characterize potential impacts of activities in other locations.

Atlantic Coast

Comparative genomics of Bacillus anthracis A and B-clades reveals genetic variation in genes responsible for spore germination

Bacillus anthracis , the causative agent of anthrax, is composed of three genetic clades (A, B, and C). Clade-A is the most common and distributed worldwide, B-clade has a narrow geographic distribution, and C-clade is rare. South Africa's Kruger National Park (KNP) has high diversity of B. anthracis , with strains from A and B clades described from its northernmost region, Pafuri. We employed whole genome sequencing to investigate the genomic diversity of B. anthracis strains isolated from animal carcasses ( n = 34) during the 2012–2015 outbreaks. Whole-genome single-nucleotide polymorphism (wgSNP) analysis assigned the 2012–2015 B. anthracis genomes to the A-clade branch, distributed across the branch's two minor sub-clades A.Br.005/006. Additionally, pan-genomic analysis distinguished the A- and B-clade genomes, identifying unique accessory genes. Notable genetic differences include the biosynthetic spore cell wall genes; long-chain fatty acid CoA ligases ( FaD13 ), Bacillus collagen-like protein of anthracis (BclA) involved in the exosporium germination, as well as a truncated murein DD-endopeptidase ( mepH ) found in the pXO2 plasmid of the B-clade strains. The tryptophan synthase subunit alpha gene ( trpA ), which results in a pseudogene in B-clade genomes separates the A- and B-clade genomes. These differences in biosynthetic cell wall genes suggest variation in adaptability or cell growth of the B-clade strains in the environment, further influenced by the truncation of the trpA gene involved in spore germination. The A.Br.005/006-clade strains in KNP exhibit higher genetic diversity, which may enhance their resilience to environmental stressors. In contrast, the KNP B-clade (B.Br.001/002) strains show limited genetic variation, potentially reducing their adaptability. This pattern is evident through whole-genome SNP analysis and pan-genomics investigating the evolution of B. anthracis .

Kruger National Park

The expectations and challenges of wildlife disease research in the era of genomics: Forecasting with a horizon scan-like exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented ‘Big Data’ tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural ‘Genomics of Disease in Wildlife’ workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed five themes considered to be the most pressing to the application of genomics in wildlife disease research: i) “Improving Communication”, ii) “Methodological and Analytical Advancements”, iii) “Translation into Practice”, iv) “Integrating Landscape Ecology and Genomics”, and v) “Emerging New Questions”. Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

Journal of Heredity

Exploring genetic variation and population structure in a threatened species, Noturus placidus, with whole-genome sequence data

The Neosho madtom ( Noturus placidus ) is a small catfish, generally less than 3 inches in length, unique to the Neosho-Spring River system within the Arkansas River Basin. It was federally listed as threatened in 1990, largely due to habitat loss. For conservation efforts, we generated whole-genome sequence data from 10 Neosho madtom individuals originating from 3 geographically separated populations to evaluate genetic diversity and population structure. A Neosho madtom genome was de novo assembled, and genome size and content were assessed. Single nucleotide polymorphisms were assessed from de Bruijn graphs, and via reference alignment with both the channel catfish ( Ictalurus punctatus) reference genome and Neosho madtom reference genome. Principal component analysis and structure analysis indicated weak population structure, suggesting fish from the 3 locations represent a single population. Using a novel method, genome-wide conservation and divergence between the Neosho madtom, channel catfish, and zebrafish ( Danio rerio ) was assessed by pairwise contig alignment, which demonstrated that genes important to embryonic development frequently had conserved sequences. This research in a threatened species with no previously published genomic resources provides novel genetic information to guide current and future conservation efforts and demonstrates that using whole-genome sequencing provides detailed information of population structure and demography using only a limited number of rare and valuable samples.

Arkansas, Kansas, Missouri, Oklahoma

The transformative impact of genomics on sage-grouse conservation and management

For over two decades, genetic studies have been used to assist in the conservation and management of both Greater Sage-grouse ( Centrocercus urophasianus ) and Gunnison Sage-grouse ( C. minimus ), addressing a wide variety of topics including taxonomy, parentage, population connectivity, and demography. The field of conservation genetics has been transformed by dramatic improvements in sequencing technology, facilitating genomic studies in many wildlife species. The quality and amount of data generated by genomic methods vastly exceed that of traditional genetic studies, allowing for increased precision in estimating genetic parameters of interest. Perhaps more importantly, genomic methods can provide insight into non-neutral evolution such as adaptive divergence. Here we recount the shift from genetic to genomic methods using two wildlife species of substantial conservation interest, focusing on the improved capabilities and advantages of genomic methods. For instance, reassessment of divergence in sage-grouse using genomic methods confirmed strong differentiation between the two species and revealed that a small population in the state of Washington was more genetically distinct than previously recognized. Further, new genomic resources and approaches have been used to identify a family of genes linked to local dietary adaptation suggesting that sage-grouse may possess digestive and metabolic adaptations that mitigate the effects of consuming plant secondary metabolites like those found in sagebrush. Genetic variation among populations in these gene regions is thought to be involved with local dietary adaptations, and therefore maintaining the tie between sage-grouse and the chemistry of local sagebrush may be an important management consideration. We posit that the integration of newly developed genomic resources combined with the vast wealth of ecological and behavioral data for sage-grouse has the potential to shed light on mechanistic relationships that ultimately are vital to the conservation and management of these species.

Book chapter

A chromosome-level genome assembly of a vernal pool specialist amphibian, the Western Spadefoot, Spea hammondii

We assembled and annotated a chromosome-level genome for the Western Spadefoot, Spea hammondii (Anura, Scaphiopodidae) representing one of only three amphibians included in the California Conservation Genomics Project (CCGP). Spea hammondii is a vernal pool breeding anuran native to California and northwestern Baja California which has undergone both range contractions and local extirpations across its distribution, primarily due to habitat loss and degradation and drought. The species is recognized by the state of California as a Species of Special Concern and is proposed for listing under the United States Endangered Species Act. Using the established CCGP pipeline, this S. hammondii genome was produced using Pacific Biosciences HiFi long-reads and Omni-C proximity ligation, resulting in a de novo genome assembly 1.14 Gb in length, distributed across 479 scaffolds (scaffold N50 = 120.8 Mb; largest scaffold = 183.6 Mb) with a BUSCO completeness score of 90.9% using a conserved tetrapod ortholog set. Our assembly shows high base accuracy (quality value [QV] = 63.7) and low frameshift error in coding regions (QV 50.42). Annotation of this genome yielded 20,434 genes with a BUSCO completeness score of 94.7%. This genome assembly, in combination with range-wide resequencing data from CCGP, will facilitate statewide population genomic assessments to delineate conservation units, quantify inbreeding and genomic load, and test for adaptive variation associated with vernal pool hydrology and drought tolerance, all of which are important considerations in the proposed federal listing.

Journal of Heredity

A chromosome-level genome assembly of a vernal pool specialist amphibian, the Western Spadefoot, Spea hammondii

We assembled and annotated a chromosome-level reference genome for the Western Spadefoot, Spea hammondii (Anura, Scaphiopodidae) representing one of only three amphibians included in the California Conservation Genomics Project (CCGP). Spea hammondii is a vernal pool breeding anuran native to California and northwestern Baja California which has undergone both range contractions and local extirpations across its distribution, primarily due to habitat loss and degradation and drought. The species is recognized by the state of California as a Species of Special Concern and is proposed for listing under the United States Endangered Species Act. Using the established CCGP pipeline, this S. hammondii genome was produced using Pacific Biosciences HiFi long-reads and Omni-C proximity ligation, resulting in a de novo genome assembly 1.14 Gb in length, distributed across 479 scaffolds (scaffold N50 = 120.8 Mb; largest scaffold = 183.6 Mb) with a BUSCO completeness score of 90.9% using a conserved tetrapod ortholog set. Our assembly shows high base accuracy (QV = 63.7) and low frameshift error in coding regions (QV 50.42). Annotation of this genome yielded 20,434 genes with a BUSCO completeness score of 94.7%. This reference genome, in combination with range-wide resequencing data from CCGP, will facilitate statewide population genomic assessments to delineate conservation units, quantify inbreeding and genomic load, and test for adaptive variation associated with vernal pool hydrology and drought tolerance, all of which are important considerations in the proposed federal listing.

BioRxiv

Species-specific responses to landscape features shaped genomic structure within Alaska galliformes

Aim Connectivity is vital to the resiliency of populations to environmental change and stochastic events, especially for cold-adapted species as Arctic and alpine tundra habitats retract as the climate warms. We examined the influence of past and current landscapes on genomic connectivity in cold-adapted galliformes as a critical first step to assess the vulnerability of Alaska ptarmigan and grouse to environmental change. We hypothesize that the mosaic of physical features and habitat within Alaska promoted the formation of genetic structure across species. Location Alaska, United States of America. Taxa Ptarmigan and Grouse (Galliformes: Tetraoninae). Methods We collected double digest restriction-site-associated DNA sequence data from six ptarmigan and grouse species ( N = 13–145/species) sampled across multiple ecosystems up to ~10 degrees of latitude. Spatial genomic structure was analysed using methods that reflect different temporal scales: (1) principal components analysis to identify major trends in the distribution of genomic variation; (2) maximum likelihood clustering analyses to test for the presence of multiple genomic groupings; (3) shared co-ancestry analyses to assess contemporary relationships and (4) effective migration surfaces to identify regions that deviate from a null model of isolation by distance. Results Levels of genomic structure varied across species (Φ ST =0.009–0.042). Three general patterns of structure emerged: (1) east-west partition located near the Yukon-Tanana uplands; (2) north-south split coinciding with the Alaska Range and (3) northern group near the Brooks Range. Species-specific patterns were observed; not all landscape features were barriers to gene flow for all ptarmigan and grouse and temporal contrasts were detected at the Brooks Range. Main conclusions Within Alaska galliformes, patterns of genomic structure coincide with physiographic features and highlight the importance of physical and ecological barriers in shaping how genomic diversity is arrayed across the landscape. Lack of concordance in spatial patterns indicates that species behaviour and habitat affinities play key roles in driving the contrasting patterns of genomic structure.

Alaska

Genomic identity of white oak species in an eastern North American syngameon

The eastern North American white oaks, a complex of approximately 16 potentially interbreeding species, have become a classic model for studying the genetic nature of species in a syngameon. Genetic work over the past two decades has demonstrated the reality of oak species, but gene flow between sympatric oaks raises the question of whether there are conserved regions of the genome that define oak species. Does gene flow homogenize the entire genome? Do the regions of the genome that distinguish a species in one part of its range differ from the regions that distinguish it in other parts of its range, where it grows in sympatry with different species? Or are there regions of the genome that are relatively conserved across species ranges? In this study, we revisit seven species of the eastern North American white oak syngameon using a set of 80 single-nucleotide polymorphisms (SNPs) selected in a previous study because they show differences among, and consistency within, the species. We test the hypothesis that there exist segments of the genome that do not become homogenized by repeated introgression, but retain distinct alleles characteristic of each species. We undertake a range-wide sampling to investigate whether SNPs that appeared to be fixed based on a relatively small sample in our previous work are fixed or nearly fixed across the range of the species. Each of the seven species remains genetically distinct across its range, given our diagnostic set of markers, with relatively few individuals exhibiting admixture of multiple species. SNPs map back to all 12 Quercus linkage groups (chromosomes) and are separated from each other by an average of 7.47 million bp (± 8.74 million bp, SD), but are significantly clustered relative to a random null distribution, suggesting that our SNP toolkit reflects genome-wide patterns of divergence while potentially being concentrated in regions of the genome that reflect a higher-than-average history of among-species divergence. This application of a DNA toolkit designed for the simple problem of identifying species in the field has two important implications. First, the eastern North American white oak syngameon is composed of entities that most taxonomists would consider “good species.” Second, and more fundamentally, species in the syngameon are genetically coherent because characteristic portions of the genome remain divergent despite a history of introgression. Understanding the conditions under which some loci diverge while others introgress is key to understanding the origins and maintenance of global tree diversity.

Estern United States

A haploid pseudo-chromosome genome assembly for a keystone sagebrush species of western North American rangelands

Increased ecological disturbances, species invasions, and climate change are creating severe conservation problems for several plant species that are widespread and foundational. Understanding the genetic diversity of these species and how it relates to adaptation to these stressors are necessary for guiding conservation and restoration efforts. This need is particularly acute for big sagebrush ( Artemisia tridentata ; Asteraceae), which was once the dominant shrub over 1,000,000 km 2 in western North America but has since retracted by half and thus has become the target of one of the largest restoration seeding efforts globally. Here, we present the first reference-quality genome assembly for an ecologically important subspecies of big sagebrush ( A. tridentata subsp. tridentata ) based on short and long reads, as well as chromatin proximity ligation data analyzed using the HiRise pipeline. The final 4.2-Gb assembly consists of 5,492 scaffolds, with nine pseudo-chromosomal scaffolds (nine scaffolds comprising at least 90% of the assembled genome; n = 9). The assembly contains an estimated 43,377 genes based on ab initio gene discovery and transcriptional data analyzed using the MAKER pipeline, with 91.37% of BUSCOs being completely assembled. The final assembly was highly repetitive, with repeat elements comprising 77.99% of the genome, making the Artemisia tridentata subsp. tridentata genome one of the most highly repetitive plant genomes to be sequenced and assembled. This genome assembly advances studies on plant adaptation to drought and heat stress and provides a valuable tool for future genomic research.

Arizona, California, Colorado, Idaho, Montana, Neb

Speciation genomics in the tiger whiptail lizards ( Aspidoscelis tigris complex)

The transition from small genetic to genome-scale datasets for studying biodiversity has revealed that genetic exchange through introgressive hybridization is a widespread phenomenon in nature. Despite this, a lack of high-quality reference genomes for most non-model species limits our understanding of the impact of this process for many taxonomic groups. This restricts the range of insights that genomic tools can provide for conservation biologists, who often hope to employ genomic datasets to accurately identify historically isolated lineages to protect and to predict their evolutionary fate in the face of environmental change. Tiger whiptail lizards ( Aspidoscelis tigris complex) are an abundant and important ecological component of ecosystems across the southwestern United States. In this study, we assembled and annotated a chromosome-level reference genome for A. t. stejnegeri from coastal California. We then used this reference genome to reconstruct patterns of speciation and admixture within the larger species complex, finding evidence that gene flow is widespread both geographically and across the genome.

northern Mexico, southwestern United States

Ecology shapes the genomic and biosynthetic diversification of Streptomyces bacteria from insectivorous bats

Streptomyces are prolific producers of secondary metabolites from which many clinically useful compounds have been derived. They inhabit diverse habitats but have rarely been reported in vertebrates. Here, we aim to determine to what extent the ecological source (bat host species and cave sites) influence the genomic and biosynthetic diversity of Streptomyces bacteria. We analysed draft genomes of 132 Streptomyces isolates sampled from 11 species of insectivorous bats from six cave sites in Arizona and New Mexico, USA. We delineated 55 species based on the genome-wide average nucleotide identity and core genome phylogenetic tree. Streptomyces isolates that colonize the same bat species or inhabit the same site exhibit greater overall genomic similarity than they do with Streptomyces from other bat species or sites. However, when considering biosynthetic gene clusters (BGCs) alone, BGC distribution is not structured by the ecological or geographical source of the Streptomyces that carry them. Each genome carried between 19–65 BGCs (median=42.5) and varied even among members of the same Streptomyces species. Nine major classes of BGCs were detected in ten of the 11 bat species and in all sites: terpene, non-ribosomal peptide synthetase, polyketide synthase, siderophore, RiPP-like, butyrolactone, lanthipeptide, ectoine, melanin. Finally, Streptomyces genomes carry multiple hybrid BGCs consisting of signature domains from two to seven distinct BGC classes. Taken together, our results bring critical insights to understanding Streptomyces -bat ecology and BGC diversity that may contribute to bat health and in augmenting current efforts in natural product discovery, especially from underexplored or overlooked environments.

Microbial Genomics