Utilizing field collected insects for next generation sequencing: Effects of sampling, storage, and DNA extraction methods

Effect of temperature and humidity on insect DNA integrity evaluated by real-time PCR

Insects collected in dry traps can degrade rapidly, especially in warm, humid environments where many biodiversity and biosecurity surveillance activities are undertaken. Degradation can severely impact diagnostics, as trap catches can become difficult to identify to species level using morphological characters or, of increasing importance, molecular approaches. This is especially problematic for biosecurity surveillance of exotic tephritid fruit flies, where diagnostics are heavily reliant on morphological characters. We tested the effects of differing temperature and humidity conditions on mock samples of tephritid fruit flies in a controlled environment and compared our results to field trap catches. DNA degradation was quantified using real-time PCR assays, including one assay newly developed and tested here. We observed a correlation between increasing DNA degradation and increasing temperature and humidity. The greatest DNA degradation occurred under combined high humidity (90% relative humidity) and constant high temperature (35 °C). Unexpectedly, fluctuating temperature did not have a significant impact on DNA. Other factors, such as trap design, time in the field, and rainfall, did not significantly correlate with DNA quality across the field samples tested. When plotted against mock samples, field samples clustered together, with no clear pattern or predictability regarding the quantity of DNA preserved, indicating other untested environmental variables may be at play. Predictably, increased exposure time was found to have a detrimental effect on DNA quality for all treatments. These findings will improve the delivery of surveillance activities through the implementation of shorter trap clearance timeframes and improved trap designs and procedures.

Vegetation structure and climate shape mountain arthropod distributions across trophic levels

Arthropods play a vital role in ecosystems; yet, their distributions remain poorly understood, particularly in mountainous regions. This study delves into the modelling of the distribution of 31 foliar arthropod genera in the French Alps, using a comprehensive approach encompassing multi-trophic sampling, community DNA metabarcoding and random forest models. The results underscore the significant importance of vegetation structure, such as herbaceous vegetation density, and forest density and heterogeneity, along with climate, in shaping the distributions of most arthropods. These responses to environmental gradients are consistent across trophic groups, with the exception of nectarivores, whose distributions are more sensitive to landscape structure and water availability. By leveraging community DNA metabarcoding, this study sheds light on the understudied drivers of arthropod distributions, emphasizing the importance of modelling across diverse trophic groups to anticipate arthropod responses to global change.

Improving the assessment of ecosystem and wildlife health: microbiome as an early indicator

Human activities are causing dramatic declines in ecosystem health, compromising the functioning of the life-support system, economic activity, and animal and human health. In this context, monitoring the health of ecosystems and wildlife populations is crucial for determining ecological dynamics and assessing management interventions. A growing body of evidence indicates that microbiome provides a meaningful early indicator of ecosystem and wildlife health. Microbiome is ubiquitous and both environmental and host-associated microbiomes rapidly reflect anthropogenic disturbances. However, we still need to overcome current limitations such as nucleic acid degradation, sequencing depth, and the establishment of baseline data to maximize the potential of microbiome studies.

Future of DNA-based insect monitoring

Insects are crucial for ecosystem health but climate change and pesticide use are driving massive insect decline. To mitigate this loss, we need new and effective monitoring techniques. Over the past decade there has been a shift to DNA-based techniques. We describe key emerging techniques for sample collection. We suggest that the selection of tools should be broadened, and that DNA-based insect monitoring data need to be integrated more rapidly into policymaking. We argue that there are four key areas for advancement, including the generation of more complete DNA barcode databases to interpret molecular data, standardisation of molecular methods, scaling up of monitoring efforts, and integrating molecular tools with other technologies that allow continuous, passive monitoring based on images and/or laser imaging, detection, and ranging (LIDAR).

A DNA Extraction Method for Insects From Sticky Traps: Targeting a Low Abundance Pest, Phthorimaea absoluta (Lepidoptera: Gelechiidae), in Mixed Species Communities

Invasive insects can cause catastrophic damage to ecosystems and cost billions of dollars each year due to management expenses and lost revenue. Rapid detection is an important step to prevent invasive insects from spreading, but improvements in detection capabilities are needed for bulk collections like those from sticky traps. Here we present a bulk DNA extraction method designed for the detection of Phthorimaea absoluta Meyrick (Lepidoptera: Gelechiidae), an invasive moth that can decimate tomato crops. We test the extraction method for insect specimens on sticky traps, subjected to different temperature and humidity conditions, and among mock insect communities left in the field for up to 21 d. We find that the extraction method yielded high success (>92%) in recovering target DNA across field and lab trials, without a decline in recovery after three weeks, across all treatments. These results may have a large impact on tomato growing regions where P. absoluta is in the early stages of invasion or not yet present. The extraction method can also be used to improve detection capabilities for other bulk insect collections, especially those using sticky traps, to the benefit of pest surveys and biodiversity studies.

A DNA barcode survey of insect biodiversity in Pakistan

