Changes in CpG Methylation of the Vitellogenin 1 Promoter in Adult Male Zebrafish after Exposure to 17α-Ethynylestradiol.

Numerous pharmaceutical and industrial chemicals are classified as endocrine-disrupting chemicals (EDCs) that interfere with hormonal homeostasis, leading to developmental disorders and other pathologies. The synthetic estrogen 17α-ethynylestradiol (EE2) is used in oral contraceptives and other hormone therapies. EE2 and other estrogens are inadvertently introduced into aquatic environments through municipal wastewater and agricultural effluents. Exposure of male fish to estrogens increases expression of the egg yolk precursor protein vitellogenin (Vtg), which is used as a molecular marker of exposure to estrogenic EDCs. The mechanisms behind Vtg induction are not fully known, and we hypothesized that it is regulated via DNA methylation. Adult zebrafish were exposed to either dimethyl sulfoxide or 20 ng/L EE2 for 14 days. Messenger RNA (mRNA) expression and DNA methylation were assessed in male zebrafish livers at 0, 0.25, 0.5, 1, 4, 7, and 14 days of exposure; and those of females were assessed at 13 days (n ≥ 4/group/time point). To test the persistence of any changes, we included a recovery group that received EE2 for 7 days and did not receive any for the following 7 days, in the total 14-day study. Methylation of DNA at the vtg1 promoter was assessed with targeted gene bisulfite sequencing in livers of adult male and female zebrafish. A significant increase in vtg1 mRNA was observed in the EE2-exposed male fish as early as 6 h. Interestingly, DNA methylation changes were observed at 4 days. Decreases in the overall methylation of the vtg1 promoter in exposed males resulted in levels comparable to those in female controls, suggesting feminization. Importantly, DNA methylation levels in males remained significantly impacted after 7 days post-EE2 removal, unlike mRNA levels. These data identify an epigenetic mark of feminization that may serve as an indicator of not only estrogenic exposure but also previous exposure to EE2. Environ Toxicol Chem 2024;43:1547-1556. © 2024 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

"Evaluating the Benefits and Limits of Multiple Displacement Amplification with Whole-Genome Oxford Nanopore Sequencing".

Multiple Displacement Amplification (MDA) outperforms conventional PCR in long fragment and whole genome amplification which makes it attractive to couple with long-read sequencing of samples with limited quantities of DNA to obtain improved genome assemblies. Here, we explore the efficacy and limits of MDA for genome sequence assembly using Oxford Nanopore Technologies (ONT) rapid library preparations and minION sequencing. We successfully generated almost complete genome sequences for all organisms examined, including Cryptosporidium meleagridis, Staphylococcus aureus, Enterococcus faecium, and Escherichia coli, with the ability to generate high-quality data from samples starting with only 0.025 ng of total DNA. Controlled sheared DNA samples exhibited a distinct pattern of size-increase after MDA, which may be associated with the amplification of long, low-abundance fragments present in the assay, as well as generating concatemeric sequences during amplification. To address concatemers, we developed a computational pipeline (CADECT: Concatemer Detection Tool) to identify and remove putative concatemeric sequences. This study highlights the efficacy of MDA in generating high-quality genome assemblies from limited amounts of input DNA. Also, the CADECT pipeline effectively mitigated the impact of concatemeric sequences, enabling the assembly of contiguous sequences even in cases where the input genomic DNA was degraded. These results have significant implications for the study of organisms that are challenging to culture in vitro, such as Cryptosporidium, and for expediting critical results in clinical settings with limited quantities of available genomic DNA.

Improved Equine Fecal Microbiome Characterization Using Target Enrichment by Hybridization Capture.

GITDs are among the most common causes of death in adult and young horses in the United States (US). Previous studies have indicated a connection between GITDs and the equine gut microbiome. However, the low taxonomic resolution of the current microbiome sequencing methods has hampered the identification of specific bacterial changes associated with GITDs in horses. Here, we have compared TEHC, a new approach for 16S rRNA gene selection and sequencing, with conventional 16S rRNA gene amplicon sequencing for the characterization of the equine fecal microbiome. Both sequencing approaches were used to determine the fecal microbiome of four adult horses and one commercial mock microbiome. Our results show that TEHC yielded significantly more operational taxonomic units (OTUs) than conventional 16S amplicon sequencing when the same number of reads were used in the analysis. This translated into a deeper and more accurate characterization of the fecal microbiome when the samples were sequenced with TEHC according to the relative abundance analysis. Alpha and beta diversity metrics corroborated these findings and demonstrated that the microbiome of the fecal samples was significantly richer when sequenced with TEHC compared to 16S amplicon sequencing. Altogether, our study suggests that the TEHC strategy provides a more extensive characterization of the fecal microbiome of horses than the current alternative based on the PCR amplification of a portion of the 16S rRNA gene.

