High-throughput sequencing of multiple amplicons for barcoding and integrative taxonomy

Environmental factors and host sex influence the skin microbiota structure of Hong Kong newt (Paramesotriton hongkongensis) in a coldspot of chytridiomycosis in subtropical East Asia

Chytridiomycosis, an infectious skin disease caused by the chytrid fungi, Batrachochytrium dendrobatidis and B. salamandrivorans, poses a significant threat to amphibian biodiversity worldwide. Antifungal bacteria found on the skin of chytrid-resistant amphibians could potentially provide defense against chytridiomycosis and lower mortality rates among resistant individuals. The Hong Kong newt (Paramesotriton hongkongensis) is native to East Asia, a region suspected to be the origin of chytrids, and has exhibited asymptomatic infection, suggesting a long-term coexistence with the chytrids. Therefore, the skin microbiota of this resistant species warrant investigation, along with other factors that can affect the microbiota. Among the 149 newts sampled in their natural habitats in Hong Kong, China, putative antifungal bacteria were found in all individuals. There were 314 amplicon sequence variants distributed over 25 genera of putative antifungal bacteria; abundant ones included Acinetobacter, Flavobacterium, and Novosphingobium spp. The skin microbiota compositions were strongly influenced by the inter-site geographical distances. Despite inter-site differences, we identified some core skin microbes across sites that could be vital to P. hongkongensis. The dominant cores included the family Comamonadaceae, family Chitinophagaceae, and class Betaproteobacteria. Moreover, habitat elevation and host sex also exhibited significant effects on skin microbiota compositions. The antifungal bacteria found on these newts offer an important resource for conservation against chytridiomycosis, such as developing probiotic treatments for susceptible species.

Scalable genotyping of microbial colonies

The sequence of the 16S region is taxonomically informative and widely used for genotyping microbes. While it is easy and inexpensive to genotype several isolates by Sanger sequencing the 16S region, this method becomes quite costly if scaled to many isolates. High-throughput sequencing provides one potential avenue for obtaining 16S sequences at scale but presents additional challenges. First, DNA purification workflows for high-throughput sample preparation are labour-intensive and expensive. Second, cost-effective multiplexing and library preparation schemes are difficult to implement for many libraries on a single sequencing run. Therefore, we implemented a scalable protocol for isolate genotyping involving colony polymerase chain reaction (PCR) with simple cell lysis as well as a four-barcode indexing scheme that enables scalable multiplexing and streamlined library preparation by amplifying with four primers simultaneously in a single reaction. We tested this protocol on 93 colonies cultured from environmental samples, and we were able to ascertain the identity of ~90% of microbial isolates.

DNA barcoding, integrative taxonomy, citizen science, and Bush Blitz surveys combine to reveal 34 new species of Apanteles (Hymenoptera, Braconidae, Microgastrinae) in Australia

Microgastrinae is a megadiverse subfamily of wasps in the family Braconidae. As parasitoids of caterpillars, members of the subfamily play important roles in regulating native caterpillar populations, and several species are used commercially as biological control agents. The genus Apanteles comprises a large portion of total microgastrine diversity, however it has not been studied in Australia for more than 30 years, with only nine described species previously known from the continent. We explore the diversity and systematics of Apanteles in Australia, using cytochrome c oxidase subunit I (COI) and Wingless (wg) DNA barcodes from more than 400 Australian Apanteles specimens. Using molecular species delimitation in combination with reduced morphological diagnoses, at least 48 distinct molecular lineages of Apanteles are confirmed in Australia, and 34 new species are formally described, all authored by Slater-Baker, Fagan-Jeffries, Fernández-Triana, Portmann & Oestmann: A.adustus, A.aeternus, A.alatomicans, A.allapsus, A.amicalis, A.apollo, A.apricus, A.artemis, A.aurantius, A.auroralis, A.banrock, A.breviflagellarius, A.brockhedgesi, A.cuprum, A.darthvaderi, A.doreenwatlerae, A.ethanbeaveri, A.fenestrinus, A.ferripulvis, A.focusalis, A.hades, A.insulanus, A.kelpiellus, A.lamingtonensis, A.ligdus, A.magicus, A.margaritarius, A.pellucidus, A.phantasmatus, A.pharusalis, A.ramsaris, A.rufiterra, A.sinusulus, and A.translucentis.

Barcode 100K Specimens: In a Single Nanopore Run

It is a global priority to better manage the biosphere, but action must be informed by comprehensive data on the abundance and distribution of species. The acquisition of such information is currently constrained by high costs. DNA barcoding can speed the registration of unknown animal species, the most diverse kingdom of eukaryotes, as the BIN system automates their recognition. However, inexpensive sequencing protocols are critical as the census of all animal species is likely to require the analysis of a billion or more specimens. Barcoding involves DNA extraction followed by PCR and sequencing with the last step dominating costs until 2017. By enabling the sequencing of highly multiplexed samples, the Sequel platforms from Pacific BioSciences slashed costs by 90%, but these instruments are only deployed in core facilities because of their expense. Sequencers from Oxford Nanopore Technologies provide an escape from high capital and service costs, but their low sequence fidelity has, until recently, constrained adoption. However, the improved performance of its latest flow cells (R10.4.1) erases this barrier. This study demonstrates that a MinION flow cell can characterise an amplicon pool derived from 100,000 specimens while a Flongle flow cell can process one derived from several thousand. At $0.01 per specimen, DNA sequencing is now the least expensive step in the barcode workflow.

Longitudinal survey of total airborne bacterial and archaeal concentrations and bacterial diversity in enriched colony housing and aviaries for laying hens

Conventional cages for laying hens will be banned in Canada as of 2036, and the egg industry is transitioning toward enriched colony housing and aviaries. While higher concentrations of particulate matter have been previously reported in aviaries and other cage-free housing systems, concentrations of total bacteria and archaea suspended in the air are still uncharacterized in Canadian enriched colonies and aviaries. The aim of the present study was to conduct a longitudinal survey of airborne total bacteria and of airborne total archaea in twelve enriched colonies and twelve aviaries in Eastern Canada during a whole laying period. High-throughput sequencing of 16S rRNA gene amplicons was used to reveal and compare bacterial diversity at the start and the end of the production cycle, and during the cold and the warm seasons. Total bacterial and archaeal concentrations were significantly higher in aviaries (p < 0.05) versus enriched colonies, and in the cold season for both housing types (p < 0.05). While flock age did not have a significant effect on total bacterial and archaeal concentrations, it did on bacterial diversity in both enriched colony houses and aviaries (p < 0.05). The 2 housing systems were significantly different in their diversity of bacteria.

Cost-efficient PCR based DNA barcoding of marine invertebrate specimens with NovaSeq amplicon sequencing

BACKGROUND

The marine environment harbors high biodiversity; however, it is poorly understood. Nucleotide sequence data of all marine organisms should be accumulated before natural and/or anthropogenic environmental changes jeopardize the marine environment. In this study, we report a cost-effective and easy DNA barcoding method. This method can be readily adopted without using library preparation kits. It includes multiplex PCR of short targets, indexing PCR, and outsourcing to a sequencing service using the NovaSeq system.

METHODS AND RESULTS

