Search

Corals’ tolerance to high temperature stress largely depends on their symbiotic microalgae (Symbiodiniaceae). However, the contributing microalgal traits and their relationships to one another are largely unclear. Here we compare the in vitro cellular profiles of seven Cladocopium C1acro microalgal strains (derived from the same ancestral strain) during a four-week exposure to 27°C or 31°C. One was an unselected wild-type strain (WT), three were selected at 31°C for nine years and shown to confer thermal tolerance on the coral host (SS+) and three others were similarly selected but did not confer tolerance (SS-). Flow cytometry was used to measure the staining intensities of cells treated with dyes representing the intracellular stress indicators reactive oxygen species (ROS), reduced glutathione (rGSH) and mitochondrial activity, as well as cell size/shape and photosynthetic pigments. Cell densities and photosynthetic efficiency (ϕPSII, Fv/Fm) were also measured. WT showed the highest levels of intracellular ROS and mitochondrial activity, lowest rGSH and largest cell sizes at both temperatures. SS+ strains had the lowest ROS and highest rGSH values and a unique pattern of correlations among parameters at 31°C. Our results support previous reports implicating the role of microalgal ROS, mitochondria and rGSH in holobiont thermal tolerance and suggest flow cytometry is a useful screening tool for identifying microalgal strains with enhanced thermal tolerance.

