Overexpression of UDP-glucuronosyltransferase and cytochrome P450 enzymes confers resistance to sulfoxaflor in field populations of the aphid, Myzus persicae.

The green peach aphid, Myzus persicae, is a highly damaging, globally distributed crop pest that has evolved multiple resistance to numerous insecticides. It is thus imperative that insecticides that are not strongly compromised by pre-existing resistance are carefully managed to maximise their effective life span. Sulfoxaflor is a sulfoximine insecticide that retains efficacy against M. persicae clones that exhibit resistance to older insecticides. In the current study we monitored the efficacy of sulfoxaflor against M. persicae populations collected in Western Australia, following reports of control failures in this region. We identified clones with low (4-23-fold across multiple independent bioassay experiments), but significant, levels of resistance to sulfoxaflor compared with a reference susceptible clone. Furthermore, we demonstrate that sulfoxaflor resistance can persist after many months of culturing in the laboratory in the absence of insecticide exposure. Resistance was not conferred by known mechanisms of resistance to neonicotinoid insecticides, that act on the same target-site as sulfoxaflor, i.e. the R81T mutation or overexpresssion of the P450 gene CYP6CY3. Rather, transcriptome profiling of multiple resistant and susceptible clones identified the P450 CYP380C40 and the UDP-glucuronosyltransferase UGT344P2 as highly overexpressed (21-76-fold and 6-33-fold respectively) in the resistant clones. Transgenic expression of these genes demonstrated that they confer, low, but significant, levels of resistance to sulfoxaflor in vivo. Taken together, our data reveal the presence of low-level resistance to sulfoxaflor in M. persicae populations in Australia and uncover two novel mechanisms conferring resistance to this compound. The findings and tools generated in this study provide a platform for the development of strategies that aim to slow, prevent or overcome the evolution of more potent resistance to sulfoxaflor.

Fitness Costs Associated with Pyrethroid Resistance in Halotydeus destructor (Tucker) (Acari: Penthaleidae) Elucidated Through Semi-field Trials.

Pyrethroid resistance in the redlegged earth mite, Halotydeus destructor (Tucker), is primarily attributed to a kdr (knockdown resistance) mutation in the parasodium channel gene. To assess fitness costs associated with this resistance, adult resistant and susceptible populations were mixed in different proportions in microcosm tubs and placed in a shade-house simulating field conditions. Three separate experiments were undertaken whereby parental mites were collected from the field and offspring were followed for two to three generations. The association between fitness costs and kdr-mediated resistance was investigated by examining differences in mite numbers and changes in resistant allele frequencies across generations. In two (of the three) experiments, the population fitness measure of mites was significantly lower in microcosms containing a higher proportion of resistant individuals compared with treatments containing susceptible mites. No differences in mite fitness were observed between treatments in the third experiment; in this instance, the starting proportion of individuals homozygous for the resistant mutation was much lower (~40%) than in the other experiments (>90%). In all three experiments, a decrease in the resistant allele frequency across mite generations was observed. These findings indicate a potential deleterious pleiotropic effect of the kdr mutation on the fitness of H. destructor and have implications for resistance management strategies aimed at this important agricultural pest. Further experiments investigating fitness costs directly in the field are warranted.

Adaptive responses of free-living and symbiotic microalgae to simulated future ocean conditions.

Marine microalgae are a diverse group of microscopic eukaryotic and prokaryotic organisms capable of photosynthesis. They are important primary producers and carbon sinks but their physiology and persistence are severely affected by global climate change. Powerful experimental evolution technologies are being used to examine the potential of microalgae to respond adaptively to current and predicted future conditions, as well as to develop resources to facilitate species conservation and restoration of ecosystem functions. This review synthesizes findings and insights from experimental evolution studies of marine microalgae in response to elevated temperature and/or pCO2 . Adaptation to these environmental conditions has been observed in many studies of marine dinoflagellates, diatoms and coccolithophores. An enhancement in traits such as growth and photo-physiological performance and an increase in upper thermal limit have been shown to be possible, although the extent and rate of change differ between microalgal taxa. Studies employing multiple monoclonal replicates showed variation in responses among replicates and revealed the stochasticity of mutations. The work to date is already providing valuable information on species' climate sensitivity or resilience to managers and policymakers but extrapolating these insights to ecosystem- and community-level impacts continues to be a challenge. We recommend future work should include in situ experiments, diurnal and seasonal fluctuations, multiple drivers and multiple starting genotypes. Fitness trade-offs, stable versus plastic responses and the genetic bases of the changes also need investigating, and the incorporation of genome resequencing into experimental designs will be invaluable.