We targeted four mitochondrial genes [cytochrome c oxidase subunit I (COI), COIII, 16S rRNA (16S), and 12S rRNA (12S)] and three nuclear genes [18S rRNA (18S), 28S rRNA (28S), internal transcribed spacer 2 (ITS2)] in 95 marine invertebrate specimens, which were primarily annelids. The primers, including adapters and indices for NovaSeq sequencing, were newly designed. Two PCR runs were conducted. The 1st PCR amplified specific loci with universal primers and the 2nd added sequencing adapters and indices to the 1st PCR products. The gene sequences obtained from the FASTQ files were subjected to BLAST search and phylogenetic analyses. One run using 95 specimens yielded sequences averaging 2816 bp per specimen for a total length of six loci. Nuclear genes were more successfully assembled compared with mitochondrial genes. A weak but significantly negative correlation was observed between the average length of each locus and success rate of the assembly. Some of the sequences were almost identical to the sequences obtained from specimens collected far from Japan, indicating the presence of potentially invasive species identified for the first time.

CONCLUSIONS

We obtained gene sequences efficiently using next-generation sequencing rather than Sanger sequencing. Although this method requires further optimization to increase the success rate for some loci, it is used as a first step to select specimens for further analyses by determining the specific loci of the targets.

Prevalence of opportunistic pathogens and anti-microbial resistance in urban aquaculture ponds

Bacterial antimicrobial resistance (AMR) has emerged as a significant concern worldwide. The microbial community profile and potential AMR level in aquaculture ponds are often undervalued and attract less attention than other aquatic environments. We used amplicon and metagenomic shotgun sequencing to study microbial communities and AMR in six freshwater polyculture ponds in rural and urban areas of Bangladesh. Amplicon sequencing revealed different community structures between rural and urban ponds, with urban ponds having a higher bacterial diversity and opportunistic pathogens including Streptococcus, Staphylococcus, and Corynebacterium. Despite proteobacterial dominance, Firmicutes was the most interactive in the community network, especially in the urban ponds. Metagenomes showed that drug resistance was the most common type of AMR found, while metal resistance was only observed in urban ponds. AMR and metal resistance genes were found mainly in beta and gamma-proteobacteria in urban ponds, while AMR was found primarily in alpha-proteobacteria in rural ponds. We identified potential pathogens with a high profile of AMR and metal resistance in urban aquaculture ponds. As these ponds provide a significant source of protein for humans, our results raise significant concerns for the environmental sustainability of this food source and the dissemination of AMR into the food chain.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Capillary-Electrophoresis-Based Species Barcoding of Big Cats: CR-mtDNA-Length Polymorphism

This study aimed to provide an overview of the methodological approach used for the species determination of big cats. The molecular system described herein employs mitochondrial DNA control region (CR-mtDNA)-length polymorphism in combination with highly sensitive and precise capillary electrophoresis. We demonstrated that the described CR-mtDNA barcoding system can be utilized for species determination where the presence of biological material from big cats is expected or used as a confirmatory test alongside Sanger or massive parallel sequencing (MPS). We have also addressed the fact that species barcoding, when based on the analysis of mtDNA targets, can be biased by nuclear inserts of the mitochondrial genome (NUMTs). The CR-mtDNA barcoding system is suitable even for problematic and challenging samples, such as hair. CR-mtDNA-length polymorphisms can also distinguish hybrids from pure breeds.

Development of an accurate and rapid method for whole genome characterization of canine parvovirus

Canine parvovirus is a highly contagious pathogen affecting domestic dogs and other carnivores globally. Monitoring CPV through continuous genomic surveillance is crucial for mapping variability and developing effective control measures. Here, we developed a method using multiplex-PCR-next-generation sequencing to obtain full-length CPV genomes directly from clinical samples. This approach utilizes tiling and tailed amplicons to amplify overlapping fragments of roughly 250 base pairs. This enables the creation of Illumina libraries by conducting two PCR reaction runs. We tested the assay in 10 fecal samples from dogs diagnosed with CPV and one CPV-2 vaccine strain. Furthermore, we applied it to a feline sample previously diagnosed with the feline panleukopenia virus. The assay provided 100 % genome coverage and high sequencing depth across all 12 samples. It successfully provided the sequence of the coding regions and the left and right non-translated regions, including tandem and terminal repeats. The assay effectively amplified viral variants from divergent evolutionary groups, including the antigenic variants (2a, 2b, and 2c) and the ancestral CPV-2 strain included in vaccine formulations. Moreover, it successfully amplified the entire genome of the feline panleukopenia virus found in cat feces. This method is cost-effective, time-efficient, and does not require lab expertise in Illumina library preparation. The multiplex-PCR-next-generation methodology facilitates large-scale genomic sequencing, expanding the limited number of complete genomes currently available in databases and enabling real-time genomic surveillance. Furthermore, the method helps identify and track emerging CPV viral variants, facilitating molecular epidemiology and control. Adopting this approach can enhance our understanding of the evolution and genetic diversity of Protoparvovirus carnivoran1.

Sex-determining region complements traditionally used in phylogenetic studies nuclear and chloroplast sequences in investigation of Aigeiros Duby and Tacamahaca Spach poplars (genus Populus L., Salicaceae)

Members of the genus Populus L. play an important role in the formation of forests in the northern hemisphere and are used in urban landscaping and timber production. Populus species of closely related sections show extensive hybridization. Therefore, the systematics of the genus is rather complicated, especially for poplars of hybrid origin. We aimed to assess the efficiency of application of the sex-determining region (SDR) in addition to the nuclear and chloroplast genome loci traditionally used in phylogenetic studies of poplars to investigate relationships in sections Aigeiros Duby and Tacamahaca Spach. Targeted deep sequencing of NTS 5S rDNA, ITS, DSH 2, DSH 5, DSH 8, DSH 12, DSH 29, 6, 15, 16, X18, trnG-psbK-psbI, rps2-rpoC2, rpoC2-rpoC1, as well as SDR and ARR17 gene was performed for 379 poplars. The SDR and ARR17 gene together with traditionally used multicopy and single-copy loci of nuclear and chloroplast DNA allowed us to obtain a clustering that is most consistent with poplar systematics based on morphological data and to shed light on several controversial hypotheses about the origin of the studied taxa (for example, the inexpediency of separating P. koreana, P. maximowiczii, and P. suaveolens into different species). We present a scheme of relationships between species and hybrids of sections Aigeiros and Tacamahaca based on molecular genetic, morphological, and geographical data. The geographical proximity of species and, therefore, the possibility of hybridization between them appear to be more important than the affiliation of species to the same section. We speculate that sections Aigeiros and Tacamahaca are distinguished primarily on an ecological principle (plain and mountain poplars) rather than on a genetic basis. Joint analysis of sequencing data for the SDR and chloroplast genome loci allowed us to determine the ancestors of P. × petrovskoe - P. laurifolia (female tree) × P. × canadensis (male tree), and P. × rasumovskoe - P. nigra (female tree) × P. suaveolens (male tree). Thus, the efficiency of using the SDR for the study of poplars of sections Aigeiros and Tacamahaca and the prospects of its use for the investigation of species of the genus Populus were shown.

Competitive interference among rhizobia reduces benefits to hosts

