USGS Science⌕ Search

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 235 records · Page 13Linked to original sources

Identification of a novel Adélie penguin circovirus at Cape Crozier (Ross Island, Antarctica)

Understanding the causes of disease in Antarctic wildlife is crucial as many of these species are already threatened by environmental changes brought about by climate change. In recent years, Antarctic penguins have been showing signs of an unknown pathology: a feather disorder characterised by missing feathers resulting in exposed skin. During the 2018-19 austral summer breeding season at Cape Crozier colony on Ross Island, Antarctica, for the first time we observed an Adlie penguin chick missing down over most of its body. A guano sample was collected from the nest of the featherless chick and, using high throughput sequencing we identified a novel circovirus. Using abutting primers, we amplified the full genome, which we cloned, and Sanger sequenced to determine the complete genome of the circovirus. The Adlie penguin guano-associated circovirus genome shares <67% genome-wide nucleotide identity to other circoviruses, representing a new species of circovirus, therefore we named it penguin circovirus (PenCV). Using the same primer pair, we screened 25 previously collected cloacal swabs taken at Cape Crozier from known-age adult Adlie penguins during the 2014-15 season displaying no clinical signs of feather-loss disorder. Three of the 25 samples (12%) were positive for a PenCV, whose genome shared >99% pairwise identity to the one identified in 2018-19. This is the first report of a circovirus associated with a penguin species and could be a possible etiological agent of the feather-loss disorder in Antarctic penguins.

Viruses↗

Mutational analyses of molecularly cloned satellite tobacco mosaic virus during serial passage in plants: Evidence for hotspots of genetic change

The high level of genetic diversity and rapid evolution of viral RNA genomes are well documented, but few studies have characterized the rate and nature of ongoing genetic change over time under controlled experimental conditions, especially in plant hosts. The RNA genome of satellite tobacco mosaic virus (STMV) was used as an effective model for such studies because of advantageous features of its genome structure and because the extant genetic heterogeneity of STMV has been characterized previously. In the present study, the process of genetic change over time was studied by monitoring multiple serial passage lines of STMV populations for changes in their consensus sequences. A total of 42 passage lines were initiated by inoculation of tobacco plants with a helper tobamovirus and one of four STMV RNA inocula that were transcribed from full-length infectious STMV clones or extracted from purified STMV type strain virions. Ten serial passages were carried out for each line and the consensus genotypes of progeny STMV populations were assessed for genetic change by RNase protection analyses of the entire 1,059-nt STMV genome. Three different types of genetic change were observed, including the fixation of novel mutations in 9 of 42 lines, mutation at the major heterogeneity site near nt 751 in 5 of the 19 lines inoculated with a single genotype, and selection of a single major genotype in 6 of the 23 lines inoculated with mixed genotypes. Sequence analyses showed that the majority of mutations were single base substitutions. The distribution of mutation sites included three clusters in which mutations occurred at or very near the same site, suggesting hot spots of genetic change in the STMV genome. The diversity of genetic changes in sibling lines is clear evidence for the important role of chance and random sampling events in the process of genetic diversification of STMV virus populations.

RNA↗

Novel avian papillomaviruses identified in a south polar skua sampled on Ross Island, Antarctica

Papillomaviruses are small circular DNA viruses that infect epithelial cells of their hosts. Avian papillomaviruses are poorly sampled/documented compared to those infecting humans. We used a viral metagenomic approach to identify viruses from the oral swab taken from a deceased south polar skua ( Stercorarius maccormicki ) found at Cape Royds, Ross Island, Antarctica in late 2024. We identified three papillomaviruses and determined their complete genomes, Stercorarius maccormicki papillomavirus (SmacPV) 1-3. SmacPV1 is the most divergent of the three SmacPVs, sharing 62% genome-wide pairwise identity to SmacPV2 and SmacPV3 and <63.5% to other avian papillomaviruses. The genomes of SmacPV2 and SmacPV3 represent two new papillomavirus types sharing 82.4% genome-wide pairwise identity with each other and <72% to other papillomaviruses. SmacPV2 and SmacPV3 phylogenetically cluster with sequences of the Rissa tridactyla papillomavirus 1 from black-legged kittiwake ( Rissa tridactyla ), Larus smithsonianus papillomavirus 1 from American herring gull ( Larus smithsonianus ) and Fratercula arctica papillomavirus 1 from Atlantic puffin ( Fratercula arctica ), and they collectively represent a new papillomavirus species. These are the first papillomaviruses to be identified in Stercorarius spp . and add to the handful of known papillomaviruses in identified avian species. We also expand the known host range of papillomaviruses in Antarctic animals, which previously included Adélie penguins ( Pygoscelis adeliae ), Weddell seals ( Leptonychotes weddellii ), Antarctic fur seals ( Arctocephalus gazella ), leopard seals ( Hydrurga leptonyx ) and emerald notothen ( Trematomus bernacchii ).

