Search

Members of the fungal genusArmillariaare necrotrophic pathogens with efficient plant biomass-degrading strategies. The genus includes some of the largest terrestrial organisms on Earth, spreading underground and causing tremendous losses in diverse ecosystems. Despite their global importance, the mechanism by whichArmillariaevolved pathogenicity in a clade of dominantly non-pathogenic wood-degraders (Agaricales) remains elusive. Here, using new genomic data, we show thatArmillariaspecies, in addition to widespread gene duplications andde novogene origins, appear to have at least 775 genes that were acquired via 101 horizontal gene transfer (HGT) events, primarily from Ascomycota. Functional and expression data suggest that HGT might have affected plant biomass-degrading and virulence abilities ofArmillaria, two pivotal traits in their lifestyle. We further assayed gene expression during root and cambium colonization, and report putative virulence factors, extensive regulation of horizontally acquired and wood-decay related genes as well as novel pathogenicity-induced small secreted proteins (PiSSPs). Two PiSSPs induced necrosis in live plants, suggesting they are potential virulence effectors conserved acrossArmillaria. Overall, this study details how evolution knitted together horizontally and vertically inherited genes in complex adaptive traits, such as plant biomass degradation and pathogenicity, paving the way for development of infection models for one of the most influential pathogens of temperate forest ecosystems.

Irradiated Queensland fruit flies (Bactrocera tryoni) used in Sterile Insect Technique (SIT) programmes are marked with fluorescent dyes to distinguish them from wild flies when recaptured in monitoring traps. However, coating sterile pupae with powdered dyes can reduce emergence rates and fly quality and can sometimes produce insufficiently certain discrimination through inadequate coating or because the dye is transferred to wild flies through contact. Here we created a phenotypically distinct B. tryoni strain that lacks typical melanisation patterns through CRISPR/Cas9-mediated mutagenesis of tandemly duplicated yellow-y genes and then assessed effects of this visible trait on fly performance. Recessive mutations are only required in one of these copies to restrict melanisation and generate a phenotype clearly distinguished from wild type. The yellow strain showed significant declines in eclosion rates and in the percentage of fliers directly after emergence. Locomotor activity was greater in the yellow strain, and these mutations did not generally affect mating probability, copula latency, or copula duration. The longevity of yellow flies was approximately 10 days shorter than wild-type flies in both sexes. Overall, replacing dyes with yellow body marker for SIT can simplify production, eliminate a step that is known to reduce fly quality, remove potentially hazardous dyes from production, enable accurate discrimination from wild flies, and improve cost-effectiveness; however, direct comparisons of the decrements in performance associated with dyes on mass-reared wild-type flies and disruption of yellow-y genes are now required to determine the relative suitability of these marking methods for B. tryoni SIT.

Invasive Tephritid fruit flies are a global threat to both agriculture and horticulture industries. Biosecurity has played a critical role in reducing their damage but becomes more and more challenging after several key chemical pesticides were banned or withdrawn for health or environmental reasons. This has led to non-chemical approaches including heat and cold treatments being broadly utilized to get rid of fruit fly infestation. However, the molecular mechanisms to kill the flies underlying these stressors are not clear yet. This knowledge will certainly help refine current post-harvest treatment strategies and develop more efficient, cost-effective and environmentally friendly approaches for fruit fly management. Previously, the molecular response of the Mediterranean fruit fly,Ceratitis capitata(Wiedemann) to heat was examined thoroughly, in which 31 key genes were identified with significant changes in expression levels and their high-resolution expression timeline was constructed across 11 timepoints. However, whether these candidate genes respond to cold in the same way was unknown. Here, a temperature bioassay was conducted and the expression profiles of these genes were investigated across the same 11 timepoints using cold treatment. The results showed that most of candidate genes exhibited divergent expression profiles compared to heat treatment, suggesting that the fly molecular response to cold may be different from those to heat. This study provides new knowledge of Tephritid fruit fly response to cold at a molecular level, which could aid in improving current fruit fly management and facilitate the development of new strategies to control this serious horticultural insect pest.

Tephritid fruit flies are highly successful invaders and some—such as the Mediterranean fruit fly, Ceratitis capitata (Wiedemann)—are able to adapt to a large range of crops. Biosecurity controls require that shipments of produce are ensured to be pest-free, which is increasingly difficult due to the ban of key pesticides. Instead, stress-based strategies including controlled atmosphere, temperature, and irradiation can be used to eradicate flies inside products. However, unlike pesticide science, we do not yet have a robust scientific approach to measure cost-effectively whether a sufficiently lethal stress has been delivered and understand what this stress does to the biology of the pest. The latter is crucial as it would enable a combination of stresses targeting multiple molecular pathways and thus allow for lower doses of each to achieve higher lethality and reduce the development of resistance. Using heat as an example, this is the first study investigating the molecular stress response to heat in Tephritidae. Using a novel setup delivering measured doses of heat on C. capitata larvae and a high-density 11 timepoint gene expression experiment, we identified key components of lethal heat-stress response. While unraveling the complete molecular mechanism of fruit fly response to lethal stress would be a long-term project, this work curates and develops 31 potential biomarkers to assess whether sufficient lethal stress has been delivered. Further, as these protocols are straightforward and less expensive than other—omic approaches, our studies and approach will assist other researchers working on stress response.

