Search

A better understanding of viral factors that contribute to influenza A virus (IAV) airborne transmission is crucial for pandemic preparedness. A limited capacity for airborne transmission was recently observed in a human A(H9N2) virus isolate (A/Anhui-Lujiang/39/2018, AL/39) that possesses a leucine (L) residue at position HA1-226 (H3 numbering), indicative of human-like receptor binding potential. To evaluate the roles of the residue at this position in virus fitness and airborne transmission, a wild-type AL/39 (AL/39-wt) and a mutant virus (AL/39-HA1-L226Q) with a single substitution at position HA1-226 from leucine to glutamine (Q), a consensus residue in avian influenza viruses, were rescued and assessed in the ferret model. The AL/39-HA1-L226Q virus lost the ability to transmit by air, although the virus had a comparable capacity for replication, induced similar levels of host innate immune responses, and was detected at comparable levels in the air surrounding the inoculated ferrets relative to AL/39-wt virus. However, ferrets showed a lower susceptibility to AL/39-HA1-L226Q virus infection compared to the AL/39-wt virus. Furthermore, the AL/39-wt and AL/39-HA1-L226Q viruses each gained dominance in different anatomic sites in the respiratory tract in a co-infection competition model in ferrets. Taken together, our findings demonstrate that the increasing dominance of HA1-L226 residue in an avian A(H9N2) virus plays multifaceted roles in virus infection and transmission in the ferret model, including improved virus fitness and infectivity., Importance: Although the capacity for human-like receptor binding is a key prerequisite for non-human origin influenza A virus (IAV) to become airborne transmissible in mammalian hosts, the underlying molecular basis is not well understood. In this study, we investigated a naturally occurring substitution (leucine to glutamine) at residue 226 in the HA of an avian-origin A(H9N2) virus and assessed the impact on virus replication and airborne transmission in the ferret model. We demonstrate that the enhanced airborne transmission associated with the HA1-L226 virus was mainly due to the increased infectivity of the virus. Interestingly, we found that, unlike most sites in the ferret respiratory tract, ferret ethmoid turbinate lined with olfactory epithelium favors replication of the AL/39-HA1-L226Q virus, suggesting that this site may serve as a unique niche for IAV with avian-like receptor binding specificity to potentially allow the virus to spread to extrapulmonary tissues and to facilitate adaptation of the virus to human hosts.

Since 2020, there has been unprecedented global spread of highly pathogenic avian influenza A(H5N1) in wild bird populations with spillover into a variety of mammalian species and sporadically humans 1 . In March 2024, clade 2.3.4.4b A(H5N1) virus was first detected in dairy cattle in the U.S., with subsequent detection in numerous states 2 , leading to over a dozen confirmed human cases 3,4 . In this study, we employed the ferret model, a well-characterized species that permits concurrent investigation of viral pathogenicity and transmissibility 5 in the evaluation of A/Texas/37/2024 (TX/37) A(H5N1) virus isolated from a dairy farm worker in Texas 6 . Here, we show that the virus has a remarkable ability for robust systemic infection in ferrets, leading to high levels of virus shedding and spread to naïve contacts. Ferrets inoculated with TX/37 rapidly exhibited a severe and fatal infection, characterized by viremia and extrapulmonary spread. The virus efficiently transmitted in a direct contact setting and was capable of indirect transmission via fomites. Airborne transmission was corroborated by the detection of infectious virus shed into the air by infected animals, albeit at lower levels compared to the highly transmissible human seasonal and swine-origin H1 subtype strains. Our results show that despite maintaining an avian-like receptor binding specificity, TX/37 displays heightened virulence, transmissibility, and airborne shedding relative to other clade 2.3.4.4b virus isolated prior to the 2024 cattle outbreaks 7 , underscoring the need for continued public health vigilance., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

Clade 2.3.4.4b highly pathogenic avian influenza A(H5N1) viruses have caused large outbreaks within avian populations on five continents, with concurrent spillover into a variety of mammalian species. Mutations associated with mammalian adaptation have been sporadically identified in avian isolates, and more frequently among mammalian isolates following infection. Reports of human infection with A(H5N1) viruses following contact with infected wildlife have been reported on multiple continents, highlighting the need for pandemic risk assessment of these viruses. In this study, the pathogenicity and transmissibility of A/Chile/25945/2023 HPAI A(H5N1) virus, a novel reassortant with four gene segments (PB1, PB2, NP, MP) from North American lineage, isolated from a severe human case in Chile, was evaluated in vitro and using the ferret model. This virus possessed a high capacity to cause fatal disease, characterized by high morbidity and extrapulmonary spread in virus-inoculated ferrets. The virus was capable of transmission to naïve contacts in a direct contact setting, with contact animals similarly exhibiting severe disease, but did not exhibit productive transmission in respiratory droplet or fomite transmission models. Our results indicate that the virus would need to acquire an airborne transmissible phenotype in mammals to potentially cause a pandemic. Nonetheless, this work warrants continuous monitoring of mammalian adaptations in avian viruses, especially in strains isolated from humans, to aid pandemic preparedness efforts.

Competing Interests: The authors have declared that no competing interests exist.

