Low Genomic Diversity Research Articles

Population genomic structure of the sea urchin Diadema africanum, a relevant species in the rocky reef systems across the Macaronesian archipelagos

The sea urchin Diadema africanum is a macro-herbivore found in the rocky reef systems of the West African region and Macaronesian archipelagos. Over several decades, high densities of this species have generated marine barrens in certain areas at the Canary Islands. In contrast, more recently, during the last few years, the species has suffered mass mortality events that continue to the present day. In this study, we used 9,109 Single Nucleotide Polymorphisms (SNPs) and a fragment of a mitochondrial gene to evaluate the species’ population structure, effects of mass mortalities on its diversity, and potential local adaptation across the Canary Islands and Cabo Verde. Our research provides compelling evidence of low genomic diversity and homogeneity across the studied area for neutral markers, along with recent demographic fluctuations. The high connectivity among distant areas likely allows a rapid recovering of the populations from local mortality events. Interestingly, we also observed genomic sub-structure from 405 SNPs identified as candidate loci under selection for seawater temperature. The lack of divergence among distant sites and the low diversity found can be attributed to the species’ divergence from a small ancestral genomic pool, followed by a contemporary demographic expansion, and ongoinggene flow.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H.equinus) and blue antelope, as well as between roan and sable antelope (H.niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Genomic characterisation of the population structure and antibiotic resistance of Salmonella enterica serovar Infantis in Chile, 2009–2022

Multidrug-resistant (MDR) Salmonella Infantis has disseminated worldwide, mainly linked to the consumption of poultry products. Evidence shows dissemination of this pathogen in Chile; however, studies are primarily limited to phenotypic data or involve few isolates. As human cases of Salmonella Infantis infections have substantially increased in recent years, this study aimed to characterise the genomic epidemiology and antimicrobial-resistance profiles of isolates obtained from different sources, aiming to inform effective surveillance and control measures. We sequenced 396 Salmonella Infantis genomes and analysed them with all publicly available genomes of this pathogen from Chile (440 genomes in total), representing isolates from environmental, food, animal, and human sources obtained from 2009 to 2022. Based on bioinformatic and phenotypic methods, we assessed the population structure, dissemination among different niches, and antimicrobial resistance (AMR) profiles of Salmonella Infantis in the country. The genomic and phylogenetic analyses showed that Salmonella Infantis from Chile comprised several clusters of highly related isolates dominated by sequence type 32. The HC20_343 cluster grouped an important proportion of all isolates. This was the only cluster associated with pESI-like megaplasmids, and up to 12 acquired AMR genes/mutations predicted to result in an MDR phenotype. Accordingly, antimicrobial-susceptibility testing revealed a strong concordance between the AMR genetic determinants and their matching phenotypic expression, indicating that a significant proportion of HC20_343 isolates produce extended-spectrum β-lactamases and have intermediate fluoroquinolone resistance. HC20_343 Salmonella Infantis were spread among environmental, animal, food, and human niches, showing a close relationship between isolates from different years and sources, and a low intra-source genomic diversity. Our findings show a widespread dissemination of MDR Salmonella Infantis from the HC20_343 cluster in Chile. The high proportion of isolates with resistance to first-line antibiotics and the evidence of active transmission between the environment, animals, food, and humans highlight the urgency of improved surveillance and control measures in the country. As HC20_343 isolates predominate in the Americas, our results suggest a high prevalence of ESBL-producing Salmonella Infantis with intermediate fluoroquinolone resistance in the continent. Partially supported by the Food and Drug Administration (FDA) of the U.S. Department of Health and Human Services as part of an award, FDU001818, with 30% percentfundedby FDA/HHS; and by Agencia de Investigación y Desarrollo de Chile (ANID) through FONDECYT de Postdoctorado Folio 3230796 and Folio 3210317, FONDECYT Regular Folio 1231082, and ANID-Millennium Science Initiative Program-ICN2021_044.

Low but highly geographically structured genomic diversity of East Asian Eurasian otters and its conservation implications.

Populations of Eurasian otters Lutra lutra, one of the most widely distributed apex predators in Eurasia, have been depleted mainly since the 1950s. However, a lack of information about their genomic diversity and how they are organized geographically in East Asia severely impedes our ability to monitor and conserve them in particular management units. Here, we re-sequenced and analyzed 20 otter genomes spanning continental East Asia, including a population at Kinmen, a small island off the Fujian coast, China. The otters form three genetic clusters (one of L. l. lutra in the north and two of L. l. chinensis in the south), which have diverged in the Holocene. These three clusters should be recognized as three conservation management units to monitor and manage independently. The heterozygosity of the East Asian otters is as low as that of the threatened carnivores sequenced. Historical effective population size trajectories inferred from genomic variations suggest that their low genomic diversity could be partially attributed to changes in the climate since the mid-Pleistocene and anthropogenic intervention since the Holocene. However, no evidence of genetic erosion, mutation load, or high level of inbreeding was detected in the presumably isolated Kinmen Island population. Any future insitu conservation efforts should consider this information for the conservation management units.

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 (Ne = 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.

Clonal expansion and rapid characterization of Klebsiella pneumoniae ST1788, an otherwise uncommon strain spreading in Wales, UK.

A multidrug-resistant strain of Klebsiella pneumoniae (Kp) sequence type (ST) 1788, an otherwise uncommon ST worldwide, was isolated from 65 patients at 11 hospitals and 11 general practices across South and West Wales, UK, between February 2019 and November 2021. A collection of 97 Kp ST1788 isolates (including 94 from Wales) was analysed to investigate the diversity and spread across Wales and to identify molecular marker(s) to aid development of a strain-specific real-time PCR. Whole genome sequencing (WGS) was performed with Illumina technology and the data were used to perform phylogenetic analyses. Pan-genome analysis of further Kp genome collections was used to identify an ST1788-specific gene target; a real-time PCR was then validated against a panel of 314 strains and 218 broth-enriched screening samples. Low genomic diversity was demonstrated amongst the 94 isolates from Wales. Evidence of spread within and across healthcare facilities was found. A yersiniabactin locus and the KL2 capsular locus were identified in 85/94 (90.4 %) and 94/94 (100 %) genomes respectively; bla SHV-232, bla TEM-1, bla CTX-M-15 and bla OXA-1 were simultaneously carried by 86/94 (91.5 %) isolates; 4/94 (4.3 %) isolates also carried bla OXA-48 carbapenemase. Aminoglycoside and fluoroquinolone resistance markers were found in 94/94 (100 %) and 86/94 (91.5 %) isolates respectively. The ST1788-specific real-time PCR was 100 % sensitive and specific. Our analyses demonstrated recent clonal expansion and spread of Kp ST1788 in the community and across healthcare facilities in South and West Wales with isolates carrying well-defined antimicrobial resistance and virulence markers. An ST1788-specific marker was also identified, enabling rapid and reliable preliminary characterization of isolates by real-time PCR. This study confirms the utility of WGS in investigating novel strains and in aiding proactive implementation of molecular tools to assist infection control specialists.

The utility of SARS-CoV-2 genomic data for informative clustering under different epidemiological scenarios and sampling.

Clustering pathogen sequence data is a common practice in epidemiology to gain insights into the genetic diversity and evolutionary relationships among pathogens. We can find groups of cases with a shared transmission history and common origin, as well as identifying transmission hotspots. Motivated by the experience of clustering SARS-CoV-2 cases using whole genome sequence data during the COVID-19 pandemic to aid with public health investigation, we investigated how differences in epidemiology and sampling can influence the composition of clusters that are identified. We performed genomic clustering on simulated SARS-CoV-2 outbreaks produced with different transmission rates and levels of genomic diversity, along with varying the proportion of cases sampled. In single outbreaks with a low transmission rate, decreasing the sampling fraction resulted in multiple, separate clusters being identified where intermediate cases in transmission chains are missed. Outbreaks simulated with a high transmission rate were more robust to changes in the sampling fraction and largely resulted in a single cluster that included all sampled outbreak cases. When considering multiple outbreaks in a sampled jurisdiction seeded by different introductions, low genomic diversity between introduced cases caused outbreaks to be merged into large clusters. If the transmission and sampling fraction, and diversity between introductions was low, a combination of the spurious break-up of outbreaks and the linking of closely related cases in different outbreaks resulted in clusters that may appear informative, but these did not reflect the true underlying population structure. Conversely, genomic clusters matched the true population structure when there was relatively high diversity between introductions and a high transmission rate. Differences in epidemiology and sampling can impact our ability to identify genomic clusters that describe the underlying population structure. These findings can help to guide recommendations for the use of pathogen clustering in public health investigations.

Challenges Associated with Investigating Salmonella Enteritidis with Low Genomic Diversity in New York State: The Impact of Adjusting Analytical Methods and Correlation with Epidemiological Data.

Defining investigation-worthy genomic clusters among strains of Salmonella Enteritidis is challenging because of their highly clonal nature. We investigated a cluster identified by core genome multilocus sequence typing (cgMLST) consisting of 265 isolates with isolation dates spanning two and a half years. This cluster experienced chaining, growing to a range of 14 alleles. The volume of isolates and broad allele range of this cluster made it difficult to ascertain whether it represented a common-source outbreak. We explored laboratory-based methods to subdivide and refine this cluster. These methods included using cgMLST with a narrower allele range, whole genome multilocus sequence typing (wgMLST) and high-quality single-nucleotide polymorphism (hqSNP) analysis. At each analysis level, epidemiologists retroactively reviewed exposures, geography, and temporality for potential commonalities. Lowering the threshold to 0 alleles using cgMLST proved an effective method to refine this analysis, resulting in this large cluster being subdivided into 34 smaller clusters. Additional analysis by wgMLST and hqSNP provided enhanced cluster resolution, with the majority of clusters being further refined. These analysis methods combined with more stringent allele thresholds and layering of epidemiologic data proved useful in helping to subdivide this large cluster into actionable subclusters.

Phylogenic position and low genomic diversity of 'Candidatus Rickettsia kotlanii' inferred by complete genome sequences of two Japanese isolates.

Many Rickettsia species of the spotted fever group (SFG) cause tick-borne diseases known as "spotted fever". One of the candidate SFG Rickettsia species is 'Candidatus Rickettsia kotlanii', which was first detected in Haemaphysalis concinna in Hungary in 2006. However, its precise phylogenetic position in the SFG is not clear because only single-gene sequence-based phylogenetic analyses were performed using very limited genes. Here, we present the complete genome sequences of two Japanese 'Ca. R. kotlanii' isolates, which differed only by a 135-bp insertion/deletion (InDel). Using these genomes and publicly available whole genome sequences of other Rickettsia species, the precise phylogenetic position of 'Ca. R. kotlanii' in Rickettsia was determined to be in a clade of the SFG. The phylogenetic relationships and average nucleotide identity of 'Ca. R. kotlanii' relative to the other species indicated that 'Ca. R. kotlanii' is an independent taxon in the SFG. Notably, although the genomes of the two isolates were almost identical, the isolates were obtained from different tick species in different regions and years, suggesting extremely low genomic diversity in 'Ca. R. kotlanii'. While the genome of 'Ca. R. kotlanii' is the smallest in the transitional group and SFG Rickettsia sequenced to date, we identified genes uniquely present or absent in 'Ca. R. kotlanii', but most were apparently degraded. Therefore, analyses of differences at the sequence (single nucleotide polymorphisms and small InDels) or gene expression level will be required to understand the functional or physiological features unique to 'Ca. R. kotlanii'. This article is protected by copyright. All rights reserved.

Investigation of hospital discharge cases and SARS-CoV-2 introduction into Lothian care homes

The first epidemic wave of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in Scotland resulted in high case numbers and mortality in care homes. In Lothian, over a third of care homes reported an outbreak while there was limited testing of hospital patients discharged to care homes. Investigate hospital discharges as a source of SARS-CoV-2 introduction into care homes during the first epidemic wave. A clinical review was performed for all discharges from hospitals to care homes starting 1st March 2020 to 31st May 2020. Episodes were ruled out based on coronavirus disease (COVID-19) test history, clinical assessment at discharge, whole genome sequencing (WGS) data and an infectious period of 14 days. Clinical samples were processed for WGS, and consensus genomes generated were used for analysis by cluster investigation and virus epidemiological tool (CIVET). Patient timelines were obtained using electronic hospital records. In total 787 hospital discharges to care homes were identified. Out of these 776 (99%) were ruled out for hospital discharge introduction. However, for 10 episodes the results were inconclusive as there was low genomic diversity in consensus genomes or no sequencing data. Only one discharge episode had a genomic, time and location link to positive cases during hospital admission leading to 10 further positive cases in the care home. Majority of hospital discharges were ruled out for introduction into Lothian care homes highlighting the importance of screening all new admissions when faced with a novel emerging virus and no vaccine available.

SARS-CoV-2 Genomic Diversity in Households Highlights the Challenges of Sequence-Based Transmission Inference.

The reliability of sequence-based inference of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission is not clear. Sequence data from infections among household members can define the expected genomic diversity of a virus along a defined transmission chain. SARS-CoV-2 cases were identified prospectively among 2,369 participants in 706 households. Specimens with a reverse transcription-PCR cycle threshold of ≤30 underwent whole-genome sequencing. Intrahost single-nucleotide variants (iSNV) were identified at a ≥5% frequency. Phylogenetic trees were used to evaluate the relationship of household and community sequences. There were 178 SARS-CoV-2 cases in 706 households. Among 147 specimens sequenced, 106 yielded a whole-genome consensus with coverage suitable for identifying iSNV. Twenty-six households had sequences from multiple cases within 14 days. Consensus sequences were indistinguishable among cases in 15 households, while 11 had ≥1 consensus sequence that differed by 1 to 2 mutations. Sequences from households and the community were often interspersed on phylogenetic trees. Identification of iSNV improved inference in 2 of 15 households with indistinguishable consensus sequences and in 6 of 11 with distinct ones. In multiple-infection households, whole-genome consensus sequences differed by 0 to 1 mutations. Identification of shared iSNV occasionally resolved linkage, but the low genomic diversity of SARS-CoV-2 limits the utility of "sequence-only" transmission inference. IMPORTANCE We performed whole-genome sequencing of SARS-CoV-2 from prospectively identified cases in three longitudinal household cohorts. In a majority of multi-infection households, SARS-CoV-2 consensus sequences were indistinguishable, and they differed by 1 to 2 mutations in the rest. Importantly, even with modest genomic surveillance of the community (3 to 5% of cases sequenced), it was not uncommon to find community sequences interspersed with household sequences on phylogenetic trees. Identification of shared minority variants only occasionally resolved these ambiguities in transmission linkage. Overall, the low genomic diversity of SARS-CoV-2 limits the utility of "sequence-only" transmission inference. Our work highlights the need to carefully consider both epidemiologic linkage and sequence data to define transmission chains in households, hospitals, and other transmission settings.

African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity

African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.

A pipeline for effectively developing highly polymorphic simple sequence repeats markers based on multi-sample genomic data.

Simple sequence repeats (SSRs) are widely used genetic markers in ecology, evolution, and conservation even in the genomics era, while a general limitation to their application is the difficulty of developing polymorphic SSR markers. Next‐generation sequencing (NGS) offers the opportunity for the rapid development of SSRs; however, previous studies developing SSRs using genomic data from only one individual need redundant experiments to test the polymorphisms of SSRs. In this study, we designed a pipeline for the rapid development of polymorphic SSR markers from multi‐sample genomic data. We used bioinformatic software to genotype multiple individuals using resequencing data, detected highly polymorphic SSRs prior to experimental validation, significantly improved the efficiency and reduced the experimental effort. The pipeline was successfully applied to a globally threatened species, the brown eared‐pheasant (Crossoptilon mantchuricum), which showed very low genomic diversity. The 20 newly developed SSR markers were highly polymorphic, the average number of alleles was much higher than the genomic average. We also evaluated the effect of the number of individuals and sequencing depth on the SSR mining results, and we found that 10 individuals and ~10X sequencing data were enough to obtain a sufficient number of polymorphic SSRs, even for species with low genetic diversity. Furthermore, the genome assembly of NGS data from the optimal number of individuals and sequencing depth can be used as an alternative reference genome if a high‐quality genome is not available. Our pipeline provided a paradigm for the application of NGS technology to mining and developing molecular markers for ecological and evolutionary studies.

Genome analysis-based studies on bacterial genetic diversity

There are a huge number of bacterial species on earth, and a huge intra-species genomic diversity are also observed in many bacteria. The high ability of bacteria to acquire foreign DNA and the presence of various mobile genetic elements contribute the generation of such genomic diversity. During the biochemical and genetic analysis of a Pseudomonas aeruginosa toxin, called cytotoxin, and its converting phage, which I first engaged in my research carrier, I became very interested in the genetic diversity of bacteria and mobile genetic elements such as bacteriophages, and realized the usefulness and power of genome analysis. Since then, I have been involved in genome analyses of various pathogenic bacteria such as enterohemorrhagic Escherichia coli (EHEC), commensal bacteria of human and other animals, and bacteria or bacterial communities in natural environments. I was so lucky that I jumped in this research field at the very begging of genome analyses and experienced a very exciting time of surprisingly rapid advance in genome sequencing technologies which revolutionized a wide range of biology. In this article, I first review the main findings which our group obtained from the genome analyses on the P. aeruginosa cytotoxin converting phage and those on the evolution and genomic diversity of EHEC and related bacteria. The results of our analyses of Rickettsiaceae family genomes, which show surprisingly very low genomic diversity, and genome sequence-based analyses of an intrahospital bacterial outbreak and within-host genomic diversity are also summarized.

Genome of the estuarine oyster provides insights into climate impact and adaptive plasticity

Understanding the roles of genetic divergence and phenotypic plasticity in adaptation is central to evolutionary biology and important for assessing adaptive potential of species under climate change. Analysis of a chromosome-level assembly and resequencing of individuals across wide latitude distribution in the estuarine oyster (Crassostrea ariakensis) revealed unexpectedly low genomic diversity and population structures shaped by historical glaciation, geological events and oceanographic forces. Strong selection signals were detected in genes responding to temperature and salinity stress, especially of the expanded solute carrier families, highlighting the importance of gene expansion in environmental adaptation. Genes exhibiting high plasticity showed strong selection in upstream regulatory regions that modulate transcription, indicating selection favoring plasticity. Our findings suggest that genomic variation and population structure in marine bivalves are heavily influenced by climate history and physical forces, and gene expansion and selection may enhance phenotypic plasticity that is critical for the adaptation to rapidly changing environments.

Genome skims analysis of betel palms (Areca spp., Arecaceae) and development of a profiling method to assess their plastome diversity.

The betel nut (Areca catechu L., Arecaceae) is a monoecious cultivated palm tree that is widespread in tropical regions. It is mainly cultivated for producing areca nuts, from which seeds are extracted and chewed by local populations principally in the Indo-Pacific region. Seeds contain alkaloids which are central nervous system stimulants and are highly addictive. Wild relatives of the betel nut are distributed in South Asia and Australasia, with ca. 40-50 Areca species currently recognized. The geographic origin(s) of the betel nut and its subsequent diffusion and diversification remains poorly documented. Here, a genome skimming approach was applied to screen nucleotidic variation in the most abundant genomic regions. Low coverage sequencing data allowed us to assemble full plastomes, mitochondrial regions (either full mitogenomes or the full set of mitochondrial genes) and the nuclear ribosomal DNA cluster for nine representatives of the Areca genus collected in the field and herbarium collections (including a 182-years old specimen collected during the Dumont d'Urville's expedition). These three genomic compartments provided similar phylogenetic signals, and revealed very low genomic diversity in our sample of cultivated betel nut. We finally developed a genotyping method targeting 34 plastid DNA microsatellites. This plastome profiling approach is useful for tracing the spread of matrilineages, and in combination with nuclear genomic data, can resolve the history of the betel nut. Our method also proves to be efficient for analyzing herbarium specimens, even those collected >100years ago.

Antimicrobial resistance and genomic characterization of Salmonella Dublin isolates in cattle from the United States.

Salmonella enterica subspecies enterica serotype Dublin is a host-adapted serotype in cattle, associated with enteritis and systemic disease. The primary clinical manifestation of Salmonella Dublin infection in cattle, especially calves, is respiratory disease. While rare in humans, it can cause severe illness, including bacteremia, with hospitalization and death. In the United States, S. Dublin has become one of the most multidrug-resistant serotypes. The objective of this study was to characterize S. Dublin isolates from sick cattle by analyzing phenotypic and genotypic antimicrobial resistance (AMR) profiles, the presence of plasmids, and phylogenetic relationships. S. Dublin isolates (n = 140) were selected from submissions to the NVSL for Salmonella serotyping (2014–2017) from 21 states. Isolates were tested for susceptibility against 14 class-representative antimicrobial drugs. Resistance profiles were determined using the ABRicate with Resfinder and NCBI databases, AMRFinder and PointFinder. Plasmids were detected using ABRicate with PlasmidFinder. Phylogeny was determined using vSNP. We found 98% of the isolates were resistant to more than 4 antimicrobials. Only 1 isolate was pan-susceptible and had no predicted AMR genes. All S. Dublin isolates were susceptible to azithromycin and meropenem. They showed 96% resistance to sulfonamides, 97% to tetracyclines, 95% to aminoglycosides and 85% to beta-lactams. The most common AMR genes were: sulf2 and tetA (98.6%), aph(6)-Id (97.9%), aph(3’’)-Ib, (97.1%), floR (94.3%), and blaCMY-2 (85.7%). All quinolone resistant isolates presented mutations in gyrA. Ten plasmid types were identified among all isolates with IncA/C2, IncX1, and IncFII(S) being the most frequent. The S. Dublin isolates show low genomic genetic diversity. This study provided antimicrobial susceptibility and genomic insight into S. Dublin clinical isolates from cattle in the U.S. Further sequence analysis integrating food and human origin S. Dublin isolates may provide valuable insight on increased virulence observed in humans.

Genome-informed approach to identify genetic determinants of Flavobacterium psychrophilum phage susceptibility.

The fish pathogen Flavobacterium psychrophilum infects farmed salmonids worldwide, and application of bacteriophages has been suggested for controlling disease outbreaks in aquaculture. Successful application of phages requires detailed knowledge about the variability in phage susceptibility of the host communities. In this study, we analysed the genetic diversity of F. psychrophilum hosts and phages from the Baltic Sea area to identify genetic determinants of phage-host interaction patterns. A host range analysis of 103 phages tested against 177 F. psychrophilum strains (18 231 phage-host interactions) identified nine phage clusters, infecting from 10% to 91% of the strain collection. The core genome-based comparison of 35 F. psychrophilum isolates revealed an extremely low overall genomic diversity (>99.5% similarity). However, a small subset of 16 ORFs, including genes involved in the type IX secretion system (T9SS), gliding motility and hypothetical cell-surface related proteins, exhibited a highly elevated genetic diversity. These specific genetic variations were linked to variability in phage infection patterns obtained from experimental studies, indicating that these genes are key determinants of phage susceptibility. These findings provide novel insights on the molecular mechanisms determining phage susceptibility in F. psychrophilum and emphasizes the importance of phages as drivers of core genomic diversity in this pathogen.

Population genomic, climatic and anthropogenic evidence suggest the role of human forces in endangerment of green peafowl (Pavo muticus).

Both anthropogenic impacts and historical climate change could contribute to population decline and species extinction, but their relative importance is still unclear. Emerging approaches based on genomic, climatic and anthropogenic data provide a promising analytical framework to address this question. This study applied such an integrative approach to examine potential drivers for the endangerment of the green peafowl (Pavo muticus). Several demographic reconstructions based on population genomes congruently retrieved a drastic population declination since the mid-Holocene. Furthermore, a comparison between historical and modern genomes suggested genetic diversity decrease during the last 50 years. However, climate-based ecological niche models predicted stationary general range during these periods and imply the little impact of climate change. Further analyses suggested that human disturbance intensities were negatively correlated with the green peafowl's effective population sizes and significantly associated with its survival status (extirpation or persistence). Archaeological and historical records corroborate the critical role of humans, leaving the footprint of low genomic diversity and high inbreeding in the survival populations. This study sheds light on the potential deep-time effects of human disturbance on species endangerment and offers a multi-evidential approach in examining underlying forces for population declines.