Background The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly’s biology and proposing alternative control methods to pesticide use. Results Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. Conclusions The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome’s organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

Multipartite viral vectors provide a simple, inexpensive and effective biotechnological tool to transiently manipulate (i.e. reduce or increase) gene expression in planta and characterise the function of genetic traits. The development of virus-induced gene regulation (VIGR) systems usually involve the targeted silencing or overexpression of genes involved in pigment biosynthesis or degradation in plastids, thereby providing rapid visual assessment of success in establishing RNA- or DNA-based VIGR systems in planta. Carotenoids pigments provide plant tissues with an array of yellow, orange, and pinkish-red colours. VIGR-induced transient manipulation of carotenoid-related gene expression has advanced our understanding of carotenoid biosynthesis, regulation, accumulation and degradation, as well as plastid signalling processes. In this review, we describe mechanisms of VIGR, the importance of carotenoids as visual markers of technology development, and knowledge gained through manipulating carotenogenesis in model plants as well as horticultural crops not always amenable to transgenic approaches. We outline how VIGR can be utilised in plants to fast-track the characterisation of gene function(s), accelerate fruit tree breeding programs, edit genomes, and biofortify plant products enriched in carotenoid micronutrients for horticultural innovation.

Three subtypes of C4 photosynthesis exist (NADP-ME, NAD-ME and PEPCK), each known to be beneficial under specific environmental conditions. However, the influence of photosynthetic subtype on transcriptomic plasticity, as well as the genes underpinning this variability, remain largely unknown. Here, we comprehensively investigate the responses of six C4 grass species, spanning all three C4 subtypes, to two controlled environmental stresses: low light (200 µmol m-2 sec-1 ) and glacial CO2 (subambient; 180 ppm). We identify a susceptibility within NADP-ME species to glacial CO2 . Notably, although glacial CO2 phenotypes could be tied to C4 subtype, biochemical and transcriptomic responses to glacial CO2 were largely species specific. Nevertheless, we were able to identify subtype specific subsets of significantly differentially expressed transcripts which link resource acquisition and allocation to NADP-ME species susceptibility to glacial CO2 . Here, low light phenotypes were comparable across species with no clear subtype response, while again, transcriptomic responses to low light were largely species specific. However, numerous functional similarities were noted within the transcriptomic responses to low light, suggesting these responses are functionally relatively conserved. Additionally, PEPCK species exhibited heightened regulation of transcripts related to metabolism in response to both stresses, likely tied to their C4 metabolic pathway. These results highlight the influence that both species and subtype can have on plant responses to abiotic stress, building on our mechanistic understanding of acclimation within C4 grasses and highlighting avenues for future crop improvements.

This review examines how next-generation nucleic acid sequencing (NGS) is being used in the diagnosis of invasive fungal disease (IFD) and how well its implementation compares to PCR-based diagnosis. The comparison allows for the assessment of the advantages and potential limitations of NGS in the diagnosis of IFD. NGS is established as a method for sequencing of microbial genomes. It is increasingly being tested as a direct diagnostic method from various sample types including samples such as plasma, formalin fixed and paraffin embedded specimens, and bronchoalveolar lavage. In these cases, the methodology was generally able to identify the causative agents of IFD. NGS-based methods will play an increasingly important role in diagnosis of IFD. The major limitation is currently cost, the need to standardise methods of nucleic acid isolation, and sequence analysis to enable broad uptake and application of the method.

Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.

Sea anemones have a wide array of toxic compounds (peptide toxins found in their venom) which have potential uses as therapeutics. To date, the majority of studies characterizing toxins in sea anemones have been restricted to species from the superfamily, Actinioidea. No highly complete draft genomes are currently available for this superfamily, however, highlighting our limited understanding of the genes encoding toxins in this important group. Here we have sequenced, assembled, and annotated a draft genome for Actinia tenebrosa. The genome is estimated to be approximately 255 megabases, with 31,556 protein‐coding genes. Quality metrics revealed that this draft genome matches the quality and completeness of other model cnidarian genomes, including Nematostella, Hydra, and Acropora. Phylogenomic analyses revealed strong conservation of the Cnidaria and Hexacorallia core‐gene set. However, we found that lineage‐specific gene families have undergone significant expansion events compared with shared gene families. Enrichment analysis performed for both gene ontologies, and protein domains revealed that genes encoding toxins contribute to a significant proportion of the lineage‐specific genes and gene families. The results make clear that the draft genome of A. tenebrosa will provide insight into the evolution of toxins and lineage‐specific genes, and provide an important resource for the discovery of novel biological compounds., To date, the majority of studies characterizing toxins in sea anemones have been restricted to species from the Actinioidea superfamily; however, no draft genomes are currently available for this superfamily. Here we have sequenced, assembled, and annotated the first Actinioidean draft genome for Actinia tenebrosa. The results make clear that the draft genome of A. tenebrosa will provide insight into the evolution of toxins, species‐specific genes, the cnidarian, and hexacorallian core‐gene set, and provide an important resource for novel biological compounds.

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.

