Metabolic disruptions and impaired reproductive fitness in wild-caught freshwater turtles (Emydura macquarii macquarii) exposed to elevated per- and polyfluoroalkyl substances (PFAS).

Per- and poly-fluoroalkyl substances (PFAS) pose a threat to organisms and ecosystems due to their persistent nature. Ecotoxicology endpoints used in regulatory guidelines may not reflect multiple, low-level but persistent stressors. This study examines the biological effects of PFAS on Eastern short-necked turtles in Queensland, Australia. In this study, blood samples were collected and analysed for PFAS, hormone levels, and functional omics endpoints. High levels of PFAS were found in turtles at the impacted site, with PFOS being the dominant constituent. The PFAS profiles of males and females differed, with males having higher PFAS concentrations. Hormone concentrations differed between impacted and reference sites in male turtles, with elevated testosterone and corticosterone indicative of stress. Further, energy utilisation, nucleotide synthesis, nitrogen metabolism, and amino acid synthesis were altered in both male and female turtles from PFAS-impacted sites. Both sexes show similar metabolic responses to environmental stressors from the PFAS-contaminated site, which may adversely affect their reproductive fitness. Purine metabolism, caffeine metabolism, and ferroptosis pathway changes in turtles can cause gout, cell death, and overall health problems. Further, the study showed that prolonged exposure to elevated PFAS levels in the wild could compromise turtle reproductive fitness by disrupting reproductive steroids and metabolic pathways.

Complex multiple introductions drive fall armyworm invasions into Asia and Australia.

The fall armyworm (FAW) Spodoptera frugiperda is thought to have undergone a rapid 'west-to-east' spread since 2016 when it was first identified in western Africa. Between 2018 and 2020, it was recorded from South Asia (SA), Southeast Asia (SEA), East Asia (EA), and Pacific/Australia (PA). Population genomic analyses enabled the understanding of pathways, population sources, and gene flow in this notorious agricultural pest species. Using neutral single nucleotide polymorphic (SNP) DNA markers, we detected genome introgression that suggested most populations in this study were overwhelmingly C- and R-strain hybrids (n = 252/262). SNP and mitochondrial DNA markers identified multiple introductions that were most parsimoniously explained by anthropogenic-assisted spread, i.e., associated with international trade of live/fresh plants and plant products, and involved 'bridgehead populations' in countries to enable successful pest establishment in neighbouring countries. Distinct population genomic signatures between Myanmar and China do not support the 'African origin spread' nor the 'Myanmar source population to China' hypotheses. Significant genetic differentiation between populations from different Australian states supported multiple pathways involving distinct SEA populations. Our study identified Asia as a biosecurity hotspot and a FAW genetic melting pot, and demonstrated the use of genome analysis to disentangle preventable human-assisted pest introductions from unpreventable natural pest spread.

Omics-based ecosurveillance for the assessment of ecosystem function, health, and resilience.

Current environmental monitoring efforts often focus on known, regulated contaminants ignoring the potential effects of unmeasured compounds and/or environmental factors. These specific, targeted approaches lack broader environmental information and understanding, hindering effective environmental management and policy. Switching to comprehensive, untargeted monitoring of contaminants, organism health, and environmental factors, such as nutrients, temperature, and pH, would provide more effective monitoring with a likely concomitant increase in environmental health. However, even this method would not capture subtle biochemical changes in organisms induced by chronic toxicant exposure. Ecosurveillance is the systematic collection, analysis, and interpretation of ecosystem health-related data that can address this knowledge gap and provide much-needed additional lines of evidence to environmental monitoring programs. Its use would therefore be of great benefit to environmental management and assessment. Unfortunately, the science of 'ecosurveillance', especially omics-based ecosurveillance is not well known. Here, we give an overview of this emerging area and show how it has been beneficially applied in a range of systems. We anticipate this review to be a starting point for further efforts to improve environmental monitoring via the integration of comprehensive chemical assessments and molecular biology-based approaches. Bringing multiple levels of omics technology-based assessment together into a systems-wide ecosurveillance approach will bring a greater understanding of the environment, particularly the microbial communities upon which we ultimately rely to remediate perturbed ecosystems.

