Arthropod diversity in the alpine tundra using metabarcoding: Spatial and temporal differences in alpha- and beta-diversity.

All ecosystems face ecological challenges in this century. Therefore, it is becoming increasingly important to understand the ecology and degree of local adaptation of functionally important Arctic-alpine biomes by looking at the most diverse taxon of metazoans: the Arthropoda. This is the first study to utilize metabarcoding in the Alpine tundra, providing insights into the effects of micro-environmental parameters on alpha- and beta-diversity of arthropods in such unique environments. To characterize arthropod diversity, pitfall traps were set at three middle-alpine sampling sites in the Scandinavian mountain range in Norway during the snow-free season in 2015. A metabarcoding approach was then used to determine the small-scale biodiversity patterns of arthropods in the Alpine tundra. All DNA was extracted directly from the preservative EtOH from 27 pitfall traps. In order to identify the controlling environmental conditions, all sampling locations were equipped with automatic data loggers for permanent measurement of the microenvironmental conditions. The variables measured were: air temperature [°C] at 15 cm height, soil temperature [°C] at 15 cm depth, and soil moisture [vol.%] at 15 cm depth. A total of 233 Arthropoda OTUs were identified. The number of unique OTUs found per sampling location (ridge, south-facing slope, and depression) was generally higher than the OTUs shared between the sampling locations, demonstrating that niche features greatly impact arthropod community structure. Our findings emphasize the fine-scale heterogeneity of arctic-alpine ecosystems and provide evidence for trait-based and niche-driven adaptation. The spatial and temporal differences in arthropod diversity were best explained by soil moisture and soil temperature at the respective locations. Furthermore, our results show that arthropod diversity is underestimated in alpine-tundra ecosystems using classical approaches and highlight the importance of integrating long-term functional environmental data and modern taxonomic techniques into biodiversity research to expand our ecological understanding of fine- and meso-scale biogeographical patterns.

Metabarcoding dietary analysis in the insectivorous bat Nyctalusleisleri and implications for conservation.