In vivo assessments of influenza A virus (IAV) pathogenicity and transmissibility in ferrets represent a crucial component of many pandemic risk assessment rubrics, but few systematic efforts to identify which data from in vivo experimentation are most useful for predicting pathogenesis and transmission outcomes have been conducted. To this aim, we aggregated viral and molecular data from 125 contemporary IAV (H1, H2, H3, H5, H7, and H9 subtypes) evaluated in ferrets under a consistent protocol. Three overarching predictive classification outcomes (lethality, morbidity, transmissibility) were constructed using machine learning (ML) techniques, employing datasets emphasizing virological and clinical parameters from inoculated ferrets, limited to viral sequence-based information, or combining both data types. Among 11 different ML algorithms tested and assessed, gradient boosting machines and random forest algorithms yielded the highest performance, with models for lethality and transmission consistently better performing than models predicting morbidity. Comparisons of feature selection among models was performed, and highest performing models were validated with results from external risk assessment studies. Our findings show that ML algorithms can be used to summarize complex in vivo experimental work into succinct summaries that inform and enhance risk assessment criteria for pandemic preparedness that take in vivo data into account., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Trypanosoma cruzi , the causative agent of Chagas disease ( American trypanosomiasis ), is a highly complex zoonosis that is present throughout South America, Central America, and Mexico. The transmission of this disease is influenced by various factors, including human activities like deforestation and land use changes, which may have altered the natural transmission cycles and their connection to the environment. In this study conducted in the Argentine Chaco region, we examined the transmission dynamics of T. cruzi by collecting blood samples from wild and domestic animals, as well as triatomine bugs from human dwellings, across five sites of varying anthropic intervention. Samples were analyzed for T. cruzi infection via qPCR, and we additionally examined triatomines for bloodmeal analysis via NGS amplicon sequencing. Our analysis revealed a 15.3% infection rate among 20 wild species (n = 123) and no T. cruzi presence in 9 species of domestic animals (n = 1359) or collected triatomines via qPCR. Additionally, we found chicken (34.28%), human (21.59%), and goat (19.36%) as the predominant bloodmeal sources across all sites. These findings suggest that anthropic intervention and other variables analyzed may have directly impacted the spillover dynamics of T. cruzi 's sylvatic cycle and potentially reduced its prevalence in human habitats.

Data from influenza A virus (IAV) infected ferrets provides invaluable information towards the study of novel and emerging viruses that pose a threat to human health. This gold standard model can recapitulate many clinical signs of infection present in IAV-infected humans, support virus replication of human, avian, swine, and other zoonotic strains without prior adaptation, and permit evaluation of virus transmissibility by multiple modes. While ferrets have been employed in risk assessment settings for >20 years, results from this work are typically reported in discrete stand-alone publications, making aggregation of raw data from this work over time nearly impossible. Here, we describe a dataset of 728 ferrets inoculated with 126 unique IAV, conducted by a single research group under a uniform experimental protocol. This collection of morbidity, mortality, and viral titer data represents the largest publicly available dataset to date of in vivo-generated IAV infection outcomes on a per-ferret level., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Viral pathogenesis and therapeutic screening studies that utilize small mammalian models rely on the accurate quantification and interpretation of morbidity measurements, such as weight and body temperature, which can vary depending on the model, agent and/or experimental design used. As a result, morbidity-related data are frequently normalized within and across screening studies to aid with their interpretation. However, such data normalization can be performed in a variety of ways, leading to differences in conclusions drawn and making comparisons between studies challenging. Here, we discuss variability in the normalization, interpretation, and presentation of morbidity measurements for four model species frequently used to study a diverse range of human viral pathogens - mice, hamsters, guinea pigs and ferrets. We also analyze findings aggregated from influenza A virus-infected ferrets to contextualize this discussion. We focus on serially collected weight and temperature data to illustrate how the conclusions drawn from this information can vary depending on how raw data are collected, normalized and measured. Taken together, this work supports continued efforts in understanding how normalization affects the interpretation of morbidity data and highlights best practices to improve the interpretation and utility of these findings for extrapolation to public health contexts., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

As use of the ferret model to study influenza A virus (IAV) pathogenicity increases, periodic assessment of data generated in this model is warranted, to identify features associated with virus replication throughout the respiratory tract and to refine future analyses. However, protocol-specific differences present between independent laboratories limit easy aggregation of virological data. We compiled viral titer and clinical data from >1,000 ferrets inoculated with 125 contemporary IAV under a consistent experimental protocol (including high- and low-pathogenicity avian, swine-origin, and human viruses, spanning H1, H2, H3, H5, H7, and H9 subtypes) and examined which meaningful and statistically supported associations were present among numerous quantitative measurements. Viral titers correlated positively between ferret nasal turbinate tissue, lung tissue, and nasal wash specimens, though the strength of the associations varied, notably regarding the particular nasal wash summary measure employed and properties of the virus itself. Use of correlation coefficients and mediation analyses further supported the interconnectedness of viral titer measurements taken at different sites throughout the respiratory tract. IAV possessing mammalian host adaptation markers in the HA and PB2 exhibited more rapid growth in the ferret upper respiratory tract early after infection, supported by quantities derived from infectious titer data to capture infection progression, compared with viruses bearing hallmarks of avian IAV. Collectively, this work identifies summary metrics most closely linked with virological and phenotypic outcomes in ferrets, supporting continued refinement of data analyzed from in vivo experimentation, notably from studies conducted to evaluate the public health risk posed by novel and emerging IAV.IMPORTANCEFerrets are frequently employed to study the pandemic potential of novel and emerging influenza A viruses. However, systematic retrospective analyses of data generated from these experiments are rarely performed, limiting our ability to identify trends in this data and explore how analyses can be refined. Using logarithmic viral titer and clinical data aggregated from one research group over 20 years, we assessed which meaningful and statistically supported associations were present among numerous quantitative measurements obtained from influenza A virus (IAV)-infected ferrets, including those capturing viral titers, infection progression, and disease severity. We identified numerous linear correlations between parameters assessing virus replication at discrete sites in vivo , including parameters capturing infection progression not frequently employed in the field, and sought to investigate the interconnected nature of these associations. This work supports continued refinement of data analyzed from in vivo experimentation, notably from studies which evaluate the public health risk posed by IAV., Competing Interests: The authors declare no conflict of interest.

