Robert S. Cornman

Critical considerations for the application of environmental DNA methods to detect aquatic species

Summary Species detection using environmental DNA (eDNA) has tremendous potential for contributing to the understanding of the ecology and conservation of aquatic species. Detecting species using eDNA methods, rather than directly sampling the organisms, can reduce impacts on sensitive species and increase the power of field surveys for rare and elusive species. The sensitivity of eDNA methods, however, requires a heightened awareness and attention to quality assurance and quality control protocols. Additionally, the interpretation of eDNA data demands careful consideration of multiple factors. As eDNA methods have grown in application, diverse approaches have been implemented to address these issues. With interest in eDNA continuing to expand, supportive guidelines for undertaking eDNA studies are greatly needed. Environmental DNA researchers from around the world have collaborated to produce this set of guidelines and considerations for implementing eDNA methods to detect aquatic macroorganisms. Critical considerations for study design include preventing contamination in the field and the laboratory, choosing appropriate sample analysis methods, validating assays, testing for sample inhibition and following minimum reporting guidelines. Critical considerations for inference include temporal and spatial processes, limits of correlation of eDNA with abundance, uncertainty of positive and negative results, and potential sources of allochthonous DNA. We present a synthesis of knowledge at this stage for application of this new and powerful detection method.

Standard methods for molecular research in Apis mellifera

Summary From studies of behaviour, chemical communication, genomics and developmental biology, among many others, honey bees have long been a key organism for fundamental breakthroughs in biology. With a genome sequence in hand, and much improved genetic tools, honey bees are now an even more appealing target for answering the major questions of evolutionary biology, population structure, and social organization. At the same time, agricultural incentives to understand how honey bees fall prey to disease, or evade and survive their many pests and pathogens, have pushed for a genetic understanding of individual and social immunity in this species. Below we describe and reference tools for using modern molecular-biology techniques to understand bee behaviour, health, and other aspects of their biology. We focus on DNA and RNA techniques, largely because techniques for assessing bee proteins are covered in detail in Hartfelder et al. (2013). We cover practical needs for bee sampling, transport, and storage, and then discuss a range of current techniques for genetic analysis. We then provide a roadmap for genomic resources and methods for studying bees, followed by specific statistical protocols for population genetics, quantitative genetics, and phylogenetics. Finally, we end with three important tools for predicting gene regulation and function in honey bees: Fluorescence in situ hybridization (FISH), RNA interference (RNAi), and the estimation of chromosomal methylation and its role in epigenetic gene regulation.

Systemic Spread and Propagation of a Plant-Pathogenic Virus in European Honeybees, Apis mellifera

ABSTRACT Emerging and reemerging diseases that result from pathogen host shifts are a threat to the health of humans and their domesticates. RNA viruses have extremely high mutation rates and thus represent a significant source of these infectious diseases. In the present study, we showed that a plant-pathogenic RNA virus, tobacco ringspot virus (TRSV), could replicate and produce virions in honeybees, Apis mellifera, resulting in infections that were found throughout the entire body. Additionally, we showed that TRSV-infected individuals were continually present in some monitored colonies. While intracellular life cycle, species-level genetic variation, and pathogenesis of the virus in honeybee hosts remain to be determined, the increasing prevalence of TRSV in conjunction with other bee viruses from spring toward winter in infected colonies was associated with gradual decline of host populations and winter colony collapse, suggesting the negative impact of the virus on colony survival. Furthermore, we showed that TRSV was also found in ectoparasitic Varroa mites that feed on bee hemolymph, but in those instances the virus was restricted to the gastric cecum of Varroa mites, suggesting that Varroa mites may facilitate the spread of TRSV in bees but do not experience systemic invasion. Finally, our phylogenetic analysis revealed that TRSV isolates from bees, bee pollen, and Varroa mites clustered together, forming a monophyletic clade. The tree topology indicated that the TRSVs from arthropod hosts shared a common ancestor with those from plant hosts and subsequently evolved as a distinct lineage after transkingdom host alteration. This study represents a unique example of viruses with host ranges spanning both the plant and animal kingdoms. IMPORTANCE Pathogen host shifts represent a major source of new infectious diseases. Here we provide evidence that a pollen-borne plant virus, tobacco ringspot virus (TRSV), also replicates in honeybees and that the virus systemically invades and replicates in different body parts. In addition, the virus was detected inside the body of parasitic Varroa mites, which consume bee hemolymph, suggesting that Varroa mites may play a role in facilitating the spread of the virus in bee colonies. This study represents the first evidence that honeybees exposed to virus-contaminated pollen could also be infected and raises awareness of potential risks of new viral disease emergence due to host shift events. About 5% of known plant viruses are pollen transmitted, and these are potential sources of future host-jumping viruses. The findings from this study showcase the need for increased surveillance for potential host-jumping events as an integrated part of insect pollinator management programs. Pathogen host shifts represent a major source of new infectious diseases. Here we provide evidence that a pollen-borne plant virus, tobacco ringspot virus (TRSV), also replicates in honeybees and that the virus systemically invades and replicates in different body parts. In addition, the virus was detected inside the body of parasitic Varroa mites, which consume bee hemolymph, suggesting that Varroa mites may play a role in facilitating the spread of the virus in bee colonies. This study represents the first evidence that honeybees exposed to virus-contaminated pollen could also be infected and raises awareness of potential risks of new viral disease emergence due to host shift events. About 5% of known plant viruses are pollen transmitted, and these are potential sources of future host-jumping viruses. The findings from this study showcase the need for increased surveillance for potential host-jumping events as an integrated part of insect pollinator management programs.