BackgroundProlonged mechanical stress (MS) causes thigmomorphogenesis, a stress acclimation response associated with increased disease resistance. What remains unclear is if; 1) plants pre-exposed to a short period of repetitive MS can prime defence responses upon subsequent challenge with necrotrophic pathogens, 2) MS mediates plant immunity via jasmonic acid (JA) signalling, and 3) a short period of repetitive MS can cause long-term changes in gene expression resembling a stress-induced memory. To address these points, 10-days old juvenileArabidopsisseedlings were mechanically stressed for 7-days using a soft brush and subsequently challenged with the necrotrophic pathogens,Alternaria brassicicola,andBotrytis cinerea. Here we assessed how MS impacted structural cell wall appositions, disease symptoms and altered gene expression in response to infection.ResultsThe MS-treated plants exhibited enhanced cell wall appositions and jasmonic acid (JA) accumulation that correlated with a reduction in disease progression compared to unstressed plants. The expression of genes involved in JA signalling, callose deposition, peroxidase and phytoalexin biosynthesis and reactive oxygen species detoxification were hyper-induced 4-days post-infection in MS-treated plants. The loss-of-function in JA signalling mediated by the JA-insensitivecoronatine-insensitive 1(coi1) mutant impaired the hyper-induction of defense gene expression and promoted pathogen proliferation in MS-treated plants subject to infection. The basal expression level ofPATHOGENESIS-RELATED GENE 1andPLANT DEFENSIN 1.2defense marker genes were constitutively upregulated in rosette leaves for 5-days post-MS, as well as in naïve cauline leaves that differentiated from the inflorescence meristem well after ceasing MS.ConclusionThis study reveals that exposure of juvenileArabidopsisplants to a short repetitive period of MS can alter gene expression and prime plant resistance upon subsequent challenge with necrotrophic pathogens via the JA-mediatedCOI1signalling pathway. MS may facilitate a stress-induced memory to modulate the plant’s response to future stress encounters. These data advance our understanding of how MS primes plant immunity against necrotrophic pathogens and how that could be utilised in sustainable agricultural practices.

Sensory neuron membrane proteins (SNMPs) play a critical role in insect chemosensory system. Previously, three SNMPs were identified, characterized and functionally investigated in a lepidopteran model insect, Bombyx mori. However, whether these results are consistent across other lepidopteran species are unknown. Here genome and transcriptome data analysis, expression profiling, quantitative real-time PCR (qRT-PCR) and the yeast hybridization system were utilized to examine snmp genes of Helicoverpa armigera, one of the most destructive lepidopteran pests in cropping areas. In silico expression and qRT-PCR analyses showed that, just as the B. mori snmp genes, H. armigera snmp1 (Harmsnmp1) is specifically expressed in adult antennae. Harmsnmp2 is broadly expressed in multiple tissues including adult antennae, tarsi, larval antennae and mouthparts. Harmsnmp3 is specifically expressed in larval midguts. Further RNAseq analysis suggested that the expression levels of Harmsnmp2 and Harmsnmp3 differed significantly depending on the plant species on which the larvae fed, indicating they may be involved in plant-feeding behaviours. Yeast hybridization results revealed a protein-protein interaction between HarmSNMP1 and the sex pheromone receptor, HarmOR13. This study demonstrated that SNMPs may share same functions and mechanisms in different lepidopteran species, which improved our understanding of insect snmp genes and their functions in lepidopterans.

Expanding knowledge of the C4 photosynthetic pathway can provide key information to aid biological improvements to crop photosynthesis and yield. While the C4 NADP-ME pathway is well characterised, there is increasing agricultural and bioengineering interest in the comparably understudied NAD-ME and PEPCK pathways. Within this study, a systematic identification of key differences across species has allowed us to investigate the evolution of C4-recruited genes in one C3 and eleven C4 grasses (Poaceae) spanning two independent origins of C4 photosynthesis. We present evidence for C4-specific paralogs of NAD-malic enzyme 2, MPC1 and MPC2 (mitochondrial pyruvate carriers) via increased transcript abundance and associated rates of evolution, implicating them as genes recruited to perform C4 photosynthesis within NAD-ME and PEPCK subtypes. We then investigate the localisation of AspAT across subtypes, using novel and published evidence to place AspAT3 in both the cytosol and peroxisome. Finally, these findings are integrated with transcript abundance of previously identified C4 genes to provide an updated model for C4 grass NAD-ME and PEPCK photosynthesis. This updated model allows us to develop on the current understanding of NAD-ME and PEPCK photosynthesis in grasses, bolstering our efforts to understand the evolutionary 'path to C4' and improve C4 photosynthesis.