Separating two tightly linked species-defining phenotypes in Bactrocera with hybrid recombinant analysis.

BACKGROUND: Bactrocera tryoni and Bactrocera neohumeralis mate asynchronously; the former mates exclusively around dusk while the latter mates during the day. The two species also differ in the colour of the post-pronotal lobe (callus), which is predominantly yellow in B. tryoni and brown in B. neohumeralis. We have examined the genetic relationship between the two characters in hybrids, backcrosses and multigeneration hybrid progeny. RESULTS: Our analysis of the mating time of the parental species revealed that while B. tryoni mate exclusively at dusk, B. neohumeralis females pair with B. neohumeralis males during the day and with B. tryoni males at dusk. We found considerable variance in mating time and callus colour among hybrid backcross individuals of both sexes but there was a strong although not invariant trend for callus colour to co-segregate with mating time in both sexes. To genetically separate these two phenotypes we allowed the interspecific F1 hybrids to propagate for 25 generations (F25) without selection for mating time or callus colour, finding that the advanced hybrid population had moved towards B. tryoni phenotypes for both traits. Selection for day mating in replicate lines at F25 resulted in significant phenotypic shifts in both traits towards B. neohumeralis phenotypes in F26. However, we were unable to completely recover the mating time profile of B. neohumeralis and relaxation of selection for day mating led to a shift back towards dusk mating, but not yellow callus colour, by F35. CONCLUSION: We conclude that the inheritance of the two major species-defining traits is separable but tightly linked and involves more than one gene in each case. It also appears that laboratory conditions select for the B. tryoni phenotypes for mating time. We discuss our findings in relation to speciation theory and the likely effects of domestication during the generation of mass release strains for sterile insect control programmes.

Climate stress resistance in male Queensland fruit fly varies among populations of diverse geographic origins and changes during domestication.

BACKGROUND: The highly polyphagous Queensland fruit fly (Bactrocera tryoni Froggatt) expanded its range substantially during the twentieth century and is now the most economically important insect pest of Australian horticulture, prompting intensive efforts to develop a Sterile Insect Technique (SIT) control program. Using a "common garden" approach, we have screened for natural genetic variation in key environmental fitness traits among populations from across the geographic range of this species and monitored changes in those traits induced during domestication. RESULTS: Significant variation was detected between the populations for heat, desiccation and starvation resistance and wing length (as a measure of body size). Desiccation resistance was correlated with both starvation resistance and wing length. Bioassay data for three resampled populations indicate that much of the variation in desiccation resistance reflects persistent, inherited differences among the populations. No latitudinal cline was detected for any of the traits and only weak correlations were found with climatic variables for heat resistance and wing length. All three stress resistance phenotypes and wing length changed significantly in certain populations with ongoing domestication but there was also a strong population by domestication interaction effect for each trait. CONCLUSIONS: Ecotypic variation in heat, starvation and desiccation resistance was detected in Australian Qfly populations, and these stress resistances diminished rapidly during domestication. Our results indicate a need to select source populations for SIT strains which have relatively high climatic stress resistance and to minimise loss of that resistance during domestication.

Genome-wide patterns of differentiation over space and time in the Queensland fruit fly.