Archives of Virology↗

Female persistence during toxicant treatment predicts survival probability of offspring in invasive brown treesnakes (Boiga irregularis)

Assessing the long-term efficacy of control methods is a critical component of invasive species management. For example, if traits related to control have significant heritability or are influenced by maternal effects, control methods may lose efficacy over time. The potential for these effects can be evaluated via parent/offspring survival analysis, which concomitantly recasts adaptive management as an evolutionary force for invasive species. However, difficulties can arise when the life history of an invasive is cryptic, precluding direct observations of familial relationships. Genomic pedigree reconstruction can facilitate such analyses by assigning offspring to parents in invasive species for which mating and reproduction are difficult to study. Here, we use genomic pedigree reconstruction to quantify parental longevity and probability of offspring survival for brown treesnakes ( Boiga irregularis ) on Guam in a landscape treated with toxic baits simulating application via an aerial delivery system (ADS). To do so, we used 398 single nucleotide polymorphisms (SNPs) to update an existing multi-generational genomic pedigree for a geographically-closed population of brown treesnakes. This facilitated assignment of parents to juveniles born during three consecutive years of toxic bait application under a simulated aerial treatment program (N = 72). We found that the offspring of dams that persisted until the end of the study displayed greater survival probability (cox proportional hazard model, P < 0.001), yet there was no such effect for sires. This sex-specific relationship between parental longevity and offspring survival indicates that heritability of traits contributing to resistance to ADS is unlikely, but it supports a role of maternal effects that could undermine ADS. Our study identifies potential risks associated with control efforts and also highlights the utility of parent-offspring survival analyses informed by genomic pedigree reconstruction as a tool for adaptive management.

Global Ecology and Conservation↗

Mercury methylation genes identified across diverse anaerobic microbial guilds in a eutrophic sulfate-enriched lake

Mercury (Hg) methylation is a microbially mediated process that converts inorganic Hg into bioaccumulative, neurotoxic methylmercury (MeHg). The metabolic activity of methylating organisms is highly dependent on biogeochemical conditions, which subsequently influences MeHg production. However, our understanding of the ecophysiology of methylators in natural ecosystems is still limited. Here, we identified potential locations of MeHg production in the anoxic, sulfidic hypolimnion of a freshwater lake. At these sites, we used shotgun metagenomics to characterize microorganisms with the Hg-methylation gene hgcA . Putative methylators were dominated by hgcA sequences divergent from those in well-studied, confirmed methylators. Using genome-resolved metagenomics, we identified organisms with hgcA (hgcA+) within the Bacteroidetes and the recently described Kiritimatiellaeota phyla. We identified hgcA+ genomes derived from sulfate-reducing bacteria, but these accounted for only 22% of hgcA+ genome coverage. The most abundant hgcA+ genomes were from fermenters, accounting for over half of the hgcA gene coverage. Many of these organisms also mediate hydrolysis of polysaccharides, likely from cyanobacterial blooms. This work highlights the distribution of the Hg-methylation genes across microbial metabolic guilds and indicate that primary degradation of polysaccharides and fermentation may play an important but unrecognized role in MeHg production in the anoxic hypolimnion of freshwater lakes.

Environmental Science and Technology↗

Toxicogenomics in regulatory ecotoxicology