Population-genomic variation within RNA viruses of the Western honey bee, Apis mellifera, inferred from deep sequencing

BackgroundDeep sequencing of viruses isolated from infected hosts is an efficient way to measure population-genetic variation and can reveal patterns of dispersal and natural selection. In this study, we mined existing Illumina sequence reads to investigate single-nucleotide polymorphisms (SNPs) within two RNA viruses of the Western honey bee (Apis mellifera), deformed wing virus (DWV) and Israel acute paralysis virus (IAPV). All viral RNA was extracted from North American samples of honey bees or, in one case, the ectoparasitic mite Varroa destructor.ResultsCoverage depth was generally lower for IAPV than DWV, and marked gaps in coverage occurred in several narrow regions (< 50 bp) of IAPV. These coverage gaps occurred across sequencing runs and were virtually unchanged when reads were re-mapped with greater permissiveness (up to 8% divergence), suggesting a recurrent sequencing artifact rather than strain divergence. Consensus sequences of DWV for each sample showed little phylogenetic divergence, low nucleotide diversity, and strongly negative values of Fu and Li’s D statistic, suggesting a recent population bottleneck and/or purifying selection. The Kakugo strain of DWV fell outside of all other DWV sequences at 100% bootstrap support. IAPV consensus sequences supported the existence of multiple clades as had been previously reported, and Fu and Li’s D was closer to neutral expectation overall, although a sliding-window analysis identified a significantly positive D within the protease region, suggesting selection maintains diversity in that region. Within-sample mean diversity was comparable between the two viruses on average, although for both viruses there was substantial variation among samples in mean diversity at third codon positions and in the number of high-diversity sites. FST values were bimodal for DWV, likely reflecting neutral divergence in two low-diversity populations, whereas IAPV had several sites that were strong outliers with very low FST.ConclusionsThis initial survey of genetic variation within honey bee RNA viruses suggests future directions for studies examining the underlying causes of population-genetic structure in these economically important pathogens.

Transcriptional Response of Honey Bee Larvae Infected with the Bacterial Pathogen Paenibacillus larvae