1AbstractThe assembly of reference-quality, chromosome-level genomes for both model and novel eukaryotic organisms is an increasingly achievable task for single research teams. However, the broad variety of sequencing technologies, assembly algorithms, and post-assembly processing tools currently available means that there is no clear consensus on a best-practice computational protocol for eukaryoticde novogenome assembly. An ever-increasing field of algorithms and packages with unique parameters, setup requirements, and environments makes it difficult for groups to pick up and test new tools, despite potential benefits. Here, we present a comprehensive Snakemake-based pipeline for eukaryotic genome assembly,Pyro, to further assist futurede novoassembly and benchmarking projects.Pyrocombines 20 assembly and eight polishing packages, comprising 30 different assembly approaches and up to 48 different polishing approaches in combination. These are available across Illumina short-read, Nanopore and PacBio CLR long-read technologies in one container, complete with data preparation, quality metric calculation and result reporting. We demonstratePyroeffectiveness by runningPyroon publicly available Illumina, Nanopore and PacBio CLR read sets forArabidopsis thaliana, producing 12 candidate assembly options with minimal initial input or configuration, each with extremely high contiguity and completeness.Pyrois highly customizable to expert needs, while also providing an accessible suggested set of tools for more casual users based on simple inputs.Pyrois available as a Singularity container suitable for execution on any Unix-compatible OS, and is freely available on GitHub (https://github.com/genomeassembler/pyro). This pipeline provides a one-stop solution for a variety ofde novoeukaryotic genome assembly needs, and will also assist in the assessment of new tools as a convenient benchmark-generating platform.

The assembly of reference-quality, chromosome-resolution genomes for both model and novel eukaryotic organisms is an increasingly achievable task for single research teams. However, the overwhelming abundance of sequencing technologies, assembly algorithms, and post-assembly processing tools currently available means that there is no clear consensus on a best-practice computational protocol for eukaryoticde novogenome assembly. Here, we provide a comprehensive benchmark of 28 state-of-the-art assembly and polishing packages, in various combinations, when assembling two eukaryotic genomes using both next-generation (Illumina HiSeq) and third-generation (Oxford Nanopore and PacBio CLR) sequencing data, at both controlled and open levels of sequencing coverage. Recommendations are made for the most effective tools for each sequencing technology and the best performing combinations of methods, evaluated against common assessment metrics such as contiguity, computational performance, gene completeness, and reference reconstruction, across both organisms and across sequencing coverage depth.

Corals’ tolerance to high temperature stress largely depends on their symbiotic microalgae (Symbiodiniaceae). However, the contributing microalgal traits are largely unclear. Here we compare thein vitrocellular profiles of sevenCladocopium C1acromicroalgal strains (derived from the same ancestral strain) during a four-week exposure to 27°C or 31°C. One was an unselected wild-type strain (WT), three were selected at 31°C for nine years and shown to confer thermal tolerance on the coral host (SS+) and three others were similarly selected but did not confer tolerance (SS-). Flow cytometry was used to measure the intracellular stress indicators reactive oxygen species (ROS), reduced glutathione (rGSH) and mitochondrial-membrane potential (ΔΨm), as well as cell size/shape and photosynthetic pigments. Cell densities and photosynthetic efficiency (ΦPSII, Fv/Fm) were also measured. WT showed the highest levels of intracellular ROS and ΔΨm, lowest rGSH and largest cell sizes at both temperatures. SS+ strains had the lowest ROS and highest rGSH values and a unique pattern of correlations among parameters at 31°C. Our results support previous reports implicating the role of microalgal ROS, ΔΨm and rGSH in holobiont thermal tolerance and suggest flow cytometry is a useful pre-screening tool for identifying microalgal strains with enhanced thermal tolerance.

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.

Incursions of the Queensland fruit fly Bactrocera tryoni (Qfly) into areas without permanent Qfly populations present serious threats to the Australian and New Zealand horticultural industries. Identifying the origins of recent incursions will help reduce future threats by enabling the targeting of problematic incursion routes for more stringent quarantine protocols. Here we present an analytical framework based on supervised and unsupervised machine learning to identify the origins and recent population history of incursion individuals. Our framework is based on a recently developed reference dataset of genome-wide markers for 35 Qfly populations from across the ranges of Qfly and the related taxon Bactrocera aquilonis (NTfly). We apply our framework to recent incursions into New Zealand, Tasmania and South Australia. Two distinct Qfly sources were identified for incursions into New Zealand (total 18 individuals), one from the east coast of Australia and one from New Caledonia. All eight recent incursion collections analysed (total 85 individuals) from South Australia and Tasmania most likely originated from just one of six clusters of populations in our reference database, Qfly from the east coast of Australia. None were found to originate from clusters containing NTfly or Qfly/NTfly hybrids in the Northern Territory or north Western Australia. Several, but not all, of the collections showed signals of small founding population size and two Tasmanian collections each included individuals apparently derived from three different sources within the east coast of Australia. In total, several more incursion events were detected than previously known, although some were founded by relatively few individuals.

Understanding the cumulative risk of chemical mixtures at environmentally realistic concentrations is a key challenge in honey bee ecotoxicology. Ecotoxicogenomics, including transcriptomics, measures responses in individual organisms at the molecular level which can provide insights into the mechanisms underlying phenotypic responses induced by one or more stressors and link impacts on individuals to populations. Here, fifth instar honey bee larvae were sampled from a previously reported field experiment exploring the phenotypic impacts of environmentally realistic chronic exposures of the pesticide imidacloprid (5 μg.kg

Reef-building corals are cnidarian animals (Anthozoa, Scleractinia) that mostly live in colonies composed of hundreds of thousands of tiny coral polyps. The polyps house dinoflagellate algal photosymbionts of the family Symbiodiniaceae inside their gastrodermal cells (Fig. 1 A ), and this mutualistic association builds the three-dimensional structure of the reefs via deposition of calcium carbonate skeletons. Scleractinian corals thus support a tremendous diversity and biomass of coral reef organisms. However, much of the biology of corals is poorly understood, partly due to a lack of advanced genetic tools for both the corals and the Symbiodiniaceae. In PNAS, Cleves et al. (1) catapult the field forward by describing an efficient protocol for applying the CRISPR technology to the coral Acropora millepora . They show that a CRISPR knockout of the gene encoding Heat Shock Transcription Factor 1 ( HSF1 ) in larvae of this coral reduces their heat shock tolerance, thus confirming a role for HSF1 in protecting corals from heat stress. Fig. 1. The scleractinian coral, Acropora millepora . ( A ) Close-up of individual polyps showing the endosymbiotic Symbiodiniaceae algae as brown pigmentation and as individual cells in the Inset . Symbiodiniaceae cells are ∼5 to 10 µm in diameter. Image credit: Katarina Damjanovic (photographer). ( B ) Bleached and unbleached adult colonies during a natural bleaching event on the Great Barrier Reef. Image credit: Ray Berkelmans (photographer). ( C ) Aposymbiotic larva. Image credit: Kate Quigley (photographer). Gaining a deeper understanding of the coral heat stress response is urgent given that climate warming, and particularly the resulting increased frequency and severity of summer heat waves, are responsible for the alarmingly rapid loss of coral reefs worldwide. Exposure to elevated temperature can result in the loss of the symbiotic algae from the coral … [↵][1]1To whom correspondence may be addressed. Email: madeleine.van{at}unimelb.edu.au. [1]: #xref-corresp-1-1

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.

Background Insights into the genetic capacities of species to adapt to future climate change can be gained by using comparative genomic and transcriptomic data to reconstruct the genetic changes associated with such adaptations in the past. Here we investigate the genetic changes associated with adaptation to arid environments, specifically climatic extremes and new cactus hosts, through such an analysis of five repleta group Drosophila species. Results We find disproportionately high rates of gene gains in internal branches in the species’ phylogeny where cactus use and subsequently cactus specialisation and high heat and desiccation tolerance evolved. The terminal branch leading to the most heat and desiccation resistant species, Drosophila aldrichi, also shows disproportionately high rates of both gene gains and positive selection. Several Gene Ontology terms related to metabolism were enriched in gene gain events in lineages where cactus use was evolving, while some regulatory and developmental genes were strongly selected in the Drosophila aldrichi branch. Transcriptomic analysis of flies subjected to sublethal heat shocks showed many more downregulation responses to the stress in a heat sensitive versus heat resistant species, confirming the existence of widespread regulatory as well as structural changes in the species’ differing adaptations. Gene Ontology terms related to metabolism were enriched in the differentially expressed genes in the resistant species while terms related to stress response were over-represented in the sensitive one. Conclusion Adaptations to new cactus hosts and hot desiccating environments were associated with periods of accelerated evolutionary change in diverse biochemistries. The hundreds of genes involved suggest adaptations of this sort would be difficult to achieve in the timeframes projected for anthropogenic climate change. Electronic supplementary material The online version of this article (10.1186/s12864-018-5413-3) contains supplementary material, which is available to authorized users.

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.

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high-quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome-wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert-adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate-continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics-restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large-scale and multifaceted genomic changes.

The cuticular layer of the insect exoskeleton contains diverse compounds that serve important biological functions, including the maintenance of homeostasis by protecting against water loss, protection from injury, pathogens and insecticides, and communication. Bactrocera tryoni (Froggatt) is the most destructive pest of fruit production in Australia, yet there are no published accounts of this species&rsquo, cuticular chemistry. We here provide a comprehensive description of B. tryoni cuticular chemistry. We used gas chromatography-mass spectrometry to identify and characterize compounds in hexane extracts of B. tryoni adults reared from larvae in naturally infested fruits. The compounds found included spiroacetals, aliphatic amides, saturated/unsaturated and methyl branched C12 to C20 chain esters and C29 to C33 normal and methyl-branched alkanes. The spiroacetals and esters were found to be specific to mature females, while the amides were found in both sexes. Normal and methyl-branched alkanes were qualitatively the same in all age and sex groups but some of the alkanes differed in amounts (as estimated from internal standard-normalized peak areas) between mature males and females, as well as between mature and immature flies. This study provides essential foundations for studies investigating the functions of cuticular chemistry in this economically important species.

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.

Coral bleaching tolerance can be enhanced through directed evolution of their microalgal symbionts., Coral reefs worldwide are suffering mass mortalities from marine heat waves. With the aim of enhancing coral bleaching tolerance, we evolved 10 clonal strains of a common coral microalgal endosymbiont at elevated temperatures (31°C) for 4 years in the laboratory. All 10 heat-evolved strains had expanded their thermal tolerance in vitro following laboratory evolution. After reintroduction into coral host larvae, 3 of the 10 heat-evolved endosymbionts also increased the holobionts’ bleaching tolerance. Although lower levels of secreted reactive oxygen species (ROS) accompanied thermal tolerance of the heat-evolved algae, reduced ROS secretion alone did not predict thermal tolerance in symbiosis. The more tolerant symbiosis exhibited additional higher constitutive expression of algal carbon fixation genes and coral heat tolerance genes. These findings demonstrate that coral stock with enhanced climate resilience can be developed through ex hospite laboratory evolution of their microalgal endosymbionts.

Burkholderia sp. strain SG-MS1 and Pseudomonas sp. strain SG-MS2 have previously been found to mineralize (+)-pinoresinol through a common catabolic pathway. Here, we used comparative genomics, proteomics, protein semipurification, and heterologous expression to identify a flavoprotein from the vanillyl alcohol oxidase/p-cresol methyl hydroxylase (VAO/PCMH) enzyme family in SG-MS2 that carries out the initial hydroxylation of (+)-pinoresinol at the benzylic carbon. The cognate gene is translationally coupled with a downstream cytochrome gene, and the cytochrome is required for activity. The flavoprotein has a unique combination of cofactor binding and cytochrome requirements for the VAO/PCMH family. The heterologously expressed enzyme has a K(m) of 1.17 μM for (+)-pinoresinol. The enzyme is overexpressed in strain SG-MS2 upon exposure to (+)-pinoresinol, along with 45 other proteins, 22 of which were found to be encoded by genes in an approximately 35.1-kb cluster also containing the flavoprotein and cytochrome genes. Homologs of 18 of these 22 genes, plus the flavoprotein and cytochrome genes, were also found in a 38.7-kb cluster in SG-MS1. The amino acid identities of four of the other proteins within the SG-MS2 cluster suggest they catalyze conversion of hydroxylated pinoresinol to protocatechuate and 2-methoxyhydroquinone. Nine other proteins upregulated in SG-MS2 on exposure to (+)-pinoresinol appear to be homologs of proteins known to comprise the protocatechuate and 2-methoxyhydroquinone catabolic pathways, but only three of the cognate genes lie within the cluster containing the flavoprotein and cytochrome genes. IMPORTANCE (+)-Pinoresinol is an important plant defense compound, a major food lignan for humans and some other animals, and the model compound used to study degradation of the β-β′ linkages in lignin. We report a gene cluster, in one strain each of Pseudomonas and Burkholderia, that is involved in the oxidative catabolism of (+)-pinoresinol. The flavoprotein component of the α-hydroxylase which heads the pathway belongs to the 4-phenol oxidizing (4PO) subgroup of the vanillyl alcohol oxidase/p-cresol methyl hydroxylase (VAO/PCMH) enzyme family but constitutes a novel combination of cofactor and electron acceptor properties for the family. It is translationally coupled with a cytochrome gene whose product is also required for activity. The work casts new light on the biology of (+)-pinoresinol and its transformation to other bioactive molecules. Potential applications of the findings include new options for deconstructing lignin into useful chemicals and the generation of new phytoestrogenic enterolactones from lignans.

Mycobacteria are major environmental microorganisms and cause many significant diseases, including tuberculosis. Mycobacteria make an unusual vitamin-like compound, F420, and use it to both persist during stress and resist antibiotic treatment. Understanding how mycobacteria make F420 is important, as this process can be targeted to create new drugs to combat infections like tuberculosis. In this study, we show that mycobacteria make F420 in a way that is different from other bacteria. We studied the molecular machinery that mycobacteria use to make F420, determining the chemical mechanism for this process and identifying a novel chemical intermediate. These findings also have clinical relevance, given that two new prodrugs for tuberculosis treatment are activated by F420., F420 is a low-potential redox cofactor used by diverse bacteria and archaea. In mycobacteria, this cofactor has multiple roles, including adaptation to redox stress, cell wall biosynthesis, and activation of the clinical antitubercular prodrugs pretomanid and delamanid. A recent biochemical study proposed a revised biosynthesis pathway for F420 in mycobacteria; it was suggested that phosphoenolpyruvate served as a metabolic precursor for this pathway, rather than 2-phospholactate as long proposed, but these findings were subsequently challenged. In this work, we combined metabolomic, genetic, and structural analyses to resolve these discrepancies and determine the basis of F420 biosynthesis in mycobacterial cells. We show that, in whole cells of Mycobacterium smegmatis, phosphoenolpyruvate rather than 2-phospholactate stimulates F420 biosynthesis. Analysis of F420 biosynthesis intermediates present in M. smegmatis cells harboring genetic deletions at each step of the biosynthetic pathway confirmed that phosphoenolpyruvate is then used to produce the novel precursor compound dehydro-F420-0. To determine the structural basis of dehydro-F420-0 production, we solved high-resolution crystal structures of the enzyme responsible (FbiA) in apo-, substrate-, and product-bound forms. These data show the essential role of a single divalent cation in coordinating the catalytic precomplex of this enzyme and demonstrate that dehydro-F420-0 synthesis occurs through a direct substrate transfer mechanism. Together, these findings resolve the biosynthetic pathway of F420 in mycobacteria and have significant implications for understanding the emergence of antitubercular prodrug resistance. IMPORTANCE Mycobacteria are major environmental microorganisms and cause many significant diseases, including tuberculosis. Mycobacteria make an unusual vitamin-like compound, F420, and use it to both persist during stress and resist antibiotic treatment. Understanding how mycobacteria make F420 is important, as this process can be targeted to create new drugs to combat infections like tuberculosis. In this study, we show that mycobacteria make F420 in a way that is different from other bacteria. We studied the molecular machinery that mycobacteria use to make F420, determining the chemical mechanism for this process and identifying a novel chemical intermediate. These findings also have clinical relevance, given that two new prodrugs for tuberculosis treatment are activated by F420.

Significance Helicoverpa armigera is a major agricultural and horticultural pest that recently spread from its historical distribution throughout much of the Old World to the Americas, where it is already causing hundreds of millions of dollars in damage every year. The species is notoriously quick to generate and disseminate pesticide resistance throughout its range and has a wider host range than the native Helicoverpa zea. Hybridization between the two species increases the opportunity for novel, agriculturally problematic ecotypes to emerge and spread through the Americas., Within the mega-pest lineage of heliothine moths are a number of polyphagous, highly mobile species for which the exchange of adaptive traits through hybridization would affect their properties as pests. The recent invasion of South America by one of the most significant agricultural pests, Helicoverpa armigera, raises concerns for the formation of novel combinations of adaptive genes following hybridization with the closely related Helicoverpa zea. To investigate the propensity for hybridization within the genus Helicoverpa, we carried out whole-genome resequencing of samples from six species, focusing in particular upon H. armigera population structure and its relationship with H. zea. We show that both H. armigera subspecies have greater genetic diversity and effective population sizes than do the other species. We find no signals for gene flow among the six species, other than between H. armigera and H. zea, with nine Brazilian individuals proving to be hybrids of those two species. Eight had largely H. armigera genomes with some introgressed DNA from H. zea scattered throughout. The ninth resembled an F1 hybrid but with stretches of homozygosity for each parental species that reflect previous hybridization. Regions homozygous for H. armigera-derived DNA in this individual included one containing a gustatory receptor and esterase genes previously associated with host range, while another encoded a cytochrome P450 that confers insecticide resistance. Our data point toward the emergence of novel hybrid ecotypes and highlight the importance of monitoring H. armigera genotypes as they spread through the Americas.

Summary Vegetable oils extracted from oilseeds are an important component of foods, but are also used in a range of high value oleochemical applications. Despite being biodegradable, nontoxic and renewable current plant oils suffer from the presence of residual polyunsaturated fatty acids that are prone to free radical formation that limit their oxidative stability, and consequently shelf life and functionality. Many decades of plant breeding have been successful in raising the oleic content to ~90%, but have come at the expense of overall field performance, including poor yields. Here, we engineer superhigh oleic (SHO) safflower producing a seed oil with 93% oleic generated from seed produced in multisite field trials spanning five generations. SHO safflower oil is the result of seed‐specific hairpin‐based RNA interference of two safflower lipid biosynthetic genes, FAD2.2 and FATB, producing seed oil containing less than 1.5% polyunsaturates and only 4% saturates but with no impact on lipid profiles of leaves and roots. Transgenic SHO events were compared to non‐GM safflower in multisite trial plots with a wide range of growing season conditions, which showed no evidence of impact on seed yield. The oxidative stability of the field‐grown SHO oil produced from various sites was 50 h at 110°C compared to 13 h for conventional ~80% oleic safflower oils. SHO safflower produces a uniquely stable vegetable oil across different field conditions that can provide the scale of production that is required for meeting the global demands for high stability oils in food and the oleochemical industry.

Carboxylesterase (CBE)-mediated metabolic resistance to organophosphate and carbamate insecticides is a major problem for the control of insect disease vectors, such as the mosquito. The most common mechanism involves overexpression of CBEs that bind to the insecticide with high affinity, thereby sequestering them before they can interact with their target. However, the absence of any structure for an insecticide-sequestering CBE limits our understanding of the molecular basis for this process. We present the first structure of a CBE involved in sequestration, Cqestβ21, from the mosquito disease vector Culex quinquefasciatus. Lysine methylation was used to obtain the crystal structure of Cqestβ21, which adopts a canonical α/β-hydrolase fold that has high similarity to the target of organophosphate and carbamate insecticides, acetylcholinesterase. Sequence similarity networks of the insect carboxyl/cholinesterase family demonstrate that CBEs associated with metabolic insecticide resistance across many spec...

The two tephritid fruit fly pests, Bactrocera tryoni and Bactrocera neohumeralis, are unusually well suited to the study of the genetics of reproductive isolating mechanisms. Sequence difference between the species is no greater than between a pair of conspecific Drosophila melanogaster populations. The two species exist in close sympatry, yet do not hybridize in the field, apparently kept separate by a strong premating isolation mechanism involving the time of day at which mating occurs. This spurred us to search for key genes for which time of day expression is regulated differently between the species. Using replicated, quantitative transcriptomes from head tissues of males of the two species, sampled in the day and night, we identified 141 transcripts whose abundance showed a significant interaction between species and time of day, indicating a difference in gene regulation. The brain transcripts showing this interaction were enriched for genes with a neurone function and 90% of these were more abundant at night than day in B. tryoni. Features of the expression patterns suggest that there may be a difference in the regulation of sleep-wake cycles between the species. In particular several genes, which in D. melanogaster are expressed in circadian pacemaker cells, are promising candidates to further explore the genetic differentiation involved in this prezygotic reproductive isolation mechanism.

Australia and Southeast Asia are hotspots of global diversity in the fruit-fly genus Bactrocera. Although a great diversity of species has been long recognised, evolutionary relationships are poorly understood, largely because previous sequencing techniques have provided insufficient historical signal for phylogenetic reconstruction. Poorly understood biogeographic history in Bactrocera has prevented a deeper understanding of migratory patterns in this economically important pest group. Using representatives from Australia and Malaysia, we tested the utility of a genome-reduction approach that generates thousands of single-nucleotide polymorphisms for phylogenetic reconstructions. This approach has high utility for species identification because of the ease of sample addition over time, and the species-level specificity able to be achieved with the markers. These data have provided a strongly supported phylogenetic tree congruent with topologies generated using more intensive sequencing approaches. In addition, our results do not support taxonomic assignments to species complex for a number of species, such as B. endiandrae in the dorsalis complex, yet find a close relationship between B. pallida and the dorsalis species. Our data have further validated non-monophyletic evolution of male response to primary attractants. We also showed at least two diversification events between Australia and Southeast Asia, indicating trans-regional dispersal in important pest species.

F420 is a low-potential redox cofactor that mediates the transformations of a wide range of complex organic compounds. Considered one of the rarest cofactors in biology, F420 is best known for its role in methanogenesis and has only been chemically identified in two phyla to date, the Euryarchaeota and Actinobacteria. In this work, we show that this cofactor is more widely distributed than previously reported. We detected the genes encoding all five known F420 biosynthesis enzymes (cofC, cofD, cofE, cofG and cofH) in at least 653 bacterial and 173 archaeal species, including members of the dominant soil phyla Proteobacteria, Chloroflexi and Firmicutes. Metagenome datamining validated that these genes were disproportionately abundant in aerated soils compared with other ecosystems. We confirmed through high-performance liquid chromatography analysis that aerobically grown stationary-phase cultures of three bacterial species, Paracoccus denitrificans, Oligotropha carboxidovorans and Thermomicrobium roseum, synthesized F420, with oligoglutamate sidechains of different lengths. To understand the evolution of F420 biosynthesis, we also analyzed the distribution, phylogeny and genetic organization of the cof genes. Our data suggest that although the Fo precursor to F420 originated in methanogens, F420 itself was first synthesized in an ancestral actinobacterium. F420 biosynthesis genes were then disseminated horizontally to archaea and other bacteria. Together, our findings suggest that the cofactor is more significant in aerobic bacterial metabolism and soil ecosystem composition than previously thought. The cofactor may confer several competitive advantages for aerobic soil bacteria by mediating their central metabolic processes and broadening the range of organic compounds they can synthesize, detoxify and mineralize.

Oligomerization has been suggested to be an important mechanism for increasing or maintaining the thermostability of proteins. Although it is evident that protein-protein contacts can result in substantial stabilization in many extant proteins, evidence for evolutionary selection for oligomerization is largely indirect and little is understood of the early steps in the evolution of oligomers. A laboratory-directed evolution experiment that selected for increased thermostability in the αE7 carboxylesterase from the Australian sheep blowfly, Lucilia cuprina, resulted in a thermostable variant, LcαE7-4a, that displayed increased levels of dimeric and tetrameric quaternary structure. A trade-off between activity and thermostability was made during the evolution of thermostability, with the higher-order oligomeric species displaying the greatest thermostability and lowest catalytic activity. Analysis of monomeric and dimeric LcαE7-4a crystal structures revealed that only one of the oligomerization-inducing mutations was located at a potential protein-protein interface. This work demonstrates that by imposing a selective pressure demanding greater thermostability, mutations can lead to increased oligomerization and stabilization, providing support for the hypothesis that oligomerization is a viable evolutionary strategy for protein stabilization.

The size of gene families associated with xenobiotic detoxification in insects may be associated with the complexity of their diets and their propensities to develop insecticide resistance. We test these hypotheses by collating the annotations of cytochrome P450, carboxyl/cholinesterase and glutathione S-transferase genes in 65 insect species with data on their host use and history of insecticide resistance. We find 2-4 fold variation across the species in the numbers of these genes and, in some orders, especially the Hymenoptera, there is a clear relationship between the numbers of genes and feeding preferences. However in other orders, in particular the Lepidoptera, no such relationship is apparent. The size of these three gene families also tend to correlate with insecticide resistance propensity but this may not be an independent effect because species with broader host ranges are more likely to be pests that are heavily sprayed with insecticides.

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted β-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.

We use annotated genomes of 14 Drosophila species covering diverse host use phenotypes to test whether 4 gene families that often have detoxification functions are associated with host shifts among species. Bark, slime flux, flower, and generalist necrotic fruit-feeding species all have similar numbers of carboxyl/cholinesterase, glutathione S-transferase, cytochrome P450, and UDP-glucuronosyltransferase genes. However, species feeding on toxic Morinda citrifolia fruit and the fresh fruit-feeding Drosophila suzukii have about 30 and 60 more, respectively. ABC transporters show a different pattern, with the flower-feeding D. elegans and the generalist necrotic fruit and cactus feeder D. hydei having about 20 and >100 more than the other species, respectively. Surprisingly, despite the complex secondary chemistry we find that 3 cactophilic specialists in the mojavensis species cluster have variably fewer genes than any of the other species across all 4 families. We also find 82 positive selection events across the 4 families, with the terminal D. suzukii and M. citrifolia-feeding D. sechellia branches again having the highest number of such events in proportion to their respective branch lengths. Many of the genes involved in these host-use-specific gene number differences or positive selection events lie in specific clades of the gene families that have been recurrently associated with detoxification. Several genes are also found to be involved in multiple duplication and/or positive selection events across the species studied regardless of their host use phenotypes; the most frequently involved are the ABC transporter CG1718, which is not in a specific clade associated with detoxification, and the α-esterase gene cluster, which is.

Over the last 40 years, many types of population genetic markers have been used to assess the population structure of the pest moth speciesHelicoverpa armigera. While this species is highly vagile, there is evidence of inter-continental population structure. Here, we examine Z-chromosome molecular markers within and between Chinese and Australian populations. Using 1352 polymorphic sites from 40 Z-linked loci, we compared two Chinese populations of moths separated by 700 km and found virtually no population structure (n= 41 andn= 54, with Dstatistic, suggesting that this species could have undergone a population expansion. Eight of the 40 loci had been examined in a previous study of Australian moths, of which six revealed very little inter-continental population structure. However, the two markers associated with theCyp303a1locus that has previously been proposed to be a target of a selective sweep, exhibited allele structuring between countries. Using a separate dataset of 19 Australian and four Chinese moths, we scanned the molecular variation distributed across the entire Z-chromosome and found distinct blocks of differentiation that include the region containingCyp303a1. We recommend some of these loci join those associated with insecticide resistance to form a set of genes best suited to analyzing population structure in this global pest.

Carbamate kinases catalyze the conversion of carbamate to carbamoyl phosphate, which is readily transformed into other compounds. Carbamate forms spontaneously from ammonia and carbon dioxide in aqueous solutions, so the kinases have potential for sequestrative utilization of the latter compounds. Here, we compare seven carbamate kinases from mesophilic, thermophilic, and hyperthermophilic sources. In addition to the known enzymes from Enterococcus faecalis and Pyrococcus furiosus , the previously unreported enzymes from the hyperthermophiles Thermococcus sibiricus and Thermococcus barophilus , the thermophiles Fervidobacterium nodosum and Thermosipho melanesiensis , and the mesophile Clostridium tetani were all expressed recombinantly, each in high yield. Only the clostridial enzyme did not show catalysis. In direct assays of carbamate kinase activity, the three hyperthermophilic enzymes display higher specific activities at elevated temperatures, greater stability, and remarkable substrate turnover at alkaline pH (9.9 to 11.4). Thermococcus barophilus and Thermococcus sibiricus carbamate kinases were found to be the most active when the enzymes were tested at 80°C, and maintained activity over broad temperature and pH ranges. These robust thermococcal enzymes therefore represent ideal candidates for biotechnological applications involving aqueous ammonia solutions, since nonbuffered 0.0001 to 1.0 M solutions have pH values of approximately 9.8 to 11.8. As proof of concept, here we also show that carbamoyl phosphate produced by the Thermococcus barophilus kinase is efficiently converted in situ to carbamoyl aspartate by aspartate transcarbamoylase from the same source organism. Using acetyl phosphate to simultaneously recycle the kinase cofactor ATP, at pH 9.9 carbamoyl aspartate is produced in high yield and directly from solutions of ammonia, carbon dioxide, and aspartate. IMPORTANCE Much of the nitrogen in animal wastes and used in fertilizers is commonly lost as ammonia in water runoff, from which it must be removed to prevent downstream pollution and evolution of nitrogenous greenhouse gases. Since carbamate kinases transform ammonia and carbon dioxide to carbamoyl phosphate via carbamate, and carbamoyl phosphate may be converted into other valuable compounds, the kinases provide a route for useful sequestration of ammonia, as well as of carbon dioxide, another greenhouse gas. At the same time, recycling the ammonia in chemical synthesis reduces the need for its energy-intensive production. However, robust catalysts are required for such biotransformations. Here we show that carbamate kinases from hyperthermophilic archaea display remarkable stability and high catalytic activity across broad ranges of pH and temperature, making them promising candidates for biotechnological applications. We also show that carbamoyl phosphate produced by the kinases may be efficiently used to produce carbamoyl aspartate.

Lignin is a complex aromatic polymer found in plant cell walls that makes up 15 to 40% of plant biomass. The degradation of lignin substructures by bacteria is of emerging interest because it could provide renewable alternative feedstocks and intermediates for chemical manufacturing industries. We have isolated a bacterium, strain SG61-1L, that rapidly degrades all of the stereoisomers of one lignin substructure, guaiacylglycerol-β-guaiacyl ether (GGE), which contains a key β- O -4 linkage found in most intermonomer linkages in lignin. In an effort to understand the rapid degradation of GGE by this bacterium, we heterologously expressed and kinetically characterized a suite of dehydrogenase candidates for the first known step of GGE degradation. We identified a clade of active GGE dehydrogenases and also several other dehydrogenases outside this clade that were all able to oxidize GGE. Several candidates exhibited stereoselectivity toward the GGE stereoisomers, while others had higher levels of catalytic performance than previously described GGE dehydrogenases for all four stereoisomers, indicating a variety of potential applications for these enzymes in the manufacture of lignin-derived commodities.

Control of Chinese Apolygus lucorum relies heavily on organophosphate insecticides. Here we describe resistance to the organophosphate chlorpyrifos in an A. lucorum strain, BZ-R, which was developed from a field-collected strain (BZ) by selection with chlorpyrifos in the laboratory. BZ-R showed 21-58 fold resistance to chlorpyrifos compared with the laboratory reference strain LSF and another susceptible strain, BZ-S, derived from BZ. BZ-R also showed several fold resistance to two other organophosphates and a carbamate. No synergism of chlorpyrifos by metabolic enzyme inhibitors nor any increase in detoxifying enzyme activities were observed in BZ-R. No sequence differences in acetylcholinesterase-2 were found to be associated with the resistance but the frequency of an alanine to serine substitution at position 216 of acetylcholinesterase-1 was 100% in BZ-R, ∼21-23% in SLF and BZ, and 0% in BZ-S. A single generation treatment of chlorpyrifos on the BZ strain also increased its frequency of the serine substitution to 64%. Recombinantly expressed acetylcholinesterase-1 carrying the serine substitution was about five fold less sensitive to inhibition by chlorpyrifos oxon than the wild-type enzyme. Quantitative real-time PCR found no differences in ace1 or ace2 expression levels among the strains tested. Thus the chlorpyrifos resistance is strongly associated with the serine substituted acetylcholinesterase-1. An equivalent substitution has been found to confer resistance to many organophosphate and carbamate insecticides in four other insect species.

Hexachlorocyclohexane (HCH), a synthetic organochloride, was first used as a broad-acre insecticide in the 1940s, and many HCH-degrading bacterial strains have been isolated from around the globe during the last 20 years. To date, the same degradation pathway (the lin pathway) has been implicated in all strains characterized, although the pathway has only been characterized intensively in two strains and for only a single HCH isomer. To further elucidate the evolution of the lin pathway, we have biochemically and genetically characterized three HCH-degrading strains from the Czech Republic and compared the genomes of these and seven other HCH-degrading bacterial strains. The three new strains each yielded a distinct set of metabolites during their degradation of HCH isomers. Variable assembly of the pathway is a common feature across the 10 genomes, eight of which (including all three Czech strains) were either missing key lin genes or containing duplicate copies of upstream lin genes (linA-F). The analysis also confirmed the important role of horizontal transfer mediated by insertion sequence IS6100 in the acquisition of the pathway, with a stronger association of IS6100 to the lin genes in the new strains. In one strain, a linA variant was identified that likely caused a novel degradation phenotype involving a shift in isomer preference. This study identifies a number of strains that are in the early stages of lin pathway acquisition and shows that the state of the pathway can explain the degradation patterns observed.

Three Australian tephritid fruit flies (Bactrocera tryoni - Q-fly, Bactrocera neohumeralis - NEO, and Bactrocera jarvisi - JAR) are promising models for genetic studies of pest status and invasiveness. The long history of ecological and physi- ological studies of the three species has been augmented by the development of a range of genetic and genomic tools, including the capacity for forced multigeneration crosses between the three species followed by selection experiments, a draft genome for Q-fly, and tissue- and stage-specific transcriptomes. The Q-fly and NEO species pair is of particular interest. The distribution of NEO is contained entirely within the wider distribution of Q-fly and the two species are ecologically extremely similar, with no known differences in pheromones, temperature tolerance, or host-fruit utilisation. However there are three clear differences be- tween them: humeral callus colour, complete pre-mating isolation based on mating time-of-day, and invasiveness. NEO is much less invasive, whereas in historical times Q-fly has invaded southeastern Australia and areas of Western Australia and the North- ern Territory. In southeastern fruit-growing regions, microsatellites suggest that some of these outbreaks might derive from ge- netically differentiated populations overwintering in or near the invaded area. Q-fly and NEO show very limited genome diffe- rentiation, so comparative genomic analyses and QTL mapping should be able to identify the regions of the genome controlling mating time and invasiveness, to assess the genetic bases for the invasive strains of Q-fly, and to facilitate a variety of improve- ments to current sterile insect control strategies for that species (Current Zoology 61 (3): 477-487, 2015). Keywords Bactrocera, Genomics, STERILE-insect-technique, Interspecific-hybridisation