No evidence of quantitative signal honesty across species of aposematic burnet moths (Lepidoptera: Zygaenidae)

This is the author accepted manuscript. The final version is available from Wiley via the DOI in this record.Many defended species use conspicuous visual warning signals to deter potential predators from attacking. Traditional theory holds that these signals should converge on similar forms, yet variation in visual traits and the levels of defensive chemicals is common, both within and between species. It is currently unclear how the strength of signals and potency of defences might be related: conflicting theories suggest that aposematic signals should be quantitatively honest, or, in contrast, that investment in one component should be prioritised over the other, while empirical tests have yielded contrasting results. Here, we advance this debate by examining the relationship between defensive chemicals and signal properties in a family of aposematic Lepidoptera, accounting for phylogenetic relationships and quantifying coloration from the perspective of relevant predators. We test for correlations between toxin levels and measures of wing colour across 14 species of day-flying burnet and forester moths (Lepidoptera: Zygaenidae), protected by highly aversive cyanogenic glucosides, and find no clear evidence of quantitative signal honesty. Significant relationships between toxin levels and coloration vary between sexes and sampling years, and several trends run contrary to expectations for signal honesty. Although toxin concentration is positively correlated with increasing luminance contrast in forewing pattern in one year, higher toxin levels are also associated with paler and less chromatically salient markings, at least in females, in another year. Our study also serves to highlight important factors, including sex-specific trends and seasonal variation, that should be accounted for in future work on signal honesty in aposematic species.The authors were funded by the BBSRC (SWBio DTP studentship, ref. 1355867) and Danish Council for Independent Research (DFF–1323-00088)

Shotgun sequencing of Yersinia enterocolitica strain W22703 (biotype 2, serotype O:9): genomic evidence for oscillation between invertebrates and mammals

<p>Abstract</p> <p>Background</p> <p><it>Yersinia enterocolitica </it>strains responsible for mild gastroenteritis in humans are very diverse with respect to their metabolic and virulence properties. Strain W22703 (biotype 2, serotype O:9) was recently identified to possess nematocidal and insecticidal activity. To better understand the relationship between pathogenicity towards insects and humans, we compared the W22703 genome with that of the highly pathogenic strain 8081 (biotype1B; serotype O:8), the only <it>Y. enterocolitica </it>strain sequenced so far.</p> <p>Results</p> <p>We used whole-genome shotgun data to assemble, annotate and analyse the sequence of strain W22703. Numerous factors assumed to contribute to enteric survival and pathogenesis, among them osmoregulated periplasmic glucan, hydrogenases, cobalamin-dependent pathways, iron uptake systems and the <it>Yersinia </it>genome island 1 (YGI-1) involved in tight adherence were identified to be common to the 8081 and W22703 genomes. However, sets of ~550 genes revealed to be specific for each of them in comparison to the other strain. The plasticity zone (PZ) of 142 kb in the W22703 genome carries an ancient flagellar cluster Flg-2 of ~40 kb, but it lacks the pathogenicity island YAPI_Ye, the secretion system <it>ysa </it>and <it>yts1</it>, and other virulence determinants of the 8081 PZ. Its composition underlines the prominent variability of this genome region and demonstrates its contribution to the higher pathogenicity of biotype 1B strains with respect to W22703. A novel type three secretion system of mosaic structure was found in the genome of W22703 that is absent in the sequenced strains of the human pathogenic <it>Yersinia </it>species, but conserved in the genomes of the apathogenic species. We identified several regions of differences in W22703 that mainly code for transporters, regulators, metabolic pathways, and defence factors.</p> <p>Conclusion</p> <p>The W22703 sequence analysis revealed a genome composition distinct from other pathogenic <it>Yersinia enterocolitica </it>strains, thus contributing novel data to the <it>Y. enterocolitica </it>pan-genome. This study also sheds further light on the strategies of this pathogen to cope with its environments.</p

Parasympathetic Activity and Blood Catecholamine Responses Following a Single Partial-Body Cryostimulation and a Whole-Body Cryostimulation