Adaptation to human-induced environmental change has the potential to profoundly influence the genomic architecture of affected species. This is particularly true in agricultural ecosystems, where anthropogenic selection pressure is strong. Heliothis virescens primarily feeds on cotton in its larval stages and US populations have been declining since the widespread planting of transgenic cotton, which endogenously expresses proteins derived from Bacillus thuringiensis (Bt). No physiological adaptation to Bt toxin has been found in the field, so adaptation in this altered environment could involve: 1) shifts in host plant selection mechanisms to avoid cotton, 2) changes in detoxification mechanisms required for cotton-feeding versus feeding on other hosts, or 3) loss of resistance to previously used management practices including insecticides. Here we begin to address whether such changes occurred in H. virescens populations between 1997-2012, as Bt cotton cultivation spread through the agricultural landscape. For our study, we produced an H. virescens genome assembly and used this in concert with a ddRAD-seq enabled genome scan to identify loci with significant allele frequency changes over the 15 year period. Genetic changes at a previously described H. virescens insecticide target of selection were detectable in our genome scan, and increased our confidence in this methodology. Additional loci were also detected as being under selection, and we quantified the selection strength required to elicit observed allele frequency changes at each locus. Potential contributions of genes near loci under selection to adaptive phenotypes in the H. virescens cotton system are discussed. This article is protected by copyright. All rights reserved.

Sensory neuron membrane proteins (SNMPs) play a critical role in the insect olfactory system but there is a deficit of functional studies beyond Drosophila. Here, we use a combination of available genome sequences, manual curation, genome and transcriptome data, phylogenetics, expression profiling and gene knockdown to investigate SNMP superfamily in various insect species with a focus on Lepidoptera. We curated 81 genes from 36 insect species and identified a novel lepidopteran SNMP gene family, SNMP3. Phylogenetic analysis shows that lepidopteran SNMP3, but not the previously annotated lepidopteran SNMP2, is the true homologue of the dipteran SNMP2. Digital expression, microarray and qPCR analyses show that the lepidopteran SNMP1 is specifically expressed in adult antennae. SNMP2 is widely expressed in multiple tissues while SNMP3 is specifically expressed in the larval midgut. Microarray analysis suggest SNMP3 may be involved in the silkworm immunity response to virus and bacterial infections. We functionally characterized SNMP1 in the silkworm using RNA interference (RNAi) and behavioral assays. Our results suggested that Bombyx mori SNMP1 is a functional orthologue of the Drosophila melanogaster SNMP1 and plays a critical role in pheromone detection. Split-ubiquitin yeast hybridization study shows that BmorSNMP1 has a protein-protein interaction with the pheromone receptor (BmorOR1), and the co-receptor (BmorOrco). Concluding, we propose a novel molecular model in which BmorOrco, BmorSNMP1 and BmorOR1 form a heteromer in the detection of the silkworm sex pheromone bombykol.

Lecanicillium psalliotae is an entomopathogenic, mycoparasitical, and nematophagous fungus known to produce antibiotic and antifungal compounds. Here, we report the first 36-Mb draft genome sequence of L. psalliotae strain HWLR35. The draft genome contains 197 scaffolds and is predicted to have 11,009 protein-coding genes.

Identifying the genomic changes that control morphological variation and understanding how they generate diversity is a major goal of evolutionary biology. In Heliconius butterflies, a small number of genes control the development of diverse wing colour patterns. Here, we used full-genome sequencing of individuals across the Heliconius erato radiation and closely related species to characterize genomic variation associated with wing pattern diversity. We show that variation around colour pattern genes is highly modular, with narrow genomic intervals associated with specific differences in colour and pattern. This modular architecture explains the diversity of colour patterns and provides a flexible mechanism for rapid morphological diversification. © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved.

Adaptation to human-induced environmental change has the potential to profoundly influence the genomic architecture of affected species. This is particularly true in agricultural ecosystems, where anthropogenic selection pressure is strong. Heliothis virescens feeds on cotton in its larval stages and US populations have been declining since the widespread planting of transgenic cotton, which endogenously express proteins derived Bacillus thuringiensis (Bt). No physiological adaptation to Bt toxin has been found in the field, so adaptation to this altered environment could involve: 1) shifts in host plant selection mechanisms to avoid cotton, 2) changes in detoxification mechanisms required for cotton-feeding versus feeding on other host plants, or 3) loss of resistance to previously used management practices including insecticides. Here we begin to address the question of whether such changes occurred in H. virescens populations between the years 1997 and 2012. As part of our study, we produced an H. virescens genome assembly and used this in concert with a ddRAD-seq enabled genome scan to identify loci with significant allele frequency changes over the 15 year period after Bt crops became widespread in the agricultural landscape. Confirmation that a known genetic change at a previously described H. virescens target of selection was detectable in our genome scan increased our confidence in this methodology. We quantified the strength of selection required to elicit the observed allele frequency changes at loci under selection. Potential contributions of genes near the loci under selection to adaptive phenotypes in the H. virescens cotton system are discussed.