The Queensland fruit fly, Bactrocera tryoni, is a major pest of Australian horticulture which has expanded its range in association with the spread of horticulture over the last ~ 150 years. Its distribution in northern Australia overlaps that of another fruit fly pest to which some authors accord full species status, Bactrocera aquilonis. We have used reduced representation genome-wide sequencing to genotype 359 individuals taken from 35 populations from across the current range of the two taxa, plus a further 73 individuals from six of those populations collected 15-22 years earlier. We find significant population differentiation along an east-west transect across northern Australia which likely reflects limited but bidirectional gene flow between the two taxa. The southward expansion of B. tryoni has led to relatively little genetic differentiation, and most of it is associated with a move into previously marginal inland habitats. Two disjunct populations elsewhere in Australia and three on Melanesian islands are each clearly differentiated from all others, with data strongly supporting establishment from relatively few founders and significant isolation subsequently. Resequencing of historical samples from one of the disjunct Australian populations shows that its genetic profile has changed little over a 15-year period, while the Melanesian data suggest a succession of 'island hopping' events with progressive reductions in genetic diversity. We discuss our results in relation to the control of B. tryoni and as a model for understanding the genetics of invasion and hybridisation processes.

An RNAi supplemented diet as a reverse genetics tool to control bluegreen aphid, a major pest of legumes.

Aphids are important agricultural pests causing major yield losses worldwide. Since aphids can rapidly develop resistance to chemical insecticides there is an urgent need to find alternative aphid pest management strategies. Despite the economic importance of bluegreen aphid (Acyrthosiphon kondoi), very few genetic resources are available to expand our current understanding and help find viable control solutions. An artificial diet is a desirable non-invasive tool to enable the functional characterisation of genes in bluegreen aphid and discover candidate target genes for future use in RNA interference (RNAi) mediated crop protection against aphids. To date no artificial diet has been developed for bluegreen aphid, so we set out to develop a suitable diet by testing and optimising existing diets. Here, we describe an artificial diet for rearing bluegreen aphid and also provide a proof of concept for the supplementation of the diet with RNAi molecules targeting the salivary gland transcript C002 and gap gene hunchback, resulting in bluegreen aphid mortality which has not yet been documented in this species. Managing this pest, for example via RNAi delivery through artificial feeding will be a major improvement to test bluegreen aphid candidate target genes for future pest control and gain significant insights into bluegreen aphid gene function.

Detoxifying enzyme complements and host use phenotypes in 160 insect species.

We use the genomes of 160 insect species to test the hypothesis that the size of detoxifying enzyme families is greater in species using more chemically diverse food resources. Phylogenetically appropriate contrasts in subsamples of the data generally support the hypothesis. We find relatively high numbers of cytochrome P450, glutathione S-transferase and carboxyl/choline esterase genes in omnivores and herbivores feeding on chemically complex tissues and relatively low numbers of these genes in specialists on relatively simple diets, including plant sap, nectar and pollen, and blood. Among Lepidoptera feeding on green plant tissue and Condylognatha feeding on sap we also find more of these genes in highly polyphagous species, many of which are major agricultural pests. These genomic signatures of food resource use are consistent with the hypothesis that some taxa are preadapted for insecticide resistance evolution.

Additive and epistatic interactions between AKR and AIN loci conferring bluegreen aphid resistance and hypersensitivity in Medicago truncatula.

Aphids, including the bluegreen aphid (BGA; Acyrthosiphon kondoi), are important pests in agriculture. Two BGA resistance genes have been identified in the model legume Medicago truncatula, namely AKR (Acyrthosiphon kondoi resistance) and AIN (Acyrthosiphon induced necrosis). In this study, progeny derived from a cross between a resistant accession named Jester and a highly susceptible accession named A20 were used to study the interaction between the AKR and AIN loci with respect to BGA performance and plant response to BGA infestation. These studies demonstrated that AKR and AIN have additive effects on the BGA resistance phenotype. However, AKR exerts dominant suppression epistasis on AIN-controlled macroscopic necrotic lesions. Nevertheless, both AKR and AIN condition production of H2O2 at the BGA feeding site. Electrical penetration graph analysis demonstrated that AKR prevents phloem sap ingestion, irrespective of the presence of AIN. Similarly, the jasmonic acid defense signaling pathway is recruited by AKR, irrespective of AIN. This research identifies an enhancement of aphid resistance through gene stacking, and insights into the interaction of distinct resistance genes against insect pests.

