Metabarcoding of insect-associated fungal communities: a comparison of internal transcribed spacer (ITS) and large-subunit (LSU) rRNA markers.

Full taxonomic characterisation of fungal communities is necessary for establishing ecological associations and early detection of pathogens and invasive species. Complex communities of fungi are regularly characterised by metabarcoding using the Internal Transcribed Spacer (ITS) and the Large-Subunit (LSU) gene of the rRNA locus, but reliance on a single short sequence fragment limits the confidence of identification. Here we link metabarcoding from the ITS2 and LSU D1-D2 regions to characterise fungal communities associated with bark beetles (Scolytinae), the likely vectors of several tree pathogens. Both markers revealed similar patterns of overall species richness and response to key variables (beetle species, forest type), but identification against the respective reference databases using various taxonomic classifiers revealed poor resolution towards lower taxonomic levels, especially the species level. Thus, Operational Taxonomic Units (OTUs) could not be linked via taxonomic classifiers across ITS and LSU fragments. However, using phylogenetic trees (focused on the epidemiologically important Sordariomycetes) we placed OTUs obtained with either marker relative to reference sequences of the entire rRNA cistron that includes both loci and demonstrated the largely similar phylogenetic distribution of ITS and LSU-derived OTUs. Sensitivity analysis of congruence in both markers suggested the biologically most defensible threshold values for OTU delimitation in Sordariomycetes to be 98% for ITS2 and 99% for LSU D1-D2. Studies of fungal communities using the canonical ITS barcode require corroboration across additional loci. Phylogenetic analysis of OTU sequences aligned to the full rRNA cistron shows higher success rate and greater accuracy of species identification compared to probabilistic taxonomic classifiers.

Coordinated RNA-Seq and peptidomics identify neuropeptides and G-protein coupled receptors (GPCRs) in the large pine weevil Hylobius abietis, a major forestry pest.

Hylobius abietis (Linnaeus), or large pine weevil (Coleoptera, Curculionidae), is a pest of European coniferous forests. In order to gain understanding of the functional physiology of this species, we have assembled a de novo transcriptome of H. abietis, from sequence data obtained by Next Generation Sequencing. In particular, we have identified genes encoding neuropeptides, peptide hormones and their putative G-protein coupled receptors (GPCRs) to gain insights into neuropeptide-modulated processes. The transcriptome was assembled de novo from pooled paired-end, sequence reads obtained from RNA from whole adults, gut and central nervous system tissue samples. Data analysis was performed on the transcripts obtained from the assembly including, annotation, gene ontology and functional assignment as well as transcriptome completeness assessment and KEGG pathway analysis. Pipelines were created using Bioinformatics tools and techniques for prediction and identification of neuropeptides and neuropeptide receptors. Peptidomic analysis was also carried out using a combination of MALDI-TOF as well as Q-Exactive Orbitrap mass spectrometry to confirm the identified neuropeptide. 41 putative neuropeptide families were identified in H. abietis, including Adipokinetic hormone (AKH), CAPA and DH31. Neuropeptide F, which has not been yet identified in the model beetle T. castaneum, was identified. Additionally, 24 putative neuropeptide and 9 leucine-rich repeat containing G protein coupled receptor-encoding transcripts were determined using both alignment as well as non-alignment methods. This information, submitted to the NCBI sequence read archive repository (SRA accession: SRP133355), can now be used to inform understanding of neuropeptide-modulated physiology and behaviour in H. abietis; and to develop specific neuropeptide-based tools for H. abietis control.

Metabarcoding of fungal communities associated with bark beetles.

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.

Responses of the two-spotted oak buprestid, Agrilus biguttatus (Coleoptera: Buprestidae), to host tree volatiles.

BACKGROUND: Agrilus bigutattus (Fabricius) is a forest pest of increasing importance in the United Kingdom. The larvae damage weakened native oaks and are thought to contribute to premature tree death. Suspected links with acute oak decline (AOD) are not yet confirmed, but AOD-predisposed trees appear to become more susceptible to A. biguttatus attack. Thus, management may be necessary for control of this insect. To explore the possibility of monitoring beetle populations by baited traps, the host tree volatiles regulating A. biguttatus-oak interactions were studied. RESULTS: Biologically active volatile organic compounds in dynamic headspace extracts of oak foliage and bark were identified initially by coupled gas chromatography-electroantennography (GC-EAG) and GC-mass spectrometry (GC-MS), and the structures were confirmed by GC coinjection with authentic compounds. Of two synthetic blends of these compounds comprising the active leaf volatiles, the simpler one containing three components evoked strongly positive behavioural responses in four-arm olfactometer tests with virgin females and males, although fresh leaf material was more efficient than the blend. The other blend, comprising a five-component mixture made up of bark volatiles, proved to be as behaviourally active for gravid females as bark tissue. CONCLUSIONS: These initial results on A. biguttatus chemical ecology reveal aspects of the role of attractive tree volatiles in the host-finding of beetles and underpin the development of semiochemically based surveillance strategies for this forest insect.

Accelerated species inventory on Madagascar using coalescent-based models of species delineation.