Long-term gut colonization with ESBL-producing Escherichia coli in participants without known risk factors from the southeastern United States.

We evaluated gut carriage of extended spectrum beta lactamase producing Enterobacteriaceae (ESBL-E) in southeastern U.S. residents without recent in-patient healthcare exposure. Study enrollment was January 2021-February 2022 in Athens, Georgia, U.S. and included a diverse population of 505 adults plus 50 child participants (age 0-5). Based on culture-based screening of stool samples, 4.5% of 555 participants carried ESBL-Es. This is slightly higher than reported in studies conducted 2012-2015, which found carriage rates of 2.5-3.9% in healthy U.S. residents. All ESBL-E confirmed isolates (n=25) were identified as Escherichia coli. Isolates belonged to 11 sequence types, with 48% classified as ST131. Ninety six percent of ESBL-E isolates carried a blaCTX-M gene. Isolated ESBL-Es frequently carried virulence genes as well as multiple classes of antibiotic resistance genes. Long-term colonization was common, with 64% of ESBL-E positive participants testing positive when rescreened three months later. One participant yielded isolates belonging to two different E. coli sequence types that carried blaCTX-M-1 genes on near-identical plasmids, suggesting intra-gut plasmid transfer. Isolation of E. coli on media without antibiotics revealed that ESBL-E. coli typically made up a minor fraction of the overall gut E. coli population, although in some cases they were the dominant strain. ESBL-E carriage was not associated with a significantly different stool microbiome composition. However, some microbial taxa were differentially abundant in ESBL-E carriers. Together, these results suggest that a small subpopulation of US residents are long-term, asymptomatic carriers of ESBL-Es, and may serve as an important reservoir for community spread of these ESBL genes.

A new chromosome-level genome assembly and annotation of Cryptosporidium meleagridis.

Cryptosporidium spp. are medically and scientifically relevant protozoan parasites that cause severe diarrheal illness in infants and immunosuppressed populations as well as animals. Although most human Cryptosporidium infections are caused by C. parvum and C. hominis, there are several other human-infecting species including C. meleagridis, which is commonly observed in developing countries. Here, we polished and annotated a long-read genome sequence assembly for C. meleagridis TU1867, a species which infects birds and humans. The genome sequence was generated using a combination of whole genome amplification (WGA) and long-read Oxford Nanopore Technologies sequencing. The assembly was then polished with Illumina data. The chromosome-level genome assembly is 9.2 Mbp with a contig N50 of 1.1 Mb. Annotation revealed 3,923 protein-coding genes. A BUSCO analysis indicates a completeness of 96.6% (n=446), including 430 (96.4%) single-copy and 1 (0.224%) duplicated apicomplexan conserved gene(s). The new C. meleagridis genome assembly is nearly gap-free and provides a valuable new resource for the Cryptosporidium community and future studies on evolution and host-specificity.

River drainage reorganization and reticulate evolution in the Two-Lined Salamander (Eurycea bislineata) species complex.

The origin and eventual loss of biogeographic barriers can create alternating periods of allopatry and secondary contact, facilitating gene flow among distinct metapopulations and generating reticulate evolutionary histories that are not adequately described by a bifurcating evolutionary tree. One such example may exist in the two-lined salamander (Eurycea bislineata) species complex, where discordance among morphological and molecular datasets has created a "vexing taxonomic challenge". Previous phylogeographic analyses of mitochondrial DNA (mtDNA) suggested that the reorganization of Miocene paleodrainages drove vicariance and dispersal, but the inherent limitations of a single-locus dataset precluded the evaluation of subsequent gene flow. Here, we generate triple-enzyme restriction site-associated DNA sequencing (3RAD) data for >100 individuals representing all major mtDNA lineages and use a suite of complementary methods to demonstrate that discordance among earlier datasets is best explained by a reticulate evolutionary history influenced by river drainage reorganization. Systematics of such groups should acknowledge these complex histories and relationships that are not strictly hierarchical.