The capacity of beneficial microbes to compete for host infection-and the ability of hosts to discriminate among them-introduces evolutionary conflict that is predicted to destabilize mutualism. We investigated fitness outcomes in associations between legumes and their symbiotic rhizobia to characterize fitness impacts of microbial competition. Diverse Bradyrhizobium strains varying in their capacity to fix nitrogen symbiotically with a common host plant, Acmispon strigosus, were tested in full-factorial coinoculation experiments involving 28 pairwise strain combinations. We analyzed the effects of interstrain competition and host discrimination on symbiotic-interaction outcomes by relativizing fitness proxies to clonally infected and uninfected controls. More than one thousand root nodules of coinoculated plants were genotyped to quantify strain occupancy, and the Bradyrhizobium strain genome sequences were analyzed to uncover the genetic bases of interstrain competition outcomes. Strikingly, interstrain competition favored a fast-growing, minimally beneficial rhizobia strain. Host benefits were significantly diminished in coinoculation treatments relative to expectations from clonally inoculated controls, consistent with competitive interference among rhizobia that reduced both nodulation and plant growth. Competition traits appear polygenic, linked with inter-strain allelopathic interactions in the rhizosphere. This study confirms that competition among strains can destabilize mutualism by favoring microbes that are superior in colonizing host tissues but provide minimal benefits to host plants. Moreover, our findings help resolve the paradox that despite efficient host control post infection, legumes nonetheless encounter rhizobia that vary in their nitrogen fixation.

Deep Mutational Scanning of an Oxygen-Independent Fluorescent Protein CreiLOV for Comprehensive Profiling of Mutational and Epistatic Effects

Oxygen-independent, flavin mononucleotide-based fluorescent proteins (FbFPs) are promising alternatives to green fluorescent protein in anaerobic contexts. Deep mutational scanning performs systematic profiling of protein sequence-function relationships but has not been applied to FbFPs. Focusing on CreiLOV from Chlamydomonas reinhardtii, we created and analyzed two comprehensive mutant collections: (1) single-residue, site-saturation mutagenesis libraries covering all 118 residues; and (2) a full combinatorial metagenesis library among 20 mutations at 15 residues, where mutation and residue selection was based on single-site mutagenesis results. Notably, the second type of library is indispensable to study higher-order epistasis but underrepresented in the literature. Using optimized FACS-seq assays, 2,185 (>92.5%) out of 2,360 possible single-site mutants and 165,428 (>89.7%) out of 184,320 possible combinatorial mutants were reliably assigned with fitness values. We constructed statistical and machine-learning models to analyze the CreiLOV data set, enabling accurate fitness prediction of higher-order mutants using lower-order mutagenesis data. In addition, we successfully isolated CreiLOV variants with improved fluorescence quantum yield and thermostability. This work provides new empirical data and design rules to engineer combinatorial protein variants.

Highly-multiplexed and efficient long-amplicon PacBio and Nanopore sequencing of hundreds of full mitochondrial genomes

BACKGROUND

Mitochondrial genome sequences have become critical to the study of biodiversity. Genome skimming and other short-read based methods are the most common approaches, but they are not well-suited to scale up to multiplexing hundreds of samples. Here, we report on a new approach to sequence hundreds to thousands of complete mitochondrial genomes in parallel using long-amplicon sequencing. We amplified the mitochondrial genome of 677 specimens in two partially overlapping amplicons and implemented an asymmetric PCR-based indexing approach to multiplex 1,159 long amplicons together on a single PacBio SMRT Sequel II cell. We also tested this method on Oxford Nanopore Technologies (ONT) MinION R9.4 to assess if this method could be applied to other long-read technologies. We implemented several optimizations that make this method significantly more efficient than alternative mitochondrial genome sequencing methods.

RESULTS

With the PacBio sequencing data we recovered at least one of the two fragments for 96% of samples (~ 80-90%) with mean coverage ~ 1,500x. The ONT data recovered less than 50% of input fragments likely due to low throughput and the design of the Barcoded Universal Primers which were optimized for PacBio sequencing. We compared a single mitochondrial gene alignment to half and full mitochondrial genomes and found, as expected, increased tree support with longer alignments, though whole mitochondrial genomes were not significantly better than half mitochondrial genomes.

CONCLUSIONS

This method can effectively capture thousands of long amplicons in a single run and be used to build more robust phylogenies quickly and effectively. We provide several recommendations for future users depending on the evolutionary scale of their system. A natural extension of this method is to collect multi-locus datasets consisting of mitochondrial genomes and several long nuclear loci at once.

The impact of nucleosome structure on CRISPR/Cas9 fidelity

The clustered regularly interspaced short palindromic repeats (CRISPR) Cas system is a powerful tool that has the potential to become a therapeutic gene editor in the near future. Cas9 is the best studied CRISPR system and has been shown to have problems that restrict its use in therapeutic applications. Chromatin structure is a known impactor of Cas9 targeting and there is a gap in knowledge on Cas9's efficacy when targeting such locations. To quantify at a single base pair resolution how chromatin inhibits on-target gene editing relative to off-target editing of exposed mismatching targets, we developed the gene editor mismatch nucleosome inhibition assay (GEMiNI-seq). GEMiNI-seq utilizes a library of nucleosome sequences to examine all target locations throughout nucleosomes in a single assay. The results from GEMiNI-seq revealed that the location of the protospacer-adjacent motif (PAM) sequence on the nucleosome edge drives the ability for Cas9 to access its target sequence. In addition, Cas9 had a higher affinity for exposed mismatched targets than on-target sequences within a nucleosome. Overall, our results show how chromatin structure impacts the fidelity of Cas9 to potential targets and highlight how targeting sequences with exposed PAMs could limit off-target gene editing, with such considerations improving Cas9 efficacy and resolving current limitations.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

Vectors as Sentinels: Rising Temperatures Increase the Risk of Xylella fastidiosa Outbreaks

Global change is expected to modify the threat posed by pathogens to plants. However, little is known regarding how a changing climate will influence the epidemiology of generalist vector-borne diseases. We developed a high-throughput screening method to test for the presence of a deadly plant pathogen, Xylella fastidiosa, in its insect vectors. Then, using data from a four-year survey in climatically distinct areas of Corsica (France), we demonstrated a positive correlation between the proportion of vectors positive to X. fastidiosa and temperature. Notably, a higher prevalence corresponded with milder winters. Our projections up to 2100 indicate an increased risk of outbreaks. While the proportion of vectors that carry the pathogen should increase, the climate conditions will remain suitable for the bacterium and its main vector, with possible range shifts towards a higher elevation. Besides calling for research efforts to limit the incidence of plant diseases in the temperate zone, this work reveals that recent molecular technologies could and should be used for massive screening of pathogens in vectors to scale-up surveillance and management efforts.

Concurrent molecular characterization of sand flies and Leishmania parasites by amplicon-based next-generation sequencing

BACKGROUND

Phlebotomine sand flies are vectors of Leishmania parasites, which are the causative agents of leishmaniasis. Herein, we developed an amplicon-based next-generation sequencing (Amp-NGS) to characterize sand flies and Leishmania parasites simultaneously targeting partial fragments of 18S rDNA and ITS1 genes, respectively.

METHODS

Our assay was optimized using reference sand fly (n = 8) and Leishmania spp. (n = 9) samples and validated using wild-caught sand flies from Palestine. The assay was highly specific, and all DNA references were successfully identified to the species level.

RESULTS