Recently, we have witnessed an explosion of different genomic approaches that, through a combination of advanced biological, instrumental, and bioinformatic techniques, can yield a previously unparalleled amount of data concerning the molecular and biochemical status of organisms. Fueled partially by large, well-publicized efforts such as the Human Genome Project, genomic research has become a rapidly growing topical area in multiple biological disciplines. Since 1999, when the term “toxicogenomics” was coined to describe the application of genomics to toxicology (1), a rapid increase in publications on the topic has occurred (Figure 1). The potential utility of toxicogenomics in toxicological research and regulatory activities has been the subject of scientific discussions and, as with any new technology, has evoked a wide range of opinion (2–6).

Environmental Science & Technology↗

Historical DNA reveals climate adaptation in an endangered songbird

To cope with climate change, species may shift their distributions or adapt in situ to changing environmental conditions. However, clear examples of genetic changes via adaptation are limited. We explore evolutionary responses to climate change in the endangered southwestern willow flycatcher ( Empidonax traillii extimus ) through whole-genome comparisons between historical specimens, collected from 1888 to 1909 near San Diego, California, United States, and contemporary individuals from across the breeding range. Genomic analyses revealed that introgression into San Diego increased adaptive potential over time and shifted genome-wide population structure towards that of neighbouring populations. In contrast, loci linked to climate (dew point temperature and precipitation) shifted away from neighbouring populations and in a direction consistent with adaptation to climate change in southern California. This research highlights the role of admixture in facilitating adaptive shifts through its impact on genome-wide genetic variation and represents one of the few studies to document climate adaptation in a wild population.

Nature Climate Change↗

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining ( MAT ) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By generating additional chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere-flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight toward the identification of a sexual cycle in Malassezia , with possible implications for pathogenicity.

PNAS↗

Identical sequence types of Yersinia ruckeri associated with lethal disease in wild-caught invasive Blue Catfish and cultured hybrid catfish (Channel Catfish ♀ × Blue Catfish ♂) from disparate aquatic ecosystems

Objective The Blue Catfish Ictalurus furcatus is commonly raised in warmwater aquaculture in the United States to produce Channel Catfish I. punctatus × Blue Catfish hybrids. It is also a prominent aquatic invasive species of concern in the mid-Atlantic region of the United States. Here, Yersina ruckeri was isolated from moribund Blue Catfish and hybrid catfish from disparate regions of the USA. The goal of the research here was to compare these Y. ruckeri strains to each other and other known strains for which adequate sequence data was available. In addition, we sought to determine if the strain from Blue Catfish was pathogenic to Rainbow Trout Oncorhynchus mykiss . Methods Moribund hybrid catfish from culture ponds in Mississippi were processed for diagnostic evaluation in March 2016. In April 2022, a moribund Blue Catfish specimen was collected from a tributary of the Nanticoke River in Maryland. Bacterial isolates were identified and characterized using biochemical tests, antimicrobial sensitivity screening, serotyping, and complete or partial genome sequencing. Disease pathology was described via histology. The isolate from Blue Catfish was used in challenge experiments to determine if it was pathogenic to Rainbow Trout. Multilocus sequencing typing was conducted using the PubMLST database. Results Biochemical testing was consistent with Y. ruckeri . A draft genome of the Y. ruckeri isolate was assembled based on Oxford Nanopore Technology sequencing and identified a single genomic replicon (3,791,418 bp) consistent in size to other Y. ruckeri genomes and a pLT plasmid (60, 933 bp). The challenge study demonstrated no significant virulence of this isolate for Rainbow Trout ( Y. ruckeri ). This isolate was most similar to other strains isolated from ictalurids. Notably, the gyrase B gene from this isolate was identical to that of archived strains isolated from moribund Mississippi hybrid catfish aquaculture during 2016 and these isolates share identical PubMLST sequence type profiles. Similarly, they shared a pLT plasmid that differed by only 6 bp. This plasmid has never been reported from trout isolates and appears to be unique to ictalurids. Conclusions Analyses here provide preliminary genetic evidence that geographically distant (Maryland and Mississippi, USA) isolates of Y. ruckeri from ictalurids are genetically similar to each other and Y. ruckeri (strain SC09) that infects ictalurids in China. This strain is not a biothreat to Rainbow Trout at typical culture temperatures.

