An aggregated dataset of serial morbidity and titer measurements from influenza A virus-infected ferrets

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.

Highly pathogenic avian influenza A(H5N1) virus of clade 2.3.4.4b isolated from a human case in Chile causes fatal disease and transmits between co-housed ferrets

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 reassortment with four gene segments (PB1, PB2, NP, MP) from North America 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.

Key considerations to improve the normalization, interpretation and reproducibility of morbidity data in mammalian models of viral disease

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.

Exploring associations between viral titer measurements and disease outcomes in ferrets inoculated with 125 contemporary influenza A viruses

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.

A naturally occurring HA-stabilizing amino acid (HA1-Y17) in an A(H9N2) low-pathogenic influenza virus contributes to airborne transmission

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.

Influenza A Virus Multicycle Replication Yields Comparable Viral Population Emergence in Human Respiratory and Ocular Cell Types

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.

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.

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.

How microclimatic variables and blood meal sources influence Rhodnius prolixus abundance and Trypanosoma cruzi infection in Attalea butyracea and Elaeis guineensis palms?

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.

Ultraconserved element bait set for trypanosomatida target enrichment and phylogenetics

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.

Microbiota of Four Tissue Types in American Alligators (Alligator mississippiensis) Following Extended Dietary Selenomethionine Exposure

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.

Co-occurrence of antibiotic, biocide, and heavy metal resistance genes in bacteria from metal and radionuclide contaminated soils at the Savannah River Site

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.

Identification and characterization of microRNAs (miRNAs) and their transposable element origins in the saltwater crocodile, Crocodylus porosus

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.

Mitochondrial, metagenomic, and phylogenetic analysis of the ground beetle Harpalus pensylvanicus (Coleoptera: Carabidae)

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%).

Population genetics of two chromatic morphs of the Chagas disease vector Rhodnius pallescens Barber, 1932 in Panamá

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.

Comparison of the ruminal and fecal microbiotas in beef calves supplemented or not with concentrate

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.

Genome comparison and transcriptome analysis of the invasive brown root rot pathogen, Phellinus noxius, from different geographic regions reveals potential enzymes associated with degradation of different wood substrates

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.

Genomic mutations after multigenerational exposure of Caenorhabditis elegans to pristine and sulfidized silver nanoparticles

Our previous study showed heritable reproductive toxicity in the nematode Caenorhabditis elegans after multigenerational exposure to AgNO₃ 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₃₀ of AgNO₃, 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₀ and F₁₀ 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₃ 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₃-induced multigenerational and transgenerational reproductive toxicity.

Regional biogeography of microbiota composition in the Chagas disease vector Rhodnius pallescens

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.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext)

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.

Adapterama III: Quadruple-indexed, double/triple-enzyme RADseq libraries (2RAD/3RAD)

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.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix)

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_Demuxy.

Analysis of the Rumen Microbiota of Beef Calves Supplemented During the Suckling Phase

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.

Generalist host species drive Trypanosoma cruzi vector infection in oil palm plantations in the Orinoco region, Colombia

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.

The Novel Evolution of the Sperm Whale Genome

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.

Conflicting Evolutionary Histories of the Mitochondrial and Nuclear Genomes in New World Myotis Bats

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.

Blood Meal Source Characterization Using Illumina Sequencing in the Chagas Disease Vector Rhodnius pallescens (Hemiptera: Reduviidae) in Panamá

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.