American foulbrood disease of honey bees is caused by the bacterium Paenibacillus larvae. Infection occurs per os in larvae and systemic infection requires a breaching of the host peritrophic matrix and midgut epithelium. Genetic variation exists for both bacterial virulence and host resistance, and a general immunity is achieved by larvae as they age, the basis of which has not been identified. To quickly identify a pool of candidate genes responsive to P. larvae infection, we sequenced transcripts from larvae inoculated with P. larvae at 12 hours post-emergence and incubated for 72 hours, and compared expression levels to a control cohort. We identified 75 genes with significantly higher expression and six genes with significantly lower expression. In addition to several antimicrobial peptides, two genes encoding peritrophic-matrix domains were also up-regulated. Extracellular matrix proteins, proteases/protease inhibitors, and members of the Osiris gene family were prevalent among differentially regulated genes. However, analysis of Drosophila homologs of differentially expressed genes revealed spatial and temporal patterns consistent with developmental asynchrony as a likely confounder of our results. We therefore used qPCR to measure the consistency of gene expression changes for a subset of differentially expressed genes. A replicate experiment sampled at both 48 and 72 hours post infection allowed further discrimination of genes likely to be involved in host response. The consistently responsive genes in our test set included a hymenopteran-specific protein tyrosine kinase, a hymenopteran specific serine endopeptidase, a cytochrome P450 (CYP9Q1), and a homolog of trynity, a zona pellucida domain protein. Of the known honey bee antimicrobial peptides, apidaecin was responsive at both time-points studied whereas hymenoptaecin was more consistent in its level of change between biological replicates and had the greatest increase in expression by RNA-seq analysis.

Transcriptome resources for the frogs Lithobates clamitans and Pseudacris regilla, emphasizing antimicrobial peptides and conserved loci for phylogenetics

We developed genetic resources for two North American frogs, Lithobates clamitans and Pseudacris regilla, widespread native amphibians that are potential indicator species of environmental health. For both species, mRNA from multiple tissues was sequenced using 454 technology. De novo assemblies with Mira3 resulted in 50 238 contigs (N50 = 687 bp) and 48 213 contigs (N50 = 686 bp) for L. clamitans and P. regilla, respectively, after clustering with CD‐Hit‐EST and purging contigs below 200 bp. We performed BLASTX similarity searches against the Xenopus tropicalis proteome and, for predicted ORFs, HMMER similarity searches against the Pfam‐A database. Because there is broad interest in amphibian immune factors, we manually annotated putative antimicrobial peptides. To identify conserved regions suitable for amplicon resequencing across a broad taxonomic range, we performed an additional assembly of public short‐read transcriptome data derived from two species of the genus Rana and identified reciprocal best TBLASTX matches among all assemblies. Although P. regilla, a hylid frog, is substantially more diverged from the ranid species, we identified 56 genes that were sufficiently conserved to allow nondegenerate primer design with Primer3. In addition to providing a foundation for comparative genomics and quantitative gene expression analysis, our results enable quick development of nuclear sequence‐based markers for phylogenetics or population genetics.

Characterization of a Novel Hepadnavirus in the White Sucker (Catostomus commersonii) from the Great Lakes Region of the United States

ABSTRACT The white sucker Catostomus commersonii is a freshwater teleost often utilized as a resident sentinel. Here, we sequenced the full genome of a hepatitis B-like virus that infects white suckers from the Great Lakes Region of the United States. Dideoxy sequencing confirmed that the white sucker hepatitis B virus (WSHBV) has a circular genome (3,542 bp) with the prototypical codon organization of hepadnaviruses. Electron microscopy demonstrated that complete virions of approximately 40 nm were present in the plasma of infected fish. Compared to avi- and orthohepadnaviruses, sequence conservation of the core, polymerase, and surface proteins was low and ranged from 16 to 27% at the amino acid level. An X protein homologue common to the orthohepadnaviruses was not present. The WSHBV genome included an atypical, presumptively noncoding region absent in previously described hepadnaviruses. Phylogenetic analyses confirmed WSHBV as distinct from previously documented hepadnaviruses. The level of divergence in protein sequences between WSHBV and other hepadnaviruses and the identification of an HBV-like sequence in an African cichlid provide evidence that a novel genus of the family Hepadnaviridae may need to be established that includes these hepatitis B-like viruses in fishes. Viral transcription was observed in 9.5% (16 of 169) of white suckers evaluated. The prevalence of hepatic tumors in these fish was 4.9%, and only 2.4% of fish were positive for both virus and hepatic tumors. These results are not sufficient to draw inferences regarding the association of WSHBV and carcinogenesis in white sucker. IMPORTANCE We report the first full-length genome of a hepadnavirus from fishes. Phylogenetic analysis of this genome indicates divergence from genomes of previously described hepadnaviruses from mammalian and avian hosts and supports the creation of a novel genus. The discovery of this novel virus may better our understanding of the evolutionary history of hepatitis B-like viruses of other hosts. In fishes, knowledge of this virus may provide insight regarding possible risk factors associated with hepatic neoplasia in the white sucker. This may also offer another model system for mechanistic research.