New T2T assembly of Cryptosporidium parvum IOWA annotated with reference genome gene identifiers.

Cryptosporidium parvum is a significant pathogen causing gastrointestinal infections in humans and animals, that is spread through the ingestion of contaminated food and water. Despite its global impact on public health, generating a C. parvum genome sequence has always been challenging due to a lack of in vitro cultivation systems and challenging sub-telomeric gene families. A gapless telomere to telomere genome assembly has been created for Cryptosporidium parvum IOWA obtained from Bunch Grass Farms, named here as CpBGF. There are 8 chromosomes that total 9,259,183 bp. The new hybrid assembly which was generated with Illumina and Oxford Nanopore resolves complex sub-telomeric regions of chromosomes 1, 7 and 8. To facilitate ease of use and consistency with the literature, whenever possible, chromosomes have been oriented and genes in this annotation have been given the same gene IDs used in the current reference genome sequence generated in 2004. The annotation of this assembly utilized considerable RNA expression evidence, thus, untranslated regions, long noncoding RNAs and antisense RNAs are annotated. The CpBGF genome assembly serves as a valuable resource for understanding the biology, pathogenesis, and transmission of C. parvum, and it facilitates the development of diagnostics, drugs, and vaccines against cryptosporidiosis.

The gut-brain axis mediates bacterial driven modulation of reward signaling.

OBJECTIVE: Our goal is to investigate if microbiota composition modulates reward signaling and assess the role of the vagus in mediating microbiota to brain communication. METHODS: Male germ-free Fisher rats were colonized with gastrointestinal contents from chow (low fat (LF) ConvLF) or HF (ConvHF) fed rats. RESULTS: Following colonization, ConvHF rats consumed significantly more food than ConvLF animals. ConvHF rats displayed lower feeding-induced extracellular DOPAC levels (a metabolite of dopamine) in the Nucleus Accumbens (NAc) as well as reduced motivation for HF foods compared to ConvLF rats. Dopamine receptor 2 (DDR2) expression levels in the NAc were also significantly lower in ConvHF animals. Similar deficits were observed in conventionally raised HF fed rats, showing that diet-driven alteration in reward can be initiated via microbiota. Selective gut to brain deafferentation restored DOPAC levels, DRD2 expression, and motivational drive in ConvHF rats. CONCLUSIONS: We concluded from these data that a HF-type microbiota is sufficient to alter appetitive feeding behavior and that bacteria to reward communication is mediated by the vagus nerve.

Phylogeography of the blacklegged tick (Ixodes scapularis) throughout the USA identifies candidate loci for differences in vectorial capacity.

The blacklegged tick (Ixodes scapularis (Journal of the Academy of Natural Sciences of Philadelphia, 1821, 2, 59)) is a vector of Borrelia burgdorferi sensu stricto (s.s.) (International Journal of Systematic Bacteriology, 1984, 34, 496), the causative bacterial agent of Lyme disease, part of a slow-moving epidemic of Lyme borreliosis spreading across the northern hemisphere. Well-known geographical differences in the vectorial capacity of these ticks are associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographical studies of these ticks have been restricted to relatively few populations or few genetic loci. Here we present the most comprehensive phylogeographical study of genome-wide markers in I. scapularis, conducted by using 3RAD (triple-enzyme restriction-site associated sequencing) and surveying 353 ticks from 33 counties throughout the species' range. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify five spatially associated genetic clusters of I. scapularis. In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme-endemic regions. Finally, we identify 10 variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome-scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex.

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M âˆ¼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

A high-quality Ixodes scapularis genome advances tick science.

Ixodes spp. and related ticks transmit prevalent infections, although knowledge of their biology and development of anti-tick measures have been hindered by the lack of a high-quality genome. In the present study, we present the assembly of a 2.23-Gb Ixodes scapularis genome by sequencing two haplotypes within one individual, complemented by chromosome-level scaffolding and full-length RNA isoform sequencing, yielding a fully reannotated genome featuring thousands of new protein-coding genes and various RNA species. Analyses of the repetitive DNA identified transposable elements, whereas the examination of tick-associated bacterial sequences yielded an improved Rickettsia buchneri genome. We demonstrate how the Ixodes genome advances tick science by contributing to new annotations, gene models and epigenetic functions, expansion of gene families, development of in-depth proteome catalogs and deciphering of genetic variations in wild ticks. Overall, we report critical genetic resources and biological insights impacting our understanding of tick biology and future interventions against tick-transmitted infections.