Journal of Aquatic Animal Health↗

Incomplete bioinformatic filtering and inadequate age and growth analysis lead to an incorrect inference of harvested-induced changes

Understanding the evolutionary impacts of harvest on fish populations is important for informing fisheries management and conservation and has become a growing research topic over the last decade. However, the dynamics of fish populations are highly complex, and phenotypes can be influenced by many biotic and abiotic factors. Therefore, it is vital to collect robust data and explore multiple alternative hypotheses before concluding that fish populations are influenced by harvest. In their recently published manuscript, Bowles et al, Evolutionary Applications, 13(6):1128 conducted age/growth and genomic analysis of walleye ( Sander vitreus ) populations sampled 13–15 years (1–2.5 generations) apart and hypothesized that observed phenotypic and genomic changes in this time period were likely due to harvest. Specifically, Bowles et al. (2020) documented differential declines in size-at-age in three exploited walleye populations compared to a separate, but presumably less-exploited, reference population. Additionally, they documented population genetic differentiation in one population pair, homogenization in another, and outlier loci putatively under selection across time points. Based on their phenotypic and genetic results, they hypothesized that selective harvest had led to fisheries-induced evolution (referred to as nascent changes) in the exploited populations in as little as 1–2.5 generations. We re-analyzed their data and found that (a) sizes declined across both exploited and reference populations during the time period studied and (b) observed genomic differentiation in their study was the result of inadequate data filtering, including retaining individuals with high amounts of missing data and retaining potentially undersplit and oversplit loci that created false signals of differentiation between time points. This re-analysis did not provide evidence for phenotypic or genetic changes attributable to harvest in any of the study populations, contrasting the hypotheses presented by Bowles et al. (2020). Our comment highlights the potential pitfalls associated with conducting age/growth analyses with low sample sizes and inadequately filtering genomic datasets.

Quebec↗

A genetic signature of the evolution of loss of flight in the Galapagos cormorant

INTRODUCTION Changes in the size and proportion of limbs and other structures have played a key role in the evolution of species. One common class of limb modification is recurrent wing reduction and loss of flight in birds. Indeed, Darwin used the occurrence of flightless birds as an argument in favor of his theory of natural selection. Loss of flight has evolved repeatedly and is found among 26 families of birds in 17 different orders. Despite the frequency of these modifications, we have a limited understanding of their underpinnings at the genetic and molecular levels. RATIONALE To better understand the evolution of changes in limb size, we studied a classic case of recent loss of flight in the Galapagos cormorant ( Phalacrocorax harrisi ). Cormorants are large water birds that live in coastal areas or near lakes, and P. harrisi is the only flightless cormorant among approximately 40 extant species. The entire population is distributed along the coastlines of Isabela and Fernandina islands in the Galapagos archipelago. P. harrisi has a pair of short wings, which are smaller than those of any other cormorant. The extreme reduction of the wings and pectoral skeleton observed in P. harrisi is an attractive model for studying the evolution of loss of flight because it occurred very recently; phylogenetic evidence suggests that P. harrisi diverged from its flighted relatives within the past 2 million years. We developed a comparative and predictive genomics approach that uses the genome sequences of P. harrisi and its flighted relatives to find candidate genetic variants that likely contributed to the evolution of loss of flight. RESULTS We sequenced and de novo assembled the whole genomes of P. harrisi and three closely related flighted cormorant species. We identified thousands of coding variants exclusive to P. harrisi and classified them according to their probability of altering protein function based on conservation. Variants most likely to alter protein function were significantly enriched in genes mutated in human skeletal ciliopathies, including Ofd1 , Evc , Wdr34 , and Ift122 . We carried out experiments in Caenorhabditis elegans to confirm that a missense variant present in the Galapagos cormorant IFT122 protein is sufficient to affect ciliary function. The primary cilium is essential for Hedgehog (Hh) signaling in vertebrates, and individuals affected by ciliopathies have small limbs and ribcages, mirroring the phenotype of P. harrisi . We also identified a 4–amino acid deletion in the regulatory domain of Cux1 , a highly conserved transcription factor that has been experimentally shown to regulate limb growth in chicken. The four missing amino acids are perfectly conserved in all birds and mammals sequenced to date. We tested the consequences of this deletion in a chondrogenic cell line and showed that it impairs the ability of CUX1 to transcriptionally up-regulate cilia-related genes (some of which contain function-altering variants in P. harrisi ) and to promote chondrogenic differentiation. Finally, we show that positive selection may have played a role in the fixation of the variants associated with loss of flight in P. harrisi . CONCLUSION Our results indicate that the combined effect of variants in genes necessary for the correct transcriptional regulation and function of the primary cilium likely contributed to the evolution of highly reduced wings and other skeletal adaptations associated with loss of flight in P. harrisi . Our approach may be generally useful for identification of variants underlying evolutionary novelty from genomes of closely related species.