Among the wild-caught sand flies (n = 187), Phlebotomus spp. represented 95% of the collected samples (177/187), including Ph. sergenti (147/187, 79%), Ph. papatasi (19/187, 10.2%), Ph. perfiliewi (3/187, 1.6%), Ph. tobbi (2/187, 1.2%) and Ph. syriacus (6/187, 3.2%). Sergentomyia spp. represented only 5% (10/187) of the collected samples and included S. dentata (n = 6), S. fallax (n = 2), S. schwetzi (n = 1) and S. ghesquiere (n = 1). The study observed strong positive correlation between sand fly identification results of the Amp-NGS and morphological identification method (r = 0.84, df = 185, P < 0.001). Some discrepancies between the two methods in the identification of closely related species (i.e. Ph. perfiliewi, Ph. tobbi and Ph. syriacus) were observed. Leishmania DNA was detected and identified as L. tropica in 14 samples (14/187, 7.5%).

CONCLUSIONS

Our assay was sensitive to detect (limit of detection was 0.0016 ng/reaction) and identify Leishmania DNA in sand flies, thus representing a new tool for studying sand flies and their associated Leishmania parasites in endemic areas.

Predicting distributions of Wolbachia strains through host ecological contact—Who's manipulating whom?

Reproductive isolation in response to divergent selection is often mediated via third‐party interactions. Under these conditions, speciation is inextricably linked to ecological context. We present a novel framework for understanding arthropod speciation as mediated by Wolbachia, a microbial endosymbiont capable of causing host cytoplasmic incompatibility (CI). We predict that sympatric host sister‐species harbor paraphyletic Wolbachia strains that provide CI, while well‐defined congeners in ecological contact and recently diverged noninteracting congeners are uninfected due to Wolbachia redundancy. We argue that Wolbachia provides an adaptive advantage when coupled with reduced hybrid fitness, facilitating assortative mating between co‐occurring divergent phenotypes—the contact contingency hypothesis. To test this, we applied a predictive algorithm to empirical pollinating fig wasp data, achieving up to 91.60% accuracy. We further postulate that observed temporal decay of Wolbachia incidence results from adaptive host purging—adaptive decay hypothesis—but implementation failed to predict systematic patterns. We then account for post‐zygotic offspring mortality during CI mating, modeling fitness clines across developmental resources—the fecundity trade‐off hypothesis. This model regularly favored CI despite fecundity losses. We demonstrate that a rules‐based algorithm accurately predicts Wolbachia infection status. This has implications among other systems where closely related sympatric species encounter adaptive disadvantage through hybridization.

A high-throughput amplicon sequencing approach for population-wide species diversity and composition survey

Management of agricultural pests requires an understanding of pest species diversity, their interactions with beneficial insects and spatial-temporal patterns of pest abundance. Invasive and agriculturally important insect pests can build up very high populations, especially in cropping landscapes. Traditionally, sampling effort for species identification involves small sample sizes and is labour intensive. Here, we describe a multi-primer high throughput sequencing (HTS) metabarcoding method and associated analytical workflow for a rapid, intensive, high-volume survey of pest species compositions. We demonstrate our method using the taxonomically challenging Bemisia pest cryptic species complex as examples. The whiteflies Bemisia including the 'tabaci' species are agriculturally important capable of vectoring diverse plant viruses that cause diseases and crop losses. Our multi-primer metabarcoding HTS amplicon approach simultaneously process high volumes of whitefly individuals, with efficiency to detect rare (i.e., 1%) test-species, while our improved whitefly primers for metabarcoding also detected beneficial hymenopteran parasitoid species from whitefly nymphs. Field-testing our redesigned Bemisia metabarcoding primer sets across the Tanzania, Uganda and Malawi cassava cultivation landscapes, we identified the sub-Saharan Africa 1 Bemisia putative species as the dominant pest species, with other cryptic Bemisia species being detected at various abundances. We also provide evidence that Bemisia species compositions can be affected by host crops and sampling techniques that target either nymphs or adults. Our multi-primer HTS metabarcoding method incorporated two over-lapping amplicons of 472bp and 518bp that spanned the entire 657bp 3' barcoding region for Bemisia, and is particularly suitable to molecular diagnostic surveys of this highly cryptic insect pest species complex that also typically exhibited high population densities in heavy crop infestation episodes. Our approach can be adopted to understand species biodiversity across landscapes, with broad implications for improving trans-boundary biosecurity preparedness, thus contributing to molecular ecological knowledge and the development of control strategies for high-density, cryptic, pest-species complexes.

Improved high throughput protocol for targeting eukaryotic symbionts in metazoan and eDNA samples

Eukaryote symbionts of animals are major drivers of ecosystems not only because of their diversity and host interactions from variable pathogenicity but also through different key roles such as commensalism and to different types of interdependence. However, molecular investigations of metazoan eukaryomes require minimising coamplification of homologous host genes. In this study we (1) identified a previously published “antimetazoan” reverse primer to theoretically enable amplification of a wider range of microeukaryotic symbionts, including more evolutionarily divergent sequence types, (2) evaluated in silico several antimetazoan primer combinations, and (3) optimised the application of the best performing primer pair for high throughput sequencing (HTS) by comparing one‐step and two‐step PCR amplification approaches, testing different annealing temperatures and evaluating the taxonomic profiles produced by HTS and data analysis. The primer combination 574*F – UNonMet_DB tested in silico showed the largest diversity of nonmetazoan sequence types in the SILVA database and was also the shortest available primer combination for broadly‐targeting antimetazoan amplification across the 18S rRNA gene V4 region. We demonstrate that the one‐step PCR approach used for library preparation produces significantly lower proportions of metazoan reads, and a more comprehensive coverage of host‐associated microeukaryote reads than the two‐step approach. Using higher PCR annealing temperatures further increased the proportion of nonmetazoan reads in all sample types tested. The resulting V4 region amplicons were taxonomically informative even when only the forward read is analysed. This region also revealed a diversity of known and putatively parasitic lineages and a wider diversity of host‐associated eukaryotes.

A culture-independent approach to understanding the role of soil fungal communities in Bromus tectorum stand failure

Cheatgrass (Bromus tectorum L.) is an invasive annual grass (Poaceae) that has colonized large portions of the Intermountain West. Cheatgrass stand failures have been observed throughout the invaded region, the cause of which may be related to the presence of several species of pathogenic fungi in the soil or surface litter. In this metabarcoding study, we compared the fungal communities between sites that have and have not experienced stand failure. Samples were taken from the soil and surface litter near Winnemucca, Nevada, and in Skull Valley, Utah. Our results show distinct fungal communities associated with stand failure based on both geography and sample type. In both the Winnemucca and Skull Valley surface litter, there was an elevated abundance of the endophyte Ramimonilia apicalis in samples that had experienced a stand failure. Winnemucca surface litter stand failure samples had an increased abundance of a potential pathogen in the genus Comoclathris. Skull Valley surface litter stand failure samples had an increased abundance of an undescribed new species in the Rustroemiaceae family which is responsible for the so‐called bleach blonde syndrome in cheatgrass, while the soils had an increased abundance of potential pathogens in the genera Olpidium and Monosporascus.

Simultaneous genotyping of snails and infecting trematode parasites using high-throughput amplicon sequencing