A draft of the genome of the Gulf Coast tick, Amblyomma maculatum.

The Gulf Coast tick, Amblyomma maculatum, inhabits the Southeastern states of the USA bordering the Gulf of Mexico, Mexico, and other Central and South American countries. More recently, its U.S. range has extended West to Arizona and Northeast to New York state and Connecticut. It is a vector of Rickettsia parkeri and Hepatozoon americanum. This tick species has become a model to study tick/Rickettsia interactions. To increase our knowledge of the basic biology of A. maculatum we report here a draft genome of this tick and an extensive functional classification of its proteome. The DNA from a single male tick was used as a genomic source, and a 10X genomics protocol determined 28,460 scaffolds having equal or more than 10 Kb, totaling 1.98 Gb. The N50 scaffold size was 19,849 Kb. The BRAKER pipeline was used to find the protein-coding gene boundaries on the assembled A. maculatum genome, discovering 237,921 CDS. After trimming and classifying the transposable elements, bacterial contaminants, and truncated genes, a set of 25,702 were annotated and classified as the core gene products. A BUSCO analysis revealed 83.4% complete BUSCOs. A hyperlinked spreadsheet is provided, allowing browsing of the individual gene products and their matches to several databases.

Tissue Distribution of Mercury in the Bodies of Wild American Alligators (Alligator mississippiensis) from a Coastal Marsh in Louisiana (USA).

Total mercury (THg) concentrations were measured in wild alligators inhabiting a coastal marsh in southern Louisiana, to determine the tissue distribution of THg among various body organs and tissue compartments. Concentrations of THg in claws and dermal tail scutes were compared to those in blood, brain, gonad, heart, kidney, liver, and skeletal muscle to determine if the former tissues, commonly available by non-lethal sampling, could be used as measures of body burdens in various internal organs. Mercury was found in all body organs and tissue compartments. However, overall, THg concentrations measured in alligators were below the FDA action level for fish consumption and were comparable to previous data reported from southwestern Louisiana. Our results suggest consumption of meat from alligators found in this region may be of little public health concern. However, the extended period of time between sampling (in this study) and the present-day highlight the need for continuous, additional, and more recent sampling to ensure consumer safety. Total mercury concentrations were highest in the kidney (3.18 ± 0.69 mg/kg dw) and liver (3.12 ± 0.76 mg/kg dw). THg levels in non-lethal samples (blood, claws, and dermal tail scutes) were positively correlated with all tissue THg concentrations (blood: R2 = 0.513-0.988; claw: R2 = 0.347-0.637, scutes: R2 = 0.333-0.649). Because THg concentrations from blood, claws, and scutes were correlated with those of the internal organs, non-lethal sampling methods may be a viable method of estimating levels of THg in other body tissues.

Estimating Movement Rates Between Eurasian and North American Birds That Are Vectors of Avian Influenza.

Avian influenza (AI) is a zoonotic disease that will likely be involved in future pandemics. Because waterbird movements are difficult to quantify, determining the host-specific risk of Eurasian-origin AI movements into North America is challenging. We estimated relative rates of movements, based on long-term evolutionary averages of gene flow, between Eurasian and North American waterbird populations to obtain bidirectional baseline rates of the intercontinental movements of these AI hosts. We used population genomics and coalescent-based demographic models to obtain these gene-flow-based movement estimates. Inferred rates of movement between these continental populations varies greatly among species. Within dabbling ducks, gene flow, relative to effective population size, varies from ∼3 to 24 individuals/generation between Eurasian and American wigeons (Mareca penelope and Mareca americana) to ∼100-300 individuals/generation between continental populations of northern pintails (Anas acuta). These are evolutionary long-term averages and provide a solid foundation for understanding the relative risks of each of these host species in potential intercontinental AI movements. We scale these values to census size for evaluation in that context. In addition to being AI hosts, many of these bird species are also important in the subsistence diets of Alaskans, increasing the risk of direct bird-to-human exposure to Eurasian-origin AI virus. We contrast species-specific rates of intercontinental movements with the importance of each species in Alaskan diets to understand the relative risk of these taxa to humans. Assuming roughly equivalent AI infection rates among ducks, greater scaup (Aythya marila), mallard (Anas platyrhynchos), and northern pintail (Anas acuta) were the top three species presenting the highest risks for intercontinental AI movement both within the natural system and through exposure to subsistence hunters. Improved data on AI infection rates in this region could further refine these relative risk assessments. These directly comparable, species-based intercontinental movement rates and relative risk rankings should help in modeling, monitoring, and mitigating the impacts of intercontinental host and AI movements.