Science↗

Detection of diazotrophy in the acetylene-fermenting anaerobe Pelobacter sp. strain SFB93

Acetylene (C2H2) is a trace constituent of the present Earth's oxidizing atmosphere, reflecting a mixture of terrestrial and marine emissions from anthropogenic, biomass-burning, and unidentified biogenic sources. Fermentation of acetylene was serendipitously discovered during C2H2 block assays of N2O reductase, and Pelobacter acetylenicus was shown to grow on C2H2 via acetylene hydratase (AH). AH is a W-containing, catabolic, low-redox-potential enzyme that, unlike nitrogenase (N2ase), is specific for acetylene. Acetylene fermentation is a rare metabolic process that is well characterized only in P. acetylenicus DSM3246 and DSM3247 and Pelobacter sp. strain SFB93. To better understand the genetic controls for AH activity, we sequenced the genomes of the three acetylene-fermenting Pelobacter strains. Genome assembly and annotation produced three novel genomes containing gene sequences for AH, with two copies being present in SFB93. In addition, gene sequences for all five compulsory genes for iron-molybdenum N2ase were also present in the three genomes, indicating the cooccurrence of two acetylene transformation pathways. Nitrogen fixation growth assays showed that DSM3426 could ferment acetylene in the absence of ammonium, but no ethylene was produced. However, SFB93 degraded acetylene and, in the absence of ammonium, produced ethylene, indicating an active N2ase. Diazotrophic growth was observed under N2 but not in experimental controls incubated under argon. SFB93 exhibits acetylene fermentation and nitrogen fixation, the only known biochemical mechanisms for acetylene transformation. Our results indicate complex interactions between N2ase and AH and suggest novel evolutionary pathways for these relic enzymes from early Earth to modern days.

Applied and Environmental Microbiology↗

Development and application of a qPCR-based genotyping assay for Ophidiomyces ophidiicola to investigate the epidemiology of ophidiomycosis

Ophidiomycosis (snake fungal disease) is an infectious disease caused by the fungus Ophidiomyces ophidiicola to which all snake species appear to be susceptible. Significant variation has been observed in clinical presentation, progression of disease, and response to treatment, which may be due to genetic variation in the causative agent. Recent phylogenetic analysis based on whole-genome sequencing identified that O . ophidiicola strains from the United States formed a clade distinct from European strains, and that multiple clonal lineages of the clade are present in the United States. The purpose of this study was to design a qPCR-based genotyping assay for O . ophidiicola , then apply that assay to swab-extracted DNA samples to investigate whether the multiple O . ophidiicola clades and clonal lineages in the United States have specific geographic, taxonomic, or temporal predilections. To this end, six full genome sequences of O . ophidiicola representing different clades and clonal lineages were aligned to identify genomic areas shared between subsets of the isolates. Eleven hydrolysis-based Taqman primer-probe sets were designed to amplify selected gene segments and produce unique amplification patterns for each isolate, each with a limit of detection of 10 or fewer copies of the target sequence and an amplification efficiency of 90–110%. The qPCR-based approach was validated using samples from strains known to belong to specific clades and applied to swab-extracted O . ophidiicola DNA samples from multiple snake species, states, and years. When compared to full-genome sequencing, the qPCR-based genotyping assay assigned 75% of samples to the same major clade (Cohen’s kappa = 0.360, 95% Confidence Interval = 0.154–0.567) with 67–77% sensitivity and 88–100% specificity, depending on clade/clonal lineage. Swab-extracted O . ophidiicola DNA samples from across the United States were assigned to six different clonal lineages, including four of the six established lineages and two newly defined groups, which likely represent recombinant strains of O . ophidiicola . Using multinomial logistic regression modeling to predict clade based on snake taxonomic group, state of origin, and year of collection, state was the most significant predictor of clonal lineage. Furthermore, clonal lineage was not associated with disease severity in the most intensely sampled species, the Lake Erie watersnake ( Nerodia sipedon insularum ). Overall, this assay represents a rapid, cost-effective genotyping method for O . ophidiicola that can be used to better understand the epidemiology of ophidiomycosis.