Motivation: Bioinformatics tools, such as assemblers and aligners, are expected to produce more accurate results when given better quality sequence data as their starting point. This expectation has led to the development of stand-alone tools whose sole purpose is to detect and remove sequencing errors. A good error-correcting tool would be a transparent component in a bioinformatics pipeline, simply taking sequence data in any of the standard formats and producing a higher quality version of the same data containing far fewer errors. It should not only be able to correct all of the types of errors found in real sequence data (substitutions, insertions, deletions and uncalled bases), but it has to be both fast enough and scalable enough to be usable on the large datasets being produced by current sequencing technologies, and work on data derived from both haploid and diploid organisms. Results: This article presents Blue, an error-correction algorithm based on k-mer consensus and context. Blue can correct substitution, deletion and insertion errors, as well as uncalled bases. It accepts both FASTQ and FASTA formats, and corrects quality scores for corrected bases. Blue also maintains the pairing of reads, both within a file and between pairs of files, making it compatible with downstream tools that depend on read pairing. Blue is memory efficient, scalable and faster than other published tools, and usable on large sequencing datasets. On the tests undertaken, Blue also proved to be generally more accurate than other published algorithms, resulting in more accurately aligned reads and the assembly of longer contigs containing fewer errors. One significant feature of Blue is that its k-mer consensus table does not have to be derived from the set of reads being corrected. This decoupling makes it possible to correct one dataset, such as small set of 454 mate-pair reads, with the consensus derived from another dataset, such as Illumina reads derived from the same DNA sample. Such cross-correction can greatly improve the quality of small (and expensive) sets of long reads, leading to even better assemblies and higher quality finished genomes. Availability and implementation: The code for Blue and its related tools are available from http://www.bioinformatics.csiro.au/Blue . These programs are written in C# and run natively under Windows and under Mono on Linux. Contact: paul.greenfield@csiro.au Supplementary information: Supplementary data are available at Bioinformatics online.

The silks of arthropods have an elementary role in the natural history of the organisms that spin them, yet they are coded by rapidly evolving genes leading some authors to speculate that silk proteins are non-homologous proteins co-opted multiple times independently for similar functions. However, some general structural patterns are emerging. In this work we identified three major silk gland proteins using a combined biochemical, proteomic, next-generation sequencing and bioinformatic approach. Biochemical characterization determined that they were phosphorylated with multiple isoforms and potentially differential phosphorylation. Structural characterization showed that their structure was more similar to silk proteins from distantly related aquatic Trichopteran species than more closely related terrestrial or aquatic Diptera. Overall, our approach is easily transferable to any non-model species and if used across a larger number of aquatic species, we will be able to better understand the processes involved in linking the secondary structure of silk proteins with their function between in an organisms and its habitat.

The Insect taste system plays a central role in feeding behaviours and co-evolution of insect-host interactions. Gustatory receptors form the interface between the insect taste system and the environment. From genome and transcriptome sequencing we identified 197 novel gustatory receptor (GR) genes from the polyphagous pest Helicoverpa armigera. These GRs include a significantly expanded bitter receptor family (180 GRs) that could be further divided into three categories based on polypeptide lengths, gene structure and amino acid sequence. Type 1 includes 29 bitter Gr genes that possess introns. Type 2 includes 13 long intronless bitter Gr genes, while Type 3 comprises 131 short intronless bitter Gr genes. Calcium imaging analysis demonstrated that three Type 3 GRs (HarmGR35, HarmGR50 and HarmGR195) can be activated by a crude extract of cotton leaves. HarmGR195, a GR specifically and selectively expressed in adult tarsi, showed a specific response to proline, an amino acid widely present in plant tissues. We hypothesise that the expansion in the H. armigera GR family may be functionally tied to its polyphagous behavior. Understanding the molecular basis of polyphagy may provide opportunities for the development of new environmentally friendly pest control strategies.

Interleukins are a group of cytokines with complex immunomodulatory functions that are important for regulating immunity in vertebrate species. Reptiles and mammals last shared a common ancestor more than 350 million years ago, so it is not surprising that low sequence identity has prevented divergent interleukin genes from being identified in the central bearded dragon lizard, Pogona vitticeps, in its genome assembly. To determine the complete nucleotide sequences of key interleukin genes, we constructed full-length transcripts, using the Trinity platform, from short paired-end read RNA sequences from stimulated spleen cells. De novo transcript reconstruction and analysis allowed us to identify interleukin genes that are missing from the published P. vitticeps assembly. Identification of key cytokines in P. vitticeps will provide insight into the essential molecular mechanisms and evolution of interleukin gene families and allow for characterization of the immune response in a lizard for comparison with mammals.

Background: Oysters have important ecological functions in their natural environment, acting as global carbon sinks and improving water quality by removing excess nutrients from the water column. During their life-time oysters are exposed to a variety of pathogens that can cause severe mortality in a range of oyster species. Environmental stressors encountered in their habitat can increase the susceptibility of oysters to these pathogens and in general have been shown to impact on oyster immunity, making immune parameters expressed in these marine animals an important research topic. Results: Paired-end Illumina high throughput sequencing of six S. glomerata tissues exposed to different environmental stressors resulted in a total of 484,121,702 paired-end reads. When reads and assembled transcripts were compared to the C. gigas genome, an overall low level of similarity at the nucleotide level, but a relatively high similarity at the protein level was observed. Examination of the tissue expression pattern showed that some transcripts coding for cathepsins, heat shock proteins and antioxidant proteins were exclusively expressed in the haemolymph of S. glomerata, suggesting a role in innate immunity. Furthermore, analysis of the S. glomerata ORFs showed a wide range of genes potentially involved in innate immunity, from pattern recognition receptors, components of the Toll-like signalling and apoptosis pathways to a complex antioxidant defence mechanism. Conclusions: This is the first large scale RNA-Seq study carried out in S. glomerata, showing the complex network of innate immune components that exist in this species. The results confirmed that many of the innate immune system components observed in mammals are also conserved in oysters; however, some, such as the TLR adaptors MAL, TRIF and TRAM are either missing or have been modified significantly. The components identified in this study could help explain the oysters' natural resilience against pathogenic microorganisms encountered in their natural environment.