Importance: Despite the accumulation of evidence showing that airborne transmissible influenza A virus (IAV) typically has a lower pH threshold for hemagglutinin (HA) fusion activation, the underlying mechanism for such a link remains unclear. In our study, by using a pair of isogenic recombinant A(H9N2) viruses with a phenotypical difference in virus airborne transmission in a ferret model due to an acid-destabilizing mutation (HA1-Y17H) in the HA, we demonstrate that an acid-stable A(H9N2) virus possesses a multitude of advantages over its less stable counterpart, including better fitness in the ferret respiratory tract, more effective aerosol emission from infected animals, and improved host susceptibility. Our study provides supporting evidence for the requirement of acid stability in efficient airborne transmission of IAV and sheds light on fundamental mechanisms for virus airborne transmission., Competing Interests: The authors declare no conflict of interest.

While primarily considered a respiratory pathogen, influenza A virus (IAV) is nonetheless capable of spreading to, and replicating in, numerous extrapulmonary tissues in humans. However, within-host assessments of genetic diversity during multicycle replication have been largely limited to respiratory tract tissues and specimens. As selective pressures can vary greatly between anatomical sites, there is a need to examine how measures of viral diversity may vary between influenza viruses exhibiting different tropisms in humans, as well as following influenza virus infection of cells derived from different organ systems. Here, we employed human primary tissue constructs emulative of the human airway or corneal surface, and we infected both with a panel of human- and avian-origin IAV, inclusive of H1 and H3 subtype human viruses and highly pathogenic H5 and H7 subtype viruses, which are associated with both respiratory disease and conjunctivitis following human infection. While both cell types supported productive replication of all viruses, airway-derived tissue constructs elicited greater induction of genes associated with antiviral responses than did corneal-derived constructs. We used next-generation sequencing to examine viral mutations and population diversity, utilizing several metrics. With few exceptions, generally comparable measures of viral diversity and mutational frequency were detected following homologous virus infection of both respiratory-origin and ocular-origin tissue constructs. Expansion of within-host assessments of genetic diversity to include IAV with atypical clinical presentations in humans or in extrapulmonary cell types can provide greater insight into understanding those features most prone to modulation in the context of viral tropism. IMPORTANCE Influenza A virus (IAV) can infect tissues both within and beyond the respiratory tract, leading to extrapulmonary complications, such as conjunctivitis or gastrointestinal disease. Selective pressures governing virus replication and induction of host responses can vary based on the anatomical site of infection, yet studies examining within-host assessments of genetic diversity are typically only conducted in cells derived from the respiratory tract. We examined the contribution of influenza virus tropism on these properties two different ways: by using IAV associated with different tropisms in humans, and by infecting human cell types from two different organ systems susceptible to IAV infection. Despite the diversity of cell types and viruses employed, we observed generally similar measures of viral diversity postinfection across all conditions tested; these findings nonetheless contribute to a greater understanding of the role tissue type contributes to the dynamics of virus evolution within a human host., Competing Interests: The authors declare no conflict of interest.

Competition assays were conducted in vitro and in vivo to examine how the Delta (B.1.617.2) variant displaced the prototype Washington/1/2020 (WA/1) strain. While WA/1 virus exhibited a moderately increased proportion compared to that in the inoculum following co-infection in human respiratory cells, Delta variant possessed a substantial in vivo fitness advantage as this virus becoming predominant in both inoculated and contact animals. This work identifies critical traits of the Delta variant that likely played a role in it becoming a dominant variant and highlights the necessities of employing multiple model systems to assess the fitness of newly emerged SARS-CoV-2 variants., (Published by Elsevier Inc.)

As influenza A viruses (IAV) continue to cross species barriers and cause human infection, the establishment of risk assessment rubrics has improved pandemic preparedness efforts. In vivo pathogenicity and transmissibility evaluations in the ferret model represent a critical component of this work. As the relative contribution of in vitro experimentation to these rubrics has not been closely examined, we sought to evaluate to what extent viral titer measurements over the course of in vitro infections are predictive or correlates of nasal wash and tissue measurements for IAV infections in vivo . We compiled data from ferrets inoculated with an extensive panel of over 50 human and zoonotic IAV (inclusive of swine-origin and high- and low-pathogenicity avian influenza viruses associated with human infection) under a consistent protocol, with all viruses concurrently tested in a human bronchial epithelial cell line (Calu-3). Viral titers in ferret nasal wash specimens and nasal turbinate tissue correlated positively with peak titer in Calu-3 cells, whereas additional phenotypic and molecular determinants of influenza virus virulence and transmissibility in ferrets varied in their association with in vitro viral titer measurements. Mathematical modeling was used to estimate more generalizable key replication kinetic parameters from raw in vitro viral titers, revealing commonalities between viral infection progression in vivo and in vitro . Meta-analyses inclusive of IAV that display a diverse range of phenotypes in ferrets, interpreted with mathematical modeling of viral kinetic parameters, can provide critical information supporting a more rigorous and appropriate contextualization of in vitro experiments toward pandemic preparedness. IMPORTANCE Both in vitro and in vivo models are employed for assessing the pandemic potential of novel and emerging influenza A viruses in laboratory settings, but systematic examinations of how well viral titer measurements obtained in vitro align with results from in vivo experimentation are not frequently performed. We show that certain viral titer measurements following infection of a human bronchial epithelial cell line are positively correlated with viral titers in specimens collected from virus-inoculated ferrets and employ mathematical modeling to identify commonalities between viral infection progression between both models. These analyses provide a necessary first step in enhanced interpretation and incorporation of in vitro -derived data in risk assessment activities and highlight the utility of employing mathematical modeling approaches to more closely examine features of virus replication not identifiable by experimental studies alone.