PLoS ONE↗

Identification of duck plague virus by polymerase chain reaction

A polymerase chain reaction (PCR) assay was developed for detecting duck plague virus. A 765-bp EcoRI fragment cloned from the genome of the duck plague vaccine (DP-VAC) virus was sequenced for PCR primer development. The fragment sequence was found by GenBank alignment searches to be similar to the 3′ ends of an undefined open reading frame and the gene for DNA polymerase protein in other herpesviruses. Three of four primer sets were found to be specific for the DP-VAC virus and 100% (7/7) of field isolates but did not amplify DNA from inclusion body disease of cranes virus. The specificity of one primer set was tested with genome templates from other avian herpesviruses, including those from a golden eagle, bald eagle, great horned owl, snowy owl, peregrine falcon, prairie falcon, pigeon, psittacine, and chicken (infectious laryngotracheitis), but amplicons were not produced. Hence, this PCR test is highly specific for duck plague virus DNA. Two primer sets were able to detect 1 fg of DNA from the duck plague vaccine strain, equivalent to five genome copies. In addition, the ratio of tissue culture infectious doses to genome copies of duck plague vaccine virus from infected duck embryo cells was determined to be 1:100, making the PCR assay 20 times more sensitive than tissue culture for detecting duck plague virus. The speed, sensitivity, and specificity of this PCR provide a greatly improved diagnostic and research tool for studying the epizootiology of duck plague.

Avian Diseases↗

A phylogenomic supertree of birds

It has long been appreciated that analyses of genomic data (e.g., whole genome sequencing or sequence capture) have the potential to reveal the tree of life, but it remains challenging to move from sequence data to a clear understanding of evolutionary history, in part due to the computational challenges of phylogenetic estimation using genome-scale data. Supertree methods solve that challenge because they facilitate a divide-and-conquer approach for large-scale phylogeny inference by integrating smaller subtrees in a computationally-efficient manner. Here, we combined information from sequence capture and whole-genome phylogenies using supertree methods. However, available phylogenomic trees had limited overlap so we used taxon-rich (but not phylogenomic) megaphylogenies to weave them together. This allowed us to construct a phylogenomic supertree, with support values, that included 707 bird species (~7% of avian species diversity). We estimated branch lengths using mitochondrial sequence data and we used this to estimate divergence times. Our time-calibrated supertree supports radiation of all three major avian clades (Palaeognathae, Galloanseres, and Neoaves) near the Cretaceous-Paleogene (K-Pg) boundary. The approach we used will permit the continued addition of taxa to this supertree as new phylogenomic data are published, and it could be applied to other taxa as well.

Diversity↗

Implications of historical and contemporary processes on genetic differentiation of a declining boreal songbird: The rusty blackbird

The arrangement of habitat features via historical or contemporary events can strongly influence genomic and demographic connectivity, and in turn affect levels of genetic diversity and resilience of populations to environmental perturbation. The rusty blackbird ( Euphagus carolinus ) is a forested wetland habitat specialist whose population size has declined sharply (78%) over recent decades. The species breeds across the expansive North American boreal forest region, which contains a mosaic of habitat conditions resulting from active natural disturbance regimes and glacial history. We used landscape genomics to evaluate how past and present landscape features have shaped patterns of genetic diversity and connectivity across the species’ breeding range. Based on reduced-representation genomic and mitochondrial DNA, genetic structure followed four broad patterns influenced by both historical and contemporary forces: (1) an east–west partition consistent with vicariance during the last glacial maximum; (2) a potential secondary contact zone between eastern and western lineages at James Bay, Ontario; (3) insular differentiation of birds on Newfoundland; and (4) restricted regional gene flow among locales within western and eastern North America. The presence of genomic structure and therefore restricted dispersal among populations may limit the species’ capacity to respond to rapid environmental change.