Several methodological issues currently hamper the study of entire trematode communities within populations of their intermediate snail hosts. Here we develop a new workflow using high-throughput amplicon sequencing to simultaneously genotype snail hosts and their infecting trematode parasites. We designed primers to amplify four snail and five trematode markers in a single multiplex PCR. While also applicable to other genera, we focused on medically and economically important snail genera within the superorder Hygrophila and targeted a broad taxonomic range of parasites within the class Trematoda. We tested the workflow using 417 Biomphalaria glabrata specimens experimentally infected with Schistosoma rodhaini, two strains of Schistosoma mansoni and combinations thereof. We evaluated the reliability of infection diagnostics, the robustness of the workflow, its specificity related to host and parasite identification, and the sensitivity to detect co-infections, immature infections and changes of parasite biomass during the infection process. Finally, we investigated its applicability in wild-caught snails of other genera naturally infected with a diverse range of trematodes. After stringent quality control the workflow allows the identification of snails to species level, and of trematodes to taxonomic levels ranging from family to strain. It is sensitive to detect immature infections and changes in parasite biomass described in previous experimental studies. Co-infections were successfully identified, opening the possibility to examine parasite-parasite interactions such as interspecific competition. Together, these results demonstrate that our workflow provides a powerful tool to analyse the processes shaping trematode communities within natural snail populations.

Microbiome of Barrier Organs in Allergy: Who Runs the World? Germs!

Over the last few decades, allergic diseases have been steadily increasing worldwide, a phenomenon that is not yet completely understood. Recent evidence, however, suggests that alterations in the microbiome may be a contributing factor. The microbiome refers to all microorganisms in a habitat including bacteria, fungi, and viruses. Using modern sequencing technologies, we are now capable of detecting and analyzing the human microbiome in more detail than ever before. Epidemiological and experimental studies have indicated that a complex intestinal microbiome supports the development of the immune system during childhood, thus providing protection from allergic diseases, including food allergy. The microbiome becomes an important part of human physiology and forms dynamic relationships with our various barrier systems. For example, bacterial dysbiosis is a hallmark of atopic eczema and correlates with disease progression. Similarly, the lung and nasopharyngeal microbiome is altered in patients with asthma and allergic rhinitis. While these results are interesting, the underlying mechanisms are still unclear and need to be investigated with functional studies. This review gives a short overview of the terminology and methods used in microbiome research before highlighting results concerning the lung, skin, and intestinal microbiome in allergic diseases.

High-throughput genotyping of high-homology mutant mouse strains by next-generation sequencing

Genotyping of knockout alleles in mice is commonly performed by end-point PCR or gene-specific/universal cassette qPCR. Both have advantages and limitations in terms of assay design and interpretation of results. As an alternative method for high-throughput genotyping, we investigated next generation sequencing (NGS) of PCR amplicons, with a focus on CRISPR-mediated exon deletions where antibiotic selection markers are not present. By multiplexing the wild type and mutant-specific PCR reactions, the genotype can be called by the relative sequence counts of each product. The system is highly scalable and can be applied to a variety of different allele types, including those produced by the International Mouse Phenotyping Consortium and associated projects. One potential challenge with any assay design is locating unique areas of the genome, especially when working with gene families or regions of high homology. These can result in misleading or ambiguous genotypes for either qPCR or end-point assays. Here, we show that genotyping by NGS can negate these issues by simple, automated filtering of undesired sequences. Analysis and genotype calls can also be fully automated, using FASTQ or FASTA input files and an in-house Perl script and SQL database.

Assessing DNA for fish identifications from reference collections: the good, bad and ugly shed light on formalin fixation and sequencing approaches

Natural history collections are repositories of biodiversity and are potentially used by molecular ecologists for comparative taxonomic, phylogenetic, biogeographic and forensic purposes. Specimens in fish collections are preserved using a combination of methods with many fixed in formalin and then preserved in ethanol for long-term storage. Formalin fixation damages DNA, thereby limiting genetic analyses. In this study, the authors compared the DNA barcoding and identification success for frozen and formalin-fixed tissues obtained from specimens in the CSIRO Australian National Fish Collection. They studied 230 samples from fishes (consisting of >160 fish species). An optimized formalin-fixed, paraffin-embedded DNA extraction method resulted in usable DNA from degraded tissues. Four mini barcoding assays of the mitochondrial DNA (mtDNA) were characterized with Sanger and Illumina amplicon sequencing. In the good quality DNA (without exposure to formalin), up to 88% of the specimens were correctly matched at the species level using the cytochrome oxidase subunit 1 (COI) mini barcodes, whereas up to 58% of the specimens exposed to formalin for less than 8 weeks were correctly identified to species. In contrast, 16S primers provided higher amplification success with formalin-exposed tissues, although the COI gene was more successful for identification. Importantly, the authors found that DNA of a certain size and quality can be amplified and sequenced despite exposure to formalin, and Illumina sequencing provided them with greater power of resolution for taxa identification even when there was little DNA present. Overall, within parameter constraints, this study highlights the possibilities of recovering DNA barcodes for identification from formalin-fixed fish specimens, and the authors provide guidelines for when successful identification could be expected.

Genomic Profiling: The Strengths and Limitations of Chloroplast Genome-Based Plant Variety Authentication

Genomic profiling is a suitable tool for variety authentication and has applications in both operational quality and regulatory raw material control. It can be used to differentiate species or varieties and to identify admixtures as well as field contaminants. To establish a molecular profile, reliable and very accurate sequence data are required. As a result of the influence of the pollinator plant, nuclear genome-based authentication is in most cases not suitable for a direct application on the fruit. Sequences must be used that come exclusively from the localized mother plant. Parts of the fruit of maternal origin, e.g., components derived from the blossom, are suitable as a basis for this. Alternatively, DNA from cell organelles that are maternally inherited, such as mitochondria or chloroplasts, can be used. The latter will be discussed in this review in closer detail. Although individual gene segments on the chloroplast genome are already used for species differentiation in barcoding studies on plants, little is known about the usefulness of the entire chloroplast genome for intraspecies differentiation in general and for differentiation between modern varieties in particular. Results from the literature as well as from our own work suggest that chloroplast genome sequences are indeed very well-suited for the differentiation of old varieties. On the other hand, they are less or not suitable for the genetic differentiation of modern cultivars, because they are often too closely related.

Efficient COI barcoding using high throughput single-end 400 bp sequencing

BACKGROUND

Over the last decade, the rapid development of high-throughput sequencing platforms has accelerated species description and assisted morphological classification through DNA barcoding. However, the current high-throughput DNA barcoding methods cannot obtain full-length barcode sequences due to read length limitations (e.g. a maximum read length of 300 bp for the Illumina's MiSeq system), or are hindered by a relatively high cost or low sequencing output (e.g. a maximum number of eight million reads per cell for the PacBio's SEQUEL II system).

RESULTS

Pooled cytochrome c oxidase subunit I (COI) barcodes from individual specimens were sequenced on the MGISEQ-2000 platform using the single-end 400 bp (SE400) module. We present a bioinformatic pipeline, HIFI-SE, that takes reads generated from the 5' and 3' ends of the COI barcode region and assembles them into full-length barcodes. HIFI-SE is written in Python and includes four function modules of filter, assign, assembly and taxonomy. We applied the HIFI-SE to a set of 845 samples (30 marine invertebrates, 815 insects) and delivered a total of 747 fully assembled COI barcodes as well as 70 Wolbachia and fungi symbionts. Compared to their corresponding Sanger sequences (72 sequences available), nearly all samples (71/72) were correctly and accurately assembled, including 46 samples that had a similarity score of 100% and 25 of ca. 99%.

CONCLUSIONS