The aim of this study was to compare the effects of a single whole-body cryostimulation (WBC) and a partial-body cryostimulation (PBC) (i.e., not exposing the head to cold) on indices of parasympathetic activity and blood catecholamines. Two groups of 15 participants were assigned either to a 3-min WBC or PBC session, while 10 participants constituted a control group (CON) not receiving any cryostimulation. Changes in thermal, physiological and subjective variables were recorded before and during the 20-min after each cryostimulation. According to a qualitative statistical analysis, an almost certain decrease in skin temperature was reported for all body regions immediately after the WBC (mean decrease±90% CL, -13.7±0.7°C) and PBC (-8.3±0.3°C), which persisted up to 20-min after the session. The tympanic temperature almost certainly decreased only after the WBC session (-0.32±0.04°C). Systolic and diastolic blood pressures were very likely increased after the WBC session, whereas these changes were trivial in the other groups. In addition, heart rate almost certainly decreased after PBC (-10.9%) and WBC (-15.2%) sessions, in a likely greater proportion for WBC compared to PBC. Resting vagal-related heart rate variability indices (the root-mean square difference of successive normal R-R intervals, RMSSD, and high frequency band, HF) were very likely increased after PBC (RMSSD: +54.4%, HF: +138%) and WBC (RMSSD: +85.2%, HF: +632%) sessions without any marked difference between groups. Plasma norepinephrine concentrations were likely to very likely increased after PBC (+57.4%) and WBC (+76.2%), respectively. Finally, cold and comfort sensations were almost certainly altered after WBC and PBC, sensation of discomfort being likely more pronounced after WBC than PBC. Both acute cryostimulation techniques effectively stimulated the autonomic nervous system (ANS), with a predominance of parasympathetic tone activation. The results of this study also suggest that a whole-body cold exposure induced a larger stimulation of the ANS compared to partial-body cold exposure

A de novo transcriptome of the Asian tiger mosquito, Aedes albopictus, to identify candidate transcripts for diapause preparation

<p>Abstract</p> <p>Background</p> <p>Many temperate insects survive the harsh conditions of winter by undergoing photoperiodic diapause, a pre-programmed developmental arrest initiated by short day lengths. Despite the well-established ecological significance of photoperiodic diapause, the molecular basis of this crucial adaptation remains largely unresolved. The Asian tiger mosquito, <it>Aedes albopictus </it>(Skuse), represents an outstanding emerging model to investigate the molecular basis of photoperiodic diapause in a well-defined ecological and evolutionary context. <it>Ae. albopictus </it>is a medically significant vector and is currently considered the most invasive mosquito in the world. Traits related to diapause appear to be important factors contributing to the rapid spread of this mosquito. To generate novel sequence information for this species, as well as to discover transcripts involved in diapause preparation, we sequenced the transcriptome of <it>Ae. albopictus </it>oocytes destined to become diapausing or non-diapausing pharate larvae.</p> <p>Results</p> <p>454 GS-FLX transcriptome sequencing yielded >1.1 million quality-filtered reads, which we assembled into 69,474 contigs (N50 = 1,009 bp). Our contig filtering approach, where we took advantage of strong sequence similarity to the fully sequenced genome of <it>Aedes aegypti</it>, as well as other reference organisms, resulted in 11,561 high-quality, conservative ESTs. Differential expression estimates based on normalized read counts revealed 57 genes with higher expression, and 257 with lower expression under diapause-inducing conditions. Analysis of expression by qPCR for 47 of these genes indicated a high correlation of expression levels between 454 sequence data and qPCR, but congruence of statistically significant differential expression was low. Seven genes identified as differentially expressed based on qPCR have putative functions that are consistent with the insect diapause syndrome; three genes have unknown function and represent novel candidates for the transcriptional basis of diapause.</p> <p>Conclusions</p> <p>Our transcriptome database provides a rich resource for the comparative genomics and functional genetics of <it>Ae. albopictus</it>, an invasive and medically important mosquito. Additionally, the identification of differentially expressed transcripts related to diapause enriches the limited knowledge base for the molecular basis of insect diapause, in particular for the preparatory stage. Finally, our analysis illustrates a useful approach that draws from a closely related reference genome to generate high-confidence ESTs in a non-model organism.</p

Comparing de novo assemblers for 454 transcriptome data

<p>Abstract</p> <p>Background</p> <p>Roche 454 pyrosequencing has become a method of choice for generating transcriptome data from non-model organisms. Once the tens to hundreds of thousands of short (250-450 base) reads have been produced, it is important to correctly assemble these to estimate the sequence of all the transcripts. Most transcriptome assembly projects use only one program for assembling 454 pyrosequencing reads, but there is no evidence that the programs used to date are optimal. We have carried out a systematic comparison of five assemblers (CAP3, MIRA, Newbler, SeqMan and CLC) to establish best practices for transcriptome assemblies, using a new dataset from the parasitic nematode <it>Litomosoides sigmodontis</it>.</p> <p>Results</p> <p>Although no single assembler performed best on all our criteria, Newbler 2.5 gave longer contigs, better alignments to some reference sequences, and was fast and easy to use. SeqMan assemblies performed best on the criterion of recapitulating known transcripts, and had more novel sequence than the other assemblers, but generated an excess of small, redundant contigs. The remaining assemblers all performed almost as well, with the exception of Newbler 2.3 (the version currently used by most assembly projects), which generated assemblies that had significantly lower total length. As different assemblers use different underlying algorithms to generate contigs, we also explored merging of assemblies and found that the merged datasets not only aligned better to reference sequences than individual assemblies, but were also more consistent in the number and size of contigs.</p> <p>Conclusions</p> <p>Transcriptome assemblies are smaller than genome assemblies and thus should be more computationally tractable, but are often harder because individual contigs can have highly variable read coverage. Comparing single assemblers, Newbler 2.5 performed best on our trial data set, but other assemblers were closely comparable. Combining differently optimal assemblies from different programs however gave a more credible final product, and this strategy is recommended.</p