Omics-based ecosurveillance uncovers the influence of estuarine macrophytes on sediment microbial function and metabolic redundancy in a tropical ecosystem.

Vertical zonation within estuarine ecosystems can strongly influence microbial diversity and function by regulating competition, predation, and environmental stability. The degree to which microbial communities exhibit horizontal patterns through an estuary has received comparatively less attention. Here, we take a multi-omics ecosurveillance approach to study environmental gradients created by the transition between dominant vegetation types along a near pristine tropical river system (Wenlock River, Far North Queensland, Australia). The study sites included intertidal mudflats fringed by saltmarsh, mangrove or mixed soft substrata habitats. Collected sediments were analyzed for eukaryotes and prokaryotes using small sub-unit (SSU) rRNA gene amplicons to profile the relative taxonomic composition. Central carbon metabolism metabolites and other associated organic polar metabolites were analyzed using established metabolomics-based approaches, coupled with total heavy metals analysis. Eukaryotic taxonomic information was found to be more informative of habitat type. Bacterial taxonomy and community composition also showed habitat-specificity, with phyla Proteobacteria and Cyanobacteria strongly linked to mangroves and saltmarshes, respectively. In contrast, metabolite profiling was critical for understanding the biochemical pathways and expressed functional outputs in these systems that were tied to predicted microbial gene function (16S rRNA). A high degree of metabolic redundancy was observed in the bacterial communities, with the metabolomics data suggesting varying degrees of metabolic criticality based on habitat type. The predicted functions of the bacterial taxa combined with annotated metabolites accounted for the conservative perspective of microbial community redundancy against the putative metabolic pathway impacts in the metabolomics data. Coupling these data demonstrates that habitat-mediated estuarine gradients drive patterns of community diversity and metabolic function and highlights the real redundancy potential of habitat microbiomes. This information is useful as a point of comparison for these sensitive ecosystems and provides a framework for identifying potentially vulnerable or at-risk systems before they are significantly degraded.

Moving to Keep Fit: Feeding Behavior and Movement of Helicoverpa armigera (Lepidoptera: Noctuidae) on Artificial Diet With Different Protein: Carbohydrate Ratios.

Insect herbivores can modify their foraging behavior to obtain a balanced food intake, and they tend to move between food sources with different nutrient values. We investigated this movement in early instar larvae of Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) using a putative optimal artificial diet (OP) and high protein (HP) and high carbohydrate (HC) artificial diets based on protein (p) and carbohydrate (c) ratios. Larvae were allowed to choose between the same kind of diet cubes (effectively no-choice), or diet cubes with different p: c ratios. In no-choice tests, we found that first instar larvae remained longest on OP diet and spent the least time on HC diet, while third instar larvae remained longest on HC diet and spent least time on OP diet. First instar larvae moved the most when provided with HC diet, while third instar larvae moved most when provided with OP diet. However, both stages moved the least when allowed to choose between diet cubes with different p: c ratios. The relative growth rate decreased when larvae increased their movement, but this influence was not evident when larvae fed on HC diet. Larvae that fed only on HC diet had the highest relative growth rate, followed by larvae with access to all diets simultaneously, indicating a behavior to mix nutrient intake. We relate these findings to behavior of this major pest species under field conditions.

Multiple recombination events between two cytochrome P450 loci contribute to global pyrethroid resistance in Helicoverpa armigera.

The cotton bollworm, Helicoverpa armigera (Hübner) is one of the most serious insect pest species to evolve resistance against many insecticides from different chemical classes. This species has evolved resistance to the pyrethroid insecticides across its native range and is becoming a truly global pest after establishing in South America and having been recently recorded in North America. A chimeric cytochrome P450 gene, CYP337B3, has been identified as a resistance mechanism for resistance to fenvalerate and cypermethrin. Here we show that this resistance mechanism is common around the world with at least eight different alleles. It is present in South America and has probably introgressed into its closely related native sibling species, Helicoverpa zea. The different alleles of CYP337B3 are likely to have arisen independently in different geographic locations from selection on existing diversity. The alleles found in Brazil are those most commonly found in Asia, suggesting a potential origin for the incursion of H. armigera into the Americas.