Despite reports of confirmed human infection following ocular exposure with both influenza A virus (IAV) and SARS-CoV-2, the dynamics of virus spread throughout oculonasal tissues and the relative capacity of virus transmission following ocular inoculation remain poorly understood. Furthermore, the impact of exposure route on subsequent release of airborne viral particles into the air has not been examined previously. To assess this, ferrets were inoculated by the ocular route with A(H1N1)pdm09 and A(H7N9) IAVs and two SARS-CoV-2 (early pandemic Washington/1 and Delta variant) viruses. Virus replication was assessed in both respiratory and ocular specimens, and transmission was evaluated in direct contact or respiratory droplet settings. Viral RNA in aerosols shed by inoculated ferrets was quantified with a two-stage cyclone aerosol sampler (National Institute for Occupational Safety and Health [NIOSH]). All IAV and SARS-CoV-2 viruses mounted a productive and transmissible infection in ferrets following ocular inoculation, with peak viral titers and release of virus-laden aerosols from ferrets indistinguishable from those from ferrets inoculated by previously characterized intranasal inoculation methods. Viral RNA was detected in ferret conjunctival washes from all viruses examined, though infectious virus in this specimen was recovered only following IAV inoculation. Low-dose ocular-only aerosol exposure or inhalation aerosol exposure of ferrets to IAV similarly led to productive infection of ferrets and shedding of aerosolized virus. Viral evolution during infection was comparable between all inoculation routes examined. These data support that both IAV and SARS-CoV-2 can establish a high-titer mammalian infection following ocular exposure that is associated with rapid detection of virus-laden aerosols shed by inoculated animals. IMPORTANCE Documented human infection with influenza viruses and SARS-CoV-2 has been reported among individuals wearing respiratory protection in the absence of eye protection, highlighting the capacity of these respiratory tract-tropic viruses to exploit nonrespiratory routes of exposure to initiate productive infection. However, comprehensive evaluations of how ocular exposure may modulate virus pathogenicity and transmissibility in mammals relative to respiratory exposure are limited and have not investigated multiple virus families side by side. Using the ferret model, we show that ocular exposure with multiple strains of either coronaviruses or influenza A viruses leads to an infection that results in shedding of detectable aerosolized virus from inoculated animals, contributing toward onward transmission of both viruses to susceptible contacts. Collectively, these studies support that the ocular surface represents a susceptible mucosal surface that, if exposed to a sufficient quantity of either virus, permits establishment of an infection which is similarly transmissible as that following respiratory exposure.

The continued spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans necessitates evaluation of variants for enhanced virulence and transmission. We used the ferret model to perform a comparative analysis of four SARS-CoV-2 strains, including an early pandemic isolate from the United States (WA1), and representatives of the Alpha, Beta, and Delta lineages. While Beta virus was not capable of pronounced replication in ferrets, WA1, Alpha, and Delta viruses productively replicated in the ferret upper respiratory tract, despite causing only mild disease with no overt histopathological changes. Strain-specific transmissibility was observed; WA1 and Delta viruses transmitted in a direct contact setting, whereas Delta virus was also capable of limited airborne transmission. Viral RNA was shed in exhaled air particles from all inoculated animals but was highest for Delta virus. Prior infection with SARS-CoV-2 offered varied protection against reinfection with either homologous or heterologous variants. Notable genomic variants in the spike protein were most frequently detected following WA1 and Delta virus infection. IMPORTANCE Continued surveillance and risk assessment of emerging SARS-CoV-2 variants are critical for pandemic response and preparedness. As such, in vivo evaluations are indispensable for early detection of variants with enhanced virulence and transmission. Here, we used the ferret model to compare the pathogenicity and transmissibility of an original SARS-CoV-2 isolate (USA-WA1/2020 [WA1]) to those of a panel of Alpha, Beta, and Delta variants, as well as to evaluate protection from homologous and heterologous reinfection. We observed strain-specific differences in replication kinetics in the ferret respiratory tract and virus load emitted into the air, revealing enhanced transmissibility of the Delta virus relative to previously detected strains. Prior infection with SARS-CoV-2 provided varied levels of protection from reinfection, with the Beta strain eliciting the lowest level of protection. Overall, we found that ferrets represent a useful model for comparative assessments of SARS-CoV-2 infection, transmission, and reinfection.

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., (© 2021 Society for Applied Microbiology and John Wiley & Sons Ltd.)

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.

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., Competing Interests: The EHS DNA lab provides oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services used in this manuscript (baddna.uga.edu). BB and AD were employed by, and thereby have financial interest in, Daicel Arbor Biosciences, who provided the in-solution capture reagents used in this work. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Beaudry, Wang, Kieran, Thomas, Bayona-Vásquez, Gao, Devault, Brunelle, Lu, Wang, Rhodes and Glenn.)

