Leon N. Court

Ecological assessment of estuarine sediments by pyrosequencing eukaryotic ribosomal DNA

Biodiversity assessment underpins our understanding of ecosystems and determines environmental management decisions on resource use and conservation priorities. Recently, a new discipline – environmental or ecological genomics (ecogenomics) – has emerged from major advances in sequencing technologies, such as pyrosequencing (a technique based on the detection of pyrophosphate during nucleotide incorporation), and enabled extraordinary progress in the way biodiversity can be assessed. Since 2008, numerous high-impact microbial metagenomic sequencing studies, which have relied on both classical and next-generation sequencing, have been published. As a result, many previously unrecognized taxa and biota have been identified, but none of these studies explored eukaryote diversity. Here, we illustrate the power of applying next-generation pyrosequencing to identify and enumerate eukaryote species assemblages in the context of assessing the impacts of human activity on ecosystems.

Impacts of inundation and drought on eukaryote biodiversity in semi-arid floodplain soils

Floodplain ecosystems are characterized by alternating wet and dry phases and periodic inundation defines their ecological character. Climate change, river regulation and the construction of levees have substantially altered natural flooding and drying regimes worldwide with uncertain effects on key biotic groups. In southern Australia, we hypothesized that soil eukaryotic communities in climate change affected areas of a semi‐arid floodplain would transition towards comprising mainly dry‐soil specialist species with increasing drought severity. Here, we used 18S rRNA amplicon pyrosequencing to measure the eukaryote community composition in soils that had been depleted of water to varying degrees to confirm that reproducible transitional changes occur in eukaryotic biodiversity on this floodplain. Interflood community structures (3 years post‐flood) were dominated by persistent rather than either aquatic or dry‐specialist organisms. Only 2% of taxa were unique to dry locations by 8 years post‐flood, and 10% were restricted to wet locations (inundated a year to 2 weeks post‐flood). Almost half (48%) of the total soil biota were detected in both these environments. The discovery of a large suite of organisms able to survive nearly a decade of drought, and up to a year submerged supports the concept of inherent resilience of Australian semi‐arid floodplain soil communities under increasing pressure from climatic induced changes in water availability.

Metabarcoding of benthic eukaryote communities predicts the ecological condition of estuaries

DNA-derived measurements of biological composition have the potential to produce data covering all of life, and provide a tantalizing proposition for researchers and managers. We used metabarcoding to compare benthic eukaryote composition from five estuaries of varying condition. In contrast to traditional studies, we found biotic richness was greatest in the most disturbed estuary, with this being due to the large volume of extraneous material (i.e. run-off from aquaculture, agriculture and other catchment activities) being deposited in the system. In addition, we found strong correlations between composition and a number of environmental variables, including nutrients, pH and turbidity. A wide range of taxa responded to these environmental gradients, providing new insights into their sensitivities to natural and anthropogenic stressors. Metabarcoding has the capacity to bolster current monitoring techniques, enabling the decisions regarding ecological condition to be based on a more holistic view of biodiversity.

Improved Inference of Taxonomic Richness from Environmental DNA

Accurate estimation of biological diversity in environmental DNA samples using high-throughput amplicon pyrosequencing must account for errors generated by PCR and sequencing. We describe a novel approach to distinguish the underlying sequence diversity in environmental DNA samples from errors that uses information on the abundance distribution of similar sequences across independent samples, as well as the frequency and diversity of sequences within individual samples. We have further refined this approach into a bioinformatics pipeline, Amplicon Pyrosequence Denoising Program (APDP) that is able to process raw sequence datasets into a set of validated sequences in formats compatible with commonly used downstream analyses packages. We demonstrate, by sequencing complex environmental samples and mock communities, that APDP is effective for removing errors from deeply sequenced datasets comprising biological and technical replicates, and can efficiently denoise single-sample datasets. APDP provides more conservative diversity estimates for complex datasets than other approaches; however, for some applications this may provide a more accurate and appropriate level of resolution, and result in greater confidence that returned sequences reflect the diversity of the underlying sample.