Zucchini yellow mosaic virus Genomic Sequences from Papua New Guinea: Lack of Genetic Connectivity with Northern Australian or East Timorese Genomes, and New Recombination Findings.

Zucchini yellow mosaic virus (ZYMV) isolates were obtained in Papua New Guinea (PNG) from cucumber (Cucumis sativus) or pumpkin (Cucurbita spp.) plants showing mosaic symptoms growing at Kongop in the Mount Hagen District, Western Highlands Province, or Zage in the Goroka District, Eastern Highlands Province. The samples were blotted onto FTA cards, which were sent to Australia, where they were subjected to high-throughput sequencing. When the coding regions of the nine new ZYMV genomic sequences found were compared with those of 64 other ZYMV sequences from elsewhere, they grouped together, forming new minor phylogroup VII within ZYMV's major phylogroup A. Genetic connectivity was lacking between ZYMV genomic sequences from PNG and its neighboring countries, Australia and East Timor; the closest match between a PNG and any other genomic sequence was a 92.8% nucleotide identity with a sequence in major phylogroup A's minor phylogroup VI from Japan. When the RDP5.2 recombination analysis program was used to compare 66 ZYMV sequences, evidence was obtained of 30 firm recombination events involving 41 sequences, and all isolates from PNG were recombinants. There were 21 sequences without recombination events in major phylogroup A, whereas there were only 4 such sequences within major phylogroup B. ZYMV's P1, Cl, N1a-Pro, P3, CP, and NIb regions contained the highest evidence of recombination breakpoints. Following removal of recombinant sequences, seven minor phylogroups were absent (I, III, IV, V, VI, VII, and VIII), leaving only minor phylogroups II and IX. By contrast, when a phylogenetic tree was constructed using recombinant sequences with their recombinationally derived tracts removed before analysis, five previous minor phylogroups remained unchanged within major phylogroup A (II, III, IV, V, and VII) while four formed two new merged phylogroups (I/VI and VIII/IX). Absence of genetic connectivity between PNG, Australian, and East Timorese ZYMV sequences, and the 92.8% nucleotide identity between a PNG sequence and the closest sequence from elsewhere, suggest that a single introduction may have occurred followed by subsequent evolution to adapt to the PNG environment. The need for enhanced biosecurity measures to protect against potentially damaging virus movements crossing the seas separating neighboring countries in this region of the world is discussed.

Genetic Connectivity Between Papaya Ringspot Virus Genomes from Papua New Guinea and Northern Australia, and New Recombination Insights.

Isolates of papaya ringspot virus (PRSV) were obtained from plants of pumpkin (Cucurbita spp.) or cucumber (Cucumis sativus) showing mosaic symptoms growing at Zage in Goroka District in the Eastern Highland Province of Papua New Guinea (PNG) or Bagl in the Mount Hagen District, Western Highlands Province. The samples were sent to Australia on FTA cards where they were subjected to High Throughput Sequencing (HTS). When the coding regions of the six new PRSV genomic sequences obtained via HTS were compared with those of 54 other complete PRSV sequences from other parts of the world, all six grouped together with the 12 northern Australian sequences within major phylogroup B minor phylogroup I, the Australian sequences coming from three widely dispersed locations spanning the north of the continent. Notably, none of the PNG isolates grouped with genomic sequences from the nearby country of East Timor in phylogroup A. The closest genetic match between Australian and PNG sequences was a nucleotide (nt) sequence identity of 96.9%, whereas between PNG and East Timorese isolates it was only 83.1%. These phylogenetic and nt identity findings demonstrate genetic connectivity between PRSV populations from PNG and Australia. Recombination analysis of the 60 PRSV sequences available revealed evidence of 26 recombination events within 18 isolates, only four of which were within major phylogroup B and none of which were from PNG or Australia. Within the recombinant genomes, the P1, Cl, NIa-Pro, NIb, 6K2, and 5'UTR regions contained the highest numbers of recombination breakpoints. After removal of nonrecombinant sequences, four minor phylogroups were lost (IV, VII, VIII, XV), only one of which was in phylogroup B. When genome regions from which recombinationally derived tracts of sequence were removed from recombinants prior to alignment with nonrecombinant genomes, seven previous minor phylogroups within major phylogroup A, and two within major phylogroup B, merged either partially or entirely forming four merged minor phylogroups. The genetic connectivity between PNG and northern Australian isolates and absence of detectable recombination within either group suggests that PRSV isolates from East Timor, rather than PNG, might pose a biosecurity threat to northern Australian agriculture should they prove more virulent than those already present.