Diversity↗

Phylogeographic genetic diversity in the white sucker hepatitis B Virus across the Great Lakes Region and Alberta, Canada

Hepatitis B viruses belong to a family of circular, double-stranded DNA viruses that infect a range of organisms, with host responses that vary from mild infection to chronic infection and cancer. The white sucker hepatitis B virus (WSHBV) was first described in the white sucker (Catostomus commersonii), a freshwater teleost, and belongs to the genus Parahepadnavirus. At present, the host range of WSHBV and its impact on fish health are unknown, and neither genetic diversity nor association with fish health have been studied in any parahepadnavirus. Given the relevance of genomic diversity to disease outcome for the orthohepadnaviruses, we sought to characterize genomic variation in WSHBV and determine how it is structured among watersheds. We identified WSHBV-positive white sucker inhabiting tributaries of Lake Michigan, Lake Superior, Lake Erie (USA), and Lake Athabasca (Canada). Copy number in plasma and in liver tissue was estimated via qPCR. Templates from 27 virus-positive fish were amplified and sequenced using a primer-specific, circular long-range amplification method coupled with amplicon sequencing on the Illumina MiSeq. Phylogenetic analysis of the WSHBV genome identified phylogeographical clustering reminiscent of that observed with human hepatitis B virus genotypes. Notably, most non-synonymous substitutions were found to cluster in the pre-S/spacer overlap region, which is relevant for both viral entry and replication. The observed predominance of p1/s3 mutations in this region is indicative of adaptive change in the polymerase open reading frame (ORF), while, at the same time, the surface ORF is under purifying selection. Although the levels of variation we observed do not meet the criteria used to define sub/genotypes of human and avian hepadnaviruses, we identified geographically associated genome variation in the pre-S and spacer domain sufficient to define five WSHBV haplotypes. This study of WSHBV genetic diversity should facilitate the development of molecular markers for future identification of genotypes and provide evidence in future investigations of possible differential disease outcomes.

Alberta, Michigan, Minnesota, Ohio, Wisconsin↗

A Lake Charr pangenome reveals highly conserved Ohnologs as drivers of phenotypic diversity

Whole-genome duplication (WGD) is hypothesized to spur evolutionary diversification by producing genome-wide duplicate gene sets (Ohnologs) that are initially functionally redundant but can diverge markedly as the effects of relaxed selection accumulate over time. However, the underlying mechanisms remain unclear, in part because genomic studies often reconstruct Ohnolog evolution over millions of years, during which subsequent mutations can obscure deep-time signals. Investigating the relationship between Ohnolog evolution and diversification on a contemporary timescale offers clearer insights. We explore this relationship in Lake Charr ( Salvelinus namaycush ), where ∼10% of genes are retained highly conserved polyploid duplicates following the Salmonid-Specific Fourth Round WGD. Using 31 chromosome-level assemblies of Lake Charr from morphologically and ecologically diverse populations, joined into a pangenome graph, we characterized 189,555 structural variants (SVs) that were significantly less likely to affect genes retained as sequence-conserved Ohnolog pairs, nuancing the hypothesis that gene redundancy, relaxed selection, and functional diversification are intertwined. However, we found that SVs affecting such conserved Ohnologs may be potent drivers of adaptive evolution. Notably, we identified a putative 938-Kb interchromosomal translocation containing 25 genes with highly conserved Ohnologs in a paralogous (but untranslocated) genomic block. This putative translocation appears to have facilitated Ohnolog divergence in ankrd11 and hp , genes putatively linked to craniofacial and lipid metabolic diversity in sympatric Lake Superior morphs. This research reveals that conserved Ohnologs previously presumed to be redundant remain a reservoir for adaptive change.

BioRxiv↗