The HIFI-SE pipeline represents an efficient way to produce standard full-length barcodes, while the reasonable cost and high sensitivity of our method can contribute considerably more DNA barcodes under the same budget. Our method thereby advances DNA-based species identification from diverse ecosystems and increases the number of relevant applications.

Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology

BACKGROUND

Amplicon sequencing of phylogenetic marker genes, e.g., 16S, 18S, or ITS ribosomal RNA sequences, is still the most commonly used method to determine the composition of microbial communities. Microbial ecologists often have expert knowledge on their biological question and data analysis in general, and most research institutes have computational infrastructures to use the bioinformatics command line tools and workflows for amplicon sequencing analysis, but requirements of bioinformatics skills often limit the efficient and up-to-date use of computational resources.

RESULTS

We present dadasnake, a user-friendly, 1-command Snakemake pipeline that wraps the preprocessing of sequencing reads and the delineation of exact sequence variants by using the favorably benchmarked and widely used DADA2 algorithm with a taxonomic classification and the post-processing of the resultant tables, including hand-off in standard formats. The suitability of the provided default configurations is demonstrated using mock community data from bacteria and archaea, as well as fungi.

CONCLUSIONS

By use of Snakemake, dadasnake makes efficient use of high-performance computing infrastructures. Easy user configuration guarantees flexibility of all steps, including the processing of data from multiple sequencing platforms. It is easy to install dadasnake via conda environments. dadasnake is available at https://github.com/a-h-b/dadasnake.

16S rRNA Gene Diversity of Bacterial Endophytes in Parasitic Cuscuta campestris and Its Helianthus annuus Host

Here, we report the results of 16S rRNA gene amplicon sequencing of bacterial endophytes from parasitized and unparasitized samples of the common sunflower (Helianthus annuus) and samples of its associated plant parasite field dodder (Cuscuta campestris), collected from one location in Fresno County, California (August 2017).

Here, we report the results of 16S rRNA gene amplicon sequencing of bacterial endophytes from parasitized and unparasitized samples of the common sunflower (Helianthus annuus) and samples of its associated plant parasite field dodder (Cuscuta campestris), collected from one location in Fresno County, California (August 2017).

Rapid Changes in Microbial Community Structures along a Meandering River

Streams and rivers convey freshwater from lands to the oceans, transporting various organic particles, minerals, and living organisms. Microbial communities are key components of freshwater food webs and take up, utilize, and transform this material. However, there are still important gaps in our understanding of the dynamic of these organisms along the river channels. Using high-throughput 16S and 18S rRNA gene sequencing and quantitative PCR on a 11-km long transect of the Saint-Charles River (Quebec, CA), starting from its main source, the Saint-Charles Lake, we show that bacterial and protist community structures in the river drifted quickly but progressively downstream of its source. The dominant Operational Taxonomic Units (OTUs) of the lake, notably related to Cyanobacteria, decreased in proportions, whereas relative proportions of other OTUs, such as a Pseudarcicella OTU, increased along the river course, becoming quickly predominant in the river system. Both prokaryotic and protist communities changed along the river transect, suggesting a strong impact of the shift from a stratified lake ecosystem to a continuously mixed river environment. This might reflect the cumulative effects of the increasing water turbulence, fluctuations of physicochemical conditions, differential predation pressure in the river, especially in the lake outlet by benthic filter feeders, or the relocation of microorganisms, through flocculation, sedimentation, resuspension, or inoculation from the watershed. Our study reveals that the transit of water in a river system can greatly impact both bacterial and micro-eukaryotic community composition, even over a short distance, and, potentially, the transformation of materials in the water column.

Takeaways from Mobile DNA Barcoding with BentoLab and MinION

Since the release of the MinION sequencer in 2014, it has been applied to great effect in the remotest and harshest of environments, and even in space. One of the most common applications of MinION is for nanopore-based DNA barcoding in situ for species identification and discovery, yet the existing sample capability is limited (n ≤ 10). Here, we assembled a portable sequencing setup comprising the BentoLab and MinION and developed a workflow capable of processing 32 samples simultaneously. We demonstrated this enhanced capability out at sea, where we collected samples and barcoded them onboard a dive vessel moored off Sisters' Islands Marine Park, Singapore. In under 9 h, we generated 105 MinION barcodes, of which 19 belonged to fresh metazoans processed immediately after collection. Our setup is thus viable and would greatly fortify existing portable DNA barcoding capabilities. We also tested the performance of the newly released R10.3 nanopore flow cell for DNA barcoding, and showed that the barcodes generated were ~99.9% accurate when compared to Illumina references. A total of 80% of the R10.3 nanopore barcodes also had zero base ambiguities, compared to 50-60% for R9.4.1, suggesting an improved homopolymer resolution and making the use of R10.3 highly recommended.

Elucidating Biodiversity Shifts in Ballast Water Tanks during a Cross-Latitudinal Transfer: Complementary Insights from Molecular Analyses

In this study, the evolution of ballast water (BW) assemblages across different trophic levels was characterized over a 21 day cross-latitudinal vessel transit using a combination of molecular methods. Triplicate BW samples were collected every second day and size-fractionated (<2.7, 10, >50 μm). Measurements of adenosine triphosphate (ATP) and metabarcoding of environmental nucleic acid (DNA and RNA) analyses, complemented by microscopy and flow cytometry, were performed on each sample. Measured ATP concentrations exhibited high variance between replicates and a strong negative trend in the large (≥50 μm) fraction over the voyage. In concert with microscopy, the metabarcoding data indicated a die-off of larger metazoans during the first week of study and gradual reductions in dinoflagellates and ochrophytes. The ATP and metabarcoding data signaled persistent or increased cellular activity of heterotrophic bacteria and protists in the BW, which was supported by flow cytometry. The metabarcoding showed the presence of active bacteria in all size fractions, suggesting that the sequential filtration approach does not ensure taxonomical differentiation, which has implications for BW quality assessment. Although our data show that ATP and metabarcoding have potential for indicative BW screening for BW compliance monitoring, further research and technological development is needed to improve representativeness of sampling and deliver the unequivocal response criteria required by the international Ballast Water Management Convention.

IonBreeders: bioinformatics plugins toward genomics-assisted breeding

Polymorphism information generated by next-generation sequencing (NGS) technologies has enabled applications of genome-wide markers assisted breeding. However, handling such large-scale data remains a challenge for experimental researchers and breeders, calling for the urgent development of a flexible and straightforward analysis tool for NGS data. We developed "IonBreeders" as bioinformatics plugins that implement general analysis steps from genotyping to genomic prediction. IonBreeders comprises three plugins, "ABH", "IMPUTATION", and "GENOMIC PREDICTION", for format conversion of genotyping data, preprocessing and imputation of genotyping data, and genomic prediction, respectively. "ABH" converts genotyping data derived from NGS into the ABH format, which is acceptable for our further plugins and with other breeding software tools, R/qtl, MapMaker, and AntMap. "IMPUTATION" filters out non-informative markers and imputes missing marker genotypes. In "GENOMIC PREDICTION", users can use four statistical methods based on their target trait, quantitative trait locus effect, and number of markers, and construct a prediction model for genomic selection. IonBreeders is operated in Torrent Suite, but can also handle genotype data in standard formats, e.g., Variant Call Format (VCF), by format conversion using free software or our provided scripts.

Microfluidic Enrichment Barcoding (MEBarcoding): a new method for high throughput plant DNA barcoding