Lack of knowledge of taxonomic biodiversity and reliable genetic markers in Trypanosomatidae limit our understanding of their phylogenetic relationships. Ultraconserved elements (UCEs) have improved phylogenetic analyses and inferences in many vertebrate and invertebrate taxa. However, it is unknown whether protozoans have these markers, their abundance, and if these could be reliably used for phylogenetics. In this study I design a target enrichment bait set for UCE loci for this group. In silico testing showed good loci recovery rates across 63 taxa and produced consistent, highly supported phylogenetic trees. This bait set adds a new resource of useful genetic markers for Trypanosomatidae phylogenetics., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Chagas disease is a zoonosis that affects several million people and is caused by the parasite Trypanosoma cruzi, which is mainly transmitted through the feces of triatomine bugs. Within triatomines, several Rhodnius species have been found inhabiting palms, and certain factors such as palm species and location have been related to the abundance and T. cruzi infection of those insects in palms. In this study, the main goal was to determine if R. prolixus abundances and infection rates in Attalea butyracea and Elaeis guineensis palms are related to ecological factors such as palm species, crown microclimate, and available blood meal sources. Triatomine sampling was performed in two municipalities of Casanare, Colombia, specifically in the intersection of riparian forests and oil palm plantations. For R. prolixus abundance per palm, the predictors showing more relationship were palm species and blood meal species identified in the palm, and for T. cruzi infection per triatomine, they were palm species and nymphal stage. Palm microclimate was very similar in both palm species and did not show a relationship with triatomine abundance. Comparing palm species, A. butyracea showed more blood meal species, including more refractory host species, than E. guineensis, but lower T. cruzi infection rate and parasitaemia. Interestingly, non-arboreal blood meal species were frequently found in the analyzed nymphs, indicating that the blood source for R. prolixus in palms corresponded to all the fauna located in the surrounded landscape and not only in the palm. These results could expose a new ecological scenario to interpret the T. cruzi transmission in sylvatic environments., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Triatominae is a subfamily of blood-sucking reduviid hemipterans of public health importance primarily in tropical and sub-tropical regions of the Americas, whose members possess various morphological adaptations closely associated to hematophagy. Despite their medical importance, the systematics of the subfamily is far from resolved, particularly within the tribe Triatomini. Here we employed mitochondrial genome DNA sequences to reconstruct the phylogenetic relationships among 19 species of the North-Central American (NCA) clade of Triatomini and to estimate the times of origin and diversification of its main clades. Twenty-nine mitogenomes were examined for representative specimens of 25 species, including the outgroup. Phylogenetic informativeness estimated for each protein-coding gene showed that cox1, cox2 and atp6 were the most informative markers, whereas atp8 and nad4 had high saturation levels. Phylogenetic analyses excluding the latter two protein-coding genes recovered an almost fully resolved topology. The NCA clade apparently originated shortly after emergence of an initial land bridge of the Panama Isthmus, ca. 15.05-20.05 Mya. An Asian/pantropical subclade with Linshcosteus costalis, Triatoma rubrofasciata and T. migrans was nested within the NCA clade, from which it diverged ca. 12.42-17.3Mya. Uncorrected cox1 and 13 protein-coding gene distances suggest the existence of additional species within the dimidiata complex. In contrast, T. phyllosoma, T. mazzottii and T. longipennis, from the phyllosoma complex, have considerably low cox1 and 13 PCG distances among them, suggesting mitochondrial introgression or conspecificity. Our study yielded a robust phylogeny for the group, which could be tested with further phylogenetic hypotheses based on nuclear genome-wide markers., Competing Interests: Declaration of Competing Interest None., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Rhodnius pallescens is the principal vector of Chagas disease in Panama. Recently a dark chromatic morph has been discovered in the highlands of Veraguas Province. Limited genetic studies have been conducted with regards to the population structure and dispersal potential of Triatominae vectors, particularly in R. pallescens. Next generation sequencing methods such as RADseq and complete mitochondrial DNA (mtDNA) genome sequencing have great potential for examining vector biology across space and time. Here we utilize a RADseq method (3RAD), along with complete mtDNA sequencing, to examine the population structure of the two chromatic morpho types of R. pallescens in Panama. We sequenced 105 R. pallescens samples from five localities in Panama. We generated a 2216 SNP dataset and 6 complete mtDNA genomes. RADseq showed significant differentiation among the five localities (F CT  = 0.695; P = .004), but most of this was between localities with the dark vs. light chromatic morphs (Veraguas vs. Panama Oeste). The mtDNA genomes showed a 97-98% similarity between dark and light chromatic morphs across all genes and a 502 bp insert in light morphs. Thus, both the RADseq and mtDNA data showed highly differentiated clades with essentially no gene flow between the dark and light chromatic morphs from Veraguas and central Panama respectively. We discuss the growing evidence showing clear distinctions between these two morpho types with the possibility that these are separate species, an area of research that requires further investigation. Finally, we discuss the cost-effectiveness of 3RAD which is a third of the cost compared to other RADseq methods used recently in Chagas disease vector research., Competing Interests: Declaration of Competing Interest Jose E. Calzada and Azael Saldana are members of the SNI (Sistema Nacional de Investigación from SENACYT of Panama). All other authors declare no competing interest., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Selenium represents an essential trace nutrient that is necessary for biological functions. Deficiencies can induce disease, but excess can induce toxicity. Selenium deficiency is a major concern in underdeveloped countries, while also posing as a toxic pollutant in waterways surrounding landfills, agricultural areas, and fossil fuel production sites. We examined the microbiome of selenomethionine (SeMet) fed American alligators (Alligator mississippiensis) at the beginning and end of a 7-week exposure experiment. Alligators were randomly divided into three groups: control and 1000 or 2000 ppm SeMet. DNA from before exposure (oral and cloaca swabs) and post-exposure (oral, cloaca, small & large intestines) sampling were extracted and amplified for bacterial 16 s rRNA. While treatment did not seem to have much effect, we observed a predominance of Fusobacteriaceae and Porpyromonodaceae across all tissue types. Cetobacterium and Clostridium are the most abundant genera as potential indicators of the aquatic and carrion feeding lifestyle of alligators.