New Isolates of Sweet potato feathery mottle virus and Sweet potato virus C: Biological and Molecular Properties, and Recombination Analysis Based on Complete Genomes.

Sweet potato feathery mottle virus (SPFMV) and Sweet potato virus C (SPVC) isolates were obtained from sweetpotato shoot or tuberous root samples from three widely separated locations in Australia's tropical north (Cairns, Darwin, and Kununurra). The samples were planted in the glasshouse and scions obtained from the plants were graft inoculated to Ipomoea setosa plants. Virus symptoms were recorded in the field in Kununurra and in glasshouse-grown sweetpotato and I. setosa plants. RNA extracts from I. setosa leaf samples were subjected to high-throughput sequencing. New complete SPFMV (n = 17) and SPVC (n = 6) genomic sequences were obtained and compared with 47 sequences from GenBank. Phylogenetic analysis revealed that the 17 new SPFMV genomes all fitted within either major phylogroup A, minor phylogroup II, formerly O; or major phylogroup B, formerly RC. Major phylogroup A's minor phylogroup I, formerly EA, only appeared when recombinants were included. Numbers of SPVC genomes were insufficient to subdivide it into phylogroups. Within phylogroup A's minor phylogroup II, the closest genetic match between an Australian and a Southeast Asian SPFMV sequence was the 97.4% nucleotide identity with an East Timorese sequence. Recombination analysis of the 43 SPFMV and 27 SPVC sequences revealed evidence of 44 recombination events, 16 of which involved interspecies sequence transfers between SPFMV and SPVC and 28 intraspecies transfers, 17 in SPFMV and 11 in SPVC. Within SPFMV, 11 intraspecies recombination events were between different major phylogroups and 6 were between members of the same major phylogroup. Phylogenetic analysis accounting for the detected recombination events within SPFMV sequences yielded evidence of minor phylogroup II and phylogroup B but the five sequences from minor phylogroup I were distributed in two separate groups among the sequences of minor phylogroup II. For the SPVC sequences, phylogenetic analysis accounting for the detected recombination events revealed three major phylogroups (A, B, and C), with major phylogroup A being further subdivided into two minor phylogroups. Within the recombinant genomes of both viruses, their PI, NIa-Pro, NIb, and CP genes contained the highest numbers of recombination breakpoints. The high frequency of interspecies and interphylogroup recombination events reflects the widespread occurrence of mixed SPVC and SPFMV infections within sweetpotato plants. The prevalence of infection in northern Australian sweetpotato samples reinforces the need for improved virus testing in healthy sweetpotato stock programs. Furthermore, evidence of genetic connectivity between Australian and East Timorese SPFMV genomes emphasizes the need for improved biosecurity measures to protect against potentially damaging international virus movements.

First Complete Genome Sequence of Cucurbit aphid-borne yellows virus from Papua New Guinea.

Analysis of an RNA-Seq library from cucumber leaf RNA revealed the first complete genome sequence of Cucurbit aphid-borne yellows virus (CABYV) from Papua New Guinea. We compared it with 36 complete CABYV genomes from other world regions. It most resembled the genome of South Korean isolate GS6.

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.

Sweet potato feathery mottle virus and Sweet potato virus C from East Timorese and Australian Sweetpotato: Biological and Molecular Properties, and Biosecurity Implications.