ABC transporter mis-splicing associated with resistance to Bt toxin Cry2Ab in laboratory- and field-selected pink bollworm.

Evolution of pest resistance threatens the benefits of genetically engineered crops that produce Bacillus thuringiensis (Bt) insecticidal proteins. Strategies intended to delay pest resistance are most effective when implemented proactively. Accordingly, researchers have selected for and analyzed resistance to Bt toxins in many laboratory strains of pests before resistance evolves in the field, but the utility of this approach depends on the largely untested assumption that laboratory- and field-selected resistance to Bt toxins are similar. Here we compared the genetic basis of resistance to Bt toxin Cry2Ab, which is widely deployed in transgenic crops, between laboratory- and field-selected populations of the pink bollworm (Pectinophora gossypiella), a global pest of cotton. We discovered that resistance to Cry2Ab is associated with mutations disrupting the same ATP-binding cassette transporter gene (PgABCA2) in a laboratory-selected strain from Arizona, USA, and in field-selected populations from India. The most common mutation, loss of exon 6 caused by alternative splicing, occurred in resistant larvae from both locations. Together with previous data, the results imply that mutations in the same gene confer Bt resistance in laboratory- and field-selected strains and suggest that focusing on ABCA2 genes may help to accelerate progress in monitoring and managing resistance to Cry2Ab.

Detection of sister-species in invasive populations of the fall armyworm Spodoptera frugiperda (Lepidoptera: Noctuidae) from Uganda.

The fall armyworm (FAW) Spodoptera frugiperda (J. E. Smith) is a species native to the Americas. This polyphagous lepidopteran pest was first reported in Nigeria and the Democratic Republic of São Tomé and Principe in 2016, but its presence in eastern Africa has not been confirmed via molecular characterisation. In this study, FAW specimens from western and central Uganda were identified based on the partial mtDNA COI gene sequences, with mtDNA COI haplotypes matching those identified in Nigeria and São Tomé. In this study, we sequence an additional partial mtDNA Cyt b gene and also the partial mtDNA COIII gene in Ugandan FAW samples. We detected identical mitochondrial DNA haplotypes for both the mtDNA Cyt b and COI partial genes, while combining the mtDNA COI/Cyt b haplotypes and mtDNA COIII haplotypes enabled a new maternal lineage in the Ugandan corn-preferred FAW samples to be identified. Our results suggested that the African incursions of S. frugiperda involved at least three maternal lineages. Recent full genome, phylogenetic and microsatellite analyses provided evidence to support S. frugiperda as likely consisted of two sympatric sister species known as the corn-preferred and rice-preferred strains. In our Ugandan FAW populations, we identified the presence of mtDNA haplotypes representative of both sister species. It is not known if both FAW sister species were originally introduced together or separately, and whether they have since spread as a single population. Further analyses of additional specimens originally collected from São Tomé, Nigeria and throughout Africa would be required to clarify this issue. Importantly, our finding showed that the genetic diversity of the African corn-preferred FAW species is higher than previously reported. This potentially contributed to the success of FAW establishment in Africa. Furthermore, with the additional maternal lineages detected, there is likely an increase in paternal lineages, thereby increasing the diversity of the African FAW population. Knowledge of the FAW genetic diversity will be needed to assess the risks of introducing Bt-resistance traits and to understand the FAW incursion pathways into the Old World and its potential onward spread. The agricultural implications of the presence of two evolutionary divergent FAW lineages (the corn and the rice lineage) in the African continent are further considered and discussed.

A genomic approach to identify and monitor a novel pyrethroid resistance mutation in the redlegged earth mite, Halotydeus destructor.