Transcriptome Sequencing and De Novo Analysis for Yesso Scallop (Patinopecten yessoensis) Using 454 GS FLX

BACKGROUND: Bivalves comprise 30,000 extant species, constituting the second largest group of mollusks. However, limited genetic research has focused on this group of animals so far, which is, in part, due to the lack of genomic resources. The advent of high-throughput sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provides a turning point for bivalve research. In the present study, we performed de novo transcriptome sequencing to first produce a comprehensive expressed sequence tag (EST) dataset for the Yesso scallop (Patinopecten yessoensis). RESULTS: In a single 454 sequencing run, 805,330 reads were produced and then assembled into 32,590 contigs, with about six-fold sequencing coverage. A total of 25,237 unique protein-coding genes were identified from a variety of developmental stages and adult tissues based on sequence similarities with known proteins. As determined by GO annotation and KEGG pathway mapping, functional annotation of the unigenes recovered diverse biological functions and processes. Transcripts putatively involved in growth, reproduction and stress/immune-response were identified. More than 49,000 single nucleotide polymorphisms (SNPs) and 2,700 simple sequence repeats (SSRs) were also detected. CONCLUSION: Our data provide the most comprehensive transcriptomic resource currently available for P. yessoensis. Candidate genes potentially involved in growth, reproduction, and stress/immunity-response were identified, and are worthy of further investigation. A large number of SNPs and SSRs were also identified and ready for marker development. This resource should lay an important foundation for future genetic or genomic studies on this species

Pyrosequencing the Bemisia tabaci Transcriptome Reveals a Highly Diverse Bacterial Community and a Robust System for Insecticide Resistance

BACKGROUND: Bemisia tabaci (Gennadius) is a phloem-feeding insect poised to become one of the major insect pests in open field and greenhouse production systems throughout the world. The high level of resistance to insecticides is a main factor that hinders continued use of insecticides for suppression of B. tabaci. Despite its prevalence, little is known about B. tabaci at the genome level. To fill this gap, an invasive B. tabaci B biotype was subjected to pyrosequencing-based transcriptome analysis to identify genes and gene networks putatively involved in various physiological and toxicological processes. METHODOLOGY AND PRINCIPAL FINDINGS: Using Roche 454 pyrosequencing, 857,205 reads containing approximately 340 megabases were obtained from the B. tabaci transcriptome. De novo assembly generated 178,669 unigenes including 30,980 from insects, 17,881 from bacteria, and 129,808 from the nohit. A total of 50,835 (28.45%) unigenes showed similarity to the non-redundant database in GenBank with a cut-off E-value of 10-5. Among them, 40,611 unigenes were assigned to one or more GO terms and 6,917 unigenes were assigned to 288 known pathways. De novo metatranscriptome analysis revealed highly diverse bacterial symbionts in B. tabaci, and demonstrated the host-symbiont cooperation in amino acid production. In-depth transcriptome analysis indentified putative molecular markers, and genes potentially involved in insecticide resistance and nutrient digestion. The utility of this transcriptome was validated by a thiamethoxam resistance study, in which annotated cytochrome P450 genes were significantly overexpressed in the resistant B. tabaci in comparison to its susceptible counterparts. CONCLUSIONS: This transcriptome/metatranscriptome analysis sheds light on the molecular understanding of symbiosis and insecticide resistance in an agriculturally important phloem-feeding insect pest, and lays the foundation for future functional genomics research of the B. tabaci complex. Moreover, current pyrosequencing effort greatly enriched the existing whitefly EST database, and makes RNAseq a viable option for future genomic analysis

De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum

<p>Abstract</p> <p>Background</p> <p>Because of their phylogenetic position and unique characteristics of their biology and life cycle, ferns represent an important lineage for studying the evolution of land plants. Large and complex genomes in ferns combined with the absence of economically important species have been a barrier to the development of genomic resources. However, high throughput sequencing technologies are now being widely applied to non-model species. We leveraged the Roche 454 GS-FLX Titanium pyrosequencing platform in sequencing the gametophyte transcriptome of bracken fern (<it>Pteridium aquilinum</it>) to develop genomic resources for evolutionary studies.</p> <p>Results</p> <p>681,722 quality and adapter trimmed reads totaling 254 Mbp were assembled <it>de novo </it>into 56,256 unique sequences (i.e. unigenes) with a mean length of 547.2 bp and a total assembly size of 30.8 Mbp with an average read-depth coverage of 7.0×. We estimate that 87% of the complete transcriptome has been sequenced and that all transcripts have been tagged. 61.8% of the unigenes had blastx hits in the NCBI nr protein database, representing 22,596 unique best hits. The longest open reading frame in 52.2% of the unigenes had positive domain matches in InterProScan searches. We assigned 46.2% of the unigenes with a GO functional annotation and 16.0% with an enzyme code annotation. Enzyme codes were used to retrieve and color KEGG pathway maps. A comparative genomics approach revealed a substantial proportion of genes expressed in bracken gametophytes to be shared across the genomes of <it>Arabidopsis</it>, <it>Selaginella </it>and <it>Physcomitrella</it>, and identified a substantial number of potentially novel fern genes. By comparing the list of <it>Arabidopsis </it>genes identified by blast with a list of gametophyte-specific <it>Arabidopsis </it>genes taken from the literature, we identified a set of potentially conserved gametophyte specific genes. We screened unigenes for repetitive sequences to identify 548 potentially-amplifiable simple sequence repeat loci and 689 expressed transposable elements.</p> <p>Conclusions</p> <p>This study is the first comprehensive transcriptome analysis for a fern and represents an important scientific resource for comparative evolutionary and functional genomics studies in land plants. We demonstrate the utility of high-throughput sequencing of a normalized cDNA library for <it>de novo </it>transcriptome characterization and gene discovery in a non-model plant.</p

Transcriptome analysis of Bupleurum chinense focusing on genes involved in the biosynthesis of saikosaponins

<p>Abstract</p> <p>Background</p> <p><it>Bupleurum chinense </it>DC. is a widely used traditional Chinese medicinal plant. Saikosaponins are the major bioactive constituents of <it>B. chinense</it>, but relatively little is known about saikosaponin biosynthesis. The 454 pyrosequencing technology provides a promising opportunity for finding novel genes that participate in plant metabolism. Consequently, this technology may help to identify the candidate genes involved in the saikosaponin biosynthetic pathway.</p> <p>Results</p> <p>One-quarter of the 454 pyrosequencing runs produced a total of 195, 088 high-quality reads, with an average read length of 356 bases (NCBI SRA accession SRA039388). A <it>de novo </it>assembly generated 24, 037 unique sequences (22, 748 contigs and 1, 289 singletons), 12, 649 (52.6%) of which were annotated against three public protein databases using a basic local alignment search tool (E-value ≤1e-10). All unique sequences were compared with NCBI expressed sequence tags (ESTs) (237) and encoding sequences (44) from the <it>Bupleurum </it>genus, and with a Sanger-sequenced EST dataset (3, 111). The 23, 173 (96.4%) unique sequences obtained in the present study represent novel <it>Bupleurum </it>genes. The ESTs of genes related to saikosaponin biosynthesis were found to encode known enzymes that catalyze the formation of the saikosaponin backbone; 246 cytochrome P450 (<it>P450</it>s) and 102 glycosyltransferases (<it>GT</it>s) unique sequences were also found in the 454 dataset. Full length cDNAs of 7 <it>P450</it>s and 7 uridine diphosphate <it>GT</it>s (<it>UGT</it>s) were verified by reverse transcriptase polymerase chain reaction or by cloning using 5' and/or 3' rapid amplification of cDNA ends. Two <it>P450</it>s and three <it>UGT</it>s were identified as the most likely candidates involved in saikosaponin biosynthesis. This finding was based on the coordinate up-regulation of their expression with <it>β-AS </it>in methyl jasmonate-treated adventitious roots and on their similar expression patterns with <it>β-AS </it>in various <it>B. chinense </it>tissues.</p> <p>Conclusions</p> <p>A collection of high-quality ESTs for <it>B. chinense </it>obtained by 454 pyrosequencing is provided here for the first time. These data should aid further research on the functional genomics of <it>B. chinense </it>and other <it>Bupleurum </it>species. The candidate genes for enzymes involved in saikosaponin biosynthesis, especially the <it>P450</it>s and <it>UGT</it>s, that were revealed provide a substantial foundation for follow-up research on the metabolism and regulation of the saikosaponins.</p