A Comparison of Honey Bee-Collected Pollen From Working Agricultural Lands Using Light Microscopy and ITS Metabarcoding.

Abstract Taxonomic identification of pollen has historically been accomplished via light microscopy but requires specialized knowledge and reference collections, particularly when identification to lower taxonomic levels is necessary. Recently, next-generation sequencing technology has been used as a cost-effective alternative for identifying bee-collected pollen; however, this novel approach has not been tested on a spatially or temporally robust number of pollen samples. Here, we compare pollen identification results derived from light microscopy and DNA sequencing techniques with samples collected from honey bee colonies embedded within a gradient of intensive agricultural landscapes in the Northern Great Plains throughout the 2010–2011 growing seasons. We demonstrate that at all taxonomic levels, DNA sequencing was able to discern a greater number of taxa, and was particularly useful for the identification of infrequently detected species. Importantly, substantial phenological overlap did occur for commonly detected taxa using either technique, suggesting that DNA sequencing is an appropriate, and enhancing, substitutive technique for accurately capturing the breadth of bee-collected species of pollen present across agricultural landscapes. We also show that honey bees located in high and low intensity agricultural settings forage on dissimilar plants, though with overlap of the most abundantly collected pollen taxa. We highlight practical applications of utilizing sequencing technology, including addressing ecological issues surrounding land use, climate change, importance of taxa relative to abundance, and evaluating the impact of conservation program habitat enhancement efforts.

Z chromosome divergence, polymorphism and relative effective population size in a genus of lekking birds

Sex chromosomes contribute disproportionately to species boundaries as they diverge faster than autosomes and often have reduced diversity. Their hemizygous nature contributes to faster divergence and reduced diversity, as do some types of selection. In birds, other factors (mating system and bottlenecks) can further decrease the effective population size of Z-linked loci and accelerate divergence (Fast-Z). We assessed Z-linked divergence and effective population sizes for two polygynous sage-grouse species and compared them to estimates from birds with various mating systems. We found lower diversity and higher FST for Z-linked loci than for autosomes, as expected. The πZ/πA ratio was 0.38 in Centrocercus minimus, 0.48 in Centrocercus urophasianus and 0.59 in a diverged, parapatric population of C. urophasianus, a broad range given the mating system among these groups is presumably equivalent. The full data set had unequal males and females across groups, so we compared an equally balanced reduced set of C. minimus and individuals pooled from both C. urophasianus subgroups recovering similar estimates: 0.54 for C. urophasianus and 0.38 for C. minimus. We provide further evidence that NeZ/NeA in birds is often lower than expected under random mating or monogamy. The lower ratio in C. minimus could be a consequence of stronger selection or drift acting on Z loci during speciation, as this species differs strongly from C. urophasianus in sexually selected characters with minimal mitochondrial divergence. As C. minimus also exhibited lower genomic diversity, it is possible that a more severe demographic history may contribute to its lower ratio.

Transcriptomic analysis of the mussel Elliptio complanata identifies candidate stress-response genes and an abundance of novel or noncoding transcripts.

Mussels are useful indicator species of environmental stress and degradation, and the global decline in freshwater mussel diversity and abundance is of conservation concern. Elliptio complanata is a common freshwater mussel of eastern North America that can serve both as an indicator and as an experimental model for understanding mussel physiology and genetics. To support genetic components of these research goals, we assembled transcriptome contigs from Illumina paired-end reads. Despite efforts to collapse similar contigs, the final assembly was in excess of 136,000 contigs with an N50 of 982 bp. Even so, comparisons to the CEGMA database of conserved eukaryotic genes indicated that ∼20% of genes remain unrepresented. However, numerous candidate stress-response genes were present, and we identified lineage-specific patterns of diversification among molluscs for cytochrome P450 detoxification genes and two saccharide-modifying enzymes: 1,3 beta-galactosyltransferase and fucosyltransferase. Less than a quarter of contigs had protein-level similarity based on modest BLAST and Hmmer3 statistical thresholds. These results add comparative genomic resources for molluscs and suggest a wealth of novel proteins and noncoding transcripts.