Resistance mechanisms are typically uncovered by identifying sequence variation in known candidate genes, however this strategy can be problematic for species with no reference data in known relatives. Here we take a genomic approach to identify resistance to pyrethroids in the redlegged earth mite, Halotydeus destructor, a member of the Penthalidae family of mites that are virtually uncharacterized genetically. Based on shallow genome sequencing followed by a genome assembly, we first identified contigs of the H. destructor parasodium channel gene. By linking variation in this gene to known resistant phenotypes, we located a single nucleotide polymorphism in resistant mites. This polymorphism results in a leucine (L) to phenylalanine (F) amino acid substitution in the II6 region (predicted) of the gene (L1024F). This novel mutation has not previously been linked to pyrethroid resistance, although other polymorphisms have been identified in the two-spotted spider mite, Tetranychus urticae at the same locus (L1024V). The sequencing approach was successful in generating a candidate polymorphism that was then validated using laboratory bioassays and field surveys. A high throughput Illumina-based sequencing diagnostic was developed to rapidly assess resistance allele frequencies in pools of mites sourced from hundreds of populations across Australia. Resistance was confirmed to be widespread in the southern wheatbelt region of Western Australia. Two different resistance mutations were identified in field populations, both resulting in the same amino acid substitution. The frequency and distribution of resistance amplicon haplotypes suggests at least two, and probably more independent origins of resistance.

CRISPR/Cas9 mediated genome editing of Helicoverpa armigera with mutations of an ABC transporter gene HaABCA2 confers resistance to Bacillus thuringiensis Cry2A toxins.

High levels of resistance to Bt toxin Cry2Ab have been identified to be genetically linked with loss of function mutations of an ABC transporter gene (ABCA2) in two lepidopteran insects, Helicoverpa armigera and Helicoverpa punctigera. To further confirm the causal relationship between the ABCA2 gene (HaABCA2) and Cry2Ab resistance in H. armigera, two HaABCA2 knockout strains were created from the susceptible SCD strain with the CRISPR/Cas9 genome editing system. One strain (SCD-A2KO1) is homozygous for a 2-bp deletion in exon 2 of HaABCA2 created by non-homologous end joining (NHEJ). The other strain (SCD-A2KO2) is homozygous for a 5-bp deletion in exon 18 of HaABCA2 made by homology-directed repair (HDR), which was produced to mimic the r2 resistance allele of a field-derived Cry2Ab-resistant strain from Australia. Both knockout strains obtained high levels of resistance to both Cry2Aa (>120-fold) and Cry2Ab (>100-fold) compared with the original SCD strain, but no or very limited resistance to Cry1Ac (<4-fold). Resistance to Cry2Ab in both knockouts is recessive, and genetic complementary tests confirmed Cry2Ab resistance alleles are at the same locus (i.e. HaABCA2) for the two strains. Brush border membrane vesicles (BBMVs) of midguts from both knockout strains lost binding with Cry2Ab, but maintained the same binding with Cry1Ac as the SCD strain. In vivo functional evidence from this study demonstrates knockout of HaABCA2 confers high levels of resistance to both Cry2Aa and Cry2Ab, confirming that HaABCA2 plays a key role in mediating toxicity of both Cry2Aa and Cry2Ab against H. armigera.

Mitochondrial DNA and trade data support multiple origins of Helicoverpa armigera (Lepidoptera, Noctuidae) in Brazil.

The Old World bollworm Helicoverpa armigera is now established in Brazil but efforts to identify incursion origin(s) and pathway(s) have met with limited success due to the patchiness of available data. Using international agricultural/horticultural commodity trade data and mitochondrial DNA (mtDNA) cytochrome oxidase I (COI) and cytochrome b (Cyt b) gene markers, we inferred the origins and incursion pathways into Brazil. We detected 20 mtDNA haplotypes from six Brazilian states, eight of which were new to our 97 global COI-Cyt b haplotype database. Direct sequence matches indicated five Brazilian haplotypes had Asian, African, and European origins. We identified 45 parsimoniously informative sites and multiple substitutions per site within the concatenated (945 bp) nucleotide dataset, implying that probabilistic phylogenetic analysis methods are needed. High diversity and signatures of uniquely shared haplotypes with diverse localities combined with the trade data suggested multiple incursions and introduction origins in Brazil. Increasing agricultural/horticultural trade activities between the Old and New Worlds represents a significant biosecurity risk factor. Identifying pest origins will enable resistance profiling that reflects countries of origin to be included when developing a resistance management strategy, while identifying incursion pathways will improve biosecurity protocols and risk analysis at biosecurity hotspots including national ports.