Sweet potato feathery mottle virus (SPFMV) and Sweet potato virus C (SPVC) isolates from sweetpotato were studied to examine genetic connectivity between viruses from Australia and Southeast Asia. East Timorese samples from sweetpotato were sent to Australia on FTA cards. Shoot and tuberous root samples were collected in Australia and planted in the glasshouse, and scions were graft inoculated to Ipomoea setosa plants. Symptoms in infected sweetpotato and I. setosa plants were recorded. RNA extracts from FTA cards and I. setosa leaf samples were subjected to high-throughput sequencing (HTS). Complete genomic sequences (CS) of SPFMV and SPVC (11 each) were obtained by HTS, and coat protein (CP) genes from them were compared with others from GenBank. SPFMV sequences clustered into two major phylogroups (A and B = RC) and two minor phylogroups (EA[I] and O[II]) within A; East Timorese sequences were in EA(I) and O(II), whereas Australian sequences were in O(II) and B(RC). With SPVC, CP trees provided sufficient diversity to distinguish major phylogroups A and B and six minor phylogroups within A (I to VI); East Timorese sequences were in minor phylogroup I, whereas Australian sequences were in minor phylogroups II and VI and in major phylogroup B. With SPFMV, Aus13B grouped with East Timorese sequence TM64B within minor phylogroup O, giving nucleotide sequence identities of 97.4% (CS) and 98.3% (CP). However, the closest match with an Australian sequence was the 97.6% (CS) and 98.7% (CP) nucleotide identity between Aus13B and an Argentinian sequence. With SPVC, closest nucleotide identity matches between Australian and East Timorese sequences were 94.1% with Aus6a and TM68A (CS) and 96.3% with Aus55-4C and TM64A (CP); however neither pair member belonged to the same minor phylogroup. Also, the closest Australian match was 99.1% (CP) nucleotide identity between Aus4C and New Zealand isolate NZ4-4. These first complete genome sequences of SPFMV and SPVC from sweetpotato plantings in the Australian continent and neighboring Southeast Asia suggest at least two (SPFMV) and three (SPVC) separate introductions to Australia since agriculture commenced more than two centuries ago. These findings have major implications for both healthy stock programs and biosecurity management in relation to pathogen entry into Australia and elsewhere.

Two Complete Genome Sequences of Squash mosaic virus from 20-Year-Old Cucurbit Leaf Samples from Australia.

We present the first complete Australian Squash mosaic virus (SqMV) genome sequences. We compared the 2 Australian genomes from 20-year-old cucurbit samples with 8 other SqMV genomes. The Australian genomes shared >99.0% nucleotide identities, and their RNA1 and RNA2 sequences most closely resembled isolates Y and Kimbe from Japan, respectively.

First Complete Squash leaf curl China virus Genomic Segment DNA-A Sequence from East Timor.

We present here the first complete Squash leaf curl China virus (SLCCV) genomic segment DNA-A sequence from East Timor. It was isolated from a pumpkin plant. When compared with 15 complete SLCCV DNA-A genome sequences from other world regions, it most resembled the Malaysian isolate MC1 sequence.

Analysis of an RNA-seq Strand-Specific Library from an East Timorese Cucumber Sample Reveals a Complete Cucurbit aphid-borne yellows virus Genome.

Analysis of an RNA-seq library from cucumber leaf RNA extracted from a fast technology for analysis of nucleic acids (FTA) card revealed the first complete genome of Cucurbit aphid-borne yellows virus (CABYV) from East Timor. We compare it with 35 complete CABYV genomes from other world regions. It most resembled the genome of the South Korean isolate HD118.

Metagenomic Analysis of Cucumber RNA from East Timor Reveals an Aphid lethal paralysis virus Genome.

We present here the first complete genomic Aphid lethal paralysis virus (ALPV) sequence isolated from cucumber plant RNA from East Timor. We compare it with two complete ALPV genome sequences from China, and one each from Israel, South Africa, and the United States. It most closely resembled the Chinese isolate LGH genome.