Unusually High Frequency of Genes Encoding Vegetative Insecticidal Proteins in an Australian Bacillus thuringiensis Collection

Of 188 Australian Bacillus thuringiensis strains screened for genes encoding soluble insecticidal proteins by polymerase chain reaction/restriction-length fragment polymorphism (RFLP) analysis, 87% showed the presence of such genes. Although 135 isolates (72%) produced an RFLP pattern identical to that expected for vip3A genes, 29 isolates possessed a novel vip-like gene. The novel vip-like gene was cloned from B. thuringiensis isolate C81, and sequence analysis demonstrated that it was 94% identical to the vip3Ba1 gene. The new gene was designated vip3Bb2. Cell-free supernatants from both the B. thuringiensis strain C81 and from Escherichia coli expressing the Vip3Bb2 protein were toxic for the cotton bollworm, Helicoverpa armigera.

DNA barcoding to support conservation: species identification, genetic structure and biogeography of fishes in the Murray)Darling River Basin, Australia

Freshwater fish stocks worldwide are under increasing threat of overfishing, disease, pollution and competition from introduced species. In the Murray—Darling Basin (MDB), the largest river system of Australia, more than half the native species are listed as rare or endangered. Active management is required to counteract reduction in population sizes, prevent local extinctions and to maintain genetic diversity. We describe the first comprehensive set of DNA barcodes able to discriminate between all 58 native and introduced species of freshwater fish recorded in the MDB. These barcodes also distinguish populations from those in adjacent basins, with estimated separation times as short as 0.1 million years ago. We demonstrate the feasibility of using DNA fingerprinting of ribosomal RNA (12S and 18S rRNA) genes and mitochondrial DNA control region (mtDNA CR) sequences to identify species from eggs, larvae, tissues and predator gut contents as well as differentiate populations, morphologically cryptic species and hybrids. The DNA barcode resource will enhance capacity in many areas of fish conservation biology that can benefit from improved knowledge of genetic provenance. These include captive breeding and restocking programs, life history studies and ecological research into the interactions between populations of native and exotic species.

Complete mitochondrial DNA genome of Bemisia tabaci cryptic pest species complex Asia I (Hemiptera: Aleyrodidae)

Abstract The complete length of the Asia I member of the Bemisia tabaci species complex mitochondrial DNA genome (mitogenome) is 15,210 bp (GenBank accession no. KJ778614) with an A-T biased nucleotide composition (A: 32.7%; T: 42.4%; G: 14.0%; C: 10.8%). The mitogenome consists of 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNA (rRNAs) and a 467 bp putative control region which also includes the A+T rich repeat region. All PCGs have an ATA (n = 8) or ATG (n = 5) start codon. Gene synteny of Asia I is overall similar to B. afer and two other members of the B. tabaci species complex Mediterranean and New World 1, and contains the tRNA-Ser2 located between the Cytb and ND1 genes found in Mediterranean and New World 1, but which is absent in B. afer. The orientation of the tRNA-Arg in Asia I is on the “plus” strand and differed from Mediterranean which is found on the “minus” strand. The Asia I mitogenome size is currently ranked the second smallest after B. afer (14,968 bp) followed by New World 1 (15,322 bp) and Mediterranean (15,632 bp).

The Trouble with MEAM2: Implications of Pseudogenes on Species Delimitation in the Globally Invasive Bemisia tabaci (Hemiptera: Aleyrodidae) Cryptic Species Complex