MicroRNAs (miRNAs) are 18-24 nucleotide regulatory RNAs. They are involved in the regulation of genetic and biological pathways through post transcriptional gene silencing and/or translational repression. Data suggests a slow evolutionary rate for the saltwater crocodile (Crocodylus porosus) over the past several million years when compared to birds, the closest extant relatives of crocodilians. Understanding gene regulation in the saltwater crocodile in the context of relatively slow genomic change thus holds potential for the investigation of genomics, evolution, and adaptation. Utilizing eleven tissue types and sixteen small RNA libraries, we report 644 miRNAs in the saltwater crocodile with >78% of miRNAs being novel to crocodilians. We also identified potential targets for the miRNAs and analyzed the relationship of the miRNA repertoire to transposable elements (TEs). Results suggest an increased association of DNA transposons with miRNAs when compared to retrotransposons. This work reports the first comprehensive analysis of miRNAs in Crocodylus porosus and addresses the potential impacts of miRNAs in regulating the genome in the saltwater crocodile. In addition, the data suggests a supporting role of TEs as a source for miRNAs, adding to the increasing evidence that TEs play a significant role in the evolution of gene regulation., Competing Interests: Declaration of competing interest None., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Contaminants such as heavy metals may contribute to the dissemination of antimicrobial resistance (AMR) by enriching resistance gene determinants via co-selection mechanisms. In the present study, a survey was performed on soils collected from four areas at the Savannah River Site (SRS), South Carolina, USA, with varying contaminant profiles: relatively pristine (Upper Three Runs), heavy metals (Ash Basins), radionuclides (Pond B) and heavy metal and radionuclides (Tim's Branch). Using 16S rRNA gene amplicon sequencing, we explored the structure and diversity of soil bacterial communities. Sites with legacies of metal and/or radionuclide contamination displayed significantly lower bacterial diversity compared to the reference site. Metagenomic analysis indicated that multidrug and vancomycin antibiotic resistance genes (ARGs) and metal resistance genes (MRGs) including those associated with copper, arsenic, iron, nickel and zinc were prominent in all soils including the reference site. However, significant differences were found in the relative abundance and diversity of certain ARGs and MRGs in soils with metal/radionuclide contaminated soils compared to the reference site. Co-occurrence patterns revealed significant ARG/MRG subtypes in predominant soil taxa including Acidobacteriaceae, Bradyrhizobium, Mycobacterium, Streptomyces, Verrumicrobium, Actinomadura and Solirubacterales. Overall, the study emphasizes the potential risk of human activities on the dissemination of AMR in the environment., (© 2020 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology.)

Harpalus pensylvanicus (Coloptera: Carabidae) is a weed seed predator common throughout the United States. While Carabidae is a very large group of beetles, limited genomic resources exist, especially mitochondrial genomes. This study expands research in this area by assembling and annotating the complete mitochondrial genome of H. pensylvanicus and performs phylogenetic analyses with closely related species. Further use of the metagenomic data was made to characterize microbial taxa and clusters of orthologous groups of proteins. The complete mitochondrial genome is 16,434 bp in length, AT rich, and consist of 13 protein coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a control region. Phylogenetic analyses were congruent with the Harpalinae and Pterostichinae clade together. Microbial classification shows a predominance of Gamma- (37.77%) and Alpha-proteobacteria (33.97%)., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Most of the research efforts involving the bovine gastrointestinal microbiota have focused on cattle's forestomach, particularly the rumen, so information concerning the bovine fecal microbiota is more scarce, especially in young beef cattle. The present study was performed to evaluate the ruminal and fecal microbiotas of beef calves as they reached the end of their nursing phase. A total of 18 Angus cow/calf pairs were selected and assigned to one of two treatment groups for the last 92 days of the calves' nursing period, as follows: 1) calves were supplemented with concentrate in a creep feeding system; or 2) control group with no supplementation of calves. After 92 days, ruminal and fecal samples were individually obtained from calves in both groups, and their microbiotas were evaluated using 16S rRNA gene sequencing. Ruminal samples were predominated by Prevotella (18 to 23% of the total bacterial abundance), regardless if calves received supplementation or not; however, in the feces, Prevotella was only the seventh most abundant genus (0.6 to 2.1% of total bacterial abundance). Both the rumen (P = 0.01) and the feces (P = 0.05) of calves that received supplementation had greater abundance of Firmicutes. In addition, calves that were supplemented had lower abundance of Fibrobacteres (P = 0.03) in their rumens. Regardless if the calves were supplemented or not, Faith's Phylogenetic Diversity index (P ≤ 0.007) and total concentration of short chain fatty acids (P < 0.001) were both greater in the rumen than in the feces of calves. In summary, the ruminal and fecal microbiotas of weanling beef calves were considerably distinct. Additionally, supplementation with creep feed caused some significant changes in the composition of the gastrointestinal microbiota of the calves, especially in the rumen, where supplementation caused an increase in Firmicutes and a decrease in abundance of Fibrobacteres., Competing Interests: The authors have declared that no competing interests exist.

Phellinus noxius is a root-decay pathogen with a pan-tropical/subtropical distribution that attacks a wide range of tree hosts. For this study, genomic sequencing was conducted on P. noxius isolate P919-02W.7 from Federated States of Micronesia (Pohnpei), and its gene expression profile was analyzed using different host wood (Acer, Pinus, Prunus, and Salix) substrates. The assembled genome was 33.92 Mbp with 2954 contigs and 9389 predicted genes. Only small differences were observed in size and gene content in comparison with two other P. noxius genome assemblies (isolates OVT-YTM/97 from Hong Kong, China and FFPRI411160 from Japan, respectively). Genome analysis of P. noxius isolate P919-02W.7 revealed 488 genes encoding proteins related to carbohydrate and lignin metabolism, many of these enzymes are associated with degradation of plant cell wall components. Most of the transcripts expressed by P. noxius isolate P919-02W.7 were similar regardless of wood substrates. This study highlights the vast suite of decomposing enzymes produced by P. noxius, which suggests potential for degrading diverse wood substrates, even from temperate host trees. This information contributes to our understanding of pathogen ecology, mechanisms of wood decomposition, and pathogenic/saprophytic lifestyle., (Copyright © 2020 British Mycological Society. All rights reserved.)