DNA barcoding is a valuable tool to support species identification with broad applications from traditional taxonomy, ecology, forensics, food analysis, and environmental science. We introduce Microfluidic Enrichment Barcoding (MEBarcoding) for plant DNA Barcoding, a cost-effective method for high-throughput DNA barcoding. MEBarcoding uses the Fluidigm Access Array to simultaneously amplify targeted regions for 48 DNA samples and hundreds of PCR primer pairs (producing up to 23,040 PCR products) during a single thermal cycling protocol. As a proof of concept, we developed a microfluidic PCR workflow using the Fluidigm Access Array and Illumina MiSeq. We tested 96 samples for each of the four primary DNA barcode loci in plants: rbcL, matK, trnH-psbA, and ITS. This workflow was used to build a reference library for 78 families and 96 genera from all major plant lineages - many currently lacking in public databases. Our results show that this technique is an efficient alternative to traditional PCR and Sanger sequencing to generate large amounts of plant DNA barcodes and build more comprehensive barcode databases.

A method to generate multilocus barcodes of pinned insect specimens using MiSeq

For molecular insect identification, amplicon sequencing methods are recommended because they offer a cost-effective approach for targeting small sets of informative genes from multiple samples. In this context, high-throughput multilocus amplicon sequencing has been achieved using the MiSeq Illumina sequencing platform. However, this approach generates short gene fragments of <500 bp, which then have to be overlapped using bioinformatics to achieve longer sequence lengths. This increases the risk of generating chimeric sequences or leads to the formation of incomplete loci. Here, we propose a modified nested amplicon sequencing method for targeting multiple loci from pinned insect specimens using the MiSeq Illumina platform. The modification exists in using a three-step nested PCR approach targeting near full-length loci in the initial PCR and subsequently amplifying short fragments of between 300 and 350 bp for high-throughput sequencing using Illumina chemistry. Using this method, we generated 407 sequences of three loci from 86% of all the specimens sequenced. Out of 103 pinned bee specimens of replicated species, 71% passed the 95% sequence similarity threshold between species replicates. This method worked best for pinned specimens aged between 0 and 5 years, with a limit of 10 years for pinned and 14 years for ethanol-preserved specimens. Hence, our method overcomes some of the challenges of amplicon sequencing using short read next generation sequencing and improves the possibility of creating high-quality multilocus barcodes from insect collections.

A reference library for Canadian invertebrates with 1.5 million barcodes, voucher specimens, and DNA samples

The reliable taxonomic identification of organisms through DNA sequence data requires a well parameterized library of curated reference sequences. However, it is estimated that just 15% of described animal species are represented in public sequence repositories. To begin to address this deficiency, we provide DNA barcodes for 1,500,003 animal specimens collected from 23 terrestrial and aquatic ecozones at sites across Canada, a nation that comprises 7% of the planet's land surface. In total, 14 phyla, 43 classes, 163 orders, 1123 families, 6186 genera, and 64,264 Barcode Index Numbers (BINs; a proxy for species) are represented. Species-level taxonomy was available for 38% of the specimens, but higher proportions were assigned to a genus (69.5%) and a family (99.9%). Voucher specimens and DNA extracts are archived at the Centre for Biodiversity Genomics where they are available for further research. The corresponding sequence and taxonomic data can be accessed through the Barcode of Life Data System, GenBank, the Global Biodiversity Information Facility, and the Global Genome Biodiversity Network Data Portal.

Cross-contamination and strong mitonuclear discordance in Empria sawflies (Hymenoptera, Tenthredinidae) in the light of phylogenomic data

In several sawfly taxa strong mitonuclear discordance has been observed, with nuclear genes supporting species assignments based on morphology, whereas the barcode region of the mitochondrial COI gene suggests different relationships. As previous studies were based on only a few nuclear genes, the causes and the degree of mitonuclear discordance remain ambiguous. Here, we obtained genomic-scale ddRAD data together with Sanger sequences of mitochondrial COI and two to three nuclear protein coding genes to investigate species limits and mitonuclear discordance in two closely related species groups of the sawfly genus Empria. As found previously based on nuclear ITS and mitochondrial COI sequences, species are in most cases supported as monophyletic based on new nuclear data reported here, but not based on mitochondrial COI. This mitonuclear discordance can be explained by occasional mitochondrial introgression with little or no nuclear gene flow, a pattern that might be common in haplodiploid taxa with slowly evolving mitochondrial genomes. Some species in the E. immersa group are not recovered as monophyletic according to either mitochondrial or nuclear data, but this could partly be because of unresolved taxonomy. Preliminary analyses of ddRAD data did not recover monophyly of E. japonica within the E. longicornis group (three Sanger sequenced nuclear genes strongly supported monophyly), but closer examination of the data and additional Sanger sequencing suggested that both specimens were substantially (possibly 10-20% of recovered loci) cross-contaminated. A reason could be specimen identification tag jumps during sequencing library preparation that in previous studies have been shown to affect up to 2.5% of the sequenced reads. We provide an R script to examine patterns of identical loci among the specimens and estimate that the cross-contamination rate is not unusually high for our ddRAD dataset as a whole (based on counting of identical sequences in the immersa and longicornis groups, which are well separated from each other and probably do not hybridise). The high rate of cross-contamination for both E. japonica specimens might be explained by the small number of recovered loci (~1000) compared to most other specimens (>10 000 in some cases) because of poor sequencing results. We caution against drawing unexpected biological conclusions when closely related specimens are pooled before sequencing and tagged only at one end of the molecule or at both ends using a unique combination of limited number of tags (less than the number of specimens).

Annual Protist Community Dynamics in a Freshwater Ecosystem Undergoing Contrasted Climatic Conditions: The Saint-Charles River (Canada)

Protists are key stone components of aquatic ecosystems, sustaining primary productivity and aquatic food webs. However, their diversity, ecology and structuring factors shaping their temporal distribution remain strongly misunderstood in freshwaters. Using high-throughput sequencing on water samples collected over 16 different months (including two summer and two winter periods), combined with geochemical measurements and climate monitoring, we comprehensively determined the pico- and nanoeukaryotic community composition and dynamics in a Canadian river undergoing prolonged ice-cover winters. Our analysis revealed a large protist diversity in this fluctuating ecosystem and clear seasonal patterns demonstrating a direct and/or indirect selective role of abiotic factors, such as water temperature or nitrogen concentrations, in structuring the eukaryotic microbial community. Nonetheless, our results also revealed that primary productivity, predatory as well as parasitism lifestyles, inferred from fine phylogenetic placements, remained potentially present over the annual cycle, despite the large seasonal fluctuations and the remodeling of the community composition under ice. In addition, potential interplays with the bacterial community composition were identified supporting a possible contribution of the bacterial community to the temporal dynamics of the protist community structure. Our results illustrate the complexity of the eukaryotic microbial community and provide a substantive and useful dataset to better understand the global freshwater ecosystem functioning.

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.

DNA metabarcoding of Neotropical ichthyoplankton: Enabling high accuracy with lower cost