In this study, we aim to uncover diet preferences for the insectivorous bat Nyctalusleisleri (Leisler's bat, the lesser noctule) and to provide recommendations for conservation of the species, based on the analysis of prey source habitats. Using a novel guano trap, we sampled bat faeces at selected roosts in a forest in Germany and tested two mitochondrial markers (COI and 16S) and three primer pairs for the metabarcoding of bat faecal pellets. We found a total of 17 arthropod prey orders comprising 358 species in N.leisleri guano. The most diverse orders were Lepidoptera (126 species), Diptera (86 species) and Coleoptera (48 species), followed by Hemiptera (28 species), Trichoptera (16 species), Neuroptera (15 species) and Ephemeroptera (10 species), with Lepidoptera species dominating in spring and Diptera in summer. Based on the ecological requirements of the most abundant arthropod species found in the bat guano, we propose some recommendations for the conservation of N.leisleri that are relevant for other insectivorous bat species.

Systematic review and meta-analysis: Water type and temperature affect environmental DNA decay.

Environmental DNA (eDNA) has been used in a variety of ecological studies and management applications. The rate at which eDNA decays has been widely studied but at present it is difficult to disentangle study-specific effects from factors that universally affect eDNA degradation. To address this, a systematic review and meta-analysis was conducted on aquatic eDNA studies. Analysis revealed eDNA decayed faster at higher temperatures and in marine environments (as opposed to freshwater). DNA type (mitochondrial or nuclear) and fragment length did not affect eDNA decay rate, although a preference for <200 bp sequences in the available literature means this relationship was not assessed with longer sequences (e.g. >800 bp). At present, factors such as ultraviolet light, pH, and microbial load lacked sufficient studies to feature in the meta-analysis. Moving forward, we advocate researching these factors to further refine our understanding of eDNA decay in aquatic environments.

Metabarcoding Malaise traps and soil eDNA reveals seasonal and local arthropod diversity shifts.

Forest habitats host enormous diversity, but little is known about the seasonal turnover of arthropod species between the above- and below ground forest layers. In this study, we used metabarcoding approaches to uncover arthropod diversity in different forest types and seasons. Our study shows that metabarcoding soil eDNA and Malaise trap bulk samples can provide valuable insights into the phenology and life cycles of arthropods. We found major differences in arthropod species diversity between soil samples and Malaise traps, with only 11.8% species overlap. Higher diversity levels were found in Malaise traps in summer whereas soil samples showed a diversity peak in winter, highlighting the seasonal habitat preferences and life strategies of arthropods. We conclude that collecting time series of bulk arthropod samples and eDNA in the same locations provides a more complete picture of local arthropod diversity and turnover rates and may provide valuable information on climate induced phenological shifts for long-term monitoring.

Whole Genome Sequencing of Hepatitis A Virus Using a PCR-Free Single-Molecule Nanopore Sequencing Approach.

Hepatitis A virus (HAV) is one of the most common causes of acute viral hepatitis in humans. Although HAV has a relatively small genome, there are several factors limiting whole genome sequencing such as PCR amplification artefacts and ambiguities in de novo assembly. The recently developed Oxford Nanopore technologies (ONT) allows single-molecule sequencing of long-size fragments of DNA or RNA using PCR-free strategies. We have sequenced the whole genome of HAV using a PCR-free approach by direct reverse-transcribed sequencing. We were able to sequence HAV cDNA and obtain reads over 7 kilobases in length containing almost the whole genome of the virus. The comparison of these raw long nanopore reads with the HAV reference wild type revealed a nucleotide sequence identity between 81.1 and 96.6%. By de novo assembly of all HAV reads we obtained a consensus sequence of 7362 bases, with a nucleotide sequence identity of 99.0% with the genome of the HAV strain pHM175/18f. When the assembly was performed using as reference the HAV strain pHM175/18f a consensus with a sequence similarity of 99.8 % was obtained. We have also used an ONT amplicon-based assay to sequence two fragments of the VP3 and VP1 regions which showed a sequence similarity of 100% with matching regions of the consensus sequence obtained using the direct cDNA sequencing approach. This study showed the applicability of ONT sequencing technologies to obtain the whole genome of HAV by direct cDNA nanopore sequencing, highlighting the utility of this PCR-free approach for HAV characterization and potentially other viruses of the Picornaviridae family.

The Application of Nanopore Sequencing Technology to the Study of Dinoflagellates: A Proof of Concept Study for Rapid Sequence-Based Discrimination of Potentially Harmful Algae.

Harmful algal blooms (HABs) are a naturally occurring global phenomena that have the potential to impact fisheries, leisure and ecosystems, as well as posing a significant hazard to animal and human health. There is significant interest in the development and application of methodologies to study all aspects of the causative organisms and toxins associated with these events. This paper reports the first application of nanopore sequencing technology for the detection of eukaryotic harmful algal bloom organisms. The MinION sequencing platform from Oxford Nanopore technologies provides long read sequencing capabilities in a compact, low cost, and portable format. In this study we used the MinION to sequence long-range PCR amplicons from multiple dinoflagellate species with a focus on the genus Alexandrium. Primers applicable to a wide range of dinoflagellates were selected, meaning that although the study was primarily focused on Alexandrium the applicability to three additional genera of toxic algae, namely; Gonyaulax, Prorocentrum, and Lingulodinium was also demonstrated. The amplicon generated here spanned approximately 3 kb of the rDNA cassette, including most of the 18S, the complete ITS1, 5.8S, ITS2 and regions D1 and D2 of the 28S. The inclusion of barcode genes as well as highly conserved regions resulted in identification of organisms to the species level. The analysis of reference cultures resulted in over 99% of all sequences being attributed to the correct species with an average identity above 95% from a reference list of over 200 species (see Supplementary Material 1). The use of mock community analysis within environmental samples highlighted that complex matrices did not prevent the ability to distinguish between phylogenetically similar species. Successful identification of causative organisms in environmental samples during natural toxic events further highlighted the potential of the assay. This study proves the suitability of nanopore sequencing technology for taxonomic identification of harmful algal bloom organisms and acquisition of data relevant to the World Health Organisations "one health" approach to marine monitoring.

Comparing diversity levels in environmental samples: DNA sequence capture and metabarcoding approaches using 18S and COI genes.

Environmental DNA studies targeting multiple taxa using metabarcoding provide remarkable insights into levels of species diversity in any habitat. The main drawbacks are the presence of primer bias and difficulty in identifying rare species. We tested a DNA sequence-capture method in parallel with the metabarcoding approach to reveal possible advantages of one method over the other. Both approaches were performed using the same eDNA samples and the same 18S and COI regions, followed by high throughput sequencing. Metabarcoded eDNA libraries were PCR amplified with one primer pair from 18S and COI genes. DNA sequence-capture libraries were enriched with 3,639 baits targeting the same gene regions. We tested amplicon sequence variants (ASVs) and operational taxonomic units (OTUs) in silico approaches for both markers and methods, using for this purpose the metabarcoding data set. ASVs methods uncovered more species for the COI gene, whereas the opposite occurred for the 18S gene, suggesting that clustering reads into OTUs could bias diversity richness especially using 18S with relaxed thresholds. Additionally, metabarcoding and DNA sequence-capture recovered 80%-90% of the control sample species. DNA sequence-capture was 8x more expensive, nonetheless it identified 1.5x more species for COI and 13x more genera for 18S than metabarcoding. Both approaches offer reliable results, sharing ca. 40% species and 72% families and retrieve more taxa when nuclear and mitochondrial markers are combined. eDNA metabarcoding is quite well established and low-cost, whereas DNA-sequence capture for biodiversity assessment is still in its infancy, is more time-consuming but provides more taxonomic assignments.

Assessing insect biodiversity with automatic light traps in Brazil: Pearls and pitfalls of metabarcoding samples in preservative ethanol.

Automated species identification based on data produced with metabarcoding offers an alternative for assessing biodiversity of bulk insect samples obtained with traps. We used a standard two-step PCR approach to amplify a 313 bp fragment of the barcoding region of the mitochondrial COI gene. The PCR products were sequenced on an Illumina MiSeq platform, and the OTUs production and taxonomic identifications were performed with a customized pipeline and database. The DNA used in the PCR procedures was extracted directly from the preservative ethanol of bulk insect samples obtained with automatic light traps in 12 sampling areas located in different biomes of Brazil, during wet and dry seasons. Agricultural field and forest edge habitats were collected for all sampling areas. A total of 119 insect OTUs and nine additional OTUs assigned to other arthropod taxa were obtained at a ≥97% sequence similarity level. The alpha and beta diversity analyses comparing biomes, habitats, and seasons were mostly inconclusive, except for a significant difference in beta diversity between biomes. In this study, we were able to metabarcode and HTS adult insects from their preservative medium. Notwithstanding, our results underrepresent the true magnitude of insect diversity expected from samples obtained with automatic light traps in Brazil. Although biological and technical factors might have impacted our results, measures to optimize and standardize eDNA HTS should be in place to improve taxonomic coverage of samples of unknown diversity and stored in suboptimal conditions, which is the case of most eDNA samples.

Spatial and temporal dynamics of Antarctic shallow soft-bottom benthic communities: ecological drivers under climate change.

BACKGROUND: Marine soft sediments are some of the most widespread habitats in the ocean, playing a vital role in global carbon cycling, but are amongst the least studied with regard to species composition and ecosystem functioning. This is particularly true of the Polar Regions, which are currently undergoing rapid climate change, the impacts of which are poorly understood. Compared to other latitudes, Polar sediment habitats also experience additional environmental drivers of strong seasonality and intense disturbance from iceberg scouring, which are major structural forces for hard substratum communities. This study compared sediment assemblages from two coves, near Rothera Point, Antarctic Peninsula, 67°S in order to understand the principal drivers of community structure, for the first time, evaluating composition across all size classes from mega- to micro-fauna. RESULTS: Morpho-taxonomy identified 77 macrofaunal species with densities of 464-16,084 individuals m-2. eDNA metabarcoding of microfauna, in summer only, identified a higher diversity, 189 metazoan amplicon sequence variants (ASVs) using the 18S ribosomal RNA and 249 metazoan ASVs using the mitochondrial COI gene. Both techniques recorded a greater taxonomic diversity in South Cove than Hangar Cove, with differences in communities between the coves, although the main taxonomic drivers varied between techniques. Morphotaxonomy identified the main differences between coves as the mollusc, Altenaeum charcoti, the cnidarian Edwardsia sp. and the polychaetes from the family cirratulidae. Metabarcoding identified greater numbers of species of nematodes, crustaceans and Platyhelminthes in South Cove, but more bivalve species in Hangar Cove. There were no detectable differences in community composition, measured through morphotaxonomy, between seasons, years or due to iceberg disturbance. CONCLUSIONS: This study found that unlike hard substratum communities the diversity of Antarctic soft sediment communities is correlated with the same factors as other latitudes. Diversity was significantly correlated with grain size and organic content, not iceberg scour. The increase in glacial sediment input as glaciers melt, may therefore be more important than increased iceberg disturbance.

Metabarcoding Marine Sediments: Preparation of Amplicon Libraries.

The accurate assessment of community composition and ultimately species identification is of utmost importance in any ecological and evolutionary study. Advances in sequencing technologies have allowed the unraveling of levels of biodiversity never imagined before when applied to large-scale environmental DNA studies (also termed metabarcoding/metagenetics/metasystematics/environmental barcoding). Here, we describe a detailed protocol to assess eukaryotic biodiversity in marine sediments, identifying key steps that should not be neglected when preparing Next-Generation Sequencing (NGS) amplicon libraries: DNA extraction, multiple PCR amplification of DNA barcode markers with index/ tag-primers, and final Illumina MiSeq sequencing library preparation.

Unravelling the Evolution of the Allatostatin-Type A, KISS and Galanin Peptide-Receptor Gene Families in Bilaterians: Insights from Anopheles Mosquitoes.

UNLABELLED: Allatostatin type A receptors (AST-ARs) are a group of G-protein coupled receptors activated by members of the FGL-amide (AST-A) peptide family that inhibit food intake and development in arthropods. Despite their physiological importance the evolution of the AST-A system is poorly described and relatively few receptors have been isolated and functionally characterised in insects. The present study provides a comprehensive analysis of the origin and comparative evolution of the AST-A system. To determine how evolution and feeding modified the function of AST-AR the duplicate receptors in Anopheles mosquitoes, were characterised. Phylogeny and gene synteny suggested that invertebrate AST-A receptors and peptide genes shared a common evolutionary origin with KISS/GAL receptors and ligands. AST-ARs and KISSR emerged from a common gene ancestor after the divergence of GALRs in the bilaterian genome. In arthropods, the AST-A system evolved through lineage-specific events and the maintenance of two receptors in the flies and mosquitoes (Diptera) was the result of a gene duplication event. Speciation of Anopheles mosquitoes affected receptor gene organisation and characterisation of AST-AR duplicates (GPRALS1 and 2) revealed that in common with other insects, the mosquito receptors were activated by insect AST-A peptides and the iCa2+-signalling pathway was stimulated. GPRALS1 and 2 were expressed mainly in mosquito midgut and ovaries and transcript abundance of both receptors was modified by feeding. A blood meal strongly up-regulated expression of both GPRALS in the midgut (p < 0.05) compared to glucose fed females. Based on the results we hypothesise that the AST-A system in insects shared a common origin with the vertebrate KISS system and may also share a common function as an integrator of metabolism and reproduction. HIGHLIGHTS: AST-A and KISS/GAL receptors and ligands shared common ancestry prior to the protostome-deuterostome divergence. Phylogeny and gene synteny revealed that AST-AR and KISSR emerged after GALR gene divergence. AST-AR genes were present in the hemichordates but were lost from the chordates. In protostomes, AST-ARs persisted and evolved through lineage-specific events and duplicated in the arthropod radiation. Diptera acquired and maintained functionally divergent duplicate AST-AR genes.

PACAP system evolution and its role in melanophore function in teleost fish skin.

Pituitary adenylate cyclase-activating polypeptide (PACAP) administered to tilapia melanophores ex-vivo causes significant pigment aggregation and this is a newly identified function for this peptide in fish. The G-protein coupled receptors (GPCRs), adcyap1r1a (encoding Pac1a) and vipr2a (encoding Vpac2a), are the only receptors in melanophores with appreciable levels of expression and are significantly (p < 0.05) down-regulated in the absence of light. Vpac2a is activated exclusively by peptide histidine isoleucine (PHI), which suggests that Pac1a mediates the melanin aggregating effect of PACAP on melanophores. Paradoxically activation of Pac1a with PACAP caused a rise in cAMP, which in fish melanophores is associated with melanin dispersion. We hypothesise that the duplicate adcyap1ra and vipr2a genes in teleosts have acquired a specific role in skin and that the melanin aggregating effect of PACAP results from the interaction of Pac1a with Ramp that attenuates cAMP-dependent PKA activity and favours the Ca(2+)/Calmodulin dependent pathway.

Environmental metabarcoding reveals heterogeneous drivers of microbial eukaryote diversity in contrasting estuarine ecosystems.

Assessing how natural environmental drivers affect biodiversity underpins our understanding of the relationships between complex biotic and ecological factors in natural ecosystems. Of all ecosystems, anthropogenically important estuaries represent a 'melting pot' of environmental stressors, typified by extreme salinity variations and associated biological complexity. Although existing models attempt to predict macroorganismal diversity over estuarine salinity gradients, attempts to model microbial biodiversity are limited for eukaryotes. Although diatoms commonly feature as bioindicator species, additional microbial eukaryotes represent a huge resource for assessing ecosystem health. Of these, meiofaunal communities may represent the optimal compromise between functional diversity that can be assessed using morphology and phenotype-environment interactions as compared with smaller life fractions. Here, using 454 Roche sequencing of the 18S nSSU barcode we investigate which of the local natural drivers are most strongly associated with microbial metazoan and sampled protist diversity across the full salinity gradient of the estuarine ecosystem. In order to investigate potential variation at the ecosystem scale, we compare two geographically proximate estuaries (Thames and Mersey, UK) with contrasting histories of anthropogenic stress. The data show that although community turnover is likely to be predictable, taxa are likely to respond to different environmental drivers and, in particular, hydrodynamics, salinity range and granulometry, according to varied life-history characteristics. At the ecosystem level, communities exhibited patterns of estuary-specific similarity within different salinity range habitats, highlighting the environmental sequencing biomonitoring potential of meiofauna, dispersal effects or both.

Feeding and the rhodopsin family g-protein coupled receptors in nematodes and arthropods.

In vertebrates, receptors of the rhodopsin G-protein coupled superfamily (GPCRs) play an important role in the regulation of feeding and energy homeostasis and are activated by peptide hormones produced in the brain-gut axis. These peptides regulate appetite and energy expenditure by promoting or inhibiting food intake. Sequence and function homologs of human GPCRs involved in feeding exist in the nematode roundworm, Caenorhabditis elegans (C. elegans), and the arthropod fruit fly, Drosophila melanogaster (D. melanogaster), suggesting that the mechanisms that regulate food intake emerged early and have been conserved during metazoan radiation. Nematodes and arthropods are the most diverse and successful animal phyla on Earth. They can survive in a vast diversity of environments and have acquired distinct life styles and feeding strategies. The aim of the present review is to investigate if this diversity has affected the evolution of invertebrate GPCRs. Homologs of the C. elegans and D. melanogaster rhodopsin receptors were characterized in the genome of other nematodes and arthropods and receptor evolution compared. With the exception of bombesin receptors (BBR) that are absent from nematodes, a similar gene complement was found. In arthropods, rhodopsin GPCR evolution is characterized by species-specific gene duplications and deletions and in nematodes by gene expansions in species with a free-living stage and gene deletions in representatives of obligate parasitic taxa. Based upon variation in GPCR gene number and potentially divergent functions within phyla we hypothesize that life style and feeding diversity practiced by nematodes and arthropods was one factor that contributed to rhodopsin GPCR gene evolution. Understanding how the regulation of food intake has evolved in invertebrates will contribute to the development of novel drugs to control nematodes and arthropods and the pests and diseases that use them as vectors.

Identification of a new cartilage-specific S100-like protein up-regulated during endo/perichondral mineralization in gilthead seabream.

Calcium ions and calcium-binding proteins play a major role in many cellular processes, in particular skeletogenesis and bone formation. We report here the discovery of a novel S100 protein in fish and the analysis of its gene expression patterns. A 648-bp full-length cDNA encoding an 86-amino acid S100-like calcium-binding protein was identified through the subtractive hybridization of a gilthead seabream (Sparus aurata) cDNA library constructed to identify genes associated with in vitro mineralization. Deduced protein lacks an identifiable signal peptide and exhibits two EF-hand motifs characteristic of S100 proteins. Phylogenetic and bioinformatic analyses of S100 sequences suggested that gilthead seabream protein represents a novel and fish-specific member of the S100 protein family. Expression of S100-like gene was up-regulated during the in vitro mineralization of bone-derived cell lines and during seabream development, from larvae throughout adulthood, reflecting skeletogenesis. Restriction of S100-like gene expression to chondrocytes of cartilaginous tissues undergoing endo/perichondral mineralization in juvenile fish further confirmed the mineralogenic role of the protein in fish and emphasized the potential of S100-like as a marker of mineralizing cartilage in developing fish.

Second-generation environmental sequencing unmasks marine metazoan biodiversity.

Biodiversity is of crucial importance for ecosystem functioning, sustainability and resilience, but the magnitude and organization of marine diversity at a range of spatial and taxonomic scales are undefined. In this paper, we use second-generation sequencing to unmask putatively diverse marine metazoan biodiversity in a Scottish temperate benthic ecosystem. We show that remarkable differences in diversity occurred at microgeographical scales and refute currently accepted ecological and taxonomic paradigms of meiofaunal identity, rank abundance and concomitant understanding of trophic dynamics. Richness estimates from the current benchmarked Operational Clustering of Taxonomic Units from Parallel UltraSequencing analyses are broadly aligned with those derived from morphological assessments. However, the slope of taxon rarefaction curves for many phyla remains incomplete, suggesting that the true alpha diversity is likely to exceed current perceptions. The approaches provide a rapid, objective and cost-effective taxonomic framework for exploring links between ecosystem structure and function of all hitherto intractable, but ecologically important, communities.

Identification of an osteopontin-like protein in fish associated with mineral formation.

Fish has been recently recognized as a suitable vertebrate model and represents a promising alternative to mammals for studying mechanisms of tissue mineralization and unravelling specific questions related to vertebrate bone formation. The recently developed Sparus aurata (gilthead seabream) osteoblast-like cell line VSa16 was used to construct a cDNA subtractive library aimed at the identification of genes associated with fish tissue mineralization. Suppression subtractive hybridization, combined with mirror orientation selection, identified 194 cDNA clones representing 20 different genes up-regulated during the mineralization of the VSa16 extracellular matrix. One of these genes accounted for 69% of the total number of clones obtained and was later identified as theS. aurata osteopontin-like gene. The 2138-bp full-length S. aurata osteopontin-like cDNA was shown to encode a 374 amino-acid protein containing domains and motifs characteristic of osteopontins, such as an integrin receptor-binding RGD motif, a negatively charged domain and numerous post-translational modifications (e.g. phosphorylations and glycosylations). The common origin of mammalian osteopontin and fish osteopontin-like proteins was indicated through an in silico analysis of available sequences showing similar gene and protein structures and was further demonstrated by their specific expression in mineralized tissues and cell cultures. Accordingly, and given its proven association with mineral formation and its characteristic protein domains, we propose that the fish osteopontin-like protein may play a role in hard tissue mineralization, in a manner similar to osteopontin in higher vertebrates.

Rapid identification of differentially expressed genes by in situ screening of bacteria.

The identification of differentially expressed genes is a key step in the understanding of specific molecular mechanisms. Various methods have been developed to search for differences in expression but most of them are time or money consuming. We present here an alternative technique that connects standard suppression subtractive hybridization with in situ screening of bacteria to isolate and identify differentially expressed transcripts. The in situ differential screening is based on the transfer of bacteria directly from cultures onto nylon membranes with no need of phenol/chloroform extraction, colony lifting, or polymerase chain reaction amplification. This improved method was successfully applied and must be seen as a simple, low-cost, time-saving, and reproducible approach to identify differentially expressed genes.