Our previous study showed heritable reproductive toxicity in the nematode Caenorhabditis elegans after multigenerational exposure to AgNO 3 and silver nanoparticles (Ag-NPs). The aim of this study was to determine whether such inheritable effects are correlated with induced germline mutations in C. elegans. Individual C. elegans lineages were exposed for 10 generations to equitoxic concentrations at EC 30 of AgNO 3 , Ag-NPs, and sulfidized Ag-NPs (sAg-NPs), a predominant environmentally transformed product of pristine Ag-NPs. The mutations were detected via whole genome DNA sequencing approach by comparing F 0 and F 10 generations. An increase in the total number of variants, though not statistically significant, was observed for all Ag treatments and the variants were mainly contributed by single nucleotide polymorphisms (SNPs). This potentially contributed towards reproductive as well as growth toxicity shown previously after ten generations of exposure in every Ag treatment. However, despite Ag-NPs and AgNO 3 inducing stronger reproductive toxicity than sAg-NPs, exposure to sAg-NPs resulted in higher mutation accumulation with significant increase in the number of transversions. Thus our results suggest that other mechanisms of inheritance, such as epigenetics, may be at play in Ag-NP- and AgNO 3 -induced multigenerational and transgenerational reproductive toxicity., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Background: Triatomine bugs are vectors of the protozoan parasite Trypanosoma cruzi, which causes Chagas disease. Rhodnius pallescens is a major vector of Chagas disease in Panama. Understanding the microbial ecology of disease vectors is important in the development of vector management strategies that target vector survival and fitness. In this study we examined the whole-body microbial composition of R. pallescens from three locations in Panama., Methods: We collected 89 R. pallescens specimens using Noireau traps in Attalea butyracea palms. We then extracted total DNA from whole-bodies of specimens and amplified bacterial microbiota using 16S rRNA metabarcoding PCR. The 16S libraries were sequenced on an Illumina MiSeq and analyzed using QIIME2 software., Results: We found Proteobacteria, Actinobacteria, Bacteroidetes and Firmicutes to be the most abundant bacterial phyla across all samples. Geographical location showed the largest difference in microbial composition with northern Veraguas Province having the most diversity and Panama Oeste Province localities being most similar to each other. Wolbachia was detected in high abundance (48-72%) at Panama Oeste area localities with a complete absence of detection in Veraguas Province. No significant differences in microbial composition were detected between triatomine age class, primary blood meal source, or T. cruzi infection status., Conclusions: We found biogeographical regions differ in microbial composition among R. pallescens populations in Panama. While overall the microbiota has bacterial taxa consistent with previous studies in triatomine microbial ecology, locality differences are an important observation for future studies. Geographical heterogeneity in microbiomes of vectors is an important consideration for future developments that leverage microbiomes for disease control.

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. The EHS DNA lab (baddna.uga.edu) and Georgia Genomics Bioinformatics Core (dna.uga.edu) provide oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers presented in this manuscript., (©2019 Glenn et al.)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I , we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr&#95;Demuxy., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. The EHS DNA lab provides oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers presented in this manuscript (baddna.uga.edu). This document has been reviewed in accordance with U.S Environmental Protection Agency policy and approved for publication. Any mention of trade names, manufacturer or products does not imply an endorsement by the United States Government or the U.S. Environmental Protection Agency. EPA and its employees do not endorse any commercial products., (©2019 Glenn et al.)