31 p.-11 fig.-2 tab.+ Erratum (2 p.) Papanikolaou, Alexie et al., Background: The Mediterranean fruit fly (medfly), Ceratitis capitata, is a major destructive insect pest due to its broad host range, which includes hundreds of fruits and vegetables. It exhibits a unique ability to invade and adapt to ecological niches throughout tropical and subtropical regions of the world, though medfly infestations have been prevented and controlled by the sterile insect technique (SIT) as part of integrated pest management programs (IPMs). The genetic analysis and manipulation of medfly has been subject to intensive study in an effort to improve SIT efficacy and other aspects of IPM control., Results: The 479 Mb medfly genome is sequenced from adult flies from lines inbred for 20 generations. A highquality assembly is achieved having a contig N50 of 45.7 kb and scaffold N50 of 4.06 Mb. In-depth curation of more than 1800 messenger RNAs shows specific gene expansions that can be related to invasiveness and host adaptation, including gene families for chemoreception, toxin and insecticide metabolism, cuticle proteins, opsins, and aquaporins. We identify genes relevant to IPM control, including those required to improve SIT., Conclusions: The medfly genome sequence provides critical insights into the biology of one of the most serious and widespread agricultural pests. This knowledge should significantly advance the means of controlling the size and invasive potential of medfly populations. Its close relationship to Drosophila, and other insect species important to agriculture and human health, will further comparative functional and structural studies of insect genomes that should broaden our understanding of gene family evolution, Support of this project was provided by the U.S. Department of Agriculture(USDA), Agricultural Research Service (ARS), Animal and Plant Health Inspection Service (APHIS), and National Institute of Food and Agriculture(NIFA)-Biotechnology Risk Assessment Grants Program (grant #2011-39211-30769 to AMH) for funding the initial phase of this project, and to the National Institutes of Health (NIH)-National Human Genome Research Institute (NHGRI) for funding the medfly genome sequencing, assembly and Maker 2.0 automated annotation as part of the i5K 30 genome pilot project (grant #U54 HG003273 to RAG). The NIH Intramural Research Program, National Library of Medicine funded the NCBI Gnomon annotation and the USDA-National Agricultural Library (NAL) provided support for the WebApollo curation website, with support for manual curation training (to MM-T) provided by NIGMS (grant #5R01GM080203),NHGRI (grant #5R01HG004483), and the U.S. Department of Energy(contract #DE-AC02-05CH11231). Support was provided for: toxin metabolism and insecticide resistance gene studies from MINECO,Spain (AGL2013-42632-R to FO and PH-C); microRNAs, horizontal gene transfer and bacterial contaminant studies from the European Social Fund and National Strategic Reference Framework-THALES (MIS375869 to KB, GT, AGH, and KM) and the U.S. National Science Foundation(DEB 1257053 to JHW); cuticle protein gene studies from USDA-NIFA(grant #2016-67012-24652 to AJR); sex-determination studies from L.R. Campania (grant 5/02, 2008 to GS); male reproduction and sexual differentiation studies from the FAO/IAEA (Technical Contract No: 16966 to GGa) and Cariplo IMPROVE (to FS); and programmed cell death gene studies and genomic data analysis (to MFS) from the Emmy Noether program, DFG(SCHE 1833/1-1) and the LOEWE Center for Insect Biotechnology & Bioresources grant of the Hessen State Ministry of Higher Education, Research and the Arts(HMWK), Germany and from the USDA-NIFA-Biotechnology Risk Assessment Grants Program (grant #2015-33522-24094 to AMH).

Motivation: Next-generation sequencing technologies have led to the widespread use of -omic applications. As a result, there is now a pronounced bioinformatic bottleneck. The general model organism database (GMOD) tool kit (http://gmod.org) has produced a number of resources aimed at addressing this issue. It lacks, however, a robust online solution that can deploy heterogeneous data and software within a Web content management system (CMS). Results: We present a bioinformatic framework for the Drupal CMS. It consists of three modules. First, GMOD-DBSF is an application programming interface module for the Drupal CMS that simplifies the programming of bioinformatic Drupal modules. Second, the Drupal Bioinformatic Software Bench (biosoftware_bench) allows for a rapid and secure deployment of bioinformatic software. An innovative graphical user interface (GUI) guides both use and administration of the software, including the secure provision of pre-publication datasets. Third, we present genes4all_experiment, which exemplifies how our work supports the wider research community. Conclusion: Given the infrastructure presented here, the Drupal CMS may become a powerful new tool set for bioinformaticians. The GMOD-DBSF base module is an expandable community resource that decreases development time of Drupal modules for bioinformatics. The biosoftware_bench module can already enhance biologists' ability to mine their own data. The genes4all_experiment module has already been responsible for archiving of more than 150 studies of RNAi from Lepidoptera, which were previously unpublished. Availability and implementation: Implemented in PHP and Perl. Freely available under the GNU Public License 2 or later from http://gmod-dbsf.googlecode.com Contact: alexie@butterflybase.org

Technological and conceptual advances of the last decade have led to an explosion of genomic data and the emergence of new research avenues. Evolutionary and ecological functional genomics, with its focus on the genes that affect ecological success and adaptation in natural populations, benefits immensely from a phylogenetically widespread sampling of biological patterns and processes. Among those organisms outside established model systems, butterflies offer exceptional opportunities for multidisciplinary research on the processes generating and maintaining variation in ecologically relevant traits. Here we highlight research on wing color pattern variation in two groups of Nymphalid butterflies, the African species Bicyclus anynana (subfamily Satyrinae) and species of the South American genus Heliconius (subfamily Heliconiinae), which are emerging as important systems for studying the nature and origins of functional diversity. Growing genomic resources including genomic and cDNA libraries, dense genetic maps, high-density gene arrays, and genetic transformation techniques are extending current gene mapping and expression profiling analysis and enabling the next generation of research questions linking genes, development, form, and fitness. Efforts to develop such resources in Bicyclus and Heliconius underscore the general challenges facing the larger research community and highlight the need for a community-wide effort to extend ongoing functional genomic research on butterflies.

Motivation: RNA interference (RNAi) technology is being developed as a weapon for pest insect control. To maximize the specificity that such an approach affords we have developed a bioinformatic web tool that searches the ever-growing arthropod transcriptome databases so that pest-specific RNAi sequences can be identified. This will help technology developers finesse the design of RNAi sequences and suggests which non-target species should be assessed in the risk assessment process. Availability and implementation: http://rnai.specifly.org. Contact: crobin@unimelb.edu.au

Background Blowflies have relevance in areas of forensic science, agriculture, and medicine, primarily due to the ability of their larvae to develop on flesh. While it is widely accepted that blowflies rely heavily on olfaction for identifying and locating hosts, there is limited research regarding the underlying molecular mechanisms. Using next generation sequencing (Illumina), this research examined the antennal transcriptome of Calliphora stygia (Fabricius) (Diptera: Calliphoridae) to identify members of the major chemosensory gene families necessary for olfaction. Results Representative proteins from all chemosensory gene families essential in insect olfaction were identified in the antennae of the blowfly C. stygia, including 50 odorant receptors, 22 ionotropic receptors, 21 gustatory receptors, 28 odorant binding proteins, 4 chemosensory proteins, and 3 sensory neuron membrane proteins. A total of 97 candidate cytochrome P450s and 39 esterases, some of which may act as odorant degrading enzymes, were also identified. Importantly, co-receptors necessary for the proper function of ligand-binding receptors were identified. Putative orthologues for the conserved antennal ionotropic receptors and candidate gustatory receptors for carbon dioxide detection were also amongst the identified proteins. Conclusions This research provides a comprehensive novel resource that will be fundamental for future studies regarding blowfly olfaction. Such information presents potential benefits to the forensic, pest control, and medical areas, and could assist in the understanding of insecticide resistance and targeted control through cross-species comparisons. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1466-8) contains supplementary material, which is available to authorized users.

We report here the first genome assembly and annotation of the human-pathogenic fungus Scedosporium aurantiacum, with a predicted 10,525 genes, and 11,661 transcripts. The strain WM 09.24 was isolated from the environment at Circular Quay, Sydney, New South Wales, Australia.

Migration is a key life history strategy for many animals and requires a suite of behavioural, morphological and physiological adaptations which together form the ‘migratory syndrome’. Genetic variation has been demonstrated for many traits that make up this syndrome, but the underlying genes involved remain elusive. Recent studies investigating migration‐associated genes have focussed on sampling migratory and nonmigratory populations from different geographic locations but have seldom explored phenotypic variation in a migratory trait. Here, we use a novel combination of tethered flight and next‐generation sequencing to determine transcriptomic differences associated with flight activity in a globally invasive moth pest, the cotton bollworm Helicoverpa armigera. By developing a state‐of‐the‐art phenotyping platform, we show that field‐collected H. armigera display continuous variation in flight performance with individuals capable of flying up to 40 km during a single night. Comparative transcriptomics of flight phenotypes drove a gene expression analysis to reveal a suite of expressed candidate genes which are clearly related to physiological adaptations required for long‐distance flight. These include genes important to the mobilization of lipids as flight fuel, the development of flight muscle structure and the regulation of hormones that influence migratory physiology. We conclude that the ability to express this complex set of pathways underlines the remarkable flexibility of facultative insect migrants to respond to deteriorating conditions in the form of migratory flight and, more broadly, the results provide novel insights into the fundamental transcriptional changes required for migration in insects and other taxa.

The African Mocker Swallowtail,Papilio dardanus, is a textbook example in evolutionary genetics. Classical breeding experiments have shown that wing pattern variation in this polymorphic Batesian mimic is determined by the polyallelicHlocus that controls a set of distinct mimetic phenotypes. Using bacterial artificial chromosome (BAC) sequencing, recombination analyses and comparative genomics, we show thatHco-segregates with an interval of less than 500 kb that is collinear with two other Lepidoptera genomes and contains 24 genes, including the transcription factor genesengrailed(en) andinvected(inv).His located in a region of conserved gene order, which argues against any role for genomic translocations in the evolution of a hypothesized multi-gene mimicry locus. Natural populations ofP. dardanusshow significant associations of specific morphs with single nucleotide polymorphisms (SNPs), centred onen. In addition, SNP variation in theHregion reveals evidence of non-neutral molecular evolution in theengene alone. We find evidence for a duplication potentially driving physical constraints on recombination in thelambornimorph. Absence of perfect linkage disequilibrium between different genes in the other morphs suggests thatHis limited to nucleotide positions in the regulatory and coding regions ofen. Our results therefore support the hypothesis that a single gene underlies wing pattern variation inP. dardanus.

Background Chemosensory receptors including olfactory receptors (ORs), gustatory receptors (GRs) and ionotropic receptors (IRs) play a central role in sensing chemical signals and guiding insect behaviours, and are potential target genes in insect pest control. The cotton bollworm Helicoverpa armigera is one of the most destructive pest species that can feed on over 200 different plant species. This diversity of host plants is likely linked to a complex chemosensory system. Here we built on previous work to characterize crucial chemosensory tissues linked to environmental interactions including larval antennae, larval mouthparts and larval fat bodies, as well as male and female adult heads, male and female adult tarsi, and female abdomens. Results Using transcriptome sequencing, Trinity RNA-seq assemblies and extensive manual curation, we identified a total of 91 candidate chemosensory receptors (60 candidate ORs, 10 GRs and 21 IRs). Thirty-five of these candidates present full-length transcripts. First, we performed in silico differential expression analysis on different sequenced tissues. Further, we created extensive expression profiles using reverse transcription (RT)-PCR on a variety of adult and larval stages. We found that the expression profile of HarmOR51 was limited to adult male antenna suggesting a role in mating that was further supported by a phylogenetic analysis clustering it into the pheromone receptor clade. HarmOR51 in calcium imaging analysis did not show responses to either of the two H. armigera sex pheromone components (Z9-16:Ald or Z11-16:Ald) inviting a future detailed study. In addition, we found four novel HarmORs (OR1, 53, 54 and 58) that appeared to be larvae-antennal specific. Finally, our expression profiling showed that four “divergent” HarmIRs (IR2, 7d.1, 7d.2 and 7d.3) were expressed in both adult and larval antennae, suggesting a functional divergence from their Drosophila homologues. Conclusions This study explored three chemoreceptor superfamily genes using a curated transcriptomic approach coupled with extensive expression profiling and a more limited functional characterization. Our results have now provided an extensive resource for investigating the chemoreceptor complement of this insect pest, and meanwhile allow for targeted experiments to identify potential molecular targets for pest control and to investigate insect-plant interactions. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-15-597) contains supplementary material, which is available to authorized users.