Estimación de las tasas de movimiento entre aves euroasiáticas y norteamericanas que son vectores de la influenza aviar. La influenza aviar es una enfermedad zoonótica que probablemente estará involucrada en futuras pandemias. Debido a que los movimientos de aves acuáticas son difíciles de cuantificar, La determinación del riesgo específico de hospedador de los movimientos de influenza aviar de origen euroasiático en América del Norte es un desafío. Se estimaron las tasas relativas de movimientos, sobre la base de promedios evolutivos a largo plazo del flujo de genes, entre las poblaciones de aves acuáticas euroasiáticas y norteamericanas para obtener tasas de referencia bidireccionales de los movimientos intercontinentales de estos huéspedes de influenza aviar. Se utilizó genómica de poblaciones y modelos demográficos basados en la teoría de la coalescencia para obtener estas estimaciones de movimiento basadas en el flujo de genes. Las tasas inferidas de movimiento entre estas poblaciones continentales varían mucho entre especies. Dentro de los patos chapuceros, el flujo de genes, en relación con el tamaño efectivo de la población, varía aproximadamente de 3 a 24 individuos/generación entre los silbones europeos y americanos (Mareca penelope y Mareca americana) hasta aproximadamente entre 100 a 300 individuos/generación entre poblaciones continentales de ánades rabudos (Anas acuta). Estos son promedios evolutivos a largo plazo y proporcionan una base sólida para comprender los riesgos relativos de cada una de estas especies hospedadoras en posibles movimientos intercontinentales de la influenza aviar. Se evaluaron estos valores al tamaño del censo para evaluarlos en ese contexto. Además de ser huéspedes de influenza aviar, muchas de estas especies de aves también son importantes en las dietas de subsistencia de los habitantes de Alaska, lo que aumenta el riesgo de exposición directa de las aves al ser humano por el virus de la influenza aviar de origen euroasiático. Se contrastaron las tasas específicas de especies de movimientos intercontinentales con la importancia de cada especie en las dietas de personas en Alaska para comprender el riesgo relativo de estos taxones para los humanos. Suponiendo tasas de infección por influenza aviar aproximadamente equivalentes entre patos, el porrón bastardo o pato boludo mayor (Aythya marila), el ánade real (Anas platyrhynchos) y el ánade rabudo eran las tres especies principales que presentaban los mayores riesgos para el movimiento de influenza aviar intercontinental tanto dentro del sistema natural como a través de la exposición a cazadores de subsistencia. La mejora de los datos sobre las tasas de infección por influenza aviar en esta región podría mejorar aún más estas evaluaciones de riesgo relativo. Estas tasas de movimiento intercontinental directamente comparables, basadas en especies, y clasificaciones de riesgo relativo deberían ayudar a modelar, monitorear y mitigar los impactos de los movimientos intercontinentales de huéspedes y de la influenza aviar.

Escaping the fate of Sisyphus: assessing resistome hybridization baits for antimicrobial resistance gene capture.

Finding, characterizing and monitoring reservoirs for antimicrobial resistance (AMR) is vital to protecting public health. Hybridization capture baits are an accurate, sensitive and cost-effective technique used to enrich and characterize DNA sequences of interest, including antimicrobial resistance genes (ARGs), in complex environmental samples. We demonstrate the continued utility of a set of 19 933 hybridization capture baits designed from the Comprehensive Antibiotic Resistance Database (CARD)v1.1.2 and Pathogenicity Island Database (PAIDB)v2.0, targeting 3565 unique nucleotide sequences that confer resistance. We demonstrate the efficiency of our bait set on a custom-made resistance mock community and complex environmental samples to increase the proportion of on-target reads as much as >200-fold. However, keeping pace with newly discovered ARGs poses a challenge when studying AMR, because novel ARGs are continually being identified and would not be included in bait sets designed prior to discovery. We provide imperative information on how our bait set performs against CARDv3.3.1, as well as a generalizable approach for deciding when and how to update hybridization capture bait sets. This research encapsulates the full life cycle of baits for hybridization capture of the resistome from design and validation (both in silico and in vitro) to utilization and forecasting updates and retirement.