Knowledge of ichthyoplankton dynamics is extremely important for conservation management as it can provide information about preferential spawning sites, reproductive period, migratory routes and recruitment success, which can be used to guide management and conservation efforts. However, identification of the eggs and larvae of Neotropical freshwater fish is a difficult task. DNA barcodes have emerged as an alternative and highly accurate approach for species identification, but DNA barcoding can be time-consuming and costly. To solve this problem, we aimed to develop a simple protocol based on DNA metabarcoding, to investigate whether it is possible to detect and quantify all species present in a pool of organisms. To do this, 230 larvae were cut in half, one half was sequenced by the Sanger technique and the other half was used to compose six arrays with a pool of larvae that were sequenced using a next-generation technique (NGS). The results of the Sanger sequencing allowed the identification of almost all larvae at species level, and the results from NGS showed high accuracy in species detection, ranging from 83% to 100%, with an average of 95% in all samples. No false positives were detected. The frequency of organisms in the two methods was positively correlated (Pearson), with low variation among species. In conclusion, this protocol represents a considerable advance in ichthyoplankton studies, allowing a rapid, cost-effective, quali-quantitative approach that improves the accuracy of identification.

A DNA barcode library for 5,200 German flies and midges (Insecta: Diptera) and its implications for metabarcoding-based biomonitoring

This study summarizes results of a DNA barcoding campaign on German Diptera, involving analysis of 45,040 specimens. The resultant DNA barcode library includes records for 2,453 named species comprising a total of 5,200 barcode index numbers (BINs), including 2,700 COI haplotype clusters without species‐level assignment, so called “dark taxa.” Overall, 88 out of 117 families (75%) recorded from Germany were covered, representing more than 50% of the 9,544 known species of German Diptera. Until now, most of these families, especially the most diverse, have been taxonomically inaccessible. By contrast, within a few years this study provided an intermediate taxonomic system for half of the German Dipteran fauna, which will provide a useful foundation for subsequent detailed, integrative taxonomic studies. Using DNA extracts derived from bulk collections made by Malaise traps, we further demonstrate that species delineation using BINs and operational taxonomic units (OTUs) constitutes an effective method for biodiversity studies using DNA metabarcoding. As the reference libraries continue to grow, and gaps in the species catalogue are filled, BIN lists assembled by metabarcoding will provide greater taxonomic resolution. The present study has three main goals: (a) to provide a DNA barcode library for 5,200 BINs of Diptera; (b) to demonstrate, based on the example of bulk extractions from a Malaise trap experiment, that DNA barcode clusters, labelled with globally unique identifiers (such as OTUs and/or BINs), provide a pragmatic, accurate solution to the “taxonomic impediment”; and (c) to demonstrate that interim names based on BINs and OTUs obtained through metabarcoding provide an effective method for studies on species‐rich groups that are usually neglected in biodiversity research projects because of their unresolved taxonomy.

CRIS.py: A Versatile and High-throughput Analysis Program for CRISPR-based Genome Editing

CRISPR-Cas9 technology allows the creation of user-defined genomic modifications in cells and whole organisms. However, quantifying editing rates in pools of cells or identifying correctly edited clones is tedious. Targeted next-generation sequencing provides a high-throughput platform for optimizing editing reagents and identifying correctly modified clones, but the large amount of data produced can be difficult to analyze. Here, we present CRIS.py, a simple and highly versatile python-based program which concurrently analyzes next-generation sequencing data for both knock-out and multiple user-specified knock-in modifications from one or many edited samples. Compared to available NGS analysis programs for CRISPR based-editing, CRIS.py has many advantages: (1) the ability to analyze from one to thousands of samples at once, (2) the capacity to check each sample for multiple sequence modifications, including those induced by base-editors, (3) an output in an easily searchable file format enabling users to quickly sort through and identify correctly targeted clones.

Application of a SSR-GBS marker system on investigation of European Hedgehog species and their hybrid zone dynamics

By applying second-generation sequencing technologies to microsatellite genotyping, sequence information is produced which can result in high-resolution population genetics analysis populations and increased replicability between runs and laboratories. In the present study, we establish an approach to study the genetic structure patterns of two European hedgehog species Erinaceaus europaeus and E. roumanicus. These species are usually associated with human settlements and are good models to study anthropogenic impacts on the genetic diversity of wild populations. The short sequence repeats genotyping by sequence (SSR-GBS) method presented uses amplicon sequences to determine genotypes for which allelic variants can be defined according to both length and single nucleotide polymorphisms (SNPs). To evaluate whether complete sequence information improved genetic structure definition, we compared this information with datasets based solely on length information. We identified a total of 42 markers which were successfully amplified in both species. Overall, genotyping based on complete sequence information resulted in a higher number of alleles, as well as greater genetic diversity and differentiation between species. Additionally, the structure patterns were slightly clearer with a division between both species and some potential hybrids. There was some degree of genetic structure within species, although only in E. roumanicus was this related to geographical distance. The statistically significant results obtained by SSR-GBS demonstrate that it is superior to electrophoresis-based methods for SSR genotyping. Moreover, the greater reproducibility and throughput with lower effort which can be obtained with SSR-GBS and the possibility to include degraded DNA into the analysis, allow for continued relevance of SSR markers during the genomic era.

RNA-sequencing of a citrus bud-feeder, Podagricomela weisei (Coleoptera: Chrysomelidae), reveals xenobiotic metabolism/core RNAi machinery-associated genes and conserved miRNAs

The citrus leaf-mining beetle, Podagricomela weisei Heikertinger, is an important citrus pest that ingests the mesophyll and new shoots. The mechanism underlying the xenobiotic metabolism of P. weisei is not well understood, in part because of a lack of available genomic and transcriptomic data, which has hampered the development of novel pest management approaches [e.g., RNA interference (RNAi)]. In this study, we completed the deep sequencing of the P. weisei transcriptome to identify factors potentially involved in xenobiotic metabolism and the core RNAi machinery. The sequencing of the P. weisei transcriptome generated >27 million clean reads, ultimately yielding 90,410 unigenes with an N50 of 1065 bp. The unigenes were used as queries to search the Nr database, which revealed that 21,847 unigenes were homologous to known genes in various species. Transcripts encoding genes involved in xenobiotic metabolism were identified, including genes encoding cytochrome P450 monooxygenase (P450, 47 unigenes), glutathione S-transferase (GST, 12 unigenes), esterase (EST, 25 unigenes), and the ATP-binding cassette transporter (ABC transporter, 32 unigenes). A parallel sequencing of small RNAs detected 30 conserved miRNAs, with the most abundant being Pwe-miR-1-3p, with an expression level reaching 517,996 reads in the prepared library, followed by Pwe-miR-8-3p (149,402 reads). Genes encoding components of the miRNA, siRNA, and piRNA pathways were also identified, and the results indicated that P. weisei possesses only one of each gene in all three pathways. In summary, this is the first detailed analysis of the transcriptome and small RNAs of P. weisei. The datasets presented herein may form the basis for future molecular characterizations of P. weisei as well as the development of enhanced pest control strategies.

Illumina midi-barcodes: quality proof and applications

DNA barcoding constitutes a supplemental genetically based characterization tool for the identification of species. Traditionally, the barcodes are generated with a length of 650 bp using standardized Sanger sequencing, but with the introduction of high-throughput sequencing (HTS) methods new opportunities for sequencing are available. To use HTS for barcode collection and identification, the amplification of shorter fragments is preferred. Reference DNA midi-barcodes of wild bees were produced using the Illumina MiSeq as well as the Sanger method. Although DNA midi-barcodes derived from Illumina were comparatively shorter (418 bp), their sequences were coherent to the morphological assignment of species. The Illumina barcodes proved to be effective and dealt better with some general limitations of DNA barcoding.