Molecular ecologists frequently use genome reduction strategies that rely upon restriction enzyme digestion of genomic DNA to sample consistent portions of the genome from many individuals (e.g., RADseq, GBS). However, researchers often find the existing methods expensive to initiate and/or difficult to implement consistently, especially because it is difficult to multiplex sufficient numbers of samples to fill entire sequencing lanes. Here, we introduce a low-cost and highly robust approach for the construction of dual-digest RADseq libraries that build on adapters and primers designed in Adapterama I . Major features of our method include: (1) minimizing the number of processing steps; (2) focusing on a single strand of sample DNA for library construction, allowing the use of a non-phosphorylated adapter on one end; (3) ligating adapters in the presence of active restriction enzymes, thereby reducing chimeras; (4) including an optional third restriction enzyme to cut apart adapter-dimers formed by the phosphorylated adapter, thus increasing the efficiency of adapter ligation to sample DNA, which is particularly effective when only low quantity/quality DNA samples are available; (5) interchangeable adapter designs; (6) incorporating variable-length internal indexes within the adapters to increase the scope of sample indexing, facilitate pooling, and increase sequence diversity; (7) maintaining compatibility with universal dual-indexed primers and thus, Illumina sequencing reagents and libraries; and, (8) easy modification for the identification of PCR duplicates. We present eight adapter designs that work with 72 restriction enzyme combinations. We demonstrate the efficiency of our approach by comparing it with existing methods, and we validate its utility through the discovery of many variable loci in a variety of non-model organisms. Our 2RAD/3RAD method is easy to perform, has low startup costs, has increased utility with low-concentration input DNA, and produces libraries that can be highly-multiplexed and pooled with other Illumina libraries., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. Kerin E. Bentley is currently affiliated with LeafWorks, Inc. The EHS DNA lab (http://baddna.uga.edu) provides oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers., (©2019 Bayona-Vásquez et al.)

A study was conducted to examine the effects of supplementing beef calves during their suckling phase (popularly known as creep feeding) with supplements that contained or did not contain the enzyme xylanase. Forty-two cow-calf pairs were divided into three groups and assigned to one of three treatments for a period of 105 days, as follows: (1) No supplemental feed for calves (control; CON); (2) Corn and soybean meal-based supplement feed for calves (positive control; PCON); and (3) Same feed regimen as PCON with xylanase added to the supplement (enzyme; ENZ). After 105 days, out of the 42 calves participating in the study, 25 male calves were randomly selected (8 from CON, 9 from PCON, and 8 from ENZ) and samples of their forestomach were collected by esophageal tubing. Immediately after this procedure, all calves were weaned, commingled, and placed in a common post-weaning diet for 4 weeks. At the end of this period, ruminal fluid was once again collected from the same 25 calves. All samples were subjected to DNA extraction and 16S rRNA gene sequencing. At weaning, most of the alpha diversity indexes were greater in CON; however, no differences ( P ≥ 0.23) in alpha diversity were observed in samples collected 4 weeks after weaning. Regardless of treatment, 2 phyla - Bacteroidetes and Firmicutes - comprised approximately 80% of the total bacterial abundance of samples collected on both days. At the genus level, an effect of diet ( P = 0.02) was observed for Prevotella in the samples collected at weaning; however, no differences were detected in the samples collected 4 weeks after weaning. Calf average daily gain (ADG) during the 105-day creep feeding trial tended ( P = 0.09) to be greater in the groups that received supplementation, with the greatest numerical value observed in ENZ. Moreover, there was a positive correlation (ρ = 0.43; P = 0.03) between ADG and abundance of Prevotella , indicating the importance of this bacterial group for ruminants. In summary, most of the significant differences found in this study were detected at weaning, and the majority of them disappeared 4 weeks after the calves were weaned and commingled.

Background: Oil palm plantation establishment in Colombia has the potential to impact Chagas disease transmission by increasing the distribution range of Rhodnius prolixus. In fact, previous studies have reported Trypanosoma cruzi natural infection in R. prolixus captured in oil palms (Elaeis guineensis) in the Orinoco region, Colombia. The aim of this study is to understand T. cruzi infection in vectors in oil palm plantations relative to community composition and host dietary specialization by analyzing vector blood meals and comparing these results to vectors captured in a native palm tree species, Attalea butyracea., Methods: Rhodnius prolixus nymphs (n = 316) were collected from A. butyracea and E. guineensis palms in Tauramena, Casanare, Colombia. Vector blood meals from these nymphs were determined by amplifying and sequencing a vertebrate-specific 12S rRNA gene fragment., Results: Eighteen vertebrate species were identified and pigs (Sus scrofa) made up the highest proportion of blood meals in both habitats, followed by house mouse (Mus musculus) and opossum (Didelphis marsupialis). Individual bugs feeding only from generalist mammal species had the highest predicted vector infection rate, suggesting that generalist mammalian species are more competent hosts for T. cruzi infection ., Conclusions: Oil palm plantations and A. butyracea palms found in altered areas provide a similar quality habitat for R. prolixus populations in terms of blood meal availability. Both habitats showed similarities in vector infection rate and potential host species, representing a single T. cruzi transmission scenario at the introduced oil palm plantation and native Attalea palm interface.

Target enrichment of conserved genomic regions facilitates collecting sequences of many orthologous loci from non-model organisms to address phylogenetic, phylogeographic, population genetic, and molecular evolution questions. Bait sets for sequence capture can simultaneously target thousands of loci, which opens new avenues of research on speciose groups. Current phylogenetic hypotheses on the >103,000 species of Hemiptera have failed to unambiguously resolve major nodes, suggesting that alternative datasets and more thorough taxon sampling may be required to resolve relationships. We use a recently designed ultraconserved element (UCE) bait set for Hemiptera, with a focus on the suborder Heteroptera, or the true bugs, to test previously proposed relationships. We present newly generated UCE data for 36 samples representing three suborders, all seven heteropteran infraorders, 23 families, and 34 genera of Hemiptera and one thysanopteran outgroup. To improve taxon sampling, we also mined additional UCE loci in silico from published hemipteran genomic and transcriptomic data. We obtained 2271 UCE loci for newly sequenced hemipteran taxa, ranging from 265 to 1696 (average 904) per sample. These were similar in number to the data mined from transcriptomes and genomes, but with fewer loci overall. The amount of missing data correlates with greater phylogenetic divergence from taxa used to design the baits. This bait set hybridizes to a wide range of hemipteran taxa and specimens of varying quality, including dried specimens as old as 1973. Our estimated phylogeny yielded topologies consistent with other studies for most nodes and was strongly-supported. We also demonstrate that UCE loci are almost exclusively from the transcribed portion of the genome, thus data can be successfully integrated with existing genomic and transcriptomic resources for more comprehensive phylogenetic sampling, an important feature in the era of phylogenomics. UCE approaches can be used by other researchers for additional studies on hemipteran evolution and other research that requires well resolved phylogenies., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

The sperm whale, made famous by Moby Dick, is one of the most fascinating of all ocean-dwelling species given their unique life history, novel physiological adaptations to hunting squid at extreme ocean depths, and their position as one of the earliest branching toothed whales (Odontoceti). We assembled the sperm whale (Physeter macrocephalus) genome and resequenced individuals from multiple ocean basins to identify new candidate genes for adaptation to an aquatic environment and infer demographic history. Genes crucial for skin integrity appeared to be particularly important in both the sperm whale and other cetaceans. We also find sperm whales experienced a steep population decline during the early Pleistocene epoch. These genomic data add new comparative insight into the evolution of whales., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Accurate blood meal identification is critical to understand hematophagous vector-host relationships. This study describes a customizable Next-Generation Sequencing (NGS) approach to identify blood meals from Rhodnius pallescens (Hemiptera: Reduviidae) triatomines using multiple barcoded primers and existing software to pick operational taxonomic units and match sequences for blood meal identification. We precisely identified all positive control samples using this method and further examined 74 wild-caught R. pallescens samples. With this novel blood meal identification method, we detected 13 vertebrate species in the blood meals, as well as single and multiple blood meals in individual bugs. Our results demonstrate the reliability and descriptive uses of our method., (© The Authors 2017. Published by Oxford University Press on behalf of Entomological Society of America. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Molecular ecologists seek to genotype hundreds to thousands of loci from hundreds to thousands of individuals at minimal cost per sample. Current methods, such as restriction-site-associated DNA sequencing (RADseq) and sequence capture, are constrained by costs associated with inefficient use of sequencing data and sample preparation. Here, we introduce RADcap, an approach that combines the major benefits of RADseq (low cost with specific start positions) with those of sequence capture (repeatable sequencing of specific loci) to significantly increase efficiency and reduce costs relative to current approaches. RADcap uses a new version of dual-digest RADseq (3RAD) to identify candidate SNP loci for capture bait design and subsequently uses custom sequence capture baits to consistently enrich candidate SNP loci across many individuals. We combined this approach with a new library preparation method for identifying and removing PCR duplicates from 3RAD libraries, which allows researchers to process RADseq data using traditional pipelines, and we tested the RADcap method by genotyping sets of 96-384 Wisteria plants. Our results demonstrate that our RADcap method: (i) methodologically reduces (to <5%) and allows computational removal of PCR duplicate reads from data, (ii) achieves 80-90% reads on target in 11 of 12 enrichments, (iii) returns consistent coverage (≥4×) across >90% of individuals at up to 99.8% of the targeted loci, (iv) produces consistently high occupancy matrices of genotypes across hundreds of individuals and (v) costs significantly less than current approaches., (© 2016 John Wiley & Sons Ltd.)