Unveiling the Gut Microbiota and Resistome of Wild Cotton Mice, Peromyscus gossypinus, from Heavy Metal- and Radionuclide-Contaminated Sites in the Southeastern United States.

The prevalence of antibiotic resistance genes (ARGs) can be driven by direct selection from antibiotic use and indirect selection from substances such as heavy metals (HMs). While significant progress has been made to characterize the influence of HMs on the enrichment and dissemination of ARGs in the environment, there is still much we do not know. To fill this knowledge gap, we present a comprehensive analysis of gut bacteria associated with wild cotton mice (Peromyscus gossypinus) trapped from several areas affected by legacies of HM and radionuclide contamination. We explore how these contaminants affect gut microbial community (GMC) composition and diversity and the enrichment of antibiotic, biocide, and metal resistance genes. Although we were able to identify that a myriad of co-occurring antimicrobial and HM resistance genes appear in mice from all areas, including those without a history of contamination, the proportions of co-occurring ARGs and metal resistance genes (MRGs) are higher in sites with radionuclide contamination. These results support those from several previous studies and enhance our understanding of the coselection process, while providing new insights into the ubiquity of antimicrobial resistance in the resistome of wild animals. IMPORTANCE Antimicrobial resistance is a serious global public health concern because of its prevalence and ubiquitous distribution. The rapid dissemination of antibiotic resistance genes is thought to be the result of the massive overuse of antibiotics in agriculture and therapeutics. However, previous studies have demonstrated that the spread of antibiotic resistance genes can also be influenced by heavy metal contamination. This coselection phenomenon, whereby different resistance determinants are genetically linked on the same genetic element (coresistance) or a single genetic element provides resistance to multiple antimicrobial agents (cross-resistance), has profound clinical and environmental implications. In contrast to antibiotics, heavy metals can persist in the environment as a selection pressure for long periods of time. Thus, it is important to understand how antibiotic resistance genes are distributed in the environment and to what extent heavy metal contaminants may be driving their selection, which we have done in one environmental setting.

Improved Microbial Community Characterization of 16S rRNA via Metagenome Hybridization Capture Enrichment.

Environmental microbial diversity is often investigated from a molecular perspective using 16S ribosomal RNA (rRNA) gene amplicons and shotgun metagenomics. While amplicon methods are fast, low-cost, and have curated reference databases, they can suffer from amplification bias and are limited in genomic scope. In contrast, shotgun metagenomic methods sample more genomic regions with fewer sequence acquisition biases, but are much more expensive (even with moderate sequencing depth) and computationally challenging. Here, we develop a set of 16S rRNA sequence capture baits that offer a potential middle ground with the advantages from both approaches for investigating microbial communities. These baits cover the diversity of all 16S rRNA sequences available in the Greengenes (v. 13.5) database, with no sequence having <78% sequence identity to at least one bait for all segments of 16S. The use of our baits provide comparable results to 16S amplicon libraries and shotgun metagenomic libraries when assigning taxonomic units from 16S sequences within the metagenomic reads. We demonstrate that 16S rRNA capture baits can be used on a range of microbial samples (i.e., mock communities and rodent fecal samples) to increase the proportion of 16S rRNA sequences (average > 400-fold) and decrease analysis time to obtain consistent community assessments. Furthermore, our study reveals that bioinformatic methods used to analyze sequencing data may have a greater influence on estimates of community composition than library preparation method used, likely due in part to the extent and curation of the reference databases considered. Thus, enriching existing aliquots of shotgun metagenomic libraries and obtaining modest numbers of reads from them offers an efficient orthogonal method for assessment of bacterial community composition.

Molecular Phylogeny and Evolution of Amazon Parrots in the Greater Antilles.

Amazon parrots (Amazona spp.) colonized the islands of the Greater Antilles from the Central American mainland, but there has not been a consensus as to how and when this happened. Today, most of the five remaining island species are listed as endangered, threatened, or vulnerable as a consequence of human activity. We sequenced and annotated full mitochondrial genomes of all the extant Amazon parrot species from the Greater Antillean (A. leucocephala (Cuba), A. agilis, A. collaria (both from Jamaica), A. ventralis (Hispaniola), and A. vittata (Puerto Rico)), A. albifrons from mainland Central America, and A. rhodocorytha from the Atlantic Forest in Brazil. The assembled and annotated mitogenome maps provide information on sequence organization, variation, population diversity, and evolutionary history for the Caribbean species including the critically endangered A. vittata. Despite the larger number of available samples from the Puerto Rican Parrot Recovery Program, the sequence diversity of the A. vittata population in Puerto Rico was the lowest among all parrot species analyzed. Our data support the stepping-stone dispersal and speciation hypothesis that has started approximately 3.47 MYA when the ancestral population arrived from mainland Central America and led to diversification across the Greater Antilles, ultimately reaching the island of Puerto Rico 0.67 MYA. The results are presented and discussed in light of the geological history of the Caribbean and in the context of recent parrot evolution, island biogeography, and conservation. This analysis contributes to understating evolutionary history and empowers subsequent assessments of sequence variation and helps design future conservation efforts in the Caribbean.