Although Pakistan has rich biodiversity, many groups are poorly known, particularly insects. To address this gap, we employed DNA barcoding to survey its insect diversity. Specimens obtained through diverse collecting methods at 1,858 sites across Pakistan from 2010-2019 were examined for sequence variation in the 658 bp barcode region of the cytochrome c oxidase 1 (COI) gene. Sequences from nearly 49,000 specimens were assigned to 6,590 Barcode Index Numbers (BINs), a proxy for species, and most (88%) also possessed a representative image on the Barcode of Life Data System (BOLD). By coupling morphological inspections with barcode matches on BOLD, every BIN was assigned to an order (19) and most (99.8%) were placed to a family (362). However, just 40% of the BINs were assigned to a genus (1,375) and 21% to a species (1,364). Five orders (Coleoptera, Diptera, Hemiptera, Hymenoptera, Lepidoptera) accounted for 92% of the specimens and BINs. More than half of the BINs (59%) are so far only known from Pakistan, but others have also been reported from Bangladesh (13%), India (12%), and China (8%). Representing the first DNA barcode survey of the insect fauna in any South Asian country, this study provides the foundation for a complete inventory of the insect fauna in Pakistan while also contributing to the global DNA barcode reference library.

Unamplified, Long-Read Metagenomic Sequencing Approach to Close Endosymbiont Genomes of Low-Biomass Insect Populations

With the current advancements in DNA sequencing technology, the limiting factor in long-read metagenomic assemblies is now the quantity and quality of input DNA. Although these requirements can be met through the use of axenic bacterial cultures or large amounts of biological material, insect systems that contain unculturable bacteria or that contain a low amount of available DNA cannot fully utilize the benefits of third-generation sequencing. The citrus greening disease insect vector Diaphorina citri is an example that exhibits both of these limitations. Although endosymbiont genomes have mostly been closed after the short-read sequencing of amplified template DNA, creating de novo long-read genomes from the unamplified DNA of an insect population may benefit communities using bioinformatics to study insect pathosystems. Here all four genomes of the infected D. citri microbiome were sequenced to closure using unamplified template DNA and two long-read sequencing technologies. Avoiding amplification bias and using long reads to assemble the bacterial genomes allowed for the circularization of the Wolbachia endosymbiont of Diaphorina citri for the first time and paralleled the annotation context of all four reference genomes without utilizing a traditional hybrid assembly. The strategies detailed here are suitable for the sequencing of other insect systems for which the input DNA, time, and cost are an issue.

Genetic and genomic architecture in eight strains of the laboratory opossum Monodelphis domestica

The gray short-tailed opossum (Monodelphis domestica) is an established laboratory-bred marsupial model for biomedical research. It is a critical species for comparative genomics research, providing the pivotal phylogenetic outgroup for studies of derived vs ancestral states of genomic/epigenomic characteristics for eutherian mammal lineages. To characterize the current genetic profile of this laboratory marsupial, we examined 79 individuals from eight established laboratory strains. Double digest restriction site-associated DNA sequencing and whole-genome resequencing experiments were performed to investigate the genetic architecture in these strains. A total of 66,640 high-quality single nucleotide polymorphisms (SNPs) were identified. We analyzed SNP density, average heterozygosity, nucleotide diversity, and population differentiation parameter Fst within and between the eight strains. Principal component and population structure analysis clearly resolve the strains at the level of their ancestral founder populations, and the genetic architecture of these strains correctly reflects their breeding history. We confirmed the successful establishment of the first inbred laboratory opossum strain LSD (inbreeding coefficient F > 0.99) and a nearly inbred strain FD2M1 (0.98 < F < 0.99), each derived from a different ancestral background. These strains are suitable for various experimental protocols requiring controlled genetic backgrounds and for intercrosses and backcrosses that can generate offspring with informative SNPs for studying a variety of genetic and epigenetic processes. Together with recent advances in reproductive manipulation and CRISPR/Cas9 techniques for Monodelphis domestica, the existence of distinctive inbred strains will enable genome editing on different genetic backgrounds, greatly expanding the utility of this marsupial model for biomedical research.

A DNA barcode-based survey of wild urban bees in the Loire Valley, France

The current decline of wild bees puts important ecosystem services such as pollination at risk. Both inventory and monitoring programs are needed to understand the causes of wild bee decline. Effective insect monitoring relies on both mass-trapping methods coupled with rapid and accurate identifications. Identifying wild bees using only morphology can be challenging, in particular, specimens from mass-trapped samples which are often in poor condition. We generated DNA barcodes for 2931 specimens representing 157 species (156 named and one unnamed species) and 28 genera. Automated cluster delineation reveals 172 BINs (Barcodes Index Numbers). A total of 36 species (22.93%) were found in highly urbanized areas. The majority of specimens, representing 96.17% of the species barcoded form reciprocally exclusive groups, allowing their unambiguous identification. This includes several closely related species notoriously difficult to identify. A total of 137 species (87.26%) show a "one-to-one" match between a named species and the BIN assignment. Fourteen species (8.92%) show deep conspecific lineages with no apparent morphological differentiation. Only two species pairs shared the same BIN making their identification with DNA barcodes alone uncertain. Therefore, our DNA barcoding reference library allows reliable identification by non-experts for the vast majority of wild bee species in the Loire Valley.