Complete Mitochondrial Genome of Helicoverpa zea (Lepidoptera: Noctuidae) and Expression Profiles of Mitochondrial-Encoded Genes in Early and Late Embryos.

The mitochondrial genome (mitogenome) of the bollworm, Helicoverpa zea (Boddie), was assembled using paired-end nucleotide sequence reads generated with a next-generation sequencing platform. Assembly resulted in a mitogenome of 15,348 bp with greater than 17,000-fold average coverage. Organization of the H. zea mitogenome (gene order and orientation) was identical to other known lepidopteran mitogenome sequences. Compared with Helicoverpa armigera (Hübner) mitogenome, there were a few differences in the lengths of gaps between genes, but the lengths of nucleotide overlaps were essentially conserved between the two species. Nucleotide composition of the H. zea mitochondrial genome was very similar to those of the related species H. armigera and Helicoverpa punctigera Wallengren. Mapping of RNA-Seq reads obtained from 2-h eggs and 48-h embryos to protein coding genes (PCG) revealed that all H. zea PCGs were processed as single mature gene transcripts except for the bicistronic atp8 + atp6 transcript. A tRNA-like sequence predicted to form a hammer-head-like secondary structure that may play a role in transcription start and mitogenome replication was identified within the control region of the H. zea mitogenome. Similar structures were also found within the control regions of several other lepidopteran species. Expression analysis revealed significant differences in levels of expression of PCGs within each developmental stage, but the pattern of variation was similar in both developmental stages analyzed in this study. Mapping of RNA-Seq reads to PCG transcripts also identified transcription termination and polyadenylation sites that differed from the sites described in other lepidopteran species.

Characterization of the complete mitochondrial genome of the Australian Heliothine moth, Australothis rubrescens (Lepidoptera: Noctuidae).

Australothis rubrescens is basal to the Helicoverpa lineage containing pests such as Helicoverpa armigera, H. assulta and H. gelotopoeon. An illumina library of DNA from A. rubrescens was constructed and shallow sequencing and assembly of the DNA was conducted. The complete mitochondrial genome was identified using similarity to the H. armigera mitochondrial genome. The mitochondrial genome of A. rubrescens is 15,382 bp in length. It contains 37 genes which are shared with the vast majority of animals: 13 protein-coding genes (PCGs), 2 ribosomal RNAs, 22 transfer RNAs and a non-coding AT-rich region (Table 1). As found in other Lepidopterans, the arrangement of all tRNAs of the A. rubrescens is identical to most insects. The complete mitochondrial genome of A. rubrescens will be an important tool in understanding the evolutionary history of the Heliothine moths.

Insect Resistance to Bacillus thuringiensis Toxin Cry2Ab Is Conferred by Mutations in an ABC Transporter Subfamily A Protein.

The use of conventional chemical insecticides and bacterial toxins to control lepidopteran pests of global agriculture has imposed significant selection pressure leading to the rapid evolution of insecticide resistance. Transgenic crops (e.g., cotton) expressing the Bt Cry toxins are now used world wide to control these pests, including the highly polyphagous and invasive cotton bollworm Helicoverpa armigera. Since 2004, the Cry2Ab toxin has become widely used for controlling H. armigera, often used in combination with Cry1Ac to delay resistance evolution. Isolation of H. armigera and H. punctigera individuals heterozygous for Cry2Ab resistance in 2002 and 2004, respectively, allowed aspects of Cry2Ab resistance (level, fitness costs, genetic dominance, complementation tests) to be characterised in both species. However, the gene identity and genetic changes conferring this resistance were unknown, as was the detailed Cry2Ab mode of action. No cross-resistance to Cry1Ac was observed in mutant lines. Biphasic linkage analysis of a Cry2Ab-resistant H. armigera family followed by exon-primed intron-crossing (EPIC) marker mapping and candidate gene sequencing identified three independent resistance-associated INDEL mutations in an ATP-Binding Cassette (ABC) transporter gene we named HaABCA2. A deletion mutation was also identified in the H. punctigera homolog from the resistant line. All mutations truncate the ABCA2 protein. Isolation of further Cry2Ab resistance alleles in the same gene from field H. armigera populations indicates unequal resistance allele frequencies and the potential for Bt resistance evolution. Identification of the gene involved in resistance as an ABC transporter of the A subfamily adds to the body of evidence on the crucial role this gene family plays in the mode of action of the Bt Cry toxins. The structural differences between the ABCA2, and that of the C subfamily required for Cry1Ac toxicity, indicate differences in the detailed mode-of-action of the two Bt Cry toxins.