Whole genome genetic variation and linkage disequilibrium in a diverse collection of Listeria monocytogenes isolates.

We performed whole-genome multi-locus sequence typing for 2554 genes in a large and heterogenous panel of 180 Listeria monocytogenes strains having diverse geographical and temporal origins. The subtyping data was used for characterizing genetic variation and evaluating patterns of linkage disequilibrium in the pan-genome of L. monocytogenes. Our analysis revealed the presence of strong linkage disequilibrium in L. monocytogenes, with ~99% of genes showing significant non-random associations with a large majority of other genes in the genome. Twenty-seven loci having lower levels of association with other genes were considered to be potential "hot spots" for horizontal gene transfer (i.e., recombination via conjugation, transduction, and/or transformation). The patterns of linkage disequilibrium in L. monocytogenes suggest limited exchange of foreign genetic material in the genome and can be used as a tool for identifying new recombinant strains. This can help understand processes contributing to the diversification and evolution of this pathogenic bacteria, thereby facilitating development of effective control measures.

An Open-Source Program (Haplo-ST) for Whole-Genome Sequence Typing Shows Extensive Diversity among Listeria monocytogenes Isolates in Outdoor Environments and Poultry Processing Plants.

A reliable and standardized classification of Listeria monocytogenes is important for accurate strain identification during outbreak investigations. Current whole-genome sequencing (WGS)-based approaches for strain characterization are either difficult to standardize, rendering them less suitable for data exchange, or are not freely available. Thus, we developed a portable and open-source tool, Haplo-ST, to improve standardization and provide maximum discriminatory potential to WGS data tied to a multilocus sequence typing (MLST) framework. Haplo-ST performs whole-genome MLST (wgMLST) for L. monocytogenes while allowing for data exchangeability worldwide. This tool takes in (i) raw WGS reads as input, (ii) cleans the raw data according to user-specified parameters, (iii) assembles genes across loci by mapping to genes from reference strains, and (iv) assigns allelic profiles to assembled genes and provides a wgMLST subtyping for each isolate. Data exchangeability relies on the tool assigning allelic profiles based on a centralized nomenclature defined by the widely used BIGSdb-Lm database. Tests of Haplo-ST's performance with simulated reads from L. monocytogenes reference strains demonstrated high sensitivity (97.5%), and coverage depths of ≥20× were found to be sufficient for wgMLST profiling. We then used Haplo-ST to characterize and differentiate between two groups of L. monocytogenes isolates derived from the natural environment and poultry processing plants. Phylogenetic reconstruction identified lineages within each group, and no lineage specificity was observed with isolate phenotypes (transient versus persistent) or origins. Genetic differentiation analyses between isolate groups identified 21 significantly differentiated loci, potentially enriched for adaptation and persistence of L. monocytogenes within poultry processing plants.IMPORTANCE We have developed an open-source tool (https://github.com/swarnalilouha/Haplo-ST) that provides allele-based subtyping of L. monocytogenes isolates at the whole-genome level. Along with allelic profiles, this tool also generates allele sequences and identifies paralogs, which is useful for phylogenetic tree reconstruction and deciphering relationships between closely related isolates. More broadly, Haplo-ST is flexible and can be adapted to characterize the genome of any haploid organism simply by installing an organism-specific gene database. Haplo-ST also allows for scalable subtyping of isolates; fewer reference genes can be used for low-resolution typing, whereas higher resolution can be achieved by increasing the number of genes used in the analysis. Our tool enabled clustering of L. monocytogenes isolates into lineages and detection of potential loci for adaptation and persistence in food processing environments. Findings from these analyses highlight the effectiveness of Haplo-ST in subtyping and evaluating relationships among isolates in studies of bacterial population genetics.