De novo assembly of RNA-seq data enables researchers to study transcriptomes without the need for a genome sequence; this approach can be usefully applied, for instance, in research on 'non-model organisms' of ecological and evolutionary importance, cancer samples or the microbiome. In this protocol we describe the use of the Trinity platform for de novo transcriptome assembly from RNA-seq data in non-model organisms. We also present Trinity-supported companion utilities for downstream applications, including RSEM for transcript abundance estimation, R/Bioconductor packages for identifying differentially expressed transcripts across samples and approaches to identify protein-coding genes. In the procedure, we provide a workflow for genome-independent transcriptome analysis leveraging the Trinity platform. The software, documentation and demonstrations are freely available from http://trinityrnaseq.sourceforge.net. The run time of this protocol is highly dependent on the size and complexity of data to be analyzed. The example data set analyzed in the procedure detailed herein can be processed in less than 5 h.

Background Heliconius butterfly wing pattern diversity offers a unique opportunity to investigate how natural genetic variation can drive the evolution of complex adaptive phenotypes. Positional cloning and candidate gene studies have identified a handful of regulatory and pigmentation genes implicated in Heliconius wing pattern variation, but little is known about the greater developmental networks within which these genes interact to pattern a wing. Here we took a large-scale transcriptomic approach to identify the network of genes involved in Heliconius wing pattern development and variation. This included applying over 140 transcriptome microarrays to assay gene expression in dissected wing pattern elements across a range of developmental stages and wing pattern morphs of Heliconius erato. Results We identified a number of putative early prepattern genes with color-pattern related expression domains. We also identified 51 genes differentially expressed in association with natural color pattern variation. Of these, the previously identified color pattern “switch gene” optix was recovered as the first transcript to show color-specific differential expression. Most differentially expressed genes were transcribed late in pupal development and have roles in cuticle formation or pigment synthesis. These include previously undescribed transporter genes associated with ommochrome pigmentation. Furthermore, we observed upregulation of melanin-repressing genes such as ebony and Dat1 in non-melanic patterns. Conclusions This study identifies many new genes implicated in butterfly wing pattern development and provides a glimpse into the number and types of genes affected by variation in genes that drive color pattern evolution.