Abstract Molecular species identification using suboptimal PCR primers can over-estimate species diversity due to coamplification of nuclear mitochondrial (NUMT) DNA/pseudogenes. For the agriculturally important whitefly Bemisia tabaci cryptic pest species complex, species identification depends primarily on characterization of the mitochondrial DNA cytochrome oxidase I (mtDNA COI) gene. The lack of robust PCR primers for the mtDNA COI gene can undermine correct species identification which in turn compromises management strategies. This problem is identified in the B. tabaci Africa/Middle East/Asia Minor clade which comprises the globally invasive Mediterranean (MED) and Middle East Asia Minor I (MEAM1) species, Middle East Asia Minor 2 (MEAM2), and the Indian Ocean (IO) species. Initially identified from the Indian Ocean island of Réunion, MEAM2 has since been reported from Japan, Peru, Turkey and Iraq. We identified MEAM2 individuals from a Peruvian population via Sanger sequencing of the mtDNA COI gene. In attempting to characterize the MEAM2 mitogenome, we instead characterized mitogenomes of MEAM1. We also report on the mitogenomes of MED, AUS, and IO thereby increasing genomic resources for members of this complex. Gene synteny (i.e., same gene composition and orientation) was observed with published B. tabaci cryptic species mitogenomes. Pseudogene fragments matching MEAM2 partial mtDNA COI gene exhibited low frequency single nucleotide polymorphisms that matched low copy number DNA fragments (<3%) of MEAM1 genomes, whereas presence of internal stop codons, loss of expected stop codons and poor primer annealing sites, all suggested MEAM2 as a pseudogene artifact and so not a real species.

Novel molecular approach to define pest species status and tritrophic interactions from historical Bemisia specimens

Museum specimens represent valuable genomic resources for understanding host-endosymbiont/parasitoid evolutionary relationships, resolving species complexes and nomenclatural problems. However, museum collections suffer DNA degradation, making them challenging for molecular-based studies. Here, the mitogenomes of a single 1912 Sri Lankan Bemisia emiliae cotype puparium, and of a 1942 Japanese Bemisia puparium are characterised using a Next-Generation Sequencing approach. Whiteflies are small sap-sucking insects including B. tabaci pest species complex. Bemisia emiliae’s draft mitogenome showed a high degree of homology with published B. tabaci mitogenomes, and exhibited 98–100% partial mitochondrial DNA Cytochrome Oxidase I (mtCOI) gene identity with the B. tabaci species known as Asia II-7. The partial mtCOI gene of the Japanese specimen shared 99% sequence identity with the Bemisia ‘JpL’ genetic group. Metagenomic analysis identified bacterial sequences in both Bemisia specimens, while hymenopteran sequences were also identified in the Japanese Bemisia puparium, including complete mtCOI and rRNA genes, and various partial mtDNA genes. At 88–90% mtCOI sequence identity to Aphelinidae wasps, we concluded that the 1942 Bemisia nymph was parasitized by an Eretmocerus parasitoid wasp. Our approach enables the characterisation of genomes and associated metagenomic communities of museum specimens using 1.5 ng gDNA, and to infer historical tritrophic relationships in Bemisia whiteflies.

Standardised molecular diagnostic tool for the identification of cryptic species within the Bemisia tabaci complex

BACKGROUND The whitefly Bemisia tabaci complex harbours over 40 cryptic species that have been placed in 11 phylogenetically distinct clades based on the molecular characterization of partial mitochondrial DNA COI (mtCOI) gene region. Four cryptic species are currently within the invasive clade, i.e. MED, MEAM1, MEAM2 and IO. Correct identification of these species is a critical step towards implementing reliable measures for plant biosecurity and border protection; however, no standardized B. tabaci-specific primers are currently available which has caused inconsistencies in the species identification processes. RESULTS We report three sets of polymerase chain reaction (PCR) primers developed to amplify the mtCOI region which can be used for genotyping MED, MEAM1 and IO species, and tested these primers on 91 MED, 35 MEAM1 and five IO individuals. PCR and sequencing of amplicons identified a total of 21, six and one haplotypes in MED, MEAM1 and IO respectively, of which six haplotypes were new to the B. tabaci database. CONCLUSION These primer pairs enabled standardization and robust molecular species identification via mtCOI screening of the targeted invasive cryptic species and will improve quarantine decisions. Use of this diagnostic tool could be extended to other species within the complex.