Genetic structure across urban and agricultural landscapes reveals evidence of resource specialization and philopatry in the Eastern carpenter bee, Xylocopa virginica L

Human activity continues to impact global ecosystems, often by altering the habitat suitability, persistence, and movement of native species. It is thus critical to examine the population genetic structure of key ecosystemservice providers across human-altered landscapes to provide insight into the forces that limit wildlife persistence and movement across multiple spatial scales. While some studies have documented declines of bee pollinators as a result of human-mediated habitat alteration, others suggest that some bee species may benefit from altered land use due to increased food or nesting resource availability; however, detailed population and dispersal studies have been lacking. We investigated the population genetic structure of the Eastern carpenter bee, Xylocopa virginica, across 14 sites spanning more than 450 km, including dense urban areas and intensive agricultural habitat. X. virginica is a large bee which constructs nests in natural and human-associated wooden substrates, and is hypothesized to disperse broadly across urbanizing areas. Using 10 microsatellite loci, we detected significant genetic isolation by geographic distance and significant isolation by land use, where urban and cultivated landscapes were most conducive to gene flow. This is one of the first population genetic analyses to provide evidence of enhanced insect dispersal in human-altered areas as compared to semi-natural landscapes. We found moderate levels of regional-scale population structure across the study system (G'ST = 0.146) and substantial co-ancestry between the sampling regions, where co-ancestry patterns align with major human transportation corridors, suggesting that human-mediated movement may be influencing regional dispersal processes. Additionally, we found a signature of strong site-level philopatry where our analyses revealed significant levels of high genetic relatedness at very fine scales (<1 km), surprising given X. virginica's large body size. These results provide unique genetic evidence that insects can simultaneously exhibit substantial regional dispersal as well as high local nesting fidelity in landscapes dominated by human activity.

Genome wide identification of QTL associated with yield and yield components in two popular wheat cultivars TAM 111 and TAM 112

Two drought-tolerant wheat cultivars, ‘TAM 111’ and ‘TAM 112’, have been widely grown in the Southern Great Plains of the U.S. and used as parents in many wheat breeding programs worldwide. This study aimed to reveal genetic control of yield and yield components in the two cultivars under both dryland and irrigated conditions. A mapping population containing 124 F_5:7 recombinant inbred lines (RILs) was developed from the cross of TAM 112/TAM 111. A set of 5,948 SNPs from the wheat 90K iSelect array and double digest restriction-site associated DNA sequencing was used to construct high-density genetic maps. Data for yield and yield components were obtained from 11 environments. QTL analyses were performed based on 11 individual environments, across all environments, within and across mega-environments. Thirty-six unique consistent QTL regions were distributed on 13 chromosomes including 1A, 1B, 1D, 2A, 2D, 3D, 4B, 4D, 6A, 6B, 6D, 7B, and 7D. Ten unique QTL with pleiotropic effects were identified on four chromosomes and eight were in common with the consistent QTL. These QTL increased dry biomass grain yield by 16.3 g m^-2, plot yield by 28.1 g m^-2, kernels spike^-1 by 0.7, spikes m^-2 by 14.8, thousand kernel weight by 0.9 g with favorable alleles from either parent. TAM 112 alleles mainly increased spikes m^-2 and thousand kernel weight while TMA 111 alleles increased kernels spike^-1, harvest index and grain yield. The saturated genetic map and markers linked to significant QTL from this study will be very useful in developing high throughput genotyping markers for tracking the desirable haplotypes of these important yield-related traits in popular parental cultivars.

Disjunction between canola distribution and the genetic structure of its recently described pest, the canola flower midge (Contarinia brassicola)

Population genomics is a useful tool to support integrated pest management as it can elucidate population dynamics, demography, and histories of invasion. Here, we use a restriction site-associated DNA sequencing approach combined with whole-genome amplification (WGA) to assess genomic population structure of a newly described pest of canola, the diminutive canola flower midge, Contarinia brassicola. Clustering analyses recovered little geographic structure across the main canola production region but differentiated several geographically disparate populations at edges of the agricultural zone. Given a lack of alternative hypotheses for this pattern, we suggest these data support alternative hosts for this species and thus our canola-centric view of this midge as a pest has limited our understanding of its biology. These results speak to the need for increased surveying efforts across multiple habitats and other potential hosts within Brassicaceae to improve both our ecological and evolutionary knowledge of this species and contribute to effective management strategies. We additionally found that use of WGA prior to library preparation was an effective method for increasing DNA quantity of these small insects prior to restriction site-associated DNA sequencing and had no discernible impact on genotyping consistency for population genetic analysis; WGA is therefore likely to be tractable for other similar studies that seek to randomly sample markers across the genome in small organisms.