High-throughput DNA sequencing has the potential to accelerate species discovery if it is able to recognize evolutionary entities from sequence data that are comparable to species. The general mixed Yule-coalescent (GMYC) model estimates the species boundary from DNA surveys by identifying independently evolving lineages as a transition from coalescent to speciation branching patterns on a phylogenetic tree. Applied here to 12 families from 4 orders of insects in Madagascar, we used the model to delineate 370 putative species from mitochondrial DNA sequence variation among 1614 individuals. These were compared with data from the nuclear genome and morphological identification and found to be highly congruent (98% and 94%). We developed a modified GMYC that allows for a variable transition from coalescent to speciation among lineages. This revised model increased the congruence with morphology (97%), suggesting that a variable threshold better reflects the clustering of sequence data into biological species. Local endemism was pronounced in all 5 insect groups. Most species (60-91%) and haplotypes (88-99%) were found at only 1 of the 5 study sites (40-1000 km apart). This pronounced endemism resulted in a 37% increase in species numbers using diagnostic nucleotides in a population aggregation analysis. Sample sizes between 7 and 10 individuals represented a threshold above which there was minimal increase in genetic diversity, broadly agreeing with coalescent theory and other empirical studies. Our results from > 1.4 Mb of empirical data suggest that the GMYC model captures species boundaries comparable to those from traditional methods without the need for prior hypotheses of population coherence. This provides a method of species discovery and biodiversity assessment using single-locus data from mixed or environmental samples while building a globally available taxonomic database for future identifications.

A comprehensive phylogenetic analysis of termites (Isoptera) illuminates key aspects of their evolutionary biology.

The first comprehensive combined molecular and morphological phylogenetic analysis of the major groups of termites is presented. This was based on the analysis of three genes (cytochrome oxidase II, 12S and 28S) and worker characters for approximately 250 species of termites. Parsimony analysis of the aligned dataset showed that the monophyly of Hodotermitidae, Kalotermitidae and Termitidae were well supported, while Termopsidae and Rhinotermitidae were both paraphyletic on the estimated cladogram. Within Termitidae, the most diverse and ecologically most important family, the monophyly of Macrotermitinae, Foraminitermitinae, Apicotermitinae, Syntermitinae and Nasutitermitinae were all broadly supported, but Termitinae was paraphyletic. The pantropical genera Termes, Amitermes and Nasutitermes were all paraphyletic on the estimated cladogram, with at least 17 genera nested within Nasutitermes, given the presently accepted generic limits. Key biological features were mapped onto the cladogram. It was not possible to reconstruct the evolution of true workers unambiguously, as it was as parsimonious to assume a basal evolution of true workers and subsequent evolution of pseudergates, as to assume a basal condition of pseudergates and subsequent evolution of true workers. However, true workers were only found in species with either separate- or intermediate-type nests, so that the mapping of nest habit and worker type onto the cladogram were perfectly correlated. Feeding group evolution, however, showed a much more complex pattern, particularly within the Termitidae, where it proved impossible to estimate unambiguously the ancestral state within the family (which is associated with the loss of worker gut flagellates). However, one biologically plausible optimization implies an initial evolution from wood-feeding to fungus-growing, proposed as the ancestral condition within the Termitidae, followed by the very early evolution of soil-feeding and subsequent re-evolution of wood-feeding in numerous lineages.

A molecular phylogenetic analysis of the Scarabaeinae (dung beetles).

The dung beetles (Scarabaeinae) include ca. 5000 species and exhibit a diverse array of morphologies and behaviors. This variation presumably reflects the adaptation to a diversity of food types and the different strategies used to avoid competition for vertebrate dung, which is the primary breeding environment for most species. The current classification gives great weight to the major behavioral types, separating the ball rollers and the tunnelers, but existing phylogenetic studies have been based on limited taxonomic or biogeographic sampling and have been contradictory. Here, we present a molecular phylogenetic analysis of 214 species of Scarabaeinae, representing all 12 traditionally recognized tribes and six biogeographical regions, using partial gene sequences from one nuclear (28S) and two mitochondrial (cox1, rrnL) genes. Length variation in 28S (588-621 bp) and rrnL (514-523 bp) was subjected to a thorough evaluation of alternative alignments, gap-coding methods, and tree searches using model-based (Bayesian and likelihood), maximum parsimony, and direct optimization analyses. The small-bodied, non-dung-feeding Sarophorus+Coptorhina were basal in all reconstructions. These were closely related to rolling Odontoloma+Dicranocara, suggesting an early acquisition of rolling behavior. Smaller tribes and most genera were monophyletic, while Canthonini and Dichotomiini each consisted of multiple paraphyletic lineages at hierarchical levels equivalent to the smaller tribes. Plasticity of rolling and tunneling was evidenced by a lack of monophyly (S-H test, p > 0.05) and several reversals within clades. The majority of previously unrecognized clades were geographical, including the well-supported Neotropical Phanaeini+Eucraniini, and a large Australian clade of rollers as well as tunneling Coptodactyla and Demarziella. Only three lineages, Gymnopleurini, Copris+Microcopris and Onthophagus, were widespread and therefore appear to be dispersive at a global scale. A reconstruction of biogeographical characters recovered 38-48 transitions between regions and an African origin for most lineages. Dispersal-vicariance analysis supported an African origin with links to all other regions and little back-migration. Our results provide a new synthesis of global-scale dung beetle evolution, demonstrating the great plasticity of behavioral and morphological traits and the importance of biogeographic distributions as the basis for a new classification.