ACR-26: a novel nicotinic receptor subunit of parasitic nematodes.

Nematode nicotinic acetylcholine receptors are the targets for many effective anthelmintics, including those recently introduced into the market. We have identified a novel nicotinic receptor subunit sequence, acr-26, that is expressed in all the animal parasitic nematodes we examined from clades III, IV and V, but is not present in the genomes of Trichinella spiralis, Caenorhabditis elegans, Pristionchus pacificus and Meloidogyne spp. In Ascaris suum, ACR-26 is expressed on muscle cells isolated from the head, but not from the mid-body region. Sequence comparisons with other vertebrate and nematode subunits suggested that ACR-26 may be capable of forming a functional homomeric receptor; when acr-26 cRNA was injected into Xenopus oocytes along with Xenopus laevis ric-3 cRNA we occasionally observed the formation of acetylcholine- and nicotine-sensitive channels. The unreliable expression of ACR-26 in vitro may suggest that additional subunits or chaperones may be required for efficient formation of the functional receptors. ACR-26 may represent a novel target for the development of cholinergic anthelmintics specific for animal parasites.

The cys-loop ligand-gated ion channel gene family of Brugia malayi and Trichinella spiralis: a comparison with Caenorhabditis elegans.

Nematode cys-loop ligand gated ion channels (CLGIC) mediate neurotransmission and are important targets for anthelmintics in parasitic nematodes. The CLGIC superfamily in nematodes includes ion channels gated by acetylcholine, gamma-amino butyric acid (GABA), glutamate, glycine and 5-HT. The macrocyclic lactones and the nicotinic agonists are important groups of anthelmintics that target the glutamate gated chloride channels and the nicotinic acetylcholine receptors, respectively. The model organism Caenorhabditis elegans has the most diverse families of cys-loop LGIC known in any organism. Many parasitic nematodes have homologues of C. elegans receptors but to date no genome wide investigations have been done. The genome sequencing projects of Brugia malayi (clade III) and Trichinella spiralis (clade I) have allowed us to characterise the CLGIC families in these species. Although the main groups of CLGICs targeted by anthelmintics are represented in both the nematode genomes investigated here, the CLGIC family is much smaller in B. malayi and T. spiralis, suggesting that care must be taken when using C. elegans as a model organism for distantly related nematodes.

Detection and measurement of benzimidazole resistance alleles in Haemonchus contortus using real-time PCR with locked nucleic acid Taqman probes.

Benzimidazole resistance is a common problem in parasitic nematodes of ruminants and early detection is vital if its spread is to be monitored and controlled. Real time PCR offers a fast and reliable method for rapid detection and measurement of resistance allele frequencies. In Haemonchus contortus a single nucleotide polymorphism at codon 200 of the beta-tubulin gene (TTC to TAC), causing a phenylalanine to tyrosine amino acid substitution, has been shown to be involved in many cases of resistance. Locked nucleic acid (LNA) Taqman probes have been used in this work to detect and measure the frequency of resistance alleles in individual and multiple H. contortus. Detection of resistant genotypes using LNA Taqman probes in individual H. contortus is simpler and more reliable than with previously described assays. Measurement of the frequency of resistant alleles in populations of H. contortus was achieved by using the cycle threshold (C(t)) values and a standard curve derived from populations with known allele frequencies. Results using the LNA probes on individual and multiple worms gave similar results to the allele specific PCR. The sensitivity of the LNA assay on multiple nematodes allowed reliable detection of > or = 10% resistance allele frequency. Using the final fluorescence method, it was possible to differentiate populations with approximately 0, 5 